Search NASA⌕ Search

SEARCH · Search NASA

Results for “Eulerian code”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗
DOE OSTI2017

High-resolution modeling of indirectly driven high-convergence layered inertial confinement fusion capsule implosions

In this study, we present the results of high-resolution simulations of the implosion of high-convergence layered indirect-drive inertial confinement fusion capsules of the type fielded on the National Ignition Facility using the xRAGE radiation-hydrodynamics code. In order to evaluate the suitability of xRAGE to model such experiments, we benchmark simulation results against available experimental data, including shock-timing, shock-velocity, and shell trajectory data, as well as hydrodynamic instability growth rates. We discuss the code improvements that were necessary in order to achieve favorable comparisons with these data. Due to its use of adaptive mesh refinement and Eulerian hydrodynamics, xRAGE is particularly well suited for high-resolution study of multi-scale engineering features such as the capsule support tent and fill tube, which are known to impact the performance of high-convergence capsule implosions. High-resolution two-dimensional (2D) simulations including accurate and well-resolved models for the capsule fill tube, support tent, drive asymmetry, and capsule surface roughness are presented. These asymmetry seeds are isolated in order to study their relative importance and the resolution of the simulations enables the observation of details that have not been previously reported. We analyze simulation results to determine how the different asymmetries affect hotspot reactivity, confinement, and confinement time and how these combine to degrade yield. Yield degradation associated with the tent occurs largely through decreased reactivity due to the escape of hot fuel mass from the hotspot. Drive asymmetries and the fill tube, however, degrade yield primarily via burn truncation, as associated instability growth accelerates the disassembly of the hotspot. Finally, modeling all of these asymmetries together in 2D leads to improved agreement with experiment but falls short of explaining the experimentally observed yield degradation, consistent with previous 2D simulations of such capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Nyx

Nyx is a highly parallel, adaptive mesh, finite-volume N-body compressible hydrodynamics solver for cosmological simulations. Nyx code solves equations of compressible hydrodynamics on an adaptive grid hierarchy coupled with an N-body treatment of dark matter. The gas dynamics in Nyx uses a finite volume methodology on a set of 3-D Eulerian grids; dark matter is represented as discrete particles moving under the influence of gravity. Particles are evolved via a particle-mesh method, using Cloud-in-Cell deposition/interpolation scheme. Both baryonic and dark matter contribute to the gravitational field. In addition, Nyx includes physics needed to accurately model the intergalactic medium: in optically thin limit and assuming ionization equilibrium, the code calculates heating and cooling processes of the primordial-composition gas in an ionizing ultraviolet background radiation field. Additional physics capabilities are under development.

Sexton, Jean [Lawrence Berkeley National Lab. (LBN↗
NASA NTRS1988

Formation of polar stratospheric clouds simulated in a two dimensional model of the atmosphere

A microphysics code has been implemented in a two dimensional model of the atmosphere to study formation of polar stratospheric clouds containing HCl or HNO3. The model range is from pole to pole in latitude and from the ground to about 20 km in altitude. Resolution in latitude is 10 deg and about 0.8 km in altitude. This is an Eulerian model with prescribed eddy diffusion coefficients and the circulation obtained from observations. The chemistry of the model follows the family approach for NO(x), Cl(x) and HO(x) while the ozone is fixed and changed seasonally. The aerosol code is based on an assigned population of condensation, coagulation and sedimentation. Aerosol growth is simulated in nine different size bins ranging between 0.01u and 2.56u. The model has been built to study aerosol layers formation in the upper troposphere and lower stratosphere and has been validated for sulfate aerosol resulting from a rather complex sulfur chemistry.

Visconti, Guido↗
DOE OSTI2016

The NETL MFiX Suite of multiphase flow models: A brief review and recent applications of MFiX-TFM to fossil energy technologies

Here, the MFiX suite of multiphase computational fluid dynamics (CFD) codes is being developed at U.S. Department of Energy's National Energy Technology Laboratory (NETL). It includes several different approaches to multiphase simulation: MFiX-TFM, a two-fluid (Eulerian–Eulerian) model; MFiX-DEM, an Eulerian fluid model with a Lagrangian Discrete Element Model for the solids phase; and MFiX-PIC, Eulerian fluid model with Lagrangian particle ‘parcels’ representing particle groups. These models are undergoing continuous development and application, with verification, validation, and uncertainty quantification (VV&UQ) as integrated activities. After a brief summary of recent progress in the verification, validation and uncertainty quantification (VV&UQ), this article highlights two recent accomplishments in the application of MFiX-TFM to fossil energy technology development. First, recent application of MFiX to the pilot-scale KBR TRIG™ Transport Gasifier located at DOE's National Carbon Capture Center (NCCC) is described. Gasifier performance over a range of operating conditions was modeled and compared to NCCC operational data to validate the ability of the model to predict parametric behavior. Second, comparison of code predictions at a detailed fundamental scale is presented studying solid sorbents for the post-combustion capture of CO 2 from flue gas. Specifically designed NETL experiments are being used to validate hydrodynamics and chemical kinetics for the sorbent-based carbon capture process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

MOSCATO Development and Integration in Fiscal Year 2024

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy evolution for molten salt equipment. In FY24, significant developments were made to the code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: 1. Implementation of tritium transport capabilities and validation with experimental data: To enable modeling of tritium and other fission gases within MSRs, we implemented gas transport within MOSCATO via inclusion of couple mass transport equations within the salt and structural alloys. Comparisons to experimental data from literature showed good agreement with respect to tritium release rates. 2. Preliminary implementation of two-phase flow models in MOSCATO: To model tritium and other gases above their solubility limits, we implemented preliminary two-phase flow models within MOSCATO to account for bubble transport. The first model adopted was the Level-Set approach, which can handle the high void fraction regime, but with a requirement for high mesh resolution thus high computational expense. In this report, we present a verification of the Level-Set method using a simple benchmark case. We also performed a demonstration of the code as applied to an experimental case involving cover gas flow through salt in an experimental vessel. The second model adopted was the Eulerian-Eulerian dispersed flow model, which is computationally cheaper but limited to low void fraction regimes, such as bubbly flow. Validation and verification have not yet been performed for the Eulerian-Eulerian approach, but a preliminary implementation was completed. 3. Validation with static corrosion experiments: Static corrosion experimental data for stainless steel coupons within molten salts was used to further validate the corrosion model in MOSCATO. To do so, we leveraged the existing models in MOSCATO and simulated the sample mass loss and mass gain phenomena. Several ion species, including Cr 2+ , Fe 2+ and H + , were simulated in salt using the PNP solver, while Cr 0 and Fe 0 were simulated with a diffusion solver in stainless steel. The mass loss of the samples was compared with experimental data, and good agreement was achieved. These combined activities served to further expand the capabilities of MOSCATO and make it more generally applicable to the full range of phenomena that can control chemistry and corrosion in molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2007

Simulating Space Capsule Water Landing with Explicit Finite Element Method

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Wang, John T.↗
DOE OSTI2017

Computation of Air Discharge into a Water Pool by Using Volume-of-Fluid Method

Suppression tanks, in-containment refueling water storage tanks (IRWSTs) and filtered containment venting systems (FCVSs) of nuclear power plants use water pools to accommodate high-pressure and high-temperature discharge of non-condensable gases and/or steam. Dynamic behavior of gas bubbles inside a water pool determines pressure and temperature loading to the system and thus we need validated method for applications. The underlying physics relevant to these problems is still a difficult issue due to the complicated nature of gas- liquid interfacial transfers of mass, momentum, and energy in a two-phase flow. Therefore, notwithstanding many studies performed over the past two or three decades, two-fluid theory is still a formidable task and direct numerical simulation (DNS) needs huge computing cost. Pellegrini et al. analyzed steam condensation and chugging in the water pool of the PPOOLEX STB-24 test by using STAR-CCM+ code implemented with two-phase interface model based on the Rayleigh-Taylor instability theory (RTI) which triggers steam depressurization. For the solution of the fine interface, they used Eulerian-Eulerian model for the fluid domain since the volume-of-fluid (VOF) model requires larger mesh points. The RTI resulted in complete steam condensation consistent with the experiment. The fluid tank was open to the atmosphere. They reported the volume fraction distributions in the pool, not the pool temperature or pressure history. For non-condensing gas discharges, however, the VOF is still a good alternative to the Eulerian two-fluid or the DNS approaches. The spatial distribution of each phase is defined in terms of volume fractions and the cell physical properties are obtained from linear combinations of those of constituent phases weighted by each volume fraction. Therefore, conditions where thermal and phase change effects become important are less likely to be simulated realistically. The VOF model is suited for the simulation of flows where each phase constitutes a large structure with a relatively small total contact area between phases. Present work is also interested in this VOF method for application to air discharge into a water pool which occurs in the early operating stages of the said facilities. Our purpose is to validate the VOF method against an existing air discharge experiment of PPOOLEX by using a commercial CFD code, STAR-CCM+ version 10.04, in terms of macroscopic flow properties. Effect of turbulence is also discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2003

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗
NASA NTRS1989

Analysis of rotary engine combustion processes based on unsteady, three-dimensional computations

A new computer code was developed for predicting the turbulent, and chemically reacting flows with sprays occurring inside of a stratified charge rotary engine. The solution procedure is based on an Eulerian Lagrangian approach where the unsteady, 3-D Navier-Stokes equations for a perfect gas mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite volume, Steger-Warming flux vector splitting scheme, and the liquid phase equations are solved in Lagrangian coordinates. Both the details of the numerical algorithm and the finite difference predictions of the combustor flow field during the opening of exhaust and/or intake, and also during fuel vaporization and combustion, are presented.

Raju, M. S.↗
NASA NTRS1990

Analysis of rotary engine combustion processes based on unsteady, three-dimensional computations

A new computer code was developed for predicting the turbulent and chemically reacting flows with sprays occurring inside of a stratified charge rotary engine. The solution procedure is based on an Eulerian Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite volume, Steger-Warming flux vector splitting scheme, and the liquid phase equations are solved in Lagrangian coordinates. Both the details of the numerical algorithm and the finite difference predictions of the combustor flow field during the opening of exhaust and/or intake, and also during fuel vaporization and combustion, are presented.

Raju, M. S.↗
NASA NTRS1998

Three-Dimensional Hydrodynamic Simulations of the Postimpact Proto-Earth

Astrophysical fluid configurations are susceptible to a variety of nonaxisymmetric instabilities under the combined effects of rotation, self-gravity, and thermal pressure. When strong enough, they can induce rapid transport of mass and angular momentum. Our own previous studies of nonaxisymmetric instabilities in model protostars and protostellar disks show that significant transport can occur on orbital timescales and that material can be ejected to large distances. In this contribution, we present three-dimensional simulations of the circumterrestrial debris belt that may have resulted from a giant impact. Our three-dimensional hydrodynamics code with self-gravity and artificial viscosity is fully second-order in space and time; the equations of hydrodynamics and the Poisson equation are solved on an Eulerian cylindrical grid. In the preliminary calculations presented here, we use a simplified EOS where the central proto-Earth is treated as an n = 1/2 polytropic fluid, surrounded by a more compressible, rapidly rotating, fluid disk that represents silicate vapor. Our initial disk parameters are generated from the endstate data of recent smoothed particle hydrodynamics giant-impact calculations. Ultimately, we wish to detennine under what conditions nonaxisymmetric instabilities grow in the postimpact disk and whether they facilitate the transport of material outside the proto-Earth's Roche Limit, leading to the formation of the Moon. In future work, we hope to include a more realistic EOS and the consequences of heating, cooling, and phase transitions.

Pickett, B. K.↗
NASA NTRS1992

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗
DOE OSTI2018

Continuum kinetic modelling of cross-separatrix plasma transport in a tokamak edge including self-consistent electric fields

Eulerian kinetic calculations are presented here for the axisymmetric cross-separatrix plasma transport at the edge of a tokamak. The simulations are performed with a high-order finite volume code COGENT, which solves the long-wavelength limit of the full-f ion gyrokinetic equation including the non-linear Fokker–Plank collision model. Self-consistent 2D electrostatic potential variations are obtained from the quasi-neutrality (∇·j = 0) equation coupled to an isothermal fluid electron response. Finally, illustrative solutions presented for parameters characteristic of the DIII-D tokamak show qualitative agreement with experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validation and application of a multiphase CFD model for hydrodynamics, temperature field and RTD simulation in a pilot-scale biomass pyrolysis vapor phase upgrading reactor

Accurate prediction of transport phenomena is critical for VPU reactor design, optimization, and scale-up. The current study focused on the validation and application of a multiphase CFD model within an open-source code MFiX for hydrodynamics, temperature field, and residence time distribution (RTD) simulation in a non-reacting circulating fluidized bed riser for biomass pyrolysis vapor phase upgrading (VPU). First, an Eulerian-Eulerian approach three-dimensional CFD model was employed to simulate the pilot-scale VPU riser on the supercomputer Joule. Excellent quantitative agreement between experimental and simulated results was achieved for pressure drops and temperature field in a range of operating conditions. Then the validated multiphase CFD model was applied to predict gas and solid residence time distributions (RTDs) since prediction and analysis of RTD is an important tool to study the complex multiphase flow behavior and mixing inside chemical reactors. The predictions show that solid mean residence time is 3.5 times the gas residence time; the solid RTD is more sensitive to the process gas flow rate than the solids circulation rate.

09 BIOMASS FUELS↗
NASA NTRS2000

Unsteady Cascade Aerodynamic Response Using a Multiphysics Simulation Code

The multiphysics code Spectrum(TM) is applied to calculate the unsteady aerodynamic pressures of oscillating cascade of airfoils representing a blade row of a turbomachinery component. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena, in the present case being between fluids and structures. Interaction constraints are enforced in a fully coupled manner using the augmented-Lagrangian method. The arbitrary Lagrangian-Eulerian method is utilized to account for deformable fluid domains resulting from blade motions. Unsteady pressures are calculated for a cascade designated as the tenth standard, and undergoing plunging and pitching oscillations. The predicted unsteady pressures are compared with those obtained from an unsteady Euler co-de refer-red in the literature. The Spectrum(TM) code predictions showed good correlation for the cases considered.

Lawrence, C.↗
DOE OSTI2019

A Review of Physical and Numerical Approaches for the Study of Gas Stirring in Ladle Metallurgy

This article presents a review of the research into gas stirring in ladle metallurgy carried out over the past few decades. Herein, the physical modeling experiments are divided into four major areas: (1) mixing and homogenization in the ladle; (2) gas bubble formation, transformation, and interactions in the plume zone; (3) inclusion behavior at the steel–slag interface and in the molten steel; and (4) open eye formation. Several industrial trials have also been carried out to optimize gas stirring and open eye formation. Approaches for selecting criteria for scaling to guarantee flow similarity between industrial trials and physical modeling experiments are discussed. To describe the bubble behavior and two-phase plume structure, four main mathematical models have been used in different research fields: (1) the quasi-single-phase model, (2) the volume of fluid (VOF) model, (3) the Eulerian multiphase (E–E) model, and (4) the Eulerian–Lagrangian (E–L) model. In recent years, the E–E model has been used to predict gas stirring conditions in the ladle, and specific models in commercial packages, as well as research codes, have been developed gradually to describe the complex physical and chemical phenomena. Furthermore, the coupling of turbulence models with multiphase models is also discussed. For physical modeling, some general empirical rules have not been analyzed sufficiently. Based on a comparison with the available experimental results, it is found that the mathematical models focusing on the mass transfer phenomenon and inclusion behaviors at the steel-slag interface, vacuum degassing at the gas–liquid interface, dissolution rate of the solid alloy at the liquid–solid interface, and the combination of fluid dynamics and thermodynamics need to be improved further. To describe industrial conditions using mathematical methods and improve numerical modeling, the results of physical modeling experiments and industrial trials must offer satisfactory validations for the improvement of numerical modeling.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Tracking of raindrops in flow over an airfoil

The splashback that occurs when raindrops impact an airfoil results in an 'ejecta fog' of small droplets around the leading edge. Acceleration of these droplets by the air flow field is a momentum sink for the air flow and has been hypothesized to contribute to the degradation of airfoil performance in heavy rain. Presented here is a one-way coupled Eulerian-Lagrangian particle tracking scheme to evaluate droplet number densities and momentum sink terms around a NACA 64-210 airfoil section for three rainfall rates. A laminar air flow field is determined with a standard CFD code and input to the particle tracking algorithm. Raindrops are assumed to be non-interacting, non-deforming, non-evaporating, and non-spinning spheres and are tracked through the curvilinear grid used by the air flow code. A simple model is used to simulate impacts and the resulting splashback on the airfoil surface.

Valentine, James R.↗