Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hydrodynamic simulations of electron-capture supernovae: progenitor and dimension dependence

ABSTRACT We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in flash with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with six progenitor models that, at the time of collapse, span a range of propagating flame front radii. For comparison, we also perform a simulation with the same set-up using the canonical, spherically symmetrical progenitor model n8.8. We found that the variations in the progenitor models inherited from simulations of stellar evolution and flame propagation do not significantly alter the global properties of the neutrino-driven ECSN explosion, such as the explosion energy (∼1.36–1.48 × 1050 erg) and the mass (∼0.017–0.018 M⊙) and composition of the ejecta. Due to aspherical perturbations induced by the 2D flame, the ejecta contains a small amount (≲1.8 × 10−3 M⊙) of low-Ye (0.35 < Ye < 0.4) component. The baryonic mass of the protoneutron star is ∼1.34 M⊙ (∼1.357 M⊙) with the new (n8.8) progenitor models when simulations end at ∼400 ms and the discrepancy is due to updated weak-interaction rates in the progenitor evolutionary simulations. Our results reflect the nature of ECSN progenitors containing a strongly degenerate oxygen–neon–magnesium (ONeMg) core and suggest a standardized ECSN explosion initialized by ONeMg core collapse. Moreover, we carry out a rudimentary three-dimensional simulation and find that the explosion properties are fairly compatible with the 2D counterpart. Our paper facilitates a more thorough understanding of ECSN explosions following the ONeMg core collapse, though more three-dimensional simulations are still needed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Hydrodynamic simulations of disrupted planetary accretion discs inside the core of an AGB star

ABSTRACT Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system, suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc. Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Hydrodynamics of nonintegrable systems from a relaxation-time approximation

We develop a general kinetic theory framework to describe the hydrodynamics of strongly interacting, nonequilibrium quantum systems in which integrability is weakly broken, leaving a few residual conserved quantities. This framework is based on a generalized relaxation-time approximation; it gives a simple, but surprisingly accurate, prescription for computing nonequilibrium transport even in strongly interacting systems. We validate the predictions of this approximation against matrix product operator calculations on chaotic quantum spin chains, finding surprisingly good agreement. Lastly, we show that despite its simplicity, our framework can capture phenomena distinctive to strongly interacting systems, such as widely separated charge and energy diffusion constants.

1-dimensional spin chains↗
DOE OSTI2023

Thermal transfer enhancement by hydrodynamic plasmons in electron bilayers

Here, we develop a theory of heat transfer induced by thermal charge fluctuations in two-dimensional electron double layers. We consider pristine systems comprised of identical layers and focus on the regime of sufficiently high temperatures and interlayer distances 𝑑, where the relevant charge fluctuations may be described using the hydrodynamic approach. In this limit heat transfer is dominated by the plasmon resonances. For systems with Galilean-invariant electron dispersion the interlayer thermal conductance 𝜘 is proportional to the kinematic viscosity of the electron liquid and decreases as 1/𝑑 4 . In the absence of Galilean invariance 𝜘∝𝜎/𝑑 3 , where 𝜎 is the intrinsic conductivity of the liquid. This strong enhancement can be traced to a drastically different broadening of plasmon resonances in systems with and without Galilean invariance.

2-dimensional systems↗
DOE OSTI2020

Pion interferometry at 200 GeV using anisotropic hydrodynamics

In this paper, we continue our phenomenological studies of heavy-ion collisions using 3 + 1 D anisotropic hydrodynamics (aHydro). In previous works, we compared quasiparticle aHydro (aHydroQP) with ALICE 2.76 TeV Pb-Pb and BNL Relativistic Heavy Ion Collider 200 GeV Au-Au collision results. At both energies, the agreement was quite good between aHydroQP and the experimental data for many observables. In this work, we present comparisons of the Hanbury Brown–Twiss (HBT) radii and their ratios determined using π + π + pairs produced in 200 GeV Au-Au collisions. We first present comparisons with STAR results for the HBT radii and their ratios. We then present comparisons with PHENIX results for the HBT radii and their ratios. In both cases, we find reasonable agreement between aHydroQP predictions and available experimental results for the ratios of HBT radii. At the level of the radii themselves, in some cases quantitative differences on the order of 10–20% remain, which deserve further research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗
DOE OSTI2021

Observation of Hydrodynamic Flows in Imploding Fusion Plasmas on the National Ignition Facility [plus supplemental information]

Inertial confinement fusion implosions designed to have minimal fluid motion at peak compression often show significant linear flows in the laboratory, attributable per simulations to percent-level imbalances in the laser drive illumination symmetry. In this work, we present experimental results which intentionally varied the mode 1 drive imbalance by up to 4% to test hydrodynamic predictions of flows and the resultant imploded core asymmetries and performance, as measured by a combination of DT neutron spectroscopy and high-resolution x-ray core imaging. Neutron yields decrease by up to 50%, and anisotropic neutron Doppler broadening increases by 20%, in agreement with simulations. Furthermore, a tracer jet from the capsule fill-tube perturbation that is entrained by the hot-spot flow confirms the average flow speeds deduced from neutron spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effective Action for Relativistic Hydrodynamics from the Crooks Fluctuation Theorem

A new effective theory framework for fluctuating hydrodynamics in the relativistic regime is derived using standard thermodynamical principles and general properties of nonequilibrium stochastic dynamics. For the first time, we establish clear and concise conditions for ensuring that the resulting effective theories are causal, stable, and well-posed within general relativity. Here, these properties are independent of spacetime foliation and are valid in the full nonlinear regime. Out-of-equilibrium fluctuations are constrained by a relativistically covariant version of the Crooks fluctuation theorem, which determines how the entropy production is distributed even when the system is driven by an external force. This leads to an emerging Z 2 symmetry responsible for imposing fluctuation-dissipation relations for n-point correlation functions, which matches the standard constraints for the Schwinger-Keldysh effective action.

Fluctuation theorems↗
DOE OSTI2025

Polarization of the 𝜙 meson in the hadronic phase with nucleon scatterings and a viscous hydrodynamic background

We extend our previous work on the spin alignment of the ϕ vector meson in the hadronic phase of the quark gluon plasma, by including effects of nucleon scatterings. The emission rates are calculated in a realistic hydrodynamic background simulated with the code FLUIDuM , for different beam energies. We find that all the effects taken into account cannot explain the out-of-plane spin alignment of the ϕ meson observed experimentally.

Polarization phenomena↗
DOE OSTI2023

COUPLING SMOOTHED PARTICLE HYDRODYNAMICS WITH FINITE ELEMENT METHOD TO SIMULATE RESIDUAL STRESSES FROM FRICTION STIR PROCESSING

Friction stir processing (FSP) is a solid-state material processing technique that locally modifies the microstructure but also induces undesirable residual stresses. A robust numerical model for the FSP can help in mitigating these residual stresses. Heat source models within a finite element method (FEM) framework suffer from inaccuracies. In contrast, smoothed particle hydrodynamics (SPH) model that explicitly captures the material flow near the tool and the associated heat generation are accurate. However, the computational expense of SPH simulations can be prohibitive. In this work, we propose a coupled SPH-FEM framework. SPH is used to model the heat generation accurately near the tool and which is then inserted into to the FEM model as a heat source. To verify this proposed coupling approach, a test case is set up with typical FSP conditions and it is modeled in both SPH and SPH-FEM. The temperatures profiles were compared after the simulations have reached steady-state temperatures. The similarity of the temperature profiles from SPH-FEA and SPH validated the proposed coupling approach. This proposed approach achieves the accuracy of the SPH method while potentially retaining the low computational expense of FEM.

smoothed particle hydrodynamics, Finite Element Me↗
DOE OSTI2020

High-performance computing in water resources hydrodynamics

In this work, we present a vision of future water resources hydrodynamics codes that can fully utilize the strengths of modern high-performance computing. The advances to computing power, formerly driven by the improvement of central processing unit processors, now focus on parallel computing and, in particular, the use of graphics processing units (GPUs). However, this shift to a parallel framework requires refactoring the code to make efficient use of the data as well as changing even the nature of the algorithm that solves the system of equations. These concepts along with other features such as the precision for the computations, dry regions management, and input/output data are analyzed in this paper. A 2D multi-GPU flood code applied to a large-scale test case is used to corroborate our statements and ascertain the new challenges for the next-generation parallel water resources codes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

HED Hydrodynamics in the Common Modeling Framework

Within the Common Modeling Framework (CMF), the HED Hydrodynamics Project has implemented a common methodology, or "Authority", for modeling high energy density (HED) experiments by extending development of the Eulerian Applications Project code Authority, "EAPA," to include relevant physics models, customization tools and templates, and in situ and post-processing capabilities. Additionally, several models for small-scale HED experiments were successfully added to the CMF and simulations reproduce the results of preexisting, benchmarked input decks. We expect use of the CMF to improve our effectiveness in developing predictive capabilities for experiments ranging from small-scale planar single-interface, single-shock to multi-interface, multi-shock configurations and up through multi-shell inertial confinement fusion (ICF) implosions. Additionally, this will strengthen our ability to design future targeted experiments. Using the CMF infrastructure and Authorities for HED modeling has three key impacts. First, it enables a version-controlled, pedigreed and archived base model for each experiment with tools that allow easy setup, execution and analysis of simulations to assess and/or expose sensitivities to various parameters. Second, having a common base set of inputs and physics definitions helps identify systematic differences between models in the codes and experimental data to support code validation, and build confidence in our predictive capabilities. Finally, designed to work in concert with multiple programs and projects including PEM, IC, OES, V&V and DSW, the CMF supports improved collaboration and integration through shared model definitions, data, simulations and results, enhancing the workflow required to underwrite the knowledge base and capabilities ultimately required for stockpile stewardship

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

GnarlyX: Eulerian multi-material hydrodynamics coupled to equation of state and hyperelastic, plastic constitutive models

GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.

42 ENGINEERING↗
DOE OSTI2026

Hydrodynamic Analysis and Optimization of Aquantis Marine Turbine: Cooperative Research and Development (Final Report)

The primary aim of this proposal is to improve the accurate prediction of hydrodynamic performance and dynamic load responses of the AQ10 floating axial-flow tidal turbine with a tri-cat mooring configuration. The validation of reduced-order modeling approaches with high-fidelity model will be implemented. Additionally, the frequency response domain, Response Amplitude Floating Wind (RAFT) toolbox plus an optimizer expanded for marine hydrokinetic turbines under the Submarine Hydrokinetic And Riverine Kilo-megawatt. Systems (SHARKS) program will be used for designing and exploring different key design parameters (platform dimension, mooring layout and its parameters) of next marine hydrokinetic (MHK) turbine generation.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Analysis of an Elasto-Hydrodynamic Seal by Using the Reynolds Equation

This paper reports numerical studies of an Elasto-Hydrodynamic (EHD) seal, which is being developed for supercritical CO2 (sCO2) turbomachinery applications. Current sCO2 turbomachinery suffers from high leakage rates, which is creating a major roadblock to the full realization of sCO2 power technology. The high leakage rates not only penalize the efficiencies but also create environmental concerns due to greenhouse effects caused by the increased CO2 discharge to the atmosphere. The proposed EHD seal needs to work at elevated pressures (10–35 MPa) and temperatures (350–700 °C) with low leakage and minimal wear. The unique mechanism of the EHD seal provides a self-regulated constriction effect to restrict the flow without substantial material contact, thereby minimizing leakage and wear. This work utilizes a physics-based modeling approach. The flow through the gradually narrowing seal clearance is modeled by the well-known Reynolds equation in EHD lubrication theory, while the deformation of the seal is modeled by using the governing equations of three-dimensional solid mechanics. As for the solution methodology, COMSOL’s Thin-Film Flow and Solid Mechanics modules were employed with their powerful capabilities. The numerical results were presented and discussed. It was observed that the Reynolds equation fully coupled with the surface deformation was able to successfully capture the constriction effect. The maximum and minimum leakages were calculated to be 2.25 g/s and 0.1 g/s at P = 5.5 MPa and P = 11 MPa for the design seal, respectively. It was interesting to observe that the seal leakage followed a quadratic trend with increasing pressure differential, which can become advantageous for high-pressure applications such as sCO2 power generation technology.

42 ENGINEERING↗