Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical hydrodynamics”

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.

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

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Diverse properties of interstellar medium embedding gamma-ray bursts at the epoch of reionization

Analysis is performed on ultra-high-resolution large-scale cosmological radiation-hydrodynamic simulations to quantify, for the first time, the physical environment of long-duration gamma-ray bursts (GRBs) at the epoch of reionization. We find that, on parsec scales, 13% of GRBs remain in high-density (≥10{sup 4} cm{sup –3}) low-temperature star-forming regions, whereas 87% of GRBs occur in low-density (∼10{sup –2.5} cm{sup –3}) high-temperature regions heated by supernovae. More importantly, the spectral properties of GRB afterglows, such as the neutral hydrogen column density, total hydrogen column density, dust column density, gas temperature, and metallicity of intervening absorbers, vary strongly from sight line to sight line. Although our model explains extant limited observationally inferred values with respect to circumburst density, metallicity, column density, and dust properties, a substantially larger sample of high-z GRB afterglows would be required to facilitate a statistically solid test of the model. Our findings indicate that any attempt to infer the physical properties (such as metallicity) of the interstellar medium (ISM) of the host galaxy based on a very small number (usually one) of sight lines would be precarious. Utilizing high-z GRBs to probe the ISM and intergalactic medium should be undertaken properly, taking into consideration the physical diversities of the ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Using statistical modeling to predict and understand fusion experiments

Over 300 cryogenic layered direct-drive inertial confinement fusion implosions have been successfully executed on the OMEGA Laser System in the last decade. However, extracting sufficient understanding from these experiments to develop new designs or to identify or mitigate degradation sources remains challenging. Recently, a statistical modeling approach was developed to successfully design and predict improved implosion experiments on OMEGA. Here, we show that one-dimensional simulations can be used to predict the outcomes of systematically perturbed three-dimensional simulations and that this statistical modeling approach can be used to identify or rule out physical mechanisms for some of the degradation sources observed on the OMEGA Laser System for direct-drive cryogenic inertial confinement fusion. In this instance, we investigate the fusion yield dependencies on the ion temperature asymmetries and laser beam size observed in experiments by comparing with trends in three-dimensional synthetic simulation databases. Using the statistical model on these systematically perturbed simulations, we find that the statistically inferred dependency on the measured ion temperature asymmetries is well explained by imposed ℓ =1 modes. However, we find that the dependency on the laser beam size is only dominated by the illumination non-uniformity for some extreme cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Galaxy assembly bias and large-scale distribution: a comparison between IllustrisTNG and a semi-analytic model

In this work, we compare large scale structure observables for stellar mass selected samples at z = 0, as predicted by two galaxy models, the hydrodynamical simulation IllustrisTNG and the Santa-Cruz semi-analytic model (SC-SAM). Although both models have been independently calibrated to match observations, rather than each other, we find good agreement between the two models for two-point clustering and galaxy assembly bias signatures. The models also show a qualitatively similar response of occupancy and clustering to secondary halo parameters other than mass, such as formation history and concentration, although with some quantitative differences. Thus, our results demonstrate that the galaxy–halo relationships in SC-SAM and TNG are quite similar to first order. However, we also find areas in which the models differ. For example, we note a strong correlation between halo gas content and environment in TNG, which is lacking in the SC-SAM, as well as differences in the occupancy predictions for low-mass haloes. Moreover, we show that higher order statistics, such as cumulants of the density field, help us to accurately describe the galaxy distribution and discriminate between models that show degenerate behaviour for two-point statistics. Our results suggest that SAMs are a promising cost-effective and intuitive method for generating mock catalogues for next generation cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

The Soft X-Ray Emission in a Large Sample of Galaxy Clusters with ROSAT PSPC

The study of soft X-ray emission of 38 X-ray selected galaxy clusters observed by ROSAT PSPC indicates that the soft excess phenomenon may be a common occurrence in galaxy clusters. Excess soft X-ray radiation, above the contribution from the hot intra-cluster medium, is evident in a large fraction of sources, and is clearly detected with large statistical significance in the deepest observations. The investigation relies on new, high resolution 21 cm HI observations. The sample selection also features analysis of infrared images, to further ensure reliability of results with respect to the characteristics of Galactic absorption. The possibility of background or calibration effects as cause of the excess emission is likewise investigated; a detailed analysis of the distribution of the excess emission with respect to detector position and Galactic HI column density shows that the excess emission is a genuine celestial phenomenon. We find evidence for a preferential distribution of the soft excess emission at distances larger than approx. 150-200 kpc from the centers of clusters; this behavior may be naturally explained in the context of a non-thermal Inverse-Compton scenario. Alternatively, we propose that the phenomenon maybe caused by thermal emission of very large-scale 'warm' filaments seen in recent hydrodynamic simulations. This new interpretation relieves the very demanding requirements of either the traditional intra-cluster 'warm' gas and the non-thermal scenarios. We also investigate the possibility of the soft excess originating from unresolved, X-ray faint cluster galaxies.

Bonamente, Massimiliano↗
DOE OSTI2013

CONDITIONS FOR SUCCESSFUL HELIUM DETONATIONS IN ASTROPHYSICAL ENVIRONMENTS

Several models for Type Ia-like supernova events rely on the production of a self-sustained detonation powered by nuclear reactions. In the absence of hydrogen, the fuel that powers these detonations typically consists of either pure helium (He) or a mixture of carbon and oxygen (C/O). Studies that systematically determine the conditions required to initiate detonations in C/O material exist, but until now no analogous investigation of He matter has been conducted. We perform one-dimensional reactive hydrodynamical simulations at a variety of initial density and temperature combinations and find critical length scales for the initiation of He detonations that range between 1 and 10{sup 10} cm. A simple estimate of the length scales over which the total consumption of fuel will occur for steady-state detonations is provided by the Chapman-Jouguet (CJ) formalism. Our initiation lengths are consistently smaller than the corresponding CJ length scales by a factor of {approx}100, providing opportunities for thermonuclear explosions in a wider range of low-mass white dwarfs (WDs) than previously thought possible. We find that virialized WDs with as little mass as 0.24 M{sub Sun} can be detonated, and that even less massive WDs can be detonated if a sizable fraction of their mass is raised to a higher adiabat. That the initiation length is exceeded by the CJ length implies that certain systems may not reach nuclear statistical equilibrium within the time it takes a detonation to traverse the object. In support of this hypothesis, we demonstrate that incomplete burning will occur in the majority of He WD detonations and that {sup 40}Ca, {sup 44}Ti, or {sup 48}Cr, rather than {sup 56}Ni, is the predominant burning product for many of these events. We anticipate that a measure of the quantity of the intermediate-mass elements and {sup 56}Ni produced in a helium-rich thermonuclear explosion can potentially be used to constrain the nature of the progenitor system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Microscopic Theory of Fluctuating Hydrodynamics in Nonlinear Lattices

The theory of fluctuating hydrodynamics has been an important tool for analyzing macroscopic behavior in nonlinear lattices. However, despite its practical success, its microscopic derivation is still incomplete. In this work, we provide the microscopic derivation of fluctuating hydrodynamics, using the coarse-graining and projection technique; the equivalence of ensembles turns out to be critical. The Green-Kubo (GK)-like formula for the bare transport coefficients are presented in a numerically computable form. Our numerical simulations show that the bare transport coefficients exist for a sufficiently large but finite coarse-graining length in the infinite lattice within the framework of the GK-like formula. Finally, this demonstrates that the bare transport coefficients uniquely exist for each physical system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2026

Subexponential Decay of Local Correlations from Diffusion-Limited Dephasing

Chaotic quantum systems at finite entropy density are expected to act as their own heat baths, rapidly dephasing local quantum superpositions. Here, we argue that in fact this dephasing is generically subexponential in one-dimensional systems with conservation laws: all local correlation functions decay as exp⁡[−𝒪⁡(𝑡 𝛼 )] with 0 ≤ 𝛼 ≤ 2/3, even when the operators are orthogonal to all hydrodynamic modes. The mechanism is diffusion-limited dephasing, in which rare low-entropy regions (“voids”) protect quantum coherences. This intrinsically quantum effect lies beyond standard hydrodynamics and disappears under extrinsic dephasing. In random charge-conserving circuits we find 𝛼 = 1/2, while in generic translation-invariant Floquet systems we bound 𝛼 ≤ 2/3. Our arguments are general, subject principally to the assumption that thermal fluctuations can create regions of zero entropy density. In systems with energy conservation, this assumption is automatically satisfied because of the third law of thermodynamics.

information scrambling↗
DOE OSTI2026

Strangeness enhancement at its extremes: multiple (multi-)strange hadron production in pp collisions at \(\sqrt{s}=5.02\) TeV

The probability to observe a specific number of strange and multi-strange hadrons (nS), denoted as P(nS), is measured by ALICE at midrapidity (|y| < 0.5) in $$\sqrt{s}=5.02$$ TeV proton-proton (pp) collisions, dividing events into several multiplicity-density classes. Exploiting, for the first time, a technique based on counting the number of strange-particle candidates event-by-event, this measurement allows one to extend the study of strangeness production beyond the mean of the distribution. This constitutes a new test bench for production mechanisms, probing events with a large imbalance between strange and non-strange content. The analysis of a large-statistics data sample makes it possible to extract P(nS) up to a maximum nS of 7 for $${\text{K}}_{\text{S}}^{0}$$, 5 for Λ and $$\overline{\Lambda }$$, 4 for Ξ− and $${\overline{\Xi } }^{+}$$, and 2 for Ω− and $${\overline{\Omega } }^{+}$$. From this, the probability of producing strange hadron multiplets per event is calculated, thereby enabling the extension of the study of strangeness enhancement to extreme situations where several strange quarks hadronize in a single event at midrapidity. Moreover, comparing hadron combinations with different u and d quark compositions and equal overall s quark content, the contribution to the enhancement pattern coming from non-strangeness related mechanisms is isolated. The results are compared with state-of-the-art phenomenological models implemented in commonly used Monte Carlo event generators, including PYTHIA 8 Monash 2013, PYTHIA 8 with QCD-based Color Reconnection and Rope Hadronization (QCD-CR + Ropes), and EPOS LHC, which incorporates both partonic interactions and hydrodynamic evolution. These comparisons show that the new approach dramatically enhances the sensitivity to the different underlying physics mechanisms modeled by each generator.

Abualrob, I J↗
DOE OSTI2024

The 3D Lyman- α forest power spectrum from eBOSS DR16

We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-α (Ly α) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of ~205 000 quasar spectra in the redshift range 2 ≤ z ≤ 4 at an effective redshift z = 2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp( i k ∙ r), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artefacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line of sight of Ly α skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly α skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range $(0.02 ≤ k ≤ 0.35)$ h Mpc -1 at the 1–2σ level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly α forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Merger Response of Halo Anisotropy Properties

Anisotropy properties—halo spin, shape, position offset, velocity offset, and orientation—are an important family of dark matter halo properties that indicate the level of directional variation of the internal structures of halos. These properties reflect the dynamical state of halos, which in turn depends on the mass assembly history. In this work, we study the evolution of anisotropy properties in response to merger activity using the IllustrisTNG simulations. We find that the response trajectories of the anisotropy properties significantly deviate from secular evolution. These trajectories have the same qualitative features and timescales across a wide range of merger and host properties. We propose explanations for the behavior of these properties and connect their evolution to the relevant stages of merger dynamics. We measure the relevant dynamical timescales. We also explore the dependence of the strength of the response on time of merger, merger ratio, and mass of the main halo. These results provide insight into the physics of halo mergers and their effects on the statistical behavior of halo properties. This study paves the way toward a physical understanding of scaling relations, particularly to how systematics in their scatter are connected to the mass assembly histories of halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Offsets between the X-ray and the Sunyaev-Zel'Dovich-effect peaks in merging galaxy clusters and their cosmological implications

Observations reveal that the peaks of the X-ray map and the Sunyaev-Zel'dovich (SZ) effect map of some galaxy clusters are offset from each other. In this paper, we perform a set of hydrodynamical simulations of mergers of two galaxy clusters to investigate the spatial offset between the maxima of the X-ray and the SZ surface brightness of the merging clusters. We find that significantly large SZ-X-ray offsets (>100 kpc) can be produced during the major mergers of galaxy clusters (with mass > 1 × 10{sup 14} M {sub ☉}). The significantly large offsets are mainly caused by a 'jump effect' that occurs between the primary and secondary pericentric passages of the two merging clusters, during which the X-ray peak may jump to the densest gas region located near the center of the small cluster, but the SZ peak remains near the center of the large one. Our simulations show that merging systems with higher masses and larger initial relative velocities may result in larger offset sizes and longer offset time durations; and only nearly head-on mergers are likely to produce significantly large offsets. We further investigate the statistical distribution of the SZ-X-ray offset sizes and find that (1) the number distribution of the offset sizes is bimodal with one peak located at low offsets ∼0 and the other at large offsets ∼350-450 h {sup –1} kpc, but the objects with intermediate offsets are scarce; and (2) the probabilities of the clusters in the mass range higher than 2 × 10{sup 14} h {sup –1} M {sub ☉} that have offsets larger than 20, 50, 200, 300, and 500 h {sup –1} kpc are 34.0%, 11.1%, 8.0%, 6.5%, and 2.0%, respectively, at z = 0.7. The probability is sensitive to the underlying pairwise velocity distribution and the merger rate of clusters. We suggest that the SZ-X-ray offsets provide a probe to the cosmic velocity fields on the cluster scale and the cluster merger rate, and future observations on the SZ-X-ray offsets for a large number of clusters may put strong constraints on them. Our simulation results suggest that the SZ-X-ray offset in the Bullet Cluster, together with the mass ratio of the two merging clusters, requires a relative velocity larger than 3000 km s{sup –1} at an initial separation 5 Mpc. The cosmic velocity distribution at the high-velocity end is expected to be crucial in determining whether there exists an incompatibility between the existence of the Bullet Cluster and the prediction of a ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Galaxy cluster profiles: a Gaussian mixture model approach to halo miscentering

Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center defined as the minimum in the potential well. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically modeled in stacked profiles of clusters with a two parameter model. We use an alternate approach in which the profiles of individual clusters are used with the corresponding likelihood computed using a Gaussian mixture model. We test the approach using halos and the corresponding subhalo profiles from the IllustrisTNG hydrodynamic simulations. We obtain significantly improved estimates of the miscentering parameters for both 3D and projected 2D profiles relevant for imaging surveys. We discuss applications to upcoming cosmological surveys. Our Python package for the Gaussian mixture model is publicly available at https://github.com/KyleMiller1/Halo-Miscentering-Mixture-Model.

Bayesian reasoning↗
DOE OSTI2020

Impact of spatiotemporal smoothing on the two-plasmon–decay instability

Higher levels of hot electrons from the two-plasmon–decay instability are observed when smoothing by spectral dispersion (SSD) is turned off in directly driven inertial confinement fusion experiments at the Omega Laser Facility. This finding is explained using a hot-spot model based on speckle statistics and simulation results from the laser–plasma simulation environment. The model accurately reproduces the relative increase in hot-electron activity at two different drive intensities, although it slightly overestimates the absolute number of hot electrons in all cases. Extrapolating from the current ≈ 360-GHz system while adhering to the logic of the hot-spot model suggests that larger SSD bandwidth should significantly mitigate hot-electron generation, and legacy 1-THz OMEGA experiments appear to support this conclusion. These results demonstrate that it is essential to account for laser speckles and spatiotemporal smoothing to obtain quantitative agreement with experiments. A compilation of hot-electron data from the past two decades reveals several other important points: (1) many prior experiments are more easily understood using recent results from multibeam absolute instability theory; and (2) experiments with ignition-scale conditions produce less hot electrons compared to OMEGA spherical experiments for a given vacuum overlapped intensity, which is a promising result for validating performance predictions based on hydrodynamic scaling relations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

X-RAY PROPERTIES EXPECTED FROM ACTIVE GALACTIC NUCLEUS FEEDBACK IN ELLIPTICAL GALAXIES

Detailed hydrodynamic simulations of active galactic nucleus feedback have been performed including the effects of radiative and mechanical momentum and energy input on the interstellar medium (ISM) of typical elliptical galaxies. We focus on the observational properties of the models in the soft and hard X-ray bands: nuclear X-ray luminosity; global X-ray luminosity and temperature of the hot ISM; and temperature and X-ray brightness profiles before, during, and after outbursts. After {approx}10 Gyr, the bolometric nuclear emission L{sub BH} is very sub-Eddington (l = L{sub BH}/L{sub Edd} {approx} 10{sup -4}), and within the range observed, though larger than typical values. Outbursts last for Almost-Equal-To 10{sup 7} yr, and the duty cycle of nuclear activity is a few Multiplication-Sign (10{sup -3} to 10{sup -2}), over the last 6 Gyr. The ISM thermal luminosity L{sub X} oscillates in phase with the nuclear luminosity, with broader peaks. This behavior helps statistically reproduce the observed large L{sub X} variation. The average gas temperature is within the observed range, in the upper half of those observed. In quiescence, the temperature profile has a negative gradient; thanks to past outbursts, the brightness profile lacks the steep shape of cooling flow models. After outbursts, disturbances are predicted in the temperature and brightness profiles (analyzed by unsharp masking). Most significantly, during major accretion episodes, a hot bubble of shocked gas is inflated at the galaxy center (within Almost-Equal-To 100 pc); the bubble would be conical in shape in real galaxies and would be radio-loud. Its detection in X-rays is within current capabilities, though it would likely remain unresolved. The ISM resumes its smooth appearance on a timescale of Almost-Equal-To 200 Myr; the duty cycle of perturbations in the ISM is of the order of 5%-10%. While showing general agreement between the models and real galaxies, this analysis indicates that additional physical input may still be required including moving to two-dimensional or three-dimensional simulations, input of relativistic jets, or allowance for a confining medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

thornado-transport: IMEX schemes for two-moment neutrino transport respecting Fermi-Dirac statistics

We develop implicit-explicit (IMEX) schemes for neutrino transport in a background material in the context of a two-moment model that evolves the angular moments of a neutrino phase-space distribution function. Considering the upper and lower bounds that are introduced by Pauli's exclusion principle on the moments, an algebraic moment closure based on Fermi-Dirac statistics and a convex-invariant time integrator both are demanded. A finite-volume/first-order discontinuous Galerkin(DG) method is used to illustrate how an algebraic moment closure based on Fermi-Dirac statistics is needed to satisfy the bounds. Several algebraic closures are compared with these bounds in mind, and the Cernohorsky and Bludman closure, which satisfies the bounds, is chosen for our IMEX schemes. For the convex-invariant time integrator, two IMEX schemes named PD-ARS have been proposed. PD-ARS denotes a convex-invariant IMEX Runge-Kutta scheme that is high-order accurate in the streaming limit, and works well in the diffusion limit. Our two PD-ARS schemes use second-and third-order, explicit, strong-stability-preserving Runge-Kutta methods as their explicit part, respectively, and therefore are second-and third-order accurate in the streaming limit, respectively. The accuracy and convex-invariance of our PD-ARS schemes are demonstrated in the numerical tests with a third-order DG method for spatial discretization and a simple Lax-Friedrichs flux. The method has been implemented in our high-order neutrino-radiation hydrodynamics (thornado) toolkit. We show preliminary results employing tabulated neutrino opacities.

Chu, Ran↗
NASA NTRS2015

Turbulent Concentration of mm-Size Particles in the Protoplanetary Nebula: Scale-Dependent Cascades

The initial accretion of primitive bodies (here, asteroids in particular) from freely-floating nebula particles remains problematic. Traditional growth-by-sticking models encounter a formidable "meter-size barrier" (or even a mm-to-cm-size barrier) in turbulent nebulae, making the preconditions for so-called "streaming instabilities" difficult to achieve even for so-called "lucky" particles. Even if growth by sticking could somehow breach the meter size barrier, turbulent nebulae present further obstacles through the 1-10km size range. On the other hand, nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids. Theoretical understanding of nebula turbulence is itself in flux; recent models of MRI (magnetically-driven) turbulence favor low-or- no-turbulence environments, but purely hydrodynamic turbulence is making a comeback, with two recently discovered mechanisms generating robust turbulence which do not rely on magnetic fields at all. An important clue regarding planetesimal formation is an apparent 100km diameter peak in the pre-depletion, pre-erosion mass distribution of asteroids; scenarios leading directly from independent nebula particulates to large objects of this size, which avoid the problematic m-km size range, could be called "leapfrog" scenarios. The leapfrog scenario we have studied in detail involves formation of dense clumps of aerodynamically selected, typically mm-size particles in turbulence, which can under certain conditions shrink inexorably on 100-1000 orbit timescales and form 10-100km diameter sandpile planetesimals. There is evidence that at least the ordinary chondrite parent bodies were initially composed entirely of a homogeneous mix of such particles. Thus, while they are arcane, turbulent concentration models acting directly on chondrule size particles are worthy of deeper study. The typical sizes of planetesimals and the rate of their formation can be estimated using a statistical model with properties inferred from large numerical simulations of turbulence. Nebula turbulence is described by its Reynolds number Re = (L/eta)(exp 4/3), where L = H alpha(exp 1/2) is the largest eddy scale, H is the nebula gas vertical scale height, alpha the turbulent viscosity parameter, and eta is the Kolmogorov or smallest scale in turbulence (typically about 1km), with eddy turnover time t(sub eta). In the nebula, Re is far larger than any numerical simulation can handle, so some physical arguments are needed to extend the results of numerical simulations to nebula conditions. In this paper, we report new physics to be incorporated into our statistical models.

mm-size particles↗