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 109 records · Page 6

Analysis of Hydrodynamic (Landau) Instability in Liquid-Propellant Combustion at Normal and Reduced Gravity

The burning of liquid propellants is a fundamental combustion problem that is applicable to various types of propulsion and energetic systems. The deflagration process is often rather complex, with vaporization and pyrolysis occurring at the liquid/gas interface and distributed combustion occurring either in the gas phase or in a spray. Nonetheless, there are realistic limiting cases in which combustion may be approximated by an overall reaction at the liquid/gas interface. In one such limit, the gas flame occurs under near-breakaway conditions, exerting little thermal or hydrodynamic influence on the burning propellant. In another such limit, distributed combustion occurs in an intrusive regime, the reaction zone lying closer to the liquid/gas interface than the length scale of any disturbance of interest. Finally, the liquid propellant may simply undergo exothermic decomposition at the surface without any significant distributed combustion, such as appears to occur in some types of HydroxylAmmonium Nitrate (HAN)-based liquid propellants at low pressures. Such limiting models have recently been formulated,thereby significantly generalizing earlier classical models that were originally introduced to study the hydrodynamic stability of a reactive liquid/gas interface. In all of these investigations, gravity appears explicitly and plays a significant role, along with surface tension, viscosity, and, in the more recent models, certain reaction-rate parameters associated with the pressure and temperature sensitivities of the reaction itself. In particular, these parameters determine the stability of the deflagration with respect to not only classical hydrodynamic disturbances, but also with respect to reactive/diffusive influences as well. Indeed, the inverse Froude number, representing the ratio of buoyant to inertial forces, appears explicitly in all of these models, and consequently, in the dispersion relation that determines the neutral stability boundaries beyond which steady, planar burning is unstable to nonsteady, and/or nonplanar (cellular) modes of burning. These instabilities thus lead to a number of interesting phenomena, such as the sloshing type of waves that have been observed in mixtures of HAN and TriEthanolAmmonium Nitrate (TEAN) with water. Although the Froude number was treated as an O(1) quantity in these studies, the limit of small inverse Froude number corresponding to the microgravity regime is increasingly of interest and can be treated explicitly, leading to various limiting forms of the models, the neutral stability boundaries, and, ultimately, the evolution equations that govern the nonlinear dynamics of the propagating reaction front. In the present work, we formally exploit this limiting parameter regime to compare some of the features of hydrodynamic instability of liquid-propellant combustion at reduced gravity with the same phenomenon at normal gravity.

Margolis, Stephen B.↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(p)(k), where A(p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for negative values of A(p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. it is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the long-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, S. B.↗

Hydrodynamic Instability in an Extended Landau/Levich Model of Liquid-Propellant Combustion at Normal and Reduced Gravity

The classical Landau/Levich models of liquid-propellant combustion, despite their relative simplicity, serve as seminal examples that correctly describe the onset of hydrodynamic instability in reactive systems. Recently, these two separate models have been combined and extended to account for a dynamic dependence, absent in the original formulations, of the local burning rate on the local pressure and temperature fields. The resulting model admits an extremely rich variety of both hydrodynamic and reactive/diffusive instabilities that can be analyzed either numerically or analytically in various limiting parameter regimes. In the present work, a formal asymptotic analysis, based on the realistic smallness of the gas-to-liquid density ratio, is developed to investigate the combined effects of gravity and other parameters on the hydrodynamic instability of the propagating liquid/gas interface. In particular, an analytical expression is derived for the neutral stability boundary A(sub p)(k), where A(sub p) is the pressure sensitivity of the burning rate and k is the wavenumber of the disturbance. The results demonstrate explicitly the stabilizing effect of gravity on long-wave disturbances, the stabilizing effect of viscosity (both liquid and gas) and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limiting case of weak gravity, it is shown that hydrodynamic instability in liquid-propellant combustion is a long-wave instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumber disturbances. It is also demonstrated that, in general, surface tension and the viscosity of both the liquid and gas phases each produce comparable stabilizing effects in the large-wavenumber regime, thereby providing important modifications to previous analyses in which one or more of these effects were neglected.

Margolis, Stephen B.↗

Hydrodynamic Properties of Planing Surfaces and Flying Boats

The study of the hydrodynamic properties of planing bottom of flying boats and seaplane floats is at the present time based exclusively on the curves of towing tests conducted in tanks. In order to provide a rational basis for the test procedure in tanks and practical design data, a theoretical study must be made of the flow at the step and relations derived that show not only qualitatively but quantitatively the inter-relations of the various factors involved. The general solution of the problem of the development of hydrodynamic forces during the motion of the seaplane float or flying boat is very difficult for it is necessary to give a three-dimensional solution, which does not always permit reducing the analysis to the form of workable computation formulas. On the other had, the problem is complicated by the fact that the object of the analysis is concerned with two fluid mediums, namely, air and water, which have a surface of density discontinuity between them. The theoretical and experimental investigations on the hydrodynamics of a ship cannot be completely carried over to the design of floats and flying-boat hulls, because of the difference in the shape of the contour lines of the bodies, and, because of the entirely different flow conditions from the hydrodynamic viewpoint.

Sokolov, N. A.↗

Hydrodynamics in a Degenerate, Strongly Attractive Fermi Gas

In summary, we use all-optical methods with evaporative cooling near a Feshbach resonance to produce a strongly interacting degenerate Fermi gas. We observe hydrodynamic behavior in the expansion dynamics. At low temperatures, collisions may not explain the expansion dynamics. We observe hydrodynamics in the trapped gas. Our observations include collisionally-damped excitation spectra at high temperature which were not discussed above. In addition, we observe weakly damped breathing modes at low temperature. The observed temperature dependence of the damping time and hydrodynamic frequency are not consistent with collisional dynamics nor with collisionless mean field interactions. These observations constitute the first evidence for superfluid hydrodynamics in a Fermi gas.

Thomas, John E.↗

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

Hot electron preheat in hydrodynamically scaled direct-drive inertial confinement fusion implosions on the NIF and OMEGA

Hot electron preheat has been quantified in warm, directly driven inertial confinement fusion implosions on OMEGA and the National Ignition Facility (NIF), to support hydrodynamic scaling studies. These CH-shell experiments were designed to be hydrodynamically equivalent, spanning a factor of 40 in laser energy and a factor of 3.4 in spatial and temporal scales, while preserving the incident laser intensity of 10 15 W/cm 2 . Experiments with similarly low levels of beam smoothing on OMEGA and NIF show a similar fraction (~0.2%) of laser energy deposited as hot electron preheat in the unablated shell on both OMEGA and NIF and similar preheat per mass (~2 kJ/mg), despite the NIF experiments generating a factor of three more hot electrons (~1.5% of laser energy) than on OMEGA (~0.5% of laser energy). This is plausibly explained by more absorption of hot electron energy in the ablated CH plasma on NIF due to larger areal density, as well as a smaller solid angle of the imploding shell as viewed from the hot electron generating region due to the hot electrons being produced at a larger standoff distance in lower-density regions by stimulated Raman scattering, in contrast to in higher-density regions by two-plasmon decay on OMEGA. Finally, the results indicate that for warm implosions at intensities of around 10 15 W/cm 2 , hydrodynamic equivalence is not violated by hot electron preheat, though for cryogenic implosions, the reduced attenuation of hot electrons in deuterium–tritium plasma will have to be considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transient Relativistic Fluid Dynamics in a General Hydrodynamic Frame

In this work, we propose a new theory of second-order viscous relativistic hydrodynamics which does not impose any frame conditions on the choice of the hydrodynamic variables. It differs from Mueller-Israel-Stewart theory by including additional transient degrees of freedom, and its first-order truncation reduces to Bemfica-Disconzi-Noronha-Kovtun theory. Conditions for causality and stability are explicitly given in the conformal regime. As an illustrative example, we consider Bjorken flow solutions to our equations and identify variables which make a hydrodynamic attractor manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis↗

Simulating Hydrodynamics in Cosmology with CRK-HACC

Abstract We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure formation simulations of the universe. The new framework couples the HACC gravitational N -body solver with a modern smoothed-particle hydrodynamics (SPH) approach called conservative reproducing kernel SPH (CRKSPH). CRKSPH utilizes smoothing functions that exactly interpolate linear fields while manifestly preserving conservation laws (momentum, mass, and energy). The CRKSPH method has been incorporated to accurately model baryonic effects in cosmology simulations—an important addition targeting the generation of precise synthetic sky predictions for upcoming observational surveys. CRK-HACC inherits the codesign strategies of the HACC solver and is built to run on modern GPU-accelerated supercomputers. In this work, we summarize the primary solver components and present a number of standard validation tests to demonstrate code accuracy, including idealized hydrodynamic and cosmological setups, as well as self-similarity measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

On the competition between mixing rate and uniformity in a coaxial hydrodynamic focusing mixer

Fast microfluidic mixers for use with line-of-sight integrating detection schemes pose unique challenges. Such detectors typically cannot discriminate signal from slow moving (e.g. near internal walls) and fast-moving portions of the fluid stream. This convolves reaction rate dynamics with fluid flow residence time dynamics. Further, the small cross sections of typical three-dimensional hydrodynamic focusing devices lead to lower detection signals. The current study focuses on achieving both small time scales of mixing and homogenous residence times. This is achieved by injecting sample through a center capillary and hydrodynamically focusing using a sheath flow within a tapered second capillary. The current design also features a third, larger coaxial capillary. The mixed stream flows into the large cross-section of this third capillary to decelerate and expand the stream by up to 14-fold to improve line-of-sight signal strength of reaction products. Hydrodynamic focusing, mixing, and expansion are studied using analytical and numerical models and also studied experimentally using a fluorescein-iodide quenching reaction. The experimentally validated models are used to explore trade-offs between mixing rate and uniformity. For the first time, this work presents detailed analysis of the Lagrangian time history of species transport during mixing inside coaxial capillaries to measure mixing nonuniformity. The mixing region enables order 100 μs mixing times and residence time widths of the same order (140 μs).

42 ENGINEERING↗

Reduced order models for Lagrangian hydrodynamics

It is reported, as a mathematical model of high-speed flow and shock wave propagation in a complex multimaterial setting, Lagrangian hydrodynamics is characterized by moving meshes, advection-dominated solutions, and moving shock fronts with sharp gradients. These challenges hinder the existing projection-based model reduction schemes from being practical. We develop several variations of projection-based reduced order model techniques for Lagrangian hydrodynamics by introducing three different reduced bases for position, velocity, and energy fields. A time-windowing approach is also developed to address the challenge imposed by the advection-dominated solutions. Lagrangian hydrodynamics is formulated as a nonlinear problem, which requires a proper hyper-reduction technique. Therefore, we apply the over-sampling DEIM and SNS approaches to reduce the complexity due to the nonlinear terms. Finally, we also present both a posteriori and a priori error bounds associated with our reduced order model. We compare the performance of the spatial and time-windowing reduced order modeling approaches in terms of accuracy and speed-up with respect to the corresponding full order model for several numerical examples, namely Sedov blast, Gresho vortices, Taylor-Green vortices, and triple-point problems.

97 MATHEMATICS AND COMPUTING↗

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗

On the study of hydrodynamic instabilities in the presence of background magnetic fields in high-energy-density plasmas

Blast-wave-driven hydrodynamic instabilities are studied in the presence of a background B-field through experiments and simulations in the high-energy-density (HED) physics regime. In experiments conducted at the Laboratoire pour l’utilisation des lasers intenses (LULI), a laser-driven shock-tube platform was used to generate a hydrodynamically unstable interface with a prescribed sinusoidal surface perturbation, and short-pulse x-ray radiography was used to characterize the instability growth with and without a 10-T B-field. The LULI experiments were modeled in FLASH using resistive and ideal magnetohydrodynamics (MHD), and comparing the experiments and simulations suggests that the Spitzer model implemented in FLASH is necessary and sufficient for modeling these planar systems. These results suggest insufficient amplification of the seed B-field, due to resistive diffusion, to alter the hydrodynamic behavior. Although the ideal-MHD simulations did not represent the experiments accurately, they suggest that similar HED systems with dynamic plasma-β (=2μ0ρv2/B2) values of less than ∼100 can reduce the growth of blast-wave-driven Rayleigh–Taylor instabilities. These findings validate the resistive-MHD FLASH modeling that is being used to design future experiments for studying B-field effects in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of a high energy density experiment to measure the suppression of hydrodynamic instability in an applied magnetic field

A transverse magnetic field can suppress the hydrodynamic instability growth of an unstable plasma interface. This effect, of interest to inertial confinement fusion and astrophysics, has mostly been studied via simulation. Here, we present the design of an experiment at the National Ignition Facility (NIF) to demonstrate this effect in a laboratory. Simulations indicate that the timescale for the diffusion of the magnetic field across the mixing region should be at least comparable to the timescale of the instability growth in order to have a measurable suppression effect. This motivates the use of lower density target materials than usual high energy density (HED) hydrodynamics experiments to permit faster hydrodynamics and higher plasma conductivities (through higher temperature), for a given laser drive and magnetic field. We discuss a target design for creating a Rayleigh–Taylor unstable HED plasma interface that uses 320 mg/cc iodine-doped carbon foam as the heavy material, 20 mg/cc carbon foam as the light material, and a 6 lm amplitude, a 120 lm wavelength ripple machined at the interface, which shows a measurable suppression effect with a nominal NIF drive and a 30 T magnetic field (the present facility limit). Here, models indicate lower density foams that may display even larger suppression effects, as the Hohlraum drive also radiatively preheats the foam to permit even higher temperatures and, hence, higher conductivities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The development of a high-resolution Eulerian radiation-hydrodynamics simulation capability for laser-driven Hohlraums

Hohlraums are hollow cylindrical cavities with high-Z material walls used to convert laser energy into uniform x-ray radiation drives for inertial confinement fusion capsule implosions and high energy density physics experiments. Credible computational modeling of hohlraums requires detailed modeling and coupling of laser physics, hydrodynamics, radiation transport, heat transport, and atomic physics. We report on improvements to Los Alamos National Laboratory's xRAGE radiation-hydrodynamics code in order to enable hohlraum modeling. xRAGE's Eulerian hydrodynamics and adaptive mesh refinement make it uniquely well suited to study the impacts of multiscale features in hohlraums. In order to provide confidence in this new modeling capability, we demonstrate xRAGE's ability to produce reasonable agreement with data from several benchmark hohlraum experiments. We also use xRAGE to perform integrated simulations of a recent layered high density carbon capsule implosion on the National Ignition Facility in order to evaluate the potential impacts of the capsule support tent, mixed cell conductivity methodologies, plasma transport, and cross-beam energy transfer (XBT). We find that XBT, seeded by plasma flows in the laser entrance hole (LEH), causes a slight decrease in energy coupling to the capsule and that all of these impact the symmetry of the x-ray drive such that they have an appreciable impact on the capsule implosion shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗