Search NASA⌕ Search

SEARCH · Search NASA

Results for “Artificial diffusion flux”

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.

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state: A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state

The development of numerical approaches to perform direct numerical simulations of compressible multiphase flows has been an active field of research for several years. Proper treatment of fluid interfaces is crucial as important physics occur in this infinitesimally small region. Furthermore, the compressibility of the fluid requires proper treatment of discontinuities. Artificial diffusivity is among a number of methods widely used for compressible flows. This study develops a general form of consistent artificial diffusion fluxes and extends the localized artificial diffusivity method for high-order central schemes to solve multiphase flows with an interface-capturing method. These fluxes ensure an oscillation-free interface for pressure, velocity, and temperature without employing a sharpening technique. Moreover, the high-order representation of all scales in the flow helps capture the wide range of instabilities inherent in these flows. The goal is to develop an approach capable of performing high-fidelity simulations supported by physics-driven validation. lastly, this is achieved by solving the five-equation model with the stiffened-gas equation of state using the proposed method for multicomponent and multiphase flows on a variety of 1D and 2D problems.

42 ENGINEERING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Scalar conservation and boundedness in simulations of compressible flow

With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The reality of artificial viscosity

Artificial viscosity is used in the computer simulation of high Reynolds number flows and is one of the oldest numerical artifices. In this work, I will describe the origin and the interpretation of artificial viscosity as a physical phenomenon. The basis of this interpretation is the finite scale theory, which describes the evolution of integral averages of the fluid solution over finite (length) scales. I will outline the derivation of finite scale Navier–Stokes equations and highlight the particular properties of the equations that depend on the finite scales. Those properties include enslavement, inviscid dissipation, and a law concerning the partition of total flux of conserved quantities into advective and diffusive components.

42 ENGINEERING↗
DOE OSTI2020

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Bi-Layer Dense/Porous Solid Electrolyte Interphase for Enhanced Lithium-Metal Stability

Due to its high theoretical capacity and low electrochemical potential, lithium metal is a highly investigated anode for next-generation high-energy batteries. However, the unstable chemical and topographical heterogeneous surface of lithium gives rise to safety and efficiency concerns that prevent it from being utilized in practical applications. Exposure to electrolyte leads to non-uniform, mixed organic-inorganic solid-electrolyte interphase (SEI) formation, starting a cascading to non-uniform flux distribution, consumption of active materials, and dendrite growth. The SEI is the key feature that will lead to harnessing the capabilities possible with lithium metal. In this work, the formation of a closed-host bi-layer solid electrolyte interphase (SEI) improves the stability of lithium anode. This is successfully realized by forming an interconnected porous LiF-rich artificial SEI in contact with Li metal, and a dense, stable in-situ formed upper layer SEI. The porous layer increases the number of Li/LiF interfaces, which reduces local volume fluctuations and improves Li+ diffusion along these interfaces. Additionally, the tortuous porous structure guides uniform Li+ flux distribution and mechanically suppresses dendrite propagation. The dense upper layer of the SEI accomplishes a closed-host design through reaction with LiNO3 additive and prevents continuous consumption of active materials seen without the additive included in the electrolyte. The duality of a dense top layer with porous bottom layer led to extended cycle life and improved rate performance, evidenced with symmetric cell testing, as well as full cell testing paired with sulfur and LiFePO4 (LFP) cathodes. This work is a good example of a rational design of the SEI, based on comprehensive consideration of various critical factors to improve Li-metal anode stability, and highlights a new pathway to improve cycling and rate performances of Li metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2018

Modeling the impact of radiation-enhanced diffusion on implanted ion profiles

We report that ion irradiations are a valuable research tool for exploring radiation effects in materials. However, it is well recognized that the implanted ions can artificially modify the radiation damage evolution, e.g., enhancing amorphization processes at low irradiation temperatures and suppressing void swelling at elevated temperatures. Therefore, the implanted ion region should be avoided for most studies of ion irradiated materials. Due to increased interest in high damage, high temperature ion irradiations studying radiation effects in materials for proposed high dose Generation IV fission and fusion energy applications, it is crucial to quantify the extent of diffusional broadening of the implanted ion profile for a variety of temperatures, irradiation fluxes, and sink strengths. The present study summarizes computational analyses of thermal and depth-dependent radiation enhanced difusion (RED) on diffusion broadening of the implanted ion profiles in Fe and Ni for a variety of irradiation conditions. For a low assumed RED coefficient (10 -20 -10 -19 m 2 /s and 10 -4 and 10 -3 peak dpa/s, respectively) or high flux, broadening of the as-implanted ion profile is very small and a suitable artifact-free midrange region for analysis exists for ion energies above 5-6 MeV at 100 peak dpa. At high RED coefficients (10 -19 -10 -18 m 2 /s and 10 -4 and 10 -3 peak dpa/s, respectively) broadening is much more significant, and no valid region for investigation exists below 6-8 MeV ion energies at any damage level >50 dpa. While increasing flux decreases irradiation time, it also increases the RED coefficient; these effects mostly offset except under recombination-dominant conditions. For typical high dose irradiation studies of void swelling that exceed ~100 displacement per atom (dpa), ion energies > 6-8 MeV must be employed in order to achieve suitable artifact-free midrange irradiation analysis regions, depending on the material system. Lastly, analysis of amorphization/precipitation require even higher energies (>10 MeV).

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2018

Ni Nanoparticle-Added Ni{sub x}/(Cu{sub 0.5}Tl{sub 0.5})Ba{sub 2}Ca{sub 2}Cu{sub 3}O{sub 10-δ} Superconductor Composites and Their Enhanced Flux Pinning Characteristics

Ni nanoparticle-added Ni{sub x}/(Cu{sub 0.5}Tl{sub 0.5})Ba{sub 2}Ca{sub 2} Cu{sub 3}O{sub 10-δ} superconductor composites are synthesized at normal pressure by two-step solid-state reaction method. The idea for adding ferromagnetic nanoparticles was to create artificial pinning centers at grain boundaries. The samples have shown orthorhombic crystal structure in which the cell parameters increase with an increase in the added concentration of Ni nanoparticles. This shows that due to the higher reactivity, Ni nanoparticles partially diffuse into the unit cell of (Cu{sub 0.5}Tl{sub 0.5})Ba{sub 2}Ca{sub 2}Cu{sub 3}O{sub 10-δ} superconductors. This also indicates as an evidence of the occupancy of these nanoparticles at grain boundaries. A metallic variation of resistivity from room temperature down to the onset of superconductivity is a typical feature of these samples with T{sub c}(R = 0) varying between 92 and 94 K. The magnitude of diamagnetism is significantly enhanced in 7 and 10% nanoparticle-added samples. The enhanced magnetization of ferromagnetic nanoparticles at reduced temperature plays a vital role in reducing the diamagnetism of nanosuperconductor composites. The softening of apical oxygen modes of type Cu(1)–O{sub A}–Cu(2) has confirmed the diffusion of Ni nanoparticles into the unit cell of the final compound. The excess conductivity analyses of conductivity data have shown that the values of B{sub c}(0), B{sub c1}, and J{sub c} are enhanced with the addition of Ni nanoparticles in the final compound. Variation of superconducting properties was attributed to pair breaking due to spin scattering across these ferromagnetic Ni nanoparticles. This shows that added Ni nanoparticles act as efficient pinning centers, which is also confirmed by the suppression of the London penetration depth.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Modeling Ring Current Proton Fluxes Using Artificial Neural Network and Van Allen Probe Measurements

Abstract Terrestrial ring current dynamics are a critical part of the near‐space environment, in that they directly drive geomagnetic field variations that control particle drifts, and define geomagnetic storms. The present study aims to specify a global and time‐varying distribution of ring current proton using geomagnetic indices and solar wind parameters with their history as input. We train an artificial neural network (ANN) model to reproduce proton fluxes measured by the Radiation Belt Storm Probes Ion Composition Experiment instrument onboard Van Allen Probes. By choosing optimal feature parameters and their history length, the model results show a high correlation and a small error between model specifications and satellite measurements. The modeled results well capture energy‐dependent proton dynamics in association with geomagnetic storms, including inward radial diffusion, acceleration and decay. Our ANN model produces proton fluxes with their corresponding 3D spatiotemporal variations, capturing the latitudinal distribution and local time asymmetry that are consistent with observations and that can further inform theory.

Li, Jinxing↗
DOE OSTI2018

Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems

The objective of this paper is to present a local bounds preserving stabilized finite element scheme for hyperbolic systems on unstructured meshes based on continuous Galerkin (CG) discretization in space. A CG semi-discrete scheme with low order artificial dissipation that satisfies the local extremum diminishing (LED) condition for systems is used to discretize a system of conservation equations in space. The low order artificial diffusion is based on approximate Riemann solvers for hyperbolic conservation laws. In this case we consider both Rusanov and Roe artificial diffusion operators. In the Rusanov case, two designs are considered, a nodal based diffusion operator and a local projection stabilization operator. The result is a discretization that is LED and has first order convergence behavior. To achieve high resolution, limited antidiffusion is added back to the semi-discrete form where the limiter is constructed from a linearity preserving local projection stabilization operator. The procedure follows the algebraic flux correction procedure usually used in flux corrected transport algorithms. To further deal with phase errors (or terracing) common in FCT type methods, high order background dissipation is added to the antidiffusive correction. The resulting stabilized semi-discrete scheme can be discretized in time using a wide variety of time integrators. Numerical examples involving nonlinear scalar Burgers equation, and several shock hydrodynamics simulations for the Euler system are considered to demonstrate the performance of the method. In conclusion, for time discretization, Crank–Nicolson scheme and backward Euler scheme are utilized.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Entropy Stable Discretization of Compressible Flows in Thermochemical Nonequilibrium

Entropy stable numerical methods for compressible flow have been demonstrated to exhibit better robustness than purely linearly stable methods and need less overall artificial dissipation for long simulations in subsonic and transonic flows. In this work we seek to extend these benefits to multicomponent, multitemperature flows in thermochemical nonequilibrium such as combustion and hypersonic flight. We first derive entropy functions that symmetrize the governing equations and allow stability proofs for such systems. The impact of diffusion model selection on provable entropy stability is considered in detail, including both rigorous models of irreversible thermodynamics and simplified models of greater practical interest. Based on the proven entropy functions we develop affordable, entropy conservative two-point flux functions for solution in conservation form. We derive entropy conservative fluxes for calorically and thermally perfect mixtures, with heat capacities described by either polynomials of the temperature or formulas from statistical thermodynamics.

42 ENGINEERING↗
DOE OSTI2021

The design of a photoionization front experiment using the Z-Machine as a driving source and estimated measurements

Radiation-driven heat fronts are present in the early universe during reionization, the circumstellar medium of supernovae, and in high-energy-density physics experiments. Dedicated experiments to observe and diagnose the behavior of these types of heat fronts can improve our understanding of these phenomena. A simulation study of photoionization fronts using the HELIOS-CR radiation hydrodynamics code provides an experimental design for the Z-Machine at Sandia National Laboratory using a measurement-calibrated input radiation flux to drive the photoionization front. The simulations use detailed atomic physics and non-diffusive radiation transport in 1D to determine an optimal gas pressure of 0.75 atm for an experiment in N gas as well as the effects of increasing the thickness of the window that seals the gas cell. Post-processing of these simulations demonstrates that ratios of atomic rate coefficients place the heat front in a physics regime where photoionization dominates the energy deposition. Here, to see the sensitivity of the simulations to changes in the model and spatial grid, this analysis performed resolution, atomic model detail, and radiation transport angular grid studies showing less than 10% deviation from the nominal model for increased complexity, when possible. An effort to emulate 3D geometric effects on the radiation flux using an artificial attenuation scheme has shown that, even for conservative estimates of the flux, simulations still produce a photoionization front. Estimations of a streaked, visible spectroscopy measurement using SPECT3D showed that line emission measurements are present early in time and that later in time thermal emission should become dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗