Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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

Nyx: A Massively Parallel AMR Code for Computational Cosmology

Nyx is a highly parallel, adaptive mesh, finite-volume N-body compressible hydrodynamics solver for cosmological simulations. It has been used to simulate different cosmological scenarios with a recent focus on the intergalactic medium and Lyman alpha forest. Together, Nyx, the compressible astrophysical simulation code, Castro, and the low Mach number code MAESTROeX, make up the AMReX-Astrophysics Suite of open-source, adaptive mesh, performance-portable astrophysical simulation codes. Other examples of cosmological simulation research codes include Enzo, Enzo-P/Cello, RAMSES, ART, FLASH, Cholla, as well as Gadget, Gasoline, Arepo, Gizmo, and SWIFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

San Juan Islands Tidal Energy Characterization

The San Juan Islands are an archipelago with multiple tidal channels that produce strong tidal currents that could be harnessed for electricity generation. At present, there is limited electrical generation on the islands, with power primarily provided by a subsea cable connection with the mainland. As part of the C-MIST program, NOAA's CO-OPS program collected ADCP data at 46 stations in the San Juan Islands and adjacent waterways, from April to August of 2017. PNNL has recently developed a 3-D tidal hydrodynamic model for tidal energy resource characterization and assessment in the Salish Sea (Yang et al. 2021). Several tidal channels in the San Juan Islands were identified as top hotspots for potential tidal energy extraction. However, due to the complex geometry and inter-connected waterways, tidal currents around San Juan Islands also exhibit strong spatial and temporal asymmetry. Therefore, it is important to analyze the NOAA ADCP data to define metrics relevant to tidal energy, and further refine PNNL's tidal hydrodynamic model to accurately simulate currents in small tidal channels and account for sharp bathymetry gradients. The outcome of this technical assistance will allow developers to identify the promising opportunities for tidal energy in the San Juan Islands.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Physics-driven Explosions of Stripped High-mass Stars: Synthetic Light Curves and Spectra of Stripped-envelope Supernovae with Broad Light Curves

Stripped-envelope supernovae (SESNe) represent a significant fraction of core-collapse supernovae, arising from massive stars that have shed their hydrogen and, in some cases, helium envelopes. The origins and explosion mechanisms of SESNe remain a topic of active investigation. In this work, we employ radiative-transfer simulations to model the light curves and spectra of a set of explosions of single, solar-metallicity, massive Wolf–Rayet stars with ejecta masses ranging from 4 to 11 M ⊙ , which were computed from a turbulence-aided and neutrino-driven explosion mechanism. We analyze these synthetic observables to explore the impact of varying ejecta mass and helium content on observable features. We find that the light curve shape of these progenitors with high ejecta masses is consistent with observed SESNe with broad light curves but not the peak luminosities. The commonly used analytic formula based on rising bolometric light curves overestimates the ejecta mass of these high-initial-mass progenitor explosions by a factor of up to 2.6. In contrast, the calibrated method by Haynie et al., which relies on late-time decay tails, reduces uncertainties to an average of 20% within the calibrated ejecta mass range. Spectroscopically, the He I 1.083 μm line remains prominent even in models with as little as 0.02 M ⊙ of helium. However, the strength of the optical He I lines is not directly proportional to the helium mass but instead depends on a complex interplay of factors such as the 56 Ni distribution, composition, and radiation field. Thus, producing realistic helium features requires detailed radiative transfer simulations for each new hydrodynamic model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2022

Evaluation of 3D Simulation Capabilities within xRAGE using a Suite of Test Problems

This presentation describes a 3D verification and validation suite for xRAGE. Three verification test problems for hydrodynamics are presented: the Kidder ball problem, the Verney shell problem, and a 5-material compression problem. There is excellent agreement between 2D and 3D XRAGE simulation results for these problems, and between the xRAGE results and the benchmark solutions. In addition, two 3D ICF test problems are presented, based on an OMEGA direct drive capsule experiment, and on a NIF indirect drive capsule experiment. These two ICF test problems are used to evaluate two hydrodynamic methods within xRAGE. In addition to hydrodynamics the simulations also include three temperature plasma physics, thermal conduction and radiation diffusion. When simulating either capsule in 3D the newer unsplit hydrodynamic method in xRAGE produces more vorticity relative to the older default method. For the indirect drive capsule the 3D simulations are in reasonable agreement with the experimental values of ion temperature and neutron production. Unstable 3D vortex rings are present near the DT ice/carbon interface which enhance the confinement of the DT ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Feasibility of the LAMMPS SPH Package for Simulating Bubble Motion in Vibrating Containers

The Smoothed Particle Hydrodynamics (SPH) package within LAMMPS is explored as a possible tool for simulating the motion of bubbles in a vibrating liquid-filled container. As an initial test case, the unphysical but computationally less intense situation of a two-dimensional single bubble rising in a quiescent liquid under the influence of gravity is considered herein. Although physically plausible behavior was obtained under certain conditions, this behavior depends strongly on the system parameters. Moreover, the large density ratio between the liquid and bubble requires extremely small timesteps, which make the simulations undesirably computationally expensive. Ultimately, it was determined that this method is not feasible for providing quantitatively accurate results for the desired application.

42 ENGINEERING↗
DOE OSTI2022

Fluid structure interaction of a subaqueous pendulum: Analyzing the effect of wake correction via large eddy simulations

The dynamic behavior of a subaqueous cylindrical pendulum and corresponding flow dynamics are investigated. The objectives were twofold: (i) to examine whether the two-dimensional model equations sufficiently capture the three-dimensional dynamics and (ii) to investigate the emerging three-dimensional vortical flow structures. Large eddy simulations with two-way coupling fluid structure interaction were carried out using the immersed boundary method to simulate the motion of the pendulum and its interactions with the initially stagnant water. The resulting pendulum motion is compared against measured data obtained in a series of experimental tests to validate the simulation results and the model equations with and without wake corrections. An analysis of the flow vorticity revealed the development of a vortex ring during the first swing and the formation of tip vortices. The evolution of the vortex rings emerging from the motion of the subaqueous cylindrical pendulum was visualized using Q-criteria showing a reasonable agreement with vortical structures observed in the experiment using particle imaging velocimetry. The hydrodynamic moments acting on the simulated pendulum and the moments calculated from the model equations are analyzed. Using the insights from these numerical simulations, a modification of the wake correction is proposed to enhance the accuracy of the rate of decay and period. The transient effect of coherent flow on pendulum dynamics, especially the added mass effect, is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Hydrodynamic conditions in laser irradiated buried layer experiments

The calculation of open shell ionization level and radiative properties of materials in Non-Local Thermal Equilibrium (NLTE) is currently still a major challenge for any atomic model. The predictions of various NLTE atomic codes at these conditions still differ significantly. In recent years, a new buried layer platform was developed at the Lawrence Livermore National Laboratory and the Laboratory for Laser Energetics. This platform is used to measure ionization distribution and emission of open L-shell, mid-Z ions and open M-shell, high-Z ions at NLTE conditions that are relevant in many laser plasma applications. These experiments offer a unique chance for benchmarking the atomic models. In order to perform these experiments, a uniform well characterized plasma source is required. In this work, we present one-dimensional (1D) and two-dimensional simulations of the experimental platform. These simulations were used for both the design and the analysis of the experiments. The simulations demonstrate the different phases of hydrodynamic evolution of the target and identify the time windows in which uniform conditions can be achieved. A 1D expansion of the target was found to be adequate to describe the target's evolution for most of the experiment duration. The fast 1D simulations were compared with recent experimental results from the Omega laser facility. The sensitivity of the results to several modeling parameters such as the electron flux limiter and laser resonant absorption is reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission

Abstract In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor or other means of deflection proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's X-ray energy deposition into the asteroid and the resulting material ablation. The relevant physics in these simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed X-ray energy deposition model developed using the radiation-hydrodynamics code Kull that can be used to initiate a nuclear mitigation mission calculation. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials at a range of porosities, two different source spectra, and a broad range of radiation fluences, source durations, and angles of incidence. Using blowoff momentum as the primary metric, the model-initiated simulation results match the full radiation-hydrodynamics results to within 10%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2023

Exploring the parameter space of MagLIF implosions using similarity scaling. III. Rise-time scaling

Magnetized liner inertial fusion (MagLIF) is a z-pinch magneto-inertial-fusion concept studied at the Z Pulsed Power Facility of Sandia National Laboratories. Two important metrics characterizing current delivery to a z-pinch load are the peak current and the current-rise time, which is roughly the time interval to reach the peak current. It is known that, when driving a z-pinch load with a longer current-rise time, the performance of the z-pinch decreases. However, a theory to understand and quantify this effect is still lacking. Here, in this paper, we utilize a framework based on similarity scaling to analytically investigate the variations in the performance of MagLIF loads when varying the current-rise time, or equivalently, the implosion timescale. To maintain similarity between the implosions, we provide scaling prescriptions of experimental input parameters defining a MagLIF load and derive the expected scaling laws for stagnation conditions and for various performance metrics. We compare predictions of the theory to 2D numerical simulations using the radiation, magneto-hydrodynamic code hydra. For several metrics, we find acceptable agreement between the theory and simulations. Our results show that the voltage φload near the MagLIF load follows a weak scaling law φload ∝$t_{φ}^{-0.12}$ with respect to the characteristic timescale t φ of the voltage source, instead of the ideal φload ∝$t_{φ}^{-1}$ scaling. This occurs because the imploding height of the MagLIF load must increase to preserve end losses. As a consequence of the longer imploding liners, the required total laser preheat energy and delivered electric energy increase. Overall, this study helps understand the trade-offs of the MagLIF design space when considering future pulsed-power generators with shorter and longer current-rise times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Differentiable programming for online training of a neural artificial viscosity function within a staggered grid Lagrangian hydrodynamics scheme

Lagrangian methods to solve the inviscid Euler equations produce numerical oscillations near shock waves. A common approach to reducing these oscillations is to add artificial viscosity (AV) to the discrete equations. The AV term acts as a dissipative mechanism that attenuates oscillations by smearing the shock across a finite number of computational cells. However, AV introduces several control parameters that are not determined by the underlying physical model, and hence, in practice are tuned to the characteristics of a given problem. We seek to improve the standard quadratic-linear AV form by replacing it with a learned neural function that reduces oscillations relative to exact solutions of the Euler equations, resulting in a hybrid numerical-neural hydrodynamic solver. Because AV is an artificial construct that exists solely to improve the numerical properties of a hydrodynamic code, there is no offline ‘viscosity data’ against which a neural network can be trained before inserting into a numerical simulation, thus requiring online training. We achieve this via differentiable programming, i.e. end-to-end backpropagation or adjoint solution through both the neural and differential equation code, using automatic differentiation of the hybrid code in the Julia programming language to calculate the necessary loss function gradients. A novel offline pre-training step accelerates training by initializing the neural network to the default numerical AV scheme, which can be learned rapidly by space-filling sampling over the AV input space. We find that online training over early time steps of simulation is sufficient to learn a neural AV function that reduces numerical oscillations in long-term hydrodynamic shock simulations. These results offer an early proof-of-principle that online differentiable training of hybrid numerical schemes with novel neural network components can improve certain performance aspects existing in purely numerical schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling framework for evaluating the impacts of hydrodynamic pressure on hydrologic exchange fluxes and residence time for a large-scale river section over a long-term period

Quantifying hydrologic exchange fluxes (HEF) at the river and subsurface interface and their residence times (RT) in subsurface are important for managing the water quality and ecosystem health in dynamic river corridor systems. In this study, a modeling framework is developed for coupling the three-dimensional (3D) multi-phase surface, subsurface flow transport, and numerical tracer model for RT in a large-scale river section over a long period. The framework is utilized to evaluate the impacts of hydrodynamic pressure on HEFs and RT for a 30 km section of the Columbia River in Washington State over a three-year period. Based on comparisons between model simulations with and without considering hydrodynamic pressure, we found that hydrodynamic pressure increases the net HEFs by 7% with river gaining water from the subsurface domain, and leads to slight reduction of RT.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Radiation hydrodynamics modelling of kilonovae with SNEC

ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Bipolar planetary nebulae from outflow collimation by common envelope evolution

ABSTRACT The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and the dense toroidal-shaped ejecta left over from common envelope (CE) evolution. Here we use the 3D hydrodynamic adaptive mesh refinement (AMR) code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the smoothed particle hydrodynamics (SPH) code phantom serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta pre-PN outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post-CE winds may be aspherical, our models show that collimation via ‘inertial confinement’ will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock-focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC-driven bipolar lobes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A positivity-preserving and conservative intersection-distribution-based remapping algorithm for staggered ALE hydrodynamics on arbitrary meshes

In this study, we introduce new intersection-distribution-based remapping tools for indirect staggered arbitrary Lagrangian-Eulerian (ALE) simulations of multi-material shock hydrodynamics on arbitrary meshes. In addition to conserving momentum and total energy, the three-stage remapper proposed in this work preserves non-negativity of the internal energy. At the first stage, we construct slope-limited piecewise-linear reconstructions of all conserved quantities on zones of the source mesh and perform intersection-based remap to obtain bound-preserving zonal quantities on the target mesh. At the second stage, we define bound-preserving nodal quantities of the staggered ALE discretization as convex combinations of corner quantities. The nodal internal energy is corrected in a way which keeps it non-negative, while providing exact conservation of total energy. At the final stage, we distribute the non-negative nodal internal energy to corners, zones and materials using non-negative weights. Proofs of positivity preservation are provided for each stage. This work is a natural extension of our paper [14] in which a similar intersection-distribution-based remapping procedure was employed. The original version used a nodal kinetic energy fix which did not provably ensure positivity preservation for the zonal internal energy after the final distribution stage. The new algorithm cures this potential drawback by using ‘coordinated’ limiters for piecewise-linear reconstructions, remapping the internal energy to nodes and correcting it before redistribution. The effectiveness of the new nodal fix is illustrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗