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 19 records

A description of multi-material hydrodynamics without reliance upon Pressure-Temperature-Equilibrium and its relevance to Eulerian codes [Slides]

This report contains "a description of the full multi-material Euler system without reliance on the assumption of the PTE closure; a discussion of the details that a numerical scheme discretizing the above system might need to have; a description of how that might play out in a PTE code after all; and, an attempt to lay out the maximally consistent set of equations."

97 MATHEMATICS AND COMPUTING↗

Plasma transport in an Eulerian AMR code

A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions to flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unit of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Immersed Finite Element Lagrangian-Eulerian Code-Coupling Framework

This report presents an assessment of immersed Eulerian-Lagrangian code-coupling techniques suitable for use in a broad range of mechanics applications. The coupling algorithm is based on an immersed finite element method that considers the Lagrangian and Eulerian overlap regions in the overall variational formulation. In this report the basic formulation details are presented followed by various aspects of the code-coupling algorithm using OpenIFEM as the Lagrangian/coupling framework. A series of representative test cases that illustrate the code-coupling algorithm are discussed. The current work provides an in-depth investigation into the immersed finite element method for the purposes of providing a rigorous coupling technique that is minimally invasive in the respective Eulerian and Lagrangian codes. A number of extensions to the base immersed finite element method have been examined. These extension include nodal and quadrature-based indicator functions, a Lagrangian volume-fraction calculation in regions of overlap, and the use of penalty constraints between the Lagrangian and Eulerian domains. A unique MPI-based coupling strategy that retains the independent MPI structure of each code has been demonstrated.

97 MATHEMATICS AND COMPUTING↗

On the Measurements of Numerical Viscosity and Resistivity in Eulerian MHD Codes

We propose a simple ansatz for estimating the value of the numerical resistivity and the numerical viscosity of any Eulerian MHD code. We test this ansatz with the help of simulations of the propagation of (magneto)sonic waves, Alfvén waves, and the tearing mode (TM) instability using the MHD code Aenus. By comparing the simulation results with analytical solutions of the resistive-viscous MHD equations and an empirical ansatz for the growth rate of TMs, we measure the numerical viscosity and resistivity of Aenus. The comparison shows that the fast magnetosonic speed and wavelength are the characteristic velocity and length, respectively, of the aforementioned (relatively simple) systems. We also determine the dependence of the numerical viscosity and resistivity on the time integration method, the spatial reconstruction scheme and (to a lesser extent) the Riemann solver employed in the simulations. From the measured results, we infer the numerical resolution (as a function of the spatial reconstruction method) required to properly resolve the growth and saturation level of the magnetic field amplified by the magnetorotational instability in the post-collapsed core of massive stars. Our results show that it is most advantageous to resort to ultra-high-order methods (e.g., the ninth-order monotonicity-preserving method) to tackle this problem properly, in particular, in three-dimensional simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Predictions for the drive capabilities of the RancheroS Flux Compression Generator into various load inductances using the Eulerian AMR Code Roxane

The Ranchero Magnetic Flux Compression Generator (FCG) has been used to create current pulses in the 10-­100 MA range for driving both “static” low inductance (0.5 nH) loads1 for generator demonstration purposes and high inductance (10-­20 nH) imploding liner loads2 for ultimate use in physics experiments at very high energy density. Simulations of the standard Ranchero generator have recently shown that it had a design issue that could lead to flux trapping in the generator, and a non-­ robust predictability in its use in high energy density experiments. A re-­examination of the design concept for the standard Ranchero generator, prompted by the possible appearance of an aneurism at the output glide plane, has led to a new generation of Ranchero generators designated the RancheroS (for swooped). This generator has removed the problematic output glide plane and replaced it with a region of constantly increasing diameter in the output end of the FCG cavity in which the armature is driven outward under the influence of an additional HE load not present in the original Ranchero. The resultant RancheroS generator, to be tested in LA43S-­L13, probably in early FY17, has a significantly increased initial inductance and may be able to drive a somewhat higher load inductance than the standard Ranchero. This report will use the Eulerian AMR code Roxane to study the ability of the new design to drive static loads, with a goal of providing a database corresponding to the load inductances for which the generator might be used and the anticipated peak currents such loads might produce in physics experiments. Such a database, combined with a simple analytic model of an ideal generator, where d(LI)/dt = 0, and supplemented by earlier estimates of losses in actual use of the standard Ranchero, scaled to estimate the increase in losses due to the longer current carrying perimeter in the RancheroS, can then be used to bound the expectations for the current drive one may apply to any load assembly in future experiments.

42 ENGINEERING↗

L4 Milestone Report for MixEOS 2016 experiments and simulations

Accurate simulations of fluid and plasma flows require accurate thermodynamic properties of the fluids or plasmas. This thermodynamic information is represented by the equations of state of the materials. For pure materials, the equations of state may be represented by analytical models for idealized circumstances, or by tabular means, such as the Sesame tables. However, when a computational cell has a mixture of two or more fluids, the equations of state are not well understood, particularly under the circumstances of high energy densities. This is a particularly difficult issue for Eulerian codes, wherein mixed cells arise simply due to the advection process. LANL Eulerian codes typically assume an “Amagat’s Law” (or Law of Partial Volumes) for the mixture in which the pressures and temperatures of fluids are at an equilibrium that is consistent with the fluids being segregated within the cell. However, for purposes of computing other EOS properties, e.g., bulk modulus, or sound speed, the fluids are considered to be fully “mixed”. LANL has also been investigating implementing instead “Dalton’s Law” in which the total pressure is considered to be the sum of the partial pressures within the cell. For ideal gases, these two laws give the same result. Other possibilities are nonpressure- temperature-equilibrated approaches in which the two fluids are not assumed to “mix” at all, and the EOS properties of the cell are computed from, say, volume-weighted averages of the individual fluid properties. The assumption of the EOS properties within a mixed cell can have a pronounced effect on the behavior of the cell, resulting in, for example, different shock speeds, pressures, temperatures and densities within the cell. There is no apparent consensus as to which approach is best under HED conditions, though we note that under typical atmospheric and near atmospheric conditions the differences may be slight.

36 MATERIALS SCIENCE↗

Numerical simulations of the ablative Rayleigh-Taylor instability in planar inertial-confinement-fusion targets using the FastRad3D code

The ablative Rayleigh-Taylor (RT) instability is a central issue in the performance of laser-accelerated inertial-confinement-fusion targets. Historically, the accurate numerical simulation of this instability has been a challenging task for many radiation hydrodynamics codes, particularly when it comes to capturing the ablatively stabilized region of the linear dispersion spectrum and modeling ab initio perturbations. Here, we present recent results from two-dimensional numerical simulations of the ablative RT instability in planar laser-ablated foils that were performed using the Eulerian code FastRad3D. Our study considers polystyrene, (cryogenic) deuterium-tritium, and beryllium target materials, quarter- and third-micron laser light, and low and high laser intensities. An initial single-mode surface perturbation is modeled in our simulations as a small modulation to the target mass density and the ablative RT growth-rate is calculated from the time history of areal-mass variations once the target reaches a steady-state acceleration. By performing a sequence of such simulations with different perturbation wavelengths, we generate a discrete dispersion spectrum for each of our examples and find that in all cases the linear RT growth-rate γ is well described by an expression of the form γ=α [kg/(1+ϵ kLm)]1/2−βkVa, where k is the perturbation wavenumber, g is the acceleration of the target, Lm is the minimum density scale-length, Va is the ablation velocity, and ϵ is either one or zero. The dimensionless coefficients α and β in the above formula depend on the particular target and laser parameters and are determined from two-dimensional simulation results through the use of a nonlinear curve-fitting procedure. While our findings are generally consistent with those of Betti et al. (Phys. Plasmas 5, 1446 (1998)), the ablative RT growth-rates predicted in this investigation are somewhat smaller than the values previously reported for the same target and laser parameters. It is speculated that differences in the equation-of-state and opacity models are largely responsible for the discrepancy. Resolution of this issue awaits the development of better experimental diagnostics capable of measuring small-wavelength (5–20 μm) perturbation growth due to the ablative RT instability in the linear regime.

Bates, J. W. (ORCID:0000000188087240)↗

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗

(U) Marker Lagrangian Simulations using Geometry from CAD Assembly Files [Slides]

The Eulerian code PAGOSA now contains an option for the forward time integration of any user selected material in the Lagrangian frame of reference. This option is accomplished using a Multifield version of Brackbill’s FLIP plus Sulsky’s MPM. FLIP+MPM is a Marker Lagrangian scheme that is not to be confused with well known Grid Lagrangian schemes. This Marker Lagrangian method allows materials to undergo completely arbitrary deformation, up to and including fracture and separation. Using this option makes PAGOSA into a mixed frame simulation tool such that fluids can be integrated in the Eulerian frame and solids integrated in the Lagrangian frame. Boundaries between the fields are said to be ’immersed’ in the computational domain, and their nature is formulated using Multifield Theory. Each Marker in a field of FLIP+MPM material tracks the thermodynamic state of a small (but finite) piece of mass. Initiation of the Markers using the geometry and material information from CAD Assembly files is the subject of the talk. The ease and speed with which the dynamics of the assembly can be studied is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Laser-driven flyer plate impact: Computational studies guided by experiments

We present a computational approach using a multimaterial, arbitrary Lagrangian–Eulerian code termed ALE3D to model the nanosecond/micrometer dynamics of the launch of 0.5–4.5 km/s laser-driven metal flyer plates and the impact with stationary targets of Pyrex and fused silica glasses, and Lexan and Plexiglas polymers, producing pressures in the target in the range of 5–20 GPa. The simulations are compared to experimental results where the flyer velocity profile and the velocity profile imparted to the target material were measured with high-speed velocimetry. The experimental flyer launch by a high-intensity pulsed laser is modeled by depositing heat into a thin vaporizable layer under the flyer plate. This model produces a flyer plate that has not been exposed to the laser pulse, allowing us to compare the properties of the real flyer to a simulated ideal flyer. The simulations of target impact are generally in good agreement with the experiment except at the highest impact velocities where the shock release process in the experiment is slower than that in the simulation. The cause of this disagreement is attributed to an inadequate description of the shock viscosity during the nanosecond unloading process.

Stekovic, Svjetlana (ORCID:0000000163878070)↗

No double detonations but core carbon ignitions in high-resolution, grid-based simulations of binary white dwarf mergers

We study the violent phase of the merger of massive binary white dwarf systems. Our aim is to characterize the conditions for explosive burning to occur, and identify a possible explosion mechanism of Type Ia supernovae. The primary components of our model systems are carbon–oxygen (C/O) white dwarfs, while the secondaries are made either of C/O or of pure helium. We account for tidal effects in the initial conditions in a self-consistent way, and consider initially well-separated systems with slow inspiral rates. We study the merger evolution using an adaptive mesh refinement, reactive, Eulerian code in three dimensions, assuming symmetry across the orbital plane. We use a corotating reference frame to minimize the effects of numerical diffusion, and solve for self-gravity using a multigrid approach. We find a novel detonation mechanism in C/O mergers with massive primaries. Here, the detonation occurs in the primary's core and relies on the combined action of tidal heating, accretion heating, and self-heating due to nuclear burning. The exploding structure is compositionally stratified, with a reverse shock formed at the surface of the dense ejecta. The existence of such a shock has not been reported elsewhere. The explosion energy (1.6 × 10 51 erg) and 56 Ni mass (0.86 M ⊙ ) are consistent with an SN Ia at the bright end of the luminosity distribution, with an approximated decline rate of Δm 15 (B) ≈ 0.99. Our study does not support double-detonation scenarios in the case of a system with a 0.6 M ⊙ helium secondary and a 0.9 M ⊙ primary. Although the accreted helium detonates, it fails to ignite carbon at the base of the boundary layer or in the primary's core.

79 ASTRONOMY AND ASTROPHYSICS↗

Return on Investment (ROI) Framework Case Study: CTH.

CTH is a Eulerian code developed at Sandia National Laboratories capable of modeling the hydrodynamic response of explosives, liquids, gases, and solids. The code solves complex multi-dimensional problems characterized by large deformations and strong shocks that are composed of various material configurations. CTH includes models for material strength, fracture, porosity, and high explosive detonation and initiation. The code is an acronym for a complex series of names relating to its origin. A full explanation can be seen in Appendix A. The software breaks penetration simulations into millions of grid-like “cells”. As a modeled projectile impacts and penetrates a target, progressively smaller blocks of cells are placed around the projectile, which show in detail deformations and breakups. Additionally, the code is uniquely suited to modeling blunt impact and blast loading leading to human body injury.

42 ENGINEERING↗

Applying the Gyrokinetic Formulation of Magnetized Plasma Dynamics to Z-pinch plasmas

This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

V&V for turbulent mixing in the intermediate asymptotic regime

Turbulent fluid flow interacts nonlinearly with a number of processes, including reactive flow (chemistry), radiation, and the transport of particles embedded in the flow. Subgrid scale modeling of these nonlinear processes is a major challenge to computational science. To this end, turbulent mixing of dissimilar fluids is a driver. Here we study validation issues for a class of turbulent mixing simulations. In previous work, we have obtained agreement between simulation and experiment. The present paper addresses a remaining issue, the role of meniscus related boundary effects for immiscible fluids in the validation process. Validation and even verification has been controversial for turbulent mixing flows, as commonly used numerical models are diffusive and the resulting numerical diffusion modifies the parameters of the problem and its solution. Additionally, numerical diffusion of concentration and temperature is well documented in the scientific literature for Eulerian simulations. Lagrangian simulations, which might appear to avoid this problem, are subject to mesh tangling, ensuing regridding and use of Arbitrary Lagrangian Eulerian codes to substitute for a pure Lagrangian algorithm. In practice, Lagrangian methods are potentially subject to the same numerical diffusion issues when used to study complex interface instabilities.

97 MATHEMATICS AND COMPUTING↗

Contaminant Entrainment in a Liquid Fuel Fire with Multi-Component Evaporation Droplet Model

A number of fire scenarios can result in the release of hazardous chemical or radiological contaminants as sub-micron solid particulates, posing a significant danger of inhalation. These particulates may be entrained in the air around the fire through a variety of mechanisms including liquid droplet formation above the boiling surface of the contaminated fuel. The formation of both film drops from the rupture of a liquid film as a boiling bubble bursts and jet drops from the collapse of the cavity of a ruptured bubble may be involved. A 2015 study by Brown and Louie focused on such a fire scenario to investigate the physics of contaminant release. This paper extends that study by supplementing the existing physics models with a multicomponent evaporation mechanism, where liquid droplets formed near the boiling fuel surface experience evaporation of their various components at different rates based on distinct material properties (density, specific heat, heat of vaporization, etc.). As in the Brown and Louie study, the Sierra low-Mach module Fuego is used as the Lagrangian-Eulerian code for simulation of this scenario. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling of direct-drive cylindrical implosion experiments with an Eulerian radiation-hydrodynamics code

Recent improvements to xRAGE, Los Alamos National Laboratory's Eulerian radiation-hydrodynamics code, have enabled the computation of laser-driven experiments relevant to inertial confinement fusion and high energy density physics. Here, previous directly driven cylindrical implosion experiments are modeled in order to benchmark xRAGE design simulations for future cylindrical implosion experiments, representing the first attempt to model such systems with an Eulerian code with adaptive mesh refinement. Simulations in 2D of transverse and axial cross-sections of the cylindrical target are performed, and the results are combined to form a 3D representation of the imploding cylinder. Synthetic radiographs are produced and analyzed from the simulation results, allowing for a direct comparison with experimentally measured quantities. The zeroth-order hydrodynamic trajectories of targets with no specified initial perturbation are well matched by the computations. Simulations of targets with a preimposed sinusoidal perturbation in the azimuthal direction show single-mode instability growth that is in agreement with the available data, but higher fidelity experimental measurements are required to enable more detailed comparisons. As a result, the mode growth observed in computations compares favorably with predictions of a linear theory for the ablative Rayleigh-Taylor instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydro-Code Implementation and Testing of a Kinetic Phase Transition Framework

In this report we describe the Kinetic Phase Transition (KPT) framework that has been worked out over the last 10 years (from around 2014) and the implementation of it into three different codes, the one-dimensional hydro- LASLO and the three-dimensional magneto-hydro- ALEGRA, Sandia codes, via subroutines in the LAMBDA Equations of State and constitutive models package, and Flag, an arbitrary Lagrangian-Eulerian multiphysics code developed within the Lagrangian Applications project (LAP) at LANL. We discuss the introduction of phase mass (and/or volume) fractions that are needed in a code for it to be ‘phase aware’, that is, not only the thermodynamic state is known in each point but also the mixture of the materials’ phases in that point. Further we point to the need of a full Equations of State for each phase in a material to achieve phase awareness and we review the equilibrium phase model, where a phase mixture is at its lowest Gibbs free energy state, to make this point clear. Contrasting the kinetic phase transition to this equilibrium model seamlessly introduce us to the KPT framework that is subsequently thoroughly discussed. While the determination of the total state and the states and mass fractions of phases in each point is a problem that can borrow many of its numerical details from Eulerian codes and mixture of materials (not phases), the update of mass fractions with time in a KPT framework needs a new set of considerations. General for any update model is that we need to prevent mass fractions from becoming unphysical (negative or their sum to be larger than one). We have solved this problem by implementing a subdivision of the hydro time step that prevents the phase from being fully present to not present at all in one subdivided time step by limiting the size of the subdivided time step. This scheme also corrects numerical problems from abrupt changes in parameter values, the so called Gibbs phenomena, that gives rise to slushing between phases in the KPT framework. Interspersed throughout the report are discussions on different thermodynamics considerations. EOS validity windows, limitations on the EOS phase space, are needed for the KPT framework and are discussed separately and exemplified. The KPT framework described in this report has been verified by code comparison, but validation is still an active area of research. There is room for improvement in the update model, both in the model for determination of rates and in how to prevent the mass fractions from becoming unphysical. In addition, the parameters in the KPT update model and the placement of the phase boundary in the EOS phase space, and interactions with other constitutive models, are closely related and interfering with each other. One possible way forward is to simultaneously develop KPT parameters, EOS, and constitutive models for each material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗