SEARCH · Engineering Papers

Results for “numerical protoplanetary disks”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Ring Formation in Protoplanetary Disks Driven by an Eccentric Instability

We find that, under certain conditions, protoplanetary disks may spontaneously generate multiple, concentric gas rings without an embedded planet through an eccentric cooling instability. Using both linear theory and nonlinear hydrodynamics simulations, we show that a variety of background states may trap a slowly processing, one-armed spiral mode that becomes unstable when a gravitationally stable disk rapidly cools. The angular momentum required to excite this spiral comes at the expense of nonuniform mass transport that generically results in multiple rings. For example, one long-term hydrodynamics simulation exhibits four long-lived, axisymmetric gas rings. We verify the instability evolution and ring-formation mechanism from first principles with our linear theory, which shows remarkable agreement with the simulation results. Dust trapped in these rings may produce observable features consistent with observed disks. Additionally, direct detection of the eccentric gas motions may be possible when the instability saturates, and any residual eccentricity left over in the rings at later times may also provide direct observational evidence of this mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2002

Vorticity Transport and Wave Emission In A Protoplanetary Disk

Higher order numerical algorithms (4th order in time, 3rd order in space) are applied to the Euler equations and are used to examine vorticity transport and wave motion in a non-self gravitating, initially isentropic Keplerian disk. In this talk we will examine the response of the disk to an isolated vortex with a circulation about equal to the rotation rate of Jupiter. The vortex is located on the 4 AU circle and the nebula is simulated from 1 to 24 AU. We show that the vortex emits pressure-supported density and Rossby-type wave packets before it decays within a few orbits. The acoustic density waves evolve into weak (non entropy preserving) shock waves that propagate over the entire disk. The Rossby waves remain in the vicinity of the initial vortex disturbance, but are rapidly damped. Temporal frequencies and spatial wavenumbers are derived from the nonlinear simulation data and correlated with analytical dispersion relations from the linearized Euler and energy equations.

Davis, S. S.↗

NASA NTRS2000

Annual Research Briefs - 2000: Center for Turbulence Research

This report contains the 2000 annual progress reports of the postdoctoral Fellows and visiting scholars of the Center for Turbulence Research (CTR). It summarizes the research efforts undertaken under the core CTR program. Last year, CTR sponsored sixteen resident Postdoctoral Fellows, nine Research Associates, and two Senior Research Fellows, hosted seven short term visitors, and supported four doctoral students. The Research Associates are supported by the Departments of Defense and Energy. The reports in this volume are divided into five groups. The first group largely consists of the new areas of interest at CTR. It includes efficient algorithms for molecular dynamics, stability in protoplanetary disks, and experimental and numerical applications of evolutionary optimization algorithms for jet flow control. The next group of reports is in experimental, theoretical, and numerical modeling efforts in turbulent combustion. As more challenging computations are attempted, the need for additional theoretical and experimental studies in combustion has emerged. A pacing item for computation of nonpremixed combustion is the prediction of extinction and re-ignition phenomena, which is currently being addressed at CTR. The third group of reports is in the development of accurate and efficient numerical methods, which has always been an important part of CTR's work. This is the tool development part of the program which supports our high fidelity numerical simulations in such areas as turbulence in complex geometries, hypersonics, and acoustics. The final two groups of reports are concerned with LES and RANS prediction methods. There has been significant progress in wall modeling for LES of high Reynolds number turbulence and in validation of the v(exp 2) - f model for industrial applications.

Source record↗

DOE OSTI2022

Evidence that the Hot Jupiter WASP-77 A b Formed Beyond Its Parent Protoplanetary Disk’s H2O Ice Line

Idealized protoplanetary disk and giant planet formation models have been interpreted to suggest that a giant planet's atmospheric abundances can be used to infer its formation location in its parent protoplanetary disk. It has recently been reported that the hot Jupiter WASP-77 A b has subsolar atmospheric carbon and oxygen abundances with a solar C/O abundance ratio. Assuming solar carbon and oxygen abundances for its host star WASP-77 A, WASP-77 A b's atmospheric carbon and oxygen abundances possibly indicate that it accreted its envelope interior to its parent protoplanetary disk's H 2 O ice line from carbon-depleted gas with little subsequent planetesimal accretion or core erosion. We show that the photospheric abundances of carbon and oxygen in WASP-77 A are supersolar with a subsolar C/O abundance ratio, implying that WASP-77 A b's atmosphere has significantly substellar carbon and oxygen abundances with a superstellar C/O ratio. Our result possibly indicates that WASP-77 A b's envelope was accreted by the planet beyond its parent protoplanetary disk's H 2 O ice line. While numerous theoretical complications to these idealized models have now been identified, the possibility of nonsolar protoplanetary disk abundance ratios confound even the most sophisticated protoplanetary disk and giant planet formation models. We therefore argue that giant planet atmospheric abundance ratios can only be meaningfully interpreted relative to the possibly nonsolar mean compositions of their parent protoplanetary disks as recorded in the photospheric abundances of their dwarf host stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Wide Dust Gaps in Protoplanetary Disks Induced by Eccentric Planets: A Mass-eccentricity Degeneracy

The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in submillimeter images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scaling of characteristic dust-gap width/gap-ring distance with respect to planet mass. Applying such scaling to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bimodal mass function due to (1) the onset of runaway gas accretion above ~20 Earth masses and (2) suppression of accretion induced by gap opening. Here, we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. Furthermore, this degeneracy may obliterate the discrepancy between the theoretical bimodal mass distribution and the observed continuous gap width distribution. Despite damping due to planet–disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in Atacama Large Millimeter/submillimeter Array images.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1996

Dust Coagulation in Protoplanetary Accretion Disks

The time evolution of dust particles in circumstellar disk-like structures around protostars and young stellar objects is discussed. In particular, we consider the coagulation of grains due to collisional aggregation. The coagulation of the particles is calculated by solving numerically the non-linear Smoluchowski equation. The different physical processes leading to relative velocities between the grains are investigated. The relative velocities may be induced by Brownian motion, turbulence and drift motion. Starting from different regimes which can be identified during the grain growth we also discuss the evolution of dust opacities. These opacities are important for both the derivation of the circumstellar dust mass from submillimeter/millimeter continuum observations and the dynamical behavior of the disks. We present results of our numerical studies of the coagulation of dust grains in a turbulent protoplanetary accretion disk described by a time-dependent one-dimensional (radial) alpha-model. For several periods and disk radii, mass distributions of coagulated grains have been calculated. From these mass spectra, we determined the corresponding Rosseland mean dust opacities. The influence of grain opacity changes due to dust coagulation on the dynamical evolution of a protostellar disk is considered. Significant changes in the thermal structure of the protoplanetary nebula are observed. A 'gap' in the accretion disk forms at the very frontier of the coagulation, i.e., behind the sublimation boundary in the region between 1 and 5 AU.

Schmitt, W.↗

NASA NTRS2003

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗

NASA NTRS2003

VLA Imaging of Protoplanetary Environments

Our program uses the techniques of millimeter interferometry to make high resolution observations of dust continuum emission to study the structure of protoplanetary disks around nearby stars. We continue to use the Very Large Array (VLA) of the National Radio Astronomy Observatories (NRAO) to make produce high angular resolution images, and we are improving our capabilities to interpret these observations with detailed numerical models. We have also made the first millimeter interferometer observations of protoplanetary disks in the southern hemisphere, using the recently upgraded Australia Telescope Compact Array (ATCA), taking advantage of the techniques we have developed in our work on northern sources.

Wilner, David J.↗

NASA NTRS1979

Dynamics of the accumulation process of the Earth group of planets: Formation of the reverse rotation of Venus

A numerical simulation of the process of formation of the terrestrial planets is carried within the framework of a new theory for the accumulation of planetary and satellite systems. The numerical simulation permitted determining the parameters of the protoplanetary disk from which Mercury, Venus and the Earth were formed as result of the evolution. The acquisition of a slow retrograde rotation for Venus was discovered during the course of the investigation, whereas Mercury and the Earth acquired direct rotation about their axes. Deviations of the semimajor axes of these three planets as well as the masses of the Earth and Venus from the true values are small as a rule (l 10%). It is shown that during the accumulation of the terrestrial planets, there existed a profound relationship between the process of formation of the orbits and masses of the planet and the process of formation of their rotation about their axes. Estimates are presented for the radii of the initial effective bodies and the time of evolution for the terrestrial accumulation zone.

Koslov, N. N.↗

NASA NTRS2001

Vorticity and Wave Motion in a Compressible Protoplanetary Disk

The impact of an isolated vortex in a compressible Keplerian disk is examined using higher order numerical solutions of the Euler and entropy-conserving Energy equations. The vortex is stretched by the background shear flow with longer lasting anticyclonic vortices persisting for about 10 vortex revolutions. Simultaneously, the vortex emits transient radial waves consisting mainly of axisymmetrical weak shock waves and a slower, nonaxisymmetric Rossby wave. These waves may contribute to certain transient events in protoplanetary disks. The vortex stretching and waves were found to have little long-term feedback on the baseline 'standard solar nebula' disk structure and confirm the extremely stable structure of non self-gravitating disks.

Davis, Sanford S.↗

DOE OSTI2024

Effects of Dust Coagulation on Streaming Instability

Abstract Streaming instability (SI) in dust has long been thought to be a promising process in triggering planetesimal formation in the protoplanetary disks (PPDs). In this study, we present the first numerical investigation that models the SI in the vertically stratified disk together with the dust coagulation process. Our simulations reveal that, even with the initially small dust sizes, because dust coagulation promotes dust size growth, SI can eventually still be triggered. Specifically, dust coagulation, limited only by dust fragmentation, broadens the parameter boundaries obtained from previous SI studies using single dust species. We describe the various stages of dust dynamics along with their size evolution and explore the impact of different dust fragmentation velocities. Implications of these results for realistic PPDs are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1988

Generation of dynamo magnetic fields in protoplanetary and other astrophysical accretion disks

A computational method for treating the generation of dynamo magnetic fields in astrophysical disks is presented. The numerical difficulty of handling the boundary condition at infinity in the cylindrical disk geometry is overcome by embedding the disk in a spherical computational space and matching the solutions to analytically tractable spherical functions in the surrounding space. The lowest lying dynamo normal modes for a 'thick' astrophysical disk are calculated. The generated modes found are all oscillatory and spatially localized. Tha potential implications of the results for the properties of dynamo magnetic fields in real astrophysical disks are discussed.

Stepinski, T. F.↗

DOE OSTI2020

Planet-induced Vortices with Dust Coagulation in Protoplanetary Disks

In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low-viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation, and drift (to the vortex center) are all quite efficient, producing dust particles ranging from 1 μm to ~1.0 cm, as well as an overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite nonuniform inside the vortex, with the mass-weighted average dust size at the vortex center (~4.0 mm) being a factor of ~10 larger than other vortex regions. Both large (~millimeter) and small (tens of microns) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. Furthermore, after the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size (~millimeter). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at the 1.33 mm band, while only distinct hot spots remain at 7.0 mm.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Fossil records of early solar irradiation and cosmolocation of the CAI factory: A reappraisal

Calcium-aluminum–rich inclusions (CAIs) in meteorites carry crucial information about the environmental conditions of the nascent Solar System prior to planet formation. Based on models of 50V– 10 Be co-production by in-situ irradiation, CAIs are considered to have formed within ~0.1 AU from the proto-Sun. Here, we present vanadium (V) and strontium (Sr) isotopic co-variations in fine- and coarse-grained CAIs and demonstrate that kinetic isotope effects during partial condensation and evaporation best explain V isotope anomalies previously attributed to solar particle irradiation. We also report initial excesses of 10 Be and argue that CV CAIs possess essentially a homogeneous level of 10 Be, inherited during their formation. Based on numerical modeling of 50V–10Be co-production by irradiation, we show that CAI formation during protoplanetary disk build-up likely occurred at greater heliocentric distances than previously considered, up to planet-forming regions (~1AU), where solar particle fluxes were sufficiently low to avoid substantial in-situ irradiation of CAIs.

58 GEOSCIENCES↗

NASA NTRS1997

Simulations of the STIS CCD clear imaging mode PSF

Detailed understanding of the properties of the point-spread function (PSF) of Space Telescope Imaging Spectrograph (STIS) imaging modes is necessary to separate real effects of geometrically complex astronomical sources, such as protoplanetary disks surrounding bright stars, from instrumental effects, such as scattering and diffraction. In order to investigate STIS imaging properties we have numerically simulated broadband stellar PSFs generated by STIS in CCD clear imaging mode, including the effects of the Lyot stop and the coronagraphic wedges. The input spectrum is a stellar model atmosphere of the appropriate spectral type, convolved with the pre-flight STIS CCD response function. The PSF modeler generates broadband PSFs by co-adding weighted monochromatic PSFs across the waveband.

Cornett, R. H.↗

NASA NTRS1998

Turbulence, Chondrules, and Planetesimals

It has been shown both numerically and experimentally that 3-D turbulence concentrates aerodynamically size-selected particles by orders of magnitude. In a previous review chapter, in "Chondrules and the protoplanetary disk" we illustrated the initial predictions of Turbulent Concentration (TC) as applied to the solar nebula. We predicted the particle size which will be most effectively concentrated by turbulence; it is the particle which has a gas drag stopping time equal to the overturn time of the smallest (Kolmogorov scale) eddy. The primary uncertainty is the level of nebula turbulence, or Reynolds number Re, which can be expressed in terms of the standard nebula eddy viscosity parameter alpha = Re(nu)(sub m)/cH, where nu(sub m) is molecular viscosity, c is sound speed, and H is vertical scale height. Several studies, and observed lifetimes of circumstellar disks, have suggested that the level of nebula turbulence can be described by alpha = 10(exp -2) - 10(exp -4). There is some recent concern about how energy is provided to maintain this turbulence, but the issue remains open. We adopt a canonical minimum mass nebula with a range of alpha > 0. We originally showed that chondrule-sized particles are selected for concentration in the terrestrial planet region if alpha = 10(exp -3) - 10(exp -4). In addition, Paque and Cuzzi found that the size distribution of chondrules is an excellent match for theoretical predictions. One then asks by what concentration factor C these particles can be concentrated; our early numerical results indicated an increase of C with alpha, and were supported by simple scaling arguments, but the extrapolation range was quite large and the predictions (C ~ 10(exp 5) - 10(exp 6) not unlikely) uncertain. The work presented here, which makes use of our recent demonstration that the particle density field is a multifractal with flow-independent properties provides a far more secure ground for such predictions. We also indicate how fine-grained dust rims on chondrules might enter into constraining the situation.

Cuzzi, Jeffrey↗

DOE OSTI2022

Titanium isotope systematics of refractory inclusions: Echoes of molecular cloud heterogeneity

Calcium-aluminum-rich inclusions (CAIs) are highly refractory objects found in different chondrite groups and represent some of the oldest known solids of the Solar System. As such, CAIs provide key information regarding the conditions prevailing in the solar protoplanetary disk as well as subsequent mixing and transport processes. Many studies have investigated CAIs for their isotopic compositions and reported nucleosynthetic isotope anomalies in numerous elements, which are typically explained by the variable incorporation of isotopically highly anomalous presolar phases. However, with the exception of 54 Cr-enriched nanospinels, the exact presolar phases responsible for the isotopic heterogeneities are yet to be identified. To address this issue, we here present in-situ Ti isotopic analyses obtained on a diverse set of CAIs from various CV3 chondrites. In this work, the in-situ measurements were performed by targeting individual mineral phases of 15 CAIs with laser-ablation mass spectrometry and indicate significant inter- and intra-CAI isotopic heterogeneity in the neutron-rich isotope 50 Ti. This is particularly pronounced for primitive fine-grained CAIs, whereas coarse-grained CAIs, which have been subject to melting, exhibit smaller degrees of Ti isotopic heterogeneity. To further investigate this Ti isotopic heterogeneity, we additionally obtained Ti isotopic compositions of sequential acid leachates from two fine-grained and two coarse-grained CAIs derived from CV3 chondrites. In contrast to potential expectations from the first part of the study, we do not observe any significant intra-CAI Ti isotopic heterogeneity between the different leaching steps. The lack of intra-CAI Ti isotopic heterogeneity in the acid leachate samples of this study likely reflects that the leaching procedure is unable to efficiently separate the carriers of isotopically anomalous Ti in CAIs. By comparing the bulk CAI Ti isotope compositions with Ti isotope data for hibonite-rich objects from the literature, we find that the range of Ti isotope compositions recorded by CAIs from various chondrite groups can be accounted for by the averaging of hibonite grains. In turn, the variable Ti isotope compositions of hibonite grains can be explained by the averaging of isotopically diverse presolar grains present in the Sun’s parental molecular cloud. This effect of averaging is statistically supported by the central limit theorem, and the concept has the potential to be useful for other isotopic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗