Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic codes”

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

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

The Advanced Simulation and Computing (ASC) Transport project’s Jayenne Implicit Monte Carlo (IMC) transport library now includes an unstructured mesh capability and is available in a Lagrangian Applications Project (LAP) code. In this presentation, we discuss recent work by the LAP and Transport projects that provides an IMC transport capability for radiation hydrodynamics in the Lagrangian frame. Verification problems testing the new capabilities have been simulated and analyzed, i.e. Marshak wave, Mach 45, Su-Olsen, picket fence, and crooked pipe, both in one and two dimensions. We also present results on two stretch goal problems: a simplified COAX high energy density physics experiment and a supernova shock simulation. Finally, we identify current limitations and future work needed to bring a full capability to the user community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1983

Closed coronal structures. V - Gasdynamic models of flaring loops and comparison with SMM observations

A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.

Peres, G.↗
DOE OSTI2024

Annual IC Progress Report

The goal of this project is to computationally investigate the progenitors of astrophysical gamma-ray bursts (GRB), the most extreme explosions in the universe across length and time scales. We have developed a cutting-edge, code-bridging approach to exploring the physics behind these luminous events to test the hypothesis that GRBs with electromagnetic radio emission originate from massive star (MS) and black hole (BH) binary systems. Our multi-scale, multi-physics theoretical framework for GRBs is based on solving well-defined and testable hydrodynamics problems associated with the dynamical evolution of the system. We solve stages of the evolution in a tractable way with a suite of well-developed state-of-the-art computational tools that include the stellar evolution code MESA, the general relativistic magnetohydrodynamics (GRMHD) code Athena++, and the binary black hole (BBH) population synthesis code COSMIC, as well as a number of self-written python post-processing tools.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Leakage predictions for Rayleigh-step, helium-purge seals

Rayleigh-step, helium purge, annular shaft seals, studied for use in liquid oxygen turbopumps, generate a hydrodynamic force that enables the seal to follow shaft perturbations. Hence, smaller clearances can be used to reduce seal leakage. FLOWCAL, a computer code developed by Mechanical Technology Incorporated, predicts gas flow rate through an annular seal with an axial pressure gradient. Analysis of a 50-mm Rayleigh-step, helium-purge, annular seal showed the flow rate increased axial pressure gradient, downstream pressure, and eccentricity ratio. Increased inlet temperature reduced leakage. Predictions made at maximum and minimum clearances (due to centrifugal and thermal growths, machining tolerances and + or - 2 percent uncertainty in the clearance measurement) placed wide boundaries on expected flow rates. The widest boundaries were set by thermal growth conditions. Predicted flow rates for a 50-mm Rayleigh-step, helium-purge, annular seal underestimated measured flow rates by three to seven times. However, the analysis did accurately predict flow rates for choked gas flow through annular seals when compared to flow rates measured in two other independent studies.

Proctor, Margaret P.↗
NASA NTRS1991

The shock breakout in SN1987A modelled with the time-dependent radiative transfer

A fully implicit high-order scheme was developed for the independent multi-group radiative transfer coupled with implicit hydrodynamics. The application of this scheme to the SN1987A explosion shows that shortly after the shock breakout a dense shell forms. A gasdynamic code describing the time dependent radiation transport in the multi-group approximation with variable Eddington factors was developed. Results and other aspects of the project are presented.

Blinnikov, S. I.↗
NASA NTRS1994

Nonlinear evolution of tidally distorted accretion disks: Two-dimensional simulations

According to a previously published linear analysis, the tidal distortion of accretion disks in binary star systems produces a local hydrodynamic instability to m = 1 internal waves, which may have arbitrarily small wavelengths in the absence of viscosity. The instability is three-dimensional and approximately incrompressible. To explore the nonlinear outcome of this instability, we develop a shearing-sheet approximation on scales comparable to the disk thickness. The large-scale azimuthal variation of the disk is represented by varying the local metric with time (local orbital phase). The hydrodynamic equations can then be posed two-dimensionally on local meridional planes. We solve these equations with a second-order gasdynamical code based on the Total-Variation-Diminishing scheme. Our simulations confirm the predicted linear growth rate. The modes saturate chaotically at velocities scaling as the product of the linear growth rate and the wavelength. If the wavelength is small compared with the disk thickness, the modes remain nearly incompressible even when nonlinear. The two-dimensional power spectrum of velocities after saturation is roughly isotropic and extends over a broad range of scales in an approximately power-law fashion. We measure the heating rate associated with the nonlinear dissipation of the modes. The dissipation implies a secular torque on the disk and a return of angular momentum to the secondary star via the tidal potential. The estimated torque is somewhat larger than the tidal torque produced by maximal disk viscosity (alpha approximately 1). At least in these two-dimensional simulations, however, there is no significant angular momentum flux within the disk.

Ryu, Dongsu↗
DOE OSTI2024

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↗
DOE OSTI2026

AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework

We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗
NASA NTRSDate not supplied

3D Modeling of Solar-Type Stars to Characterize Stellar Jitter

Detection of Earth-mass planets requires measurements of radial velocity with extreme precision. To capture the tiny disturbances caused by a planet's motion, it is necessary to understand and characterize the host star’s turbulent dynamics in order to apply proper filtering to the observational data. We take advantage of current computational and technological capabilities to develop 3D realistic models of the stellar subsurface convection and atmospheres and thereby estimate the photospheric jitter. We have identified an initial set of target stars, obtained initial conditions using the MESA code, and obtained initial 3D radiative models of the stellar surfaces and atmospheres with a spatial resolution of 50km. We present initial 3D radiative hydrodynamic model results of the planet-hosting star HD209458.

SMD↗
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

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↗
NASA NTRS2005

Forecasting of Storm Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite element based codes. It affords a capability for simulating tidal circulation and storm surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate.

Valenti, Elizabeth↗
NASA NTRS2006

Forecasting of Storm-Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm-surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm-surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm-surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and, most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite-element-based codes. It affords a capability for simulating tidal circulation and storm-surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm-surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate. To keep the computational load manageable, this work was conducted using only the wind stress, calculated by using historical data from Hurricane Camille, as the input condition for the model. Hurricane storm-surge simulations were performed on an eight-node Linux computer cluster. Each node contained dual 2-GHz processors, 2GB of memory, and a 40GB hard drive. The digital elevation model (DEM) for this region was specified using a combination of Navy data (over water), NOAA data (for the coastline), and optimized Interferometric Synthetic Aperture Radar data (over land). This high-resolution topographical data of the Mississippi coastal region provided the ADCIRC model with improved input with which to calculate improved storm-surge forecasts.

Valenti, Elizabeth↗
DOE OSTI2023

A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations

ABSTRACT We perform two-dimensional (2D) numerical simulations of core convection for zero-age main-sequence stars covering a mass range from 3 to 20 M⊙. The simulations are performed with the fully compressible time-implicit code music. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 M⊙ models. The simulations cover hundreds to thousands of convective turnover time-scales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance (dov) scales with the stellar luminosity (L) and the convective core radius (rconv). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$, which is implemented in a one-dimensional stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses ≳10 M⊙. Our 2D simulations show the formation of a nearly adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent three-dimensional simulations of convection. The most luminous models show a growth in size with time of the nearly adiabatic layer. This growth seems to slow down as the upper edge of the nearly adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive time-scale in the overshooting layer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)

The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.

20 FOSSIL-FUELED POWER PLANTS↗