Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic limit”

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 19 records

Mass fractionation of noble gases in diffusion-limited hydrodynamic hydrogen escape

The theory of mass fractionation by hydrogen is presently extended to atmospheres in which hydrogen is not the major constituent. This theoretical framework is applied to three different cases. In the first, it is shown that the fractionation of terrestrial atmospheric neon with respect to mantle neon is explainable as a consequence of diffusion-limited hydrogen escape from a steam atmosphere toward the end of the accretion process. In the second, the anomalously high Ar-38/Ar-36 ratio of Mars is shown to be due to hydrodynamic fractionation by a vigorously escaping and very pure hydrogen wind. In the last case, it is speculated that the currently high Martian D/H ratio emerged during the hydrodynamic escape phase which fractionated Ar.

Zahnle, Kevin↗

Early Hydrodynamic Escape Limits Rocky Planets to Less Than or Equal to 1.6 Earth Radii

In the past decade thousands of exoplanet candidates and hundreds of confirmed exoplanets have been found. For sub-Neptune-sized planets, those less than approx. 10 Earth masses, we can separate planets into two broad categories: predominantly rocky planets, and gaseous planets with thick volatile sheaths. Observations and subsequent analysis of these planets show that rocky planets are only found with radii less than approx. 1.6 Earth radii. No rocky planet has yet been found that violates this limit. We propose that hydrodynamic escape of hydrogen rich protoatmospheres, accreted by forming planets, explains the limit in rocky planet size. Following the hydrodynamic escape model employed by Luger et al. (2015), we modelled the XUV driven escape from young planets (less than approx.100 Myr in age) around a Sun-like star. With a simple, first-order model we found that the rocky planet radii limit occurs consistently at approx. 1.6 Earth radii across a wide range of plausible parameter spaces. Our model shows that hydrodynamic escape can explain the observed cutoff between rocky and gaseous planets. Fig. 1 shows the results of our model for rocky planets between 0.5 and 10 Earth masses that accrete 3 wt. % H2/He during formation. The simulation was run for 100 Myr, after that time the XUV flux drops off exponentially and hydrodynamic escape drops with it. A cutoff between rocky planets and gaseous ones is clearly seen at approx. 1.5-1.6 Earth radii. We are only interested in the upper size limit for rocky planets. As such, we assumed pure hydrogen atmospheres and the highest possible isothermal atmospheric temperatures, which will produce an upper limit on the hydrodynamic loss rate. Previous work shows that a reasonable approximation for an upper temperature limit in a hydrogen rich protoatmosphere is 2000-3000 K, consistent with our assumptions. From these results, we propose that the observed dichotomy between mini-Neptunes and rocky worlds is simply explained by an early episode of thermally-driven hydrodynamic escape when host stars have saturated XUV fluxes.

Lehmer, O. R.↗

Shock compression of cemented tungsten carbides to 100 GPa: Structure of shock waves, Hugoniot relations, and strength

Plate impact experiments are conducted on cemented tungsten carbides (WC) with a 3.7 and 6.0 wt. % cobalt binder to better understand their dynamic, high-pressure response to 100 GPa. The measured wave profiles show propagation of steady structured waves. Standard impedance matching procedures are used to determine the Hugoniot relations in the shock velocity–particle velocity (U s –v p ) and Hugoniot stress–specific volume (P–V/V o ) planes. The Hugoniot elastic limit of the samples is controlled by ductility of the Co binder and is determined to be 4.45 ± 0.29 GPa for cemented WC with 3.7 wt. % cobalt and 3.72 ± 0.24 GPa for a 6.0 wt. % cobalt binder. Both grades show a non-linear U s –v p relationship depending on whether the particle velocity is in the strength dominated or hydrodynamic regime. In the strength dominated regime, a non-linear decrease in U s is observed as v p increases from ambient to the material’s hydrodynamic limit. In the hydrodynamic regime, the U s –v p Hugoniot is linear and is determined to be U s = 4.97(±0.006)+1.446(±0.018)v p km/s for WC with 3.7 wt. % Co and U s = 4.93(±0.006)+1.463(±0.017)v p km/s for 6 wt. % Co. Both WC grades indicate shear-stress hardening with mean stress immediately after yield, followed by pressure softening, and then a sharp fall in stress carrying capacity as the mean stress is increased to ≈70 GPa (hydrodynamic limit) and beyond. This behavior is in contrast to pure WC ceramics, which show continued shear-stress hardening with mean stress to ≈80 GPa.

36 MATERIALS SCIENCE↗

Multiscale Control of Generic Second Order Traffic Models by Driver-Assist Vehicles

We study the derivation of generic high order macroscopic traffic models from a follow-the-leader particle description via a kinetic approach. First, we recover a third order traffic model as the hydrodynamic limit of an Enskog-type kinetic equation. Next, we introduce in the vehicle interactions a binary control modeling the automatic feedback provided by driver-assist vehicles and we upscale such a new particle description by means of another Enskog-based hydrodynamic limit. The resulting macroscopic model is now a generic second order model (GSOM), which contains in turn a control term inherited from the microscopic interactions. We show that such a control may be chosen so as to optimize global traffic trends, such as the vehicle flux or the road congestion, constrained by the GSOM dynamics. By means of numerical simulations, we investigate the effect of this control hierarchy in some specific case studies, which exemplify the multiscale path from the vehiclewise implementation of a driver-assist control to its optimal hydrodynamic design.

GSOM↗

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗

Verification study of xRAGE’s multi-ion viscosity model

This report on the adequacy of the plasma viscosity model implemented in LANL’s ASC program code, xRAGE is organized as follows. Section I provides a brief introduction on the kinetic theory which underpins plasma viscosity, and Section II describes how plasma kinetic theory can be used to calculate plasma viscosity coefficients accurate within the hydrodynamic limit. Section III describes the existing model for calculating the plasma viscosity coefficient in xRAGE. Next, Section IV shows tests of the xRAGE viscosity model against established hydrodynamic theory in the case of binary and trinary plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flux-limited diffusion in hydrodynamics

The radiation transport equation, correct to order v/c, is used to derive a flux limiter valid in high-velocity hydrodynamic flows. In the limit v = 0 the widely used result of Levermore and Pomraning (LP) (1982) is recovered. In the case of homologous flow (v varies as R) the LP result remains valid, but in the presence of shock waves and other nonhomologous flows differences may be large. We calculate the Eddington factors, which are significantly altered, and discuss the importance of having them possess the correct form in the presence of a shock wave, e.g., on the neutrino transport during the formation of a supernova.

Cooperstein, J.↗

Non-conformal attractor in boost-invariant plasmas

We study the dissipative evolution of (0+1)-dimensionally expanding media with Bjorken symmetry using the Boltzmann equation for massive particles in relaxation-time approximation. Breaking conformal symmetry by a mass induces a non-zero bulk viscous pressure in the medium. It is shown that even a small mass (in units of the local temperature) drastically modifies the well-known attractor for the shear Reynolds number previously observed in massless systems. For generic nonzero particle mass, neither the shear nor the bulk viscous pressure relax quickly to a non-equilibrium attractor; they approach the hydrodynamic limit only late, at small values of the inverse Reynolds numbers. Only the longitudinal pressure, which is a combination of thermal, shear and bulk viscous pressures, continues to show early approach to a far-off-equilibrium attractor, driven by the rapid longitudinal expansion at early times. Second-order dissipative hydrodynamics based on a gradient expansion around locally isotropic thermal equilibrium fails to reproduce this attractor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hydrodynamic nonlinear response of interacting integrable systems

In this work, we develop a formalism for computing the nonlinear response of interacting integrable systems. Our results are asymptotically exact in the hydrodynamic limit where perturbing fields vary sufficiently slowly in space and time. We show that spatially resolved nonlinear response distinguishes interacting integrable systems from noninteracting ones, exemplifying this for the Lieb–Liniger gas. We give a prescription for computing finite-temperature Drude weights of arbitrary order, which is in excellent agreement with numerical evaluation of the third-order response of the XXZ spin chain. We identify intrinsically nonperturbative regimes of the nonlinear response of integrable systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Amphiphilic Block Copolymers for Flocculation and Hydrophobization of Legacy Waste Suspensions in Flotation Driven Dewatering Operations

There exists a degrading legacy fuel pond at the Sellafield site (UK), known as the First-Generation Magnox Storage Pond (FGMSP). After being placed into a passive care and maintenance regime resulting in a long storage period, the fuel rods, primarily their cladding, corroded in the pond forming magnesium hydroxide based sludges and suspended material. These ponds are a grave concern to the British government in terms of hazard and risk reduction. Therefore, decommissioning of these ponds is a top priority beginning with the contaminant retrievals process [1]. The pond has accumulated significant quantities of waste materials amongst the skips of fuel, including but not limited to: large inventories of corroded Magnox sludge, fuel rod fragments, metal fragments (from fuel skips), concrete degradation products (from the pond infrastructure), wind-blown sand, and other materials such as bird guano and animal remains. This challenge requires a chemically robust technology to complete sludge retrievals [2]. As nuclear is different, a more stringent process operation criteria is required. From the criteria, flotation was selected as a viable technology. Flotation involves the application of collector molecules which modify the hydrophobicity of suspended particles allowing them to adsorb to rising bubbles. As some particles lack inertia for flotation due to their size, dual flocculant/collector agents can be deployed for greater particle recovery, in this case, amphiphilic diblock copolymers. Two copolymers of Poly(acrylic acid)-b-poly(n-butyl acrylate) (or PAA-b-PnBA), of different hydrophobic chain lengths (PBA) were synthesized for flotation campaigns. Flotation has shown promise to be a valuable rapid dewatering strategy for decommissioning of legacy waste ponds whilst upholding the required operational criteria. This research has shown: PAA-b-PnBA copolymers promote flocculation of Mg(OH){sub 2} particulates. Diblock copolymers retain more fluid than traditional surfactant based collectors wrt. relative particulate recovery performance- where longer PnBA chain length polymers performed best. Performance appears to be hydrodynamically limited due to over-flocculation of Mg(OH){sub 2}. Potential to combine with sedimentation for high particle recovery with low fluid carryover.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Amphiphilic Block Copolymers for Flocculation and Hydrophobization of Legacy Waste Suspensions in Flotation Driven Dewatering Operations - 20162

Flotation has been effectively demonstrated as a low footprint, high efficiency separation process when utilising simple anionic surface modifying agents such as sodium dodecyl sulphate (SDS). These agents increase hydrophobicity to remove suspended cationic surface charged particulates such as Mg(OH){sub 2} from waste suspensions similar to those at British nuclear fuel management sites- where corrosion processes have degraded fuel cladding. The technology is advantageous as it could effectively be retrofitted into existing waste management facilities at nuclear waste management sites such Sellafield and Hanford as the simplicity (no moving parts) and size of this technology coupled with low cost of construction could be rapidly deployed to aid in risk and hazard reduction decommissioning operations with minimum impact to secondary waste generation. Previously, sodium dodecyl sulphate has been deployed as a collector agent to increase the hydrophobicity of Mg(OH){sub 2} suspended particles, achieving a maximum particulate extraction percentage of 93%. To achieve further particulate extraction, flotation performance has been determined to be a function of the hydrodynamic limitations extracting fine material (>20 μm). Fine particulates have been shown in the minerals industry to have an inability to overcome slipstreams created from rising bubbles in flotation cells, likely due to the low mass and thus momentum of the colloidal material, preventing interaction between the hydrophobized particles and air-water interfaces (bubbles). To combat this mass transfer resistance poly(acrylic acid)-block-poly(n-butyl acrylate) amphiphilic diblock copolymers were synthesised via reversible addition-fragmentation chain-transfer (RAFT) polymerisation. The anticipated mode of action of these polymers involves the hydrophilic poly(acrylic acid) block polymer acting as a flocculation agent to increase the mass of the aggregates, and the hydrophobic poly(n-butyl acrylate) modifies the hydrophobicity of the resulting flocs. The concentration of this polymer (termed a macro-collector) was varied and the change in particle size distribution of the particulate suspension was analysed using static light scattering. The flotation performance was then characterised by several methods including percentage particulate recovery, volume reduction factor and residual bulk particulate concentration. This was then benchmarked against the traditional surfactant collector agent SDS. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A window into atmospheric escape on early Mars provided by argon isotopes in a Martian meteorite

Novel analytical approaches to determining the trapped 38 Ar/ 36 Ar ratio of gases contained within Martian meteorites are presented and applied to the Martian regolith breccia Northwest Africa (NWA) 7034 and paired stone NWA 11220. The resulting data indicate that extensive mass-dependent fractionation of atmospheric Ar may have occurred within 150 million years of planetary formation, ostensibly as a result of diffusion-limited hydrodynamic escape. The inferred fractional loss of Ar and lighter atmospheric constituents exceeds 50%. In conclusion, these data suggest that volatiles derived from planetary outgassing and/or impactors may dominate the present abundance of atmospheric Ar.

07 ISOTOPE AND RADIATION SOURCES↗

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Self-interacting dark scalar spikes around black holes via relativistic Bondi accretion

For this work, we consider the spike mass density profile in a dark halo by self-consistently solving the relativistic Bondi accretion of dark matter onto a non-spining black hole of mass M. We assume that the dominant component of the dark matter in the halo is a Standard model gauge-singlet scalar. Its mass m ≃ 10 -5 eV and quartic self-coupling λ ≲ 10 19 are constrained to be compatible with the properties of galactic dark halos. In the hydrodynamic limit, we find that the accretion rate is bounded from below, M min = 96πG 2 M 2 m 4 /λ$\hbar$ 3 . Therefore, for M = 10 6 M ⊙ we have M min ≃ 1.41 x 10 -9 M ⊙ yr -1 , which is subdominant compared to the Eddington accretion of baryons. The spike density profile ρ0(r) within the self-gravitating regime cannot be fitted well by a single-power law but a double-power one. Despite that, we can fit ρ 0 (r) piecewise and find that ρ 0 (r) ∝ r -1.20 near the sound horizon, ρ 0 (r) ∝ r -1.00 towards the Bondi radius and ρ 0 (r) ∝ r -1.08 for the region in between. This contrasts with more cuspy ρ 0 (r) ∝ r -1.75 for dark matter with Coulomb-like self-interaction.

79 ASTRONOMY AND ASTROPHYSICS↗

Static and dynamic properties of multi-ionic plasma mixtures

Complex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multi-component systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion and viscosity for ternary mixtures at extreme conditions. These quantities can be produced from first principles using orbital free molecular dynamics at the computational expense of very intensive simulations to reach good statistics. Utilizing the first-principles results as reference data, we assess the merit of a global analytic model for transport coefficients, "Pseudo-Ions in Jellium" (PIJ), based on an iso-electronic assumption (iso-n e ). With a multi-component hypernetted-chain integral equation, we verify the quality of the iso-n e prescription for describing the static structure of the mixtures. This semi-analytical modeling compares well with the simulation results and allows one to consider plasma mixtures not accessible to simulations. Further, applications are given for the mix of materials in ICF experiments. A reduction of a multicomponent mixture to an effective binary mixture is also established in the hydrodynamic limit and compared with PIJ estimations for ICF relevant mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗