Engineering topics

Dobrovolskis, Anthony R.

Publications and source records attributed to Dobrovolskis, Anthony R..

At least 19 records

NASA NTRS2003

Blowing in the Wind: II. Creation and Redistribution of Refractory Inclusions in a Turbulent Protoplanetary Nebula

Ca-Al rich refractory mineral inclusions (CAIs) found at 1-6% mass fraction in primitive chondrites appear to be 1-3 million years older than the dominant (chondrule) components which were accreted into the same parent bodies. A prevalent concern is that it is difficult to retain CAIs for this long against gas-drag-induced radial drift into the sun. We reassess the situation in terms of a hot inner (turbulent) nebula context for CAI formation, using analytical models of nebula evolution and particle diffusion. We show that outward radial diffusion in a weakly turbulent nebula can overwhelm inward drift, and prevent significant numbers of CAI-size particles from being lost into the sun for times on the order of 10(exp 6) years. CAIs can form early, when the inner nebula was hot, and persist in sufficient abundance to be incorporated into primitive planetesimals at a much later time. Small (less than or approx. equal to 0.1 mm diameter) CAIs persist for longer times than large (greater than or approx. equal to 5mm diameter ones). To obtain a quantitative match t o the observed volume fractions of CAIs in chondrites, another process must be allowed for: a substantial enhancement of the inner hot nebula in silicate-forming material, which we suggest was caused by rapid inward drift of meter-sized objects. This early in nebula history, the drifting rubble would have a carbon content probably an order of magnitude larger than even the most primitive (CI) carbonaceous chondrites. Abundant carbon in the evaporating material would help keep the nebula oxygen fugacity low, plausably solar, as inferred for the formation environment of CAIs. The associated production of a larger than canonical amount of CO2 might also play a role in mass-independent fractionation of oxygen isotopes, leaving the gas rich in l60 as inferred from CAIs and other high temperature condensates.

Cuzzi, Jeffrey N.↗

NASA NTRS2003

Blowing in the Wind: II. Creation and Redistribution of Refractory Inclusions in a Turbulent Protoplanetary Nebula

Ca-A1 rich refractory mineral inclusions (CAIs) found at 1-6% mass fraction in primitive chondrites appear to be 1-3 million years older than the dominant (chondrule) components which were accreted into the same parent bodies. A prevalent concern is that it is difficult to retain CAIs for this long against gas-drag-induced radial drift into the sun. We reassess the situation in terms of a hot inner (turbulent) nebula context for CAI formation, using analytical models of nebula evolution and particle diffusion. We show that outward radial diffusion in a weakly turbulent nebula can overwhelm inward drift, and prevent significant numbers of CAI-size particles from being lost into the sun for times on the order of 10(exp 6) years. CAIs can form early, when the inner nebula was hot, and persist in sufficient abundance to be incorporated into primitive planetesimals at a much later time. Small (less than or approximately 0.1 mm diameter) CAIs persist for longer times than large (greater than or approximately 5mm diameter ones. To obtain a quantitative match to the observed volume fractions of CAIs in chondrites, another process must be allowed for: a substantial enhancement of the inner hot nebula in silicate-forming material, which we suggest was caused by rapid inward drift of meter-sized objects. This early in nebula history, the drifting rubble would have a carbon content probably an order of magnitude larger than even the most primitive (CI) carbonaceous chondrites. Abundant carbon in the evaporating material would help keep the nebula oxygen fugacity low, plausibly solar; as inferred for the formation environment of CAIs. The associated production of a larger than canonical amount of CO2 might also play a role in mass-independent fractionation of oxygen isotopes, leaving the gas rich in O-16 as inferred from CAIs and other high temperature condensates.

Cuzzi, Jeffrey N.↗

NASA NTRS2002

The Fate of Ejecta from Hyperion

Ejecta from Saturn's moon Hyperion is subject to powerful perturbations from nearby Titan, which control its ultimate fate. We have performed numerical integrations to simulate a simplified system consisting of Saturn (including oblateness), Tethys, Dione, Titan, Hyperion, Iapetus, and the Sun (treated simply as a massive satellite). In addition, 1050 massless particles were ejected from Hyperion at five different points in its orbit. These particles started more or less evenly distributed over latitude and longitude, 1 km above Hyperion's mean radius, and were ejected radially outward at speeds 10\% faster than its escape speed. Only about 4\% of the particles survived for the 100,000-year course of the integration, while $\sim$8/% escaped from the Saturnian system. Titan accreted $\sim$77\% of all the particles, while Hyperion reaccreted only $\sim$5\%. This may help to account for Hyperion's rugged shape. Three particles hit Rhea and 2 hit Dione, but $\sim$5\% of the particles were removed when they penetrated within 150,000 km of Saturn. Most removals occurred within the first few thousand years. In general, ejecta from Hyperion are much more widely scattered than previously thought, and cross the orbits of all of the other classical satellites.

Lissauer, Jack J.↗

NASA NTRS2000

Size-Selective Concentration of Chondrules and Other Small Particles in Protoplanetary Nebula Turbulence

Size-selective concentration of particles in a weakly turbulent protoplanetary nebula may be responsible for the initial collection of chondrules and other constituents into primitive bodY precursors. This paper presents the main elements of this process of turbulent concentration. In the terrestrial planet region. both the characteristic size and size distribution of chondrules are explained. "Fluffier" particles would be concentrated in nebula regions which were at a lower density and/or more intensely turbulent. The spatial distribution of concentrated particle density obeys multifractal scaling, suggesting a close tie to the turbulent cascade process. This scaling behavior allows predictions of the probability distributions for concentration in the protoplanetary nebula to be made. Large concentration factors (less than 10(exp 5)) are readily obtained, implying that numerous zones of particle density significantly exceeding the gas density could exist. If most of the available solids were actually in chondrule sized particles, the ensuing particle mass density would become so large that the feedback effects on gas turbulence due to mass loading could no longer be neglected. This paper describes the process, presenting its basic elements and some implications, without including the effects of mass loading.

Cuzzi, Jeffrey N.↗

NASA NTRS1999

Turbulent Concentration of Chondrules: Size Distribution and Multifractal Scaling

Size-selective concentration of particles in 3D turbulence may be related to collection of chondrules and other constituents into primitive bodies in a weakly turbulent protoplanetary nebula. In the terrestrial planet region, both the characteristic size and narrow size distribution of chondrules are explained, whereas "fluffier" particles would be concentrated in lower density, or more intensely turbulent, regions of the nebula. The spatial distribution of concentrated particle density obeys multifractal scaling, suggesting a dose tie to the turbulent cascade process. This scaling behavior allows predictions of the concentration probabilities to be made in the protoplanetary nebula, which are so large (> 10(exp 3) - 10(exp 4)) that further studies must be made of the role of mass loading.

Cuzzi, Jeffrey N.↗

NASA NTRS1998

Scaling Properties of Particle Density Fields Formed in Simulated Turbulent Flows

Direct numerical simulations (DNS) of particle concentrations in fully developed 3D turbulence were carried out in order to study the nonuniform structure of the particle density field. Three steady-state turbulent fluid fields with Taylor microscale Reynolds numbers (Re(sub lambda)) of 40, 80 and 140 were generated by solving the Navier-Stokes equations with pseudospectral methods. Large scale forcing was used to drive the turbulence and maintain temporal stationarity. The response of the particles to the fluid was parameterized by the particle Stokes number St, defined as the ratio of the particle's stopping time to the mean period of eddies on the Kolmogorov scale (eta). In this paper, we consider only passive particles optimally coupled to these eddies (St approx. = 1) because of their tendency to concentrate more than particles with lesser or greater St values. The trajectories of up to 70 million particles were tracked in the equilibrated turbulent flows until the particle concentration field reached a statistically stationary state. The nonuniform structure of the concentration fields was characterized by the multifractal singularity spectrum, f(alpha), derived from measures obtained after binning particles into cells ranging from 2(eta) to 15(eta) in size. We observed strong systematic variations of f(alpha) across this scale range in all three simulations and conclude that the particle concentration field is not statistically self similar across the scale range explored. However, spectra obtained at the 2(eta), 4(eta), and 8(eta) scales of each flow case were found to be qualitatively similar. This result suggests that the local structure of the particle concentration field may be flow-Independent. The singularity spectra found for 2n-sized cells were used to predict concentration distributions in good agreement with those obtained directly from the particle data. This Singularity spectrum has a shape similar to the analogous spectrum derived for the inertial-range energy dissipation fields of experimental turbulent flows at Re(sub lambda) = 110 and 1100. Based on this agreement, and the expectation that both dissipation and particle concentration are controlled by the same cascade process, we hypothesize that singularity spectra similar to the ones found in this work provide a good characterization of the spatially averaged statistical properties of preferentially concentrated particles in higher Re(sub lambda) turbulent flows.

Hogan, Robert C.↗

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↗

NASA NTRS1996

From Chondrules to Planetesimals: Solids and Turbulence in the Solar Nebula

The usual theory of planetesimal formation is untenable because turbulence inhibits gravitational instability. However, turbulence can actually concentrate chondrule-sized particles by factors up to a million near stagnation points. The implications for accretion may be profound.

Dobrovolskis, Anthony R.↗

NASA NTRS1995

A Numerical Turbulence Model for Multiphase Flows in the Protoplanetary Nebula

It is thought that planets form from solid particles in a flattened, rotating, 99% gaseous nebula. These grains gradually coagulate into millimeter-to-meter sized aggregates which settle toward the midplane of the nebula. It is widely believed that the resulting dense layer eventually becomes gravitationally unstable and collapses into 'planetesimals.' A new numerical model is presented to simulate the predominant processes (gravitation, vertical convection, and shear-driven turbulence) during the stage while the particulate material is still dispersed about the midplane of the nebula. In our previous work, particles were assumed to be spheres of a single radius; in the present work, particles are spheres of different radii. Results indicate that neither a broad nor a narrow distribution of particle sizes is likely to become gravitationally unstable.

Champney, Joelle M.↗

NASA NTRS1995

Chaotic Rotation of Nereid

The shape and spin of Neptune's outermost satellite Nereid are still unknown. Ground-based photometry indicates large brightness variations, but different observers report very different lightcurve amplitudes and periods. On the contrary, Voyager 2 images spanning 12 days show no evidence of variations greater than 0.1 mag. The latter suggest either that Nereid is nearly spherical, or that it is rotating slowly. We propose that tides have already despun Nereid's rotation to a period of a few weeks, during the time before the capture of Triton when Nereid was closer to Neptune. Since Nereid reached its present orbit, tides have further despun Nereid to a period on the order of a month. For Nereid's orbital eccentricity of 0.75, tidal evolution ceases when the spin period is still approximately 1/8 of the orbital period. Furthermore, the synchronous resonance becomes quite weak for such high eccentricities, along with other low-order spin orbit commensurabilities. In contrast, high-order resonances become very strong particularly the 6:1, 6.5:1, 7:1, 7.5:1, and 8:1 spin states. If Nereid departs by more than approximately 1% from a sphere, however, these resonances overlap, generating chaos. Our simulations show that Nereid is likely to be in chaotic rotation for any spin period longer than about 2 weeks.

Dobrovolskis, Anthony R.↗

NASA NTRS1994

Turbulent Diffusion and Concentration of "Chondrules" in the Protoplanetary Nebula

Turbulence is known to possess structure on many scales. The largest or integral scale L is usually taken to be the largest relevant dimension of the system; for the nebula as a whole, it would be on the order of a vertical gas scale height H approximately 0.1R where R is the distance from the sun. Turbulent kinetic energy cascades towards smaller scales; for fully developed, homogeneous turbulence, the turbulent kinetic energy peaks at the integral scale and is distributed within its "inertial range" according to the Kolmogorov energy spectrum dE(k) = E(sub 0)(kL)(sup -5/3)dk where k is the wavenumber or inverse wavelength. The smallest scale in a turbulent regime is the Kolmogorov scale (n)(sometimes referred to as the inner scale). On this scale, molecular viscosity can dissipate turbulent kinetic energy . Turbulent flows are characterized by their Reynolds number Re = UL/v, where U is a characteristic velocity scale and v is molecular viscosity. More energetic (higher Re) flows drive turbulence to smaller Kolmogorov scales: n = LRe(sup -3/4). Two different turbulent regimes are likely to be of importance; early stage convective turbulence in a hot nebula probably extends over the entire nebula scale height, and probably has typical eddy frequency comparable to the orbit frequency. At a later stage, midplane shear layer turbulence can be driven by the presence of a differentially rotating, densely settled particle layer, with different length and time scales.

Cuzzi, Jeffrey N.↗

NASA NTRS1993

The Laplace Planes of Uranus and Pluto

Satellite orbits close to an oblate planet precess about its equatorial plane, while distant satellites precess around the plane of the planet's heliocentric orbit. In between, satellites in nearly circular orbits precess about a warped intermediate surface called the Laplace 'plane.' Herein we derive general formulas for locating the Laplace plane. Because Uranus and Pluto have high obliquities, their Laplace planes are severely warped. We present maps of these Laplace planes, of interest in telescopic searches for new satellites. The Laplace plane of the Solar System as a whole is similarly distorted, but comets in the inner Oort cloud precess too slowly to sense the Laplace plane.

Dobrovolskis, Anthony R.↗

NASA NTRS1993

Atmospheric tides on Venus. IV - Topographic winds and sediment transport

A novel theory is presented for the Venus boundary layer which encompasses the effects of topography and uses the mixing-length hypothesis to preclude the unknown eddy viscosity. The maps of mass-flux and erosion/deposition rate presented are based on Pioneer Venus orbiter relief measurements. The typically 19 cm/sec friction speeds associated with the present theory are several times greater than those estimated on the basis of Venera 9 and 10 anemometry, and mean aeolian transport is generally away from the equator, contrary to Magellan orbiter windstreak directions.

Dobrovolskis, Anthony R.↗

NASA NTRS1993

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗

NASA NTRS1991

Relationships between topographic roughness and aeolian processes

The interaction between winds and desert surfaces has important implications for sediment transport on Earth, Mars, and Venus, and for understanding the relationships between radar backscatter and aerodynamic roughness as part of the NASA Shuttle Imaging radar (SIR-C) Mission. Here, researchers report results from measurements of boundary layer wind profiles and surface roughness at sites in Death Valley and discuss their implications. The sites included a flat to undulating gravel and sand reg, alluvial fans, and a playa. Estimates of average particle size composition of Death Valley sites and arithmetic mean values of aerodynamic roughness are given in tabular form.

Greeley, Ronald↗

NASA NTRS1991

Particle-gas dynamics in the protoplanetary nebula

In the past year we made significant progress in improving our fundamental understanding of the physics of particle-gas dynamics in the protoplanetary nebula. Having brought our code to a state of fairly robust functionality, we devoted significant effort to optimizing it for running long cases. We optimized the code for vectorization to the extent that it now runs eight times faster than before. The following subject areas are covered: physical improvements to the model; numerical results; Reynolds averaging of fluid equations; and modeling of turbulence and viscosity.

Cuzzi, Jeffrey N.↗

NASA NTRS1991

Planetary ring studies

The following topics are covered: (1) characterization of the fine scale structure in Saturn's A and B rings; (2) ballistic transport modeling and evolution of fine ring structure; (3) faint features in the rings of Saturn; (4) the Encke moonlet; (5) dynamics in ringmoon systems; (6) a nonclassical radiative transfer model; and (7) particle properties from stellar occultation data.

Cuzzi, Jeffrey N.↗

NASA NTRS1991

Venus - Concentrations of radar-reflective minerals by wind

The effectiveness of wind in concentrating minerals with high radar reflectiveness on the surface of Venus is investigated experimentally in the Venus Wind Tunnel (Greeley et al., 1984) under CO2 densities typical of Venusian conditions. Density sorting of sand particles during the formation of microdunes is demonstrated, and calculations show that wind-blown deposits of dense conductive material such as ilmenite need to be only a few cm thick to account for the local enhancements of radar reflectivity observed by Pioneer Venus at wavelength 17 cm.

Greeley, Ronald↗