Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sensitivity of He Flames in X-Ray Bursts to Nuclear Physics

Through the use of axisymmetric 2D hydrodynamic simulations, we further investigate laterally propagating flames in X-ray bursts (XRBs). Our aim is to understand the sensitivity of a propagating helium flame to different nuclear physics. Using the Castro simulation code, we confirm the phenomenon of enhanced energy generation shortly after a flame is established by adding 12 C(p, γ) 13 N(α, p) 16 O to the network, in agreement with the past literature. This sudden outburst of energy leads to a short accelerating phase, causing a drastic alteration in the overall dynamics of the flame in XRBs. Furthermore, we investigate the influence of different plasma screening routines on the propagation of the XRB flame. We finally examine the performance of simplified spectral deferred correction, a novel approach to hydrodynamics and reaction coupling incorporated in Castro, as an alternative to operator splitting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Dynamical chemistry: non-equilibrium effective actions for reactive fluids

Abstract We present two approaches for describing chemical reactions taking place in fluid phase. The first method mirrors the usual derivation of the hydrodynamic equations of motion by relating conserved—or to account for chemical reactions, non-conserved—currents to local-equilibrium parameters. The second method involves a higher-brow approach in which we attack the same problem from the perspective of non-equilibrium effective field theory (EFT). Non-equilibrium effective actions are defined using the in–in formalism on the Schwinger–Keldysh contour and are therefore capable of describing thermal fluctuations and dissipation as well as quantum effects. The non-equilibrium EFT approach is especially powerful as all terms in the action are fully specified by the symmetries of the system; in particular the second law of thermodynamics does not need to be included by hand, but is instead derived from the action itself. We find that the equations of motion generated by both methods agree, but the EFT approach yields certain advantages. To demonstrate some of these advantages we construct a quadratic action that is valid to very small distance scales—much smaller than the scales at which ordinary hydrodynamic theories break down. Such an action captures the full thermodynamic and quantum behavior of reactions and diffusion at quadratic order. Finally, taking the low-frequency and low-wavenumber limit, we reproduce the linearized version of the well-known reaction–diffusion equations as a final coherence check.

Mechanics↗
DOE OSTI2021

Meshfree simulation and experimental validation of extreme thermomechanical conditions in friction stir extrusion

Abstract Friction stir extrusion (FSE) is a novel solid-phase processing technique that consolidates and extrudes metal powders, flakes, chips, or billets into high-performance parts by plastic deformation, which has the potential to save substantial processing time and energy. Currently, most studies on FSE are experimental and only a few numerical models have been developed to explain and predict the complex physics of the process. In this work, a meshfree simulation framework based on smoothed particle hydrodynamics (SPH) was developed for FSE. Unlike traditional grid-based methods, SPH is a Lagrangian particle-based method that can handle severe material deformations, capture moving interfaces and surfaces, and monitor the field variable histories explicitly without complicated tracking schemes. These aspects of SPH make it attractive for the FSE process, where in situ evolution of field variables is difficult to observe experimentally. To this end, a 3-D, fully thermomechanically coupled SPH model was developed to simulate the FSE of aluminum wires. The developed model was thoroughly validated by comparing the numerically predicted material flow, strain, temperature history, and extrusion force with experimental results for a certain set of process parameters. The validated SPH model can serve as an effective tool to predict and better understand the extreme thermomechanical conditions during the FSE process.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Intersection-distribution-based remapping between arbitrary meshes for staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics

In this paper, we present a new intersection-distribution-based remapping method between arbitrary polygonal meshes for indirect staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics. All cell-centered material quantities are conservatively remapped using intersections between the Lagrangian (old, source) mesh and the rezoned (new, target) mesh. The new nodal masses are obtained by conservative distribution of all material masses in each new cell to the cell's corners and then collecting those corner masses at new nodes. This distribution is done using a local constrained optimization approach for each cell in the new mesh. In order to remap nodal momentum we first define cell-centered momentum for each cell in the old mesh, conservatively remap this to the new mesh and then conservatively distribute the new zonal momentum to each cell's bounding nodes, again using local constrained optimization. Our method also conserves total energy by applying a new nodal kinetic energy correction that relies on a process similar to that used for remapping nodal mass and momentum. Cell-centered kinetic energy is computed, conservatively remapped and then distributed to nodes. The discrepancy between this conservatively remapped and actual nodal kinetic energy is then conservatively distributed to the internal energies of the materials in the cells surrounding each node. Unlike conventional cell-based corrections of this type, this new nodal kinetic energy correction has not been observed to drive material internal energy negative in any of our testing. Unlike flux based remapping, our new intersection-distribution method can be applied to remapping between source and target meshes that are arbitrarily different, which provides superior flexibility in the rezoning strategy. Our method is accurate, essentially conservative and essentially bounds preserving.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Measurement of the small-scale 3D Lyman- α forest power spectrum

Small-scale correlations measured in the Lyman-α (Lyα) forest encode information about the intergalactic medium and the primordial matter power spectrum. In this article, we present and implement a simple method to measure the 3-dimensional power spectrum, P 3D , of the Lyα forest at wavenumbers k corresponding to small, ~ Mpc scales. In order to estimate P 3D from sparsely and unevenly distributed data samples, we rely on averaging 1-dimensional Fourier Transforms, as previously carried out to estimate the 1-dimensional power spectrum of the Lyα forest, P 1D . Further, this methodology exhibits a very low computational cost. We confirm the validity of this approach through its application to Nyx cosmological hydrodynamical simulations. Subsequently, we apply our method to the eBOSS DR16 Lyα forest sample, providing as a proof of principle, a first P 3D measurement averaged over two redshift bins z = 2.2 and z = 2.4. This work highlights the potential for forthcoming P 3D measurements, from upcoming large spectroscopic surveys, to untangle degeneracies in the cosmological interpretation of P 1D .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Design optimization for Richtmyer–Meshkov instability suppression at shock-compressed material interfaces

We report the Richtmyer–Meshkov instability (RMI) is a phenomenon that occurs at the interface of two substances of different densities due to an impulsive acceleration, such as a shock wave passing through this interface. Under these conditions, the instability can be seen as interface perturbations begin to grow into narrow jets or spikes of one substance that propagate into the other. In some cases, this interface may involve an elastic–plastic material, which can play a significant role in the development and behavior of the RMI. The ability to effectively control RMI jetting and spike growth is one major limiting factor in technological challenges, such as inertial confinement fusion, that involve using high-pressure shock waves to implode a fuel target. The propagation of RMI growth can lead to increased asymmetry in this implosion process and significantly reduce the obtained energy yield. We use hydrodynamics simulations of impactor shock-compression experiments and methods based in design optimization to suppress RMI spike growth by altering the geometry and other properties of a shock-compressed elastic–plastic material target that shares an interface with atmospheric air. These hydrodynamics simulations use an arbitrary Lagrangian–Eulerian method with a high-order finite element approach. Our results demonstrate that RMI suppression can be achieved by intentionally creating a separate upstream interface instability to counteract the growth of long narrow RMI spikes at an interface with initial perturbations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Review of recent WEC-Sim (v6.1) advanced features

WEC-Sim (Wave Energy Converter SIMulator) is an opensource software for modeling the motions, loads and power generation of wave energy converters. WEC-Sim performs simulations in the time domain using hydrodynamic coefficients calculated by boundary element method (BEM) frequency-domain potential flow solvers such as WAMIT, NEMOH, Capytaine, or Ansys AQWA. WEC-Sim development is ongoing, including various features, applications, and example cases used to demonstrate potential use-cases to meet the needs of the growing marine energy industry. Through input from a broad user base and an extensive team of developers and collaborators, new features of WEC-Sim are developed to expand the software’s use cases and improve overall functionality. Three new WEC-Sim features highlighted in this paper include updating WEC-Sim to be compatible with MoorDyn Version 2, incorporation of second order excitation loads (quadratic transfer functions) and allowing for dynamically changing hydrodynamics.

13 HYDRO ENERGY↗
DOE OSTI2022

A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors

Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers ($\mathscr{M}$ son ≲ 10 –2 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant–Friedrichs–Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume SEVEN-LEAGUE HYDRO (SLH) code, and it makes use of a low-Mach version of the five-wave Harten–Lax–van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit–explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at $\mathscr{M}$ son ~ 10 –3 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A polar direct drive liquid deuterium–tritium wetted foam target concept for inertial confinement fusion

We propose a new approach to inertial confinement fusion (ICF) that could potentially lead to ignition and propagating thermonuclear burn at the National Ignition Facility (NIF). The proposal is based upon a combination of two concepts, referred to as polar direct drive and liquid deuterium–tritium wetted foam capsules. With this new concept, 2D radiation hydrodynamic simulations indicate that ICF ignition and propagating thermonuclear burn are possible with the laser power and energy capabilities available today on the NIF.

3D printing↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

DOE WPTO FY2021 Seedling End Report: Coupling Soil-Foundation Models to MHK Device Dynamic Models to Tighten the Design Envelope

The purpose of this Seedling project is to couple a marine renewable energy (MRE) dynamics simulation software with the soil-foundation models in the OC6 Phase II project [Bergua et al., 2021] and evaluate the software’s performance. This is a first step to accurately evaluating soil-foundation impacts on other types of MRE, like wave or current energy converters (WECs, CECs). OC6 Phase II compares offshore wind turbine (OWT) simulations using several different soil-foundation models to identify and fill key gaps in soil-foundation analyses. WEC-Sim was chosen to model the OC6 Phase II offshore wind turbine and various load cases because of its adaptability, accuracy of hydrodynamic loads, and ability to apply an arbitrary wind loading. Of the four methods used in OC6, the apparent fixity soil-foundation method was coupled with WEC-Sim. Technical challenges with flexible hydrodynamic bodies, added mass and external function libraries inhibited the ability to compare the WEC-Sim results to other OC6 participants. These challenges required that the WEC-Sim model of the OC6 OWT use a combination of rigid and flexible bodies to ensure a numerically stable solution. The rigid monopile creates a more stiff system and causes smaller amplitude motion under hydrodynamic loading and higher dominant frequency of motion under wind loading. These discrepancies are expected based on the increased stiffness of the WEC-Sim case.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Tracking the environmental impacts of ecological engineering on coastal wetlands with numerical modeling and remote sensing

Coastal wetlands are the most valuable ecosystems on the earth but facing severe degradation and losses owing to climate change and anthropogenic activities. Many ecological engineering projects (EEP) have been conducted to mitigate the degradation of coastal wetlands. However, the geomorphological impacts of EEP on coastal wetlands have not been well documented. In this study, a method employed a process-based hydrodynamic model and remote sensing (RS) was developed to evaluate the impacts of EEP on the geomorphological change of a prototype Ramsar site. Results demonstrated that RS has great potential in improving the quality of bathymetry data for the numerical model with a decrease of RMSE from 0.52 m to 0.3 m. It also showed good capacity in trend detection of geomorphological change spatially. Results showed the Chongming Dongtan wetland experienced erosion with an annual rate of -0.035 m/yr from 2013 to 2016 after the implementation of EEP. The deposition rate changed significantly in the area within 200 m of the EEP. It is found that the EEP modified the composition of vegetation, sediment transportation, as well as substrate stability, affecting the geomorphological change of coastal wetlands. It is suggested that the EEP with moderate anthropogenic disturbance is a direct and effective way to recover the coastal habitats for waterbirds. However, the modification of the coastal wetland ecosystem by EEP will lead to the potential vulnerability to global climate change. Therefore, how to mitigate the advantages and disadvantages of the EEP is needed to be further studied to find a more sustainable way for coastal management.

geomorphological change, coastal wetlands, Ecologi↗
DOE OSTI2025

Grain Boundary-Limited Thermal Transport in Suspended Thin Graphite across an Unexplored Thickness Regime

Here, we present systematic thermal conductivity (κ) measurements of suspended thin graphite ribbons, 234–527 nm thick, using a four-probe 3ω method. Unlike recent reports of phonon hydrodynamics and exceptionally high κ in micrometer-thick graphite ( Science, 2020), we observe significantly lower κ and no signatures of collective phonon flow in this intermediate thickness regime. Instead, our measured κ lies between few-layer graphene and bulk graphite. These results agree with a first-principles-informed Peierls–Boltzmann transport model with spatially resolved Monte Carlo sampling. Additionally, the temperature for the peak κ shifts lower with increasing thickness, due to the interplay of phonon-boundary and phonon-isotope scattering. Incorporating grain boundary scattering into simulations is necessary to replicate the experimental trends. These findings delineate the boundary between ballistic, hydrodynamic, and diffusive transport regimes in graphite and underscore the dominant role of disorder and geometry in phonon transport in quasi-two-dimensional materials, offering insights for nanoscale thermal management.

Boltzmann transport model↗
DOE OSTI2024

Filaments in and between galaxy clusters at low and mid-frequency with the SKA

Context. Understanding the magnetised Universe is a major challenge in modern astrophysics, and cosmic magnetism has been acknowledged as one of the key scientific drivers of the most ambitious radio instrument ever planned, the Square Kilometre Array. Aims. With this work, we aim to investigate the potential of the Square Kilometre Array and its precursors and pathfinders in the study of magnetic fields in galaxy clusters and filaments through diffuse synchrotron radio emission. Galaxy clusters and filaments of the cosmic web are indeed unique laboratories in which to investigate turbulent fluid motions and large-scale magnetic fields in action, and much of what is known about magnetic fields in galaxy clusters comes from sensitive radio observations. Methods. Based on cosmological magneto-hydrodynamic simulations, we predict radio properties (total intensity and polarisation) of a pair of galaxy clusters connected by a cosmic-web filament. Results. We use our theoretical expectations to explore the potential of polarimetric observations to study large-scale structure magnetic fields in the frequency ranges 50–350 MHz and 950–1760 MHz. We also present predictions for galaxy cluster polarimetric observations with the Square Kilometre Array precursors and pathfinders, such as the LOw frequency ARray 2.0 and the MeerKAT+ telescope. Conclusions. Our findings point out that polarisation observations are particularly powerful for the study of large-scale magnetic fields, since they are not significantly affected by confusion noise. The unprecedented sensitivity and spatial resolution of the intermediatefrequency radio telescopes make them the favourite instruments for the study of these sources through polarimetric data, potentially allowing us to understand if the energy density of relativistic electrons is in equipartition with the magnetic field or rather coupled with the thermal gas density. Our results show that low-frequency instruments also represent a precious tool to study diffuse synchrotron emission in total intensity and polarisation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Towards a quantum fluid theory of correlated many-fermion systems from first principles

Correlated many-fermion systems emerge in a broad range of phenomena in warm dense matter, plasmonics, and ultracold atoms. Quantum hydrodynamics (QHD) complements first-principles methods for many-fermion systems at larger scales. We illustrate the failure of the standard Bohm potential central to QHD for strong perturbations when the density perturbation is larger than about 10^{-3} 10 − 3 of the mean density. We then extend QHD to this regime via the many-fermion Bohm potential from first-principles. This may lead to more accurate QHD simulations beyond their common application domain in the presence of strong perturbations at scales unattainable with first-principles methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗