Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic efficiency”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effective field theory of stochastic diffusion from gravity

Planar black holes in AdS have long-lived quasinormal modes which capture the physics of charge and momentum diffusion in the dual field theory. How should we characterize the effective dynamics of a probe system coupled to the conserved currents of the dual field theory? Specifically, how would such a probe record the long-lived memory of the black hole and its Hawking fluctuations? We address this question by exhibiting a universal gauge invariant framework which captures the physics of stochastic diffusion in holography: a designer scalar with a gravitational coupling governed by a single parameter, the Markovianity index. We argue that the physics of gauge and gravitational perturbations of a planar Schwarzschild-AdS black hole can be efficiently captured by such designer scalars. We demonstrate that this framework allows one to decouple, at the quadratic order, the long-lived quasinormal and Hawking modes from the short-lived ones. It furthermore provides a template for analyzing fluctuating open quantum field theories with memory. In particular, we use this set-up to analyze the diffusive Hawking photons and gravitons about a planar Schwarzschild-AdS black hole and derive the quadratic effective action that governs fluctuating hydrodynamics of the dual CFT. Along the way we also derive results relevant for probes of hyperscaling violating backgrounds at finite temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Automatic Time Step Control to Resolve Hydromechanically Driven Fault Reactivation, Spontaneous Nucleation, and Seismic Arrest

Abstract A physical understanding of the progression from flow‐driven (quasi‐static) poromechanical deformation to dynamic fault rupture is critical to the resilient operations of several engineering systems. These processes are bridged by a progression from fault reactivation to the spontaneous nucleation of unstable sliding. Toward addressing this challenge, novel automatic time step size control methods are developed to enable accurate and efficient simulation of these dynamics and transitions from the first principles. The controllers combine local models for discretization error and Coulomb failure conditions to automatically adjust the time step size across several orders of magnitude. The methods do not require additional empirical or theoretical input and can resolve the pre‐rupture, interseismic, and seismic periods to the allowed accuracy. The computational results reveal that the proposed methods automatically capture the onset of reactivation and nucleation for homogeneous and heterogeneous fields. Hydrodynamic and structural heterogeneity lead to disparate critical nucleation sizes compared to those predicted by theory. The results highlight its potential in predicting induced seismicity in realistic subsurface engineering systems and at practical scales.

Environmental Sciences & Ecology↗
NASA NTRS1994

Shock propagation and the generation of magnetohydrodynamic wave fields in inhomogeneous molecular clouds

We develop a simple one-dimensional model for the interaction of a steady, thin, planar shock wave with a nonrigid cloud which may be in motion relative to the surrounding medium, and we apply the model to shocks impinging on, and propagating through, molecular clouds. Both 'adiabatic' (gamma = 5/3) and radiative (gamma = 1) shocks are considered and we allow for the presence of a uniform magnetic field directed either parallel or perpendicular to the shock normal. The former field orientation is equivalent to the hydrodynamic case, and the latter involves only fast MHD shocks. We focus on the manner in which such shocks can generate internal kinetic motions in the cloud on a range of size and density scales through the direct acceleration of cores and clumps by shocks transmitted into them and through the generation of an MHD wavefield via the reflection of the incident shock at clump boundaries. We find that stronger incident Mach numbers and smaller density contrasts lead to more efficient cloud acceleration, as do isothermal intercloud shocks and small intercloud magnetic field strengths. The acceleration efficiency is insensitive to the adiabatic index and the magnetic field strength in the cloud itself. For typical parameter choices, the direct acceleration of clouds and clumps by strong shocks is found to be substantial and could at least in part account for their observed velocity dispersions. If the shocks are moderately weak, the final velocity of the cloud is linearly related to its initial velocity, with higher acceleration giving shallower slopes (i.e., final velocity distributions which are less sensitive to the initial distribution). Compared to the kinetic energy of the postshock cloud, the energy given to the wavefield at each encounter is small, and the heating of the interclump medium by the dissipation of this wavefield is found to be insufficient to balance the cooling rate in the cloud as a whole (although it may be important in particular regions), even if this medium is warm, unless it is also extremely tenuous (n approximately less than 0.1/cu cm). Nevertheless, the correction for the velocity imparted to the cloud leads to a substantial increase in the critical incident Mach numner for wave emission over that reported by Spitzer for the rigid case. The implications of our model for shock-induced star formation are discussed briefly.

Miesch, Mark S.↗
DOE OSTI2019

From chiral anomaly to two-fluid hydrodynamics for electronic vortices

Highlights: • Chiral transformations with account for chiral anomaly were applied to the CDW model. • Two-fluid hydrodynamics was constructed of the order parameter and the normal liquid. • Topologically nontrivial dynamics appears under applied fields or charge injection. • Spontaneous spacetime vorticity under the driving current was modeled numerically. • Periodic sequence of amplitude nodes with 2π phase increments was demonstrated. - Abstract: Many recent experiments addressed manifestations of electronic crystals, particularly the charge density waves, in nano-junctions, under electric field effect, at high magnetic fields, together with real space visualizations by STM and micro X-ray diffraction. This activity returns the interest to stationary or transient states with static and dynamic topologically nontrivial configurations: electronic vortices as dislocations, instantons as phase slip centers, and ensembles of microscopic solitons. Describing and modeling these states and processes calls for an efficient phenomenological theory which should take into account the degenerate order parameter, various kinds of normal carriers and the electric field. Here we notice that the commonly employed time-depend Ginzburg–Landau approach suffers with violation of the charge conservation law resulting in unphysical generation of particles which is particularly strong for nucleating or moving electronic vortices. We present a consistent theory which exploits the chiral transformations taking into account the principle contribution of the fermionic chiral anomaly to the effective action. The resulting equations clarify partitions of charges, currents and rigidity among subsystems of the condensate and normal carriers. On this basis we perform the numerical modeling of a spontaneously generated coherent sequence of phase slips – the spacetime vortices – serving for the conversion among the injected normal current and the collective one.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

The dual action gas thrust bearing - A new high load bearing concept

The principle of utilizing hydrodynamic effects in diverging films for improving load capacity in gas thrust bearings is discussed. A new concept of dual action bearing based on that principle is described and analyzed. The potential of the new bearing is demonstrated both analytically for an infinitely long slider and by numerical solution for a flat sector shaped thrust bearing. It is shown that the dual action bearing can extend substantially the range of load carrying capacity in gas lubricated thrust bearings and improve their efficiency.

Etsion, I.↗
NASA NTRS2013

New Evidence for Efficient Collisionless Heating of Electrons at the Reverse Shock of a Young Supernova Remnant

Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the shock front, is one of the main open issues in shock physics. Here, we present the first clear evidence for efficient collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe K diagnostics using high-quality X-ray data obtained by the Suzaku satellite. We detect K beta (3p yields 1s) fluorescence emission from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a smaller radius than Fe K alpha (2p yields 1s) emission dominated by a relatively highly ionized component. Comparisons with our hydrodynamical simulations imply instantaneous electron heating to a temperature 1000 times higher than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the physical mechanism responsible for this heating and favors a cross-shock potential created by charge deflection at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% smaller than the previous measurement using the Fe K alpha morphology from the Chandra observations. Since strong Fe K beta fluorescence is expected only from low-ionization plasma where Fe ions still have many 3p electrons, this feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR.

Fe K diagnostics↗
DOE OSTI2016

An efficient direct solver for rarefied gas flows with arbitrary statistics

A new numerical methodology associated with a unified treatment is presented to solve the Boltzmann–BGK equation of gas dynamics for the classical and quantum gases described by the Bose–Einstein and Fermi–Dirac statistics. Utilizing a class of globally-stiffly-accurate implicit–explicit Runge–Kutta scheme for the temporal evolution, associated with the discrete ordinate method for the quadratures in the momentum space and the weighted essentially non-oscillatory method for the spatial discretization, the proposed scheme is asymptotic-preserving and imposes no non-linear solver or requires the knowledge of fugacity and temperature to capture the flow structures in the hydrodynamic (Euler) limit. The proposed treatment overcomes the limitations found in the work by Yang and Muljadi (2011) [33] due to the non-linear nature of quantum relations, and can be applied in studying the dynamics of a gas with internal degrees of freedom with correct values of the ratio of specific heat for the flow regimes for all Knudsen numbers and energy wave lengths. The present methodology is numerically validated with the unified treatment by the one-dimensional shock tube problem and the two-dimensional Riemann problems for gases of arbitrary statistics. Descriptions of ideal quantum gases including rotational degrees of freedom have been successfully achieved under the proposed methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A thermodynamic cycle for the solar cell

A solar cell is a heat engine, but textbook treatments are not wholly satisfactory from a thermodynamic standpoint, since they present solar cells as directly converting the energy of light into electricity, and the current in the circuit as maintained by an electrostatic potential. Here we propose a thermodynamic cycle in which the gas of electrons in the ρ phase serves as the working substance. The interface between the ρ and n phases acts as a self-oscillating piston that modulates the absorption of heat from the photons so that it may perform a net positive work during a complete cycle of its motion, in accordance with the laws of thermodynamics. We draw a simple hydrodynamical analogy between this model and the “putt-putt” engine of toy boats, in which the interface between the water’s liquid and gas phases serves as the piston. We point out some testable consequences of this model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Early hydrodynamic evolution of a stellar collision

The early phase of the hydrodynamic evolution following the collision of two stars is analyzed. Two strong shocks propagate from the contact surface and move toward the center of each star at a velocity that is a small fraction of the velocity of the approaching stars. The shocked region near the contact surface has a planar symmetry and a uniform pressure. The density vanishes at the (Lagrangian) surface of contact, and the speed of sound diverges there. The temperature, however, reaches a finite value, since as the density vanishes, the finite pressure is radiation dominated. For carbon-oxygen white dwarf (CO WD) collisions, this temperature is too low for any appreciable nuclear burning shortly after the collision, which allows for a significant fraction of the mass to be highly compressed to the density required for efficient {sup 56}Ni production in the detonation wave that follows. This property is crucial for the viability of collisions of typical CO WD as progenitors of type Ia supernovae, since otherwise only massive (>0.9 M {sub ☉}) CO WDs would have led to such explosions (as required by all other progenitor models). The divergence of the speed of sound limits numerical studies of stellar collisions, as it makes convergence tests exceedingly expensive unless dedicated schemes are used. We provide a new one-dimensional Lagrangian numerical scheme to achieve this. A self-similar planar solution is derived for zero-impact parameter collisions between two identical stars, under some simplifying assumptions (including a power-law density profile), which is the planar version of previous piston problems that were studied in cylindrical and spherical symmetries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Detailed simulations of the first deuterium-tritium-filled double shell implosions on the National Ignition Facility

The first indirectly driven, liquid DT-filled double shell inertial confinement fusion (ICF) implosions have recently been successfully performed on the National Ignition Facility (NIF). Double shells are a class of alternative designs that use a low-Z outer shell to compress a foam cushion that accelerates a high-Z inner shell to efficiently compress a liquid DT core. Double shells are challenging to fabricate, field, and model. Important engineering features enabling double shell fabrication include a fill-tube penetrating all shells and a carefully designed and very narrow (few μm) step-joint in the ablator. Due to the higher density materials involved, high Atwood number instabilities are also important at many material interfaces. In this paper, numerical simulations of double shell implosions using the Los Alamos National Laboratory multi-physics radiation-hydrodynamics code xRAGE will be discussed. An extensive effort has been under way for several years to develop the code capabilities for ICF simulations in a common modeling framework to allow ease of simulation setup and standardization of the computational methodology. This paper will present a wide range of simulation results capturing, quantifying, and comparing the impact of all these degradation mechanisms on implosion performance. Brief comparisons with recent experimental results and suggestions for future improvements will also be discussed. Our results suggest that capsule surface roughness and the step-joint gap have the largest impact on implosion performance. Initial experimental data may suggest that the sensitivity to the step-joint gap could provide the dominant explanation for DT-filled double shell experiments that have been fielded on NIF thus far.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

EXPLAINING THE MOST ENERGETIC SUPERNOVAE WITH AN INEFFICIENT JET-FEEDBACK MECHANISM

We suggest that the energetic radiation from core-collapse super-energetic supernovae (SESNe) is due to a long-lasting accretion process onto the newly born neutron star (NS), resulting from an inefficient operation of the jet-feedback mechanism (JFM). The jets that are launched by the accreting NS or black hole maintain their axis due to a rapidly rotating pre-collapse core and do not manage to eject core material from near the equatorial plane. The jets are able to eject material from the core along the polar directions and reduce the gravity near the equatorial plane. The equatorial gas expands, and part of it falls back over a timescale of minutes to days to prolong the jet-launching episode. According to the model for SESNe proposed in the present paper, the principal parameter that distinguishes between the different cases of core-collapse supernova (CCSN) explosions, such as between normal CCSNe and SESNe, is the efficiency of the JFM. This efficiency, in turn, depends on the pre-collapse core mass, envelope mass, core convection, and, most of all, the angular momentum profile in the core. One prediction of the inefficient JFM for SESNe is the formation of a slow equatorial outflow in the explosion. The typical velocity and mass of this outflow are estimated to be v {sub eq} ≈ 1000 km s{sup −1} and M {sub eq} ≳ 1 M {sub ⊙}, respectively, though quantitative values will have to be checked in future hydrodynamic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

The rate of collisions due to Brownian or gravitational motion of small drops

Quantitative predictions of the collision rate of two spherical drops undergoing Brownian diffusion or gravitational sedimentation are presented. The diffusion equation for relative Brownian motion of two drops is derived, and the relative motion of pairs of drops in gravitational sedimentation is traced via a trajectory analysis in order to develop theoretical models to determine the collision efficiencies, both with and without interparticle forces applied between the drops. It is concluded that finite collision rates between nondeforming fluid drops are possible for Brownian diffusion or gravitational sedimentation in the absence of attractive forces, in stark contrast to the prediction that lubrication forces prevent rigid spheres from contacting each other unless an attractive force that becomes infinite as the separation approaches zero is applied. Collision rates are shown to increase as the viscosity of the drop-phase decreases. In general, hydrodynamic interactions reduce the collision rates more for gravitational collisions than for Brownian collisions.

Zhang, Xiaoguang↗
DOE OSTI2014

Accelerating dissipative particle dynamics simulations for soft matter systems

Dissipative particle dynamics (DPD) is a coarse-grained particle-based simulation method that offers microscopic-scale insights into soft matter systems. We present an efficient implementation of a DPD model for graphical processing units (GPUs). As implemented in the LAMMPS molecular dynamics package, it can run effectively on current-generation supercomputers which often have hybrid nodes containing multi-core CPUs and (one or more) GPUs. Using efficient communication of information between the CPUs and GPUs, DPD interactions can be computed on the GPU while other portions of a full simulation model (boundary conditions, constraints, bonded interactions, diagnostic calculations, etc.) can be performed on the CPU. Our GPU-enhanced runs show a speedup of up to 9.5× versus many-core CPU simulations, and can run scalably across thousands of compute nodes. We briefly discuss how the new GPU implementation was validated against the CPU version for thermodynamics, diffusion, and hydrodynamic behavior. We also highlight large-scale models which the faster DPD implementation has enabled, for studies of monolayer self-assembly and thin-film instabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Axisymmetric electrostatic magnetohydrodynamic oscillations in tokamaks with general cross-sections and toroidal flow

The frequency spectrum and mode structure of axisymmetric electrostatic oscillations [the zonal flow (ZF), sound waves (SW), geodesic acoustic modes (GAM), and electrostatic mean flows (EMF)] in tokamaks with general cross-sections and toroidal flows are studied analytically using the electrostatic approximation for magnetohydrodynamic modes. These modes constitute the “electrostatic continua.” Starting from the energy principle for a tokamak plasma with toroidal rotation, we showed that these modes are completely stable. The ZF, the SW, and the EMF could all be viewed as special cases of the general GAM. The Euler equations for the general GAM are obtained and are solved analytically for both the low and high range of Mach numbers. The solution consists of the usual countable infinite set of eigen-modes with discrete eigen-frequencies, and two modes with lower frequencies. The countable infinite set is identified with the regular GAM. The lower frequency mode, which is also divergence free as the plasma rotation tends to zero, is identified as the ZF. The other lower (zero) frequency mode is a pure geodesic E×B flow and not divergence free is identified as the EMF. The frequency of the EMF is shown to be exactly 0 independent of plasma cross-section or its flow Mach number. We also show that in general, sound waves with no geodesic components are (almost) completely lost in tokamaks with a general cross-sectional shape. The exception is the special case of strict up-down symmetry. In this case, half of the GAMs would have no geodesic displacements. They are identified as the SW. Present day tokamaks, although not strictly up-down symmetric, usually are only slightly up-down asymmetric. They are expected to share the property with the up-down symmetric tokamak in that half of the GAMs would be more sound wave-like, i.e., have much weaker coupling to the geodesic components than the other half of non-sound-wave-like modes with stronger coupling to the geodesic displacements. Based on the general notion that the geodesic component of the GAM is more effective in tearing up the eddies in the electrostatic turbulence, it is important to preferentially excite the GAMs that are non-sound-wave like to maximize the efficiency on turbulence suppression through external means. Finally, approximate formulae for the frequencies of the EMF, ZF, SW, and the GAM for a large aspect ratio circular tokamak rotating at low Mach numbers are also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neutrino signatures of 100 2D Axisymmetric Core-Collapse Supernova Simulations

ABSTRACT We present in this paper a public data release of an unprecedentedly large set of core-collapse supernova (CCSN) neutrino emission models, comprising 100 detailed 2D axisymmetric radiation-hydrodynamic simulations evolved out to as late as ∼5 s post-bounce and spanning an extensive range of massive-star progenitors. The motivation for this paper is to provide a physically and numerically uniform benchmark data set to the broader neutrino detection community to help it characterize and optimize subsurface facilities for what is likely to be a once-in-a-lifetime galactic supernova burst event. With this release, we hope to (1) help the international experiment and modelling communities more efficiently optimize the retrieval of physical information about the next galactic CCSN, (2) facilitate the better understanding of core-collapse theory and modelling among interested experimentalists, and (3) help further integrate the broader supernova neutrino community.

Astronomy & Astrophysics↗
DOE OSTI2026

CHEQUP v0.1

CHEQUP (Castro-based Hofi Expansion with QUasineutral Plasma) is a simulation code for modeling the formation of hydrodynamic optical-field-ionized (HOFI) plasma channels, which are used as waveguides in laser-plasma acceleration experiments. This includes experiments performed at LBNL's BELLA facility as well as other laser facilities across the world. CHEQUP extends the open-source Castro hydrodynamics framework with physics modules tailored for modeling HOFI plasma channels -- including multi-species ionization and three-body recombination for mixtures of hydrogen, nitrogen, helium, and argon ; a two-temperature model tracking electron and heavy-species temperatures separately ; and coupling with other codes of the BLAST ecosystem (https://blast.lbl.gov/) such as WarpX, via the openPMD standard. CHEQUP inherits from Castro the ability to run on modern GPU architectures (NVIDIA CUDA, AMD HIP) and supports adaptive mesh refinement (AMR) for efficient multi-scale resolution. Compared to existing tools, CHEQUP would be, to our knowledge, the first open-source code implementing the full HOFI channel formation physics, and the first implementation capable of running on GPUs. This enables significantly faster, large-scale parameter scans critical for the design of next-generation LPA-based accelerators and light sources.

Lehe, Remi [Lawrence Berkeley National Laboratory ↗
DOE OSTI2025

Terrestrial atmospheric ion implantation occurred in the nearside lunar regolith during the history of Earth’s dynamo

Light volatile elements in lunar regolith are thought to be a mixture of the solar wind and Earth’s atmosphere, the latter sourced in the absence of geomagnetic field. However, the extent to which both the current and primitive geodynamo influence the transport of terres trial ions still remains unclear, and this uncertainty is further complicated by the enigmatic composition and poorly constrained location of the Eoarchean exosphere. Here we use 3-D MHD numerical simulations with contemporary magnetized and Archean unmagnetized atmospheres to investigate how Earth’s intrinsic magnetic field affects this transfer, aiming to constrain how and when the lunar isotopic signature provides a record of Earth’s pale-oatmosphere. We find that atmospheric transfer is efficient only when the Moon is within Earth’s magnetotail. The non-solar contribution to the lunar soil is best explained by implantation during the long history of the geodynamo with present-day solar wind conditions, rather than any short, putatively unmagnetized epoch of early Archean Earth. This further suggests the history of the terrestrial atmosphere, spanning billions of years, could be preserved in buried lunar soils. Our results indicate that the elemental abundances of Apollo samples are very sensitive to Earth’s hydrodynamic exobase altitude, which, at the time of ion implantation, was never smaller than 190 km.

Atmospheric dynamics↗
NASA NTRS1976

Dual-action gas thrust bearing for improving load capacity

The principle of utilizing hydrodynamic effects in diverging films to improve the load carrying capacity in gas thrust bearings is discussed. A new concept of a dual action bearing based on that principle is described and analyzed. The potential of the new bearing is demonstrated both analytically for an infinitely long slider and by numerical solution for a flat sector shaped thrust bearing. It is shown that the dual action bearing can extend substantially the range of load carrying capacity in gas lubricated thrust bearings and can improve their efficiency.

Etsion, I.↗