Engineering topics

Lyra, Wladimir

Publications and source records attributed to Lyra, Wladimir.

9 records

NASA NTRS2018

The Interplay Between Radiation Pressure and the Photoelectric Instability in Optically Thin Disks of Gas and Dust

In optically thin disks, dust grains are photoelectrically stripped of electrons by starlight, heating nearby gas and possibly creating a dust clumping instability-the photoelectric instability (PeI)-that significantly alters global disk structure. In the current work, we use the Pencil Code to perform the first numerical models of the PeI that include stellar radiation pressure on dust grains in order to explore the parameter regime in which the instability operates. In some models with low gas and dust surface densities, we see a variety of dust structures, including sharp concentric rings. In the most gas- and dust-rich models, nonaxisymmetric clumps, arcs, and spiral arms emerge that represent dust surface density enhancements of factors of ∼5-20. In one high gas surface density model, we include a large, low-order gas viscosity and find that it observably smooths the structures that form in the gas and dust, suggesting that resolved images of a given disk may be useful for deriving constraints on the effective viscosity of its gas. Our models show that radiation pressure does not preclude the formation of complex structure from the PeI, but the qualitative manifestation of the PeI depends strongly on the parameters of the system. The PeI may provide an explanation for unusual disk morphologies, such as the moving blobs of the AU Mic disk, the asymmetric dust distribution of the 49Ceti disk, and the rings and arcs found in the HD 141569A disk.

photoelectric instability (PeI)↗

DOE OSTI2017

Orbital Advection with Magnetohydrodynamics and Vector Potential

Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear terms, and we find that we do not need an advective gauge since solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited to non-magnetic problems. Benchmarked results of (hydrodynamical) planet–disk interaction and of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high-frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficient for straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

ON SHOCKS DRIVEN BY HIGH-MASS PLANETS IN RADIATIVELY INEFFICIENT DISKS. II. THREE-DIMENSIONAL GLOBAL DISK SIMULATIONS

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted by models of disk–planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for these signatures. However, such interpretation is not free of problems. The observed spirals have large pitch angles, and in at least one case (HD 100546) it appears effectively unpolarized, implying thermal emission of the order of 1000 K (465 ± 40 K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. Here we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5M{sub J} planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient gas. The gas is accelerated vertically away from the midplane to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf. This turbulence, although localized, has high α values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We find evidence that the disk regions heated up by the shocks become superadiabatic, generating convection far from the planet’s orbit.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

On shocks driven by high-mass planets in radiatively inefficient disks. I. Two-dimensional global disk simulations

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet–disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1–10 M{sub J} for a 1 M{sub ⊙} star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets (M≳5 M{sub J}) produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a 10 M{sub J} planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet–disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

PARTICLE TRAPPING AND STREAMING INSTABILITY IN VORTICES IN PROTOPLANETARY DISKS

We analyze the concentration of solid particles in vortices created and sustained by radial buoyancy in protoplanetary disks, e.g., baroclinic vortex growth. Besides the gas drag acting on particles, we also allow for back-reaction from dust onto the gas. This becomes important when the local dust-to-gas ratio approaches unity. In our two-dimensional, local, shearing sheet simulations, we see high concentrations of grains inside the vortices for a broad range of Stokes numbers, St. An initial dust-to-gas ratio of 1:100 can easily be reversed to 100:1 for St = 1.0. The increased dust-to-gas ratio triggers the streaming instability, thus counter-intuitively limiting the maximal achievable overdensities. We find that particle trapping inside vortices opens the possibility for gravity assisted planetesimal formation even for small particles (St = 0.01) and a low initial dust-to-gas ratio of 1:10{sup 4}, e.g., much smaller than in the previously studied magnetohydrodynamic zonal flow case.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

CSI 2264: Characterizing Young Stars in NGC 2264 With Short-Duration Periodic Flux Dips in Their Light Curves

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical CoRoT light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian corotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3 < P < 11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard \disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous Spitzer IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel- ow accretion columns arising near the inner disk wall.

CSI 2264↗

DOE OSTI2014

Convective overstability in accretion disks: Three-dimensional linear analysis and nonlinear saturation

Recently, Klahr and Hubbard claimed that a hydrodynamical linear overstability exists in protoplanetary disks, powered by buoyancy in the presence of thermal relaxation. We analyze this claim, confirming it through rigorous compressible linear analysis. We model the system numerically, reproducing the linear growth rate for all cases studied. We also study the saturated properties of the overstability in the shearing box, finding that the saturated state produces finite amplitude fluctuations strong enough to trigger the subcritical baroclinic instability (SBI). Saturation leads to a fast burst of enstrophy in the box, and a large-scale vortex develops in the course of the next ≈100 orbits. The amount of angular momentum transport achieved is of the order of α ≈ 10{sup –3}, as in compressible SBI models. For the first time, a self-sustained three-dimensional vortex is produced from linear amplitude perturbation of a quiescent base state.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

STEADY STATE DUST DISTRIBUTIONS IN DISK VORTICES: OBSERVATIONAL PREDICTIONS AND APPLICATIONS TO TRANSITIONAL DISKS

The Atacama Large Millimeter Array has returned images of transitional disks in which large asymmetries are seen in the distribution of millimeter sized dust in the outer disk. The explanation in vogue borrows from the vortex literature and suggests that these asymmetries are the result of dust trapping in giant vortices, excited via Rossby wave instabilities at planetary gap edges. Due to the drag force, dust trapped in vortices will accumulate in the center and diffusion is needed to maintain a steady state over the lifetime of the disk. While previous work derived semi-analytical models of the process, in this paper we provide analytical steady-steady solutions. Exact solutions exist for certain vortex models. The solution is determined by the vortex rotation profile, the gas scale height, the vortex aspect ratio, and the ratio of dust diffusion to gas-dust friction. In principle, all of these quantities can be derived from observations, which would validate the model and also provide constrains on the strength of the turbulence inside the vortex core. Based on our solution, we derive quantities such as the gas-dust contrast, the trapped dust mass, and the dust contrast at the same orbital location. We apply our model to the recently imaged Oph IRS 48 system, finding values within the range of the observational uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

A PARAMETER STUDY FOR BAROCLINIC VORTEX AMPLIFICATION

Recent studies have shown that baroclinic vortex amplification is strongly dependent on certain factors, namely, the global entropy gradient, the efficiency of thermal diffusion and/or relaxation as well as numerical resolution. We conduct a comprehensive study of a broad range and combination of various entropy gradients, thermal diffusion and thermal relaxation timescales via local shearing sheet simulations covering the parameter space relevant for protoplanetary disks. We measure the Reynolds stresses as a function of our control parameters and see that there is angular momentum transport even for entropy gradients as low as {beta} = -dln s/dln r = 1/2. Values we expect in protoplanetary disks are between {beta} = 0.5-2.0 The amplification-rate of the perturbations, {Gamma}, appears to be proportional to {beta}{sup 2} and thus proportional to the square of the Brunt-Vaeisaelae frequency ({Gamma}{proportional_to}{beta}{sup 2}{proportional_to}N {sup 2}). The saturation level of Reynolds stresses, on the other hand, seems to be proportional to {beta}{sup 1/2}. This highlights the importance of baroclinic effects even for the low entropy gradients expected in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗