Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conservation equations”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Next Generation Solvent Vapor Pressure Testing

Previous work by Savannah River National Laboratory (SRNL) indicated that the actual Next Generation Solvent vapor pressure would be higher than the Original Caustic Side Solvent Extraction (CSSX) solvent vapor pressure, but below a bounding vapor pressure using Raoult’s Law. Since solvent vapor pressure has a significant impact on Composite Lower Flammability Limits (CLFL) and attendant accident analyses, obtaining an additional margin from the bounding NGS vapor pressure would be valuable. In order to quantify how much margin might be gained from the NGS bounding solvent, SRNL researchers were requested to perform vapor pressure testing with the Next Generation Solvent (NGS) by Savannah River Mission Completion (SRMC). The vapor pressure curve for the NGS formulation set to be deployed at the Salt Waste Processing Facility (SWPF) has been determined by SRNL up to 55°C (131°F) using headspace Gas Chromatography (GC). It was expected that NGS would have a higher vapor pressure than the Original CSSX solvent; however, experimental results indicate a lower vapor pressure. At this time, it is uncertain if this difference is due to the 7x increase in concentration of the large calixarene in the solvent (0.007M BOBCalix in Original CSSX solvent vs. 0.05M MaxCalix in NGS) or to minor batch-to-batch variations in Isopar-L constituents. The NGS and the Original CSSX solvent vapor pressure data were fitted to the Antoine Equation. The Antoine equation gives a more accurate representation of the expected vapor pressure of the solvents outside of the temperature range tested. The Antoine equation fitting for NGS is given below: $P=10^{6.364⁻\frac{1788.4}{T+219.3}}$. Where, p is the vapor pressure (partial pressure) of NGS in mmHg and T is the temperature in °C. It is recommended that SRMC either continue using the more conservative Isopar-L vapor pressure calculations at SWPF with the equation developed for the Original CSSX solvent, or the equation presented above for the NGS solvent. Additionally, it is recommended to study the variability in Isopar-L vapor pressure between lots.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗
DOE OSTI2021

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
DOE OSTI2024

Structure preserving transport stabilized compatible finite element methods for magnetohydrodynamics

Here, we present compatible finite element space discretizations for the ideal compressible magnetohydrodynamic equations. The magnetic field is considered both in div- and curl-conforming spaces, leading to a strongly or weakly preserved zero-divergence condition, respectively. The equations are discretized in space such that transfers between the kinetic, internal, and magnetic energies are consistent, leading to a preserved total energy. We also discuss further adjustments to the discretization required to additionally achieve magnetic helicity preservation. Finally, we describe new transport stabilization methods for the magnetic field equation which maintain the zero-divergence and energy conservation properties, including one method which also preserves magnetic helicity. The methods' preservation and improved stability properties are confirmed numerically using a steady state and a magnetic dynamo test case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hanford Double Shell Waste Tank Corrosion Studies (Final Report FY2019)

During fiscal year (FY) 2019, the Savannah River National Laboratory (SRNL) focused on four primary tasks for Hanford Double Shell Tanks (DSTs) corrosion studies. The first task, New Limits, was a continuation from the electrochemical work started in FY16, to expand electrochemical testing to elevated temperatures and elevated hydroxide concentrations and evaluate how these factors will influence pitting corrosion. Using the pitting factor (PF) equation, these new conditions were validated to provide a conservative estimate for the susceptibility of pitting corrosion on legacy carbon steel. The second task, Secondary Liner Corrosion, focused on Vapor Space Corrosion (VSC) and immersion testing studies using two commercially available vapor corrosion inhibitors (VCIs). For this FY, the VCI strategy was applied mid-experiment to determine performance on weathered coupons. The third task, Long-Term Open Circuit Potential (OCP) Drift, was performed using simplified simulated chemistries to evaluate change in pitting corrosion risk due to evolution of corrosion potential for mill-scale and freshly polished 600-grit surfaces. Finally, the fourth task, Microbiologically-Influenced Corrosion (MIC) studies, was focused on the study of leak detection pit (LDP) water sent from the Hanford tank farm for determination of bacteria that can be conducive to corrosion using commercially available kits: BART TM and MICkit® 5. In addition, this FY, a secondary task was added which was a continuation of the study of the Quartz-enhanced Photoacoustic Spectroscopy (QEPAS) system. Work was started but was ultimately not pursued. A summary of each task performed is presented below.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Adaptive Uncertainty Quantification for Stochastic Hyperbolic Conservation Laws

Here, we propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. Furthermore, we augment the existing SFV theory with a priori convergence results for statistical quantities, in particular push-forward densities, which we demonstrate through numerical experiments. By linking refinement indicators to regions of the physical and stochastic spaces, we drive anisotropic refinements of the discretizations, introducing new degrees of freedom where deemed profitable. To illustrate our proposed method, we consider a series of numerical examples for nonlinear hyperbolic PDEs based on Burgers’ and Euler’s equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Curvature of the chiral phase transition line from the magnetic equation of state of ( 2 + 1 )-flavor QCD

We analyze the dependence of the chiral phase transition temperature on baryon number and strangeness chemical potentials by calculating the leading order curvature coefficients in the light and strange quark flavor basis as well as in the conserved charge ( B , S ) basis. Making use of scaling properties of the magnetic equation of state (MEoS) and including diagonal as well as off-diagonal contributions in the expansion of the energylike scaling variable that enters the parametrization of the MEoS, allows to explore the variation of T c ( μ B , μ S ) = T c ( 1 − ( κ 2 B μ ^ B 2 + κ 2 S μ ^ S 2 + 2 κ 11 B S μ ^ B μ ^ S ) ) along different lines in the ( μ B , μ S ) plane. On lattices with fixed cutoff in units of temperature, a T = 1 / 8 , we find κ 2 B = 0.015 ( 1 ) , κ 2 S = 0.0124 ( 5 ) and κ 11 B S = − 0.0050 ( 7 ) . We show that the chemical potential dependence along the line of vanishing strangeness chemical potential is about 10% larger than along the strangeness neutral line. The latter differs only by about 3% from the curvature on a line of vanishing strange quark chemical potential, μ s = 0 . We also show that close to the chiral limit the strange quark mass contributes like an energylike variable in scaling relations for pseudocritical temperatures. The chiral phase transition temperature decreases with decreasing strange quark mass, T c ( m s ) = T c ( m s phy ) ( 1 − 0.097 ( 2 ) ( m s − m s phys ) / m s phy + O ( ( Δ m s ) 2 ) . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ab initio equation of state for vanadium

Accurate materials’ equations of state (EOS) are essential for understanding materials properties as well as for use in multiphysics simulations. In particular, hydrodynamics simulations are based on three fundamental conservation laws (mass, momentum and energy) that form an under-determined system of equations. The equation of state serves as an additional closure relation between thermodynamic variables for a given material that enables numerical hydrodynamics simulation. In this report, we focus on the development of an ab initio EOS for the body centered cubic (BCC) phase of Vanadium (V) for eventual integration into a multiphase EOS in the OpenSesame EOS database.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Solving Newton’s equations of motion with large timesteps using recurrent neural networks based operators

Classical molecular dynamics simulations are based on solving Newton’s equations of motion. Using a small timestep, numerical integrators such as Verlet generate trajectories of particles as solutions to Newton’s equations. We introduce operators derived using recurrent neural networks that accurately solve Newton’s equations utilizing sequences of past trajectory data, and produce energy-conserving dynamics of particles using timesteps up to 4000 times larger compared to the Verlet timestep. We demonstrate significant speedup in many example problems including 3D systems of up to 16 particles.

Newton’s equations↗
DOE OSTI2021

On Two-Phase Pressure and Temperature Equilibration with Mie-Grüneisen Equations of State

In the course of reactive flow modeling, the phases of high explosive (HE) reactants and products introduce additional degrees of freedom to the underlying hydrodynamic equations of compressible reactive fluid mechanics. Hydrodynamic simulations track the bulk (mixture) conservation properties of mass, momentum and energy. When multiple phases are introduced, one must decide on appropriate closure rules to determine the relative phase volumes and energy distribution between the phases. One typical closure is to assume that the phases are in pressure and temperature equilibrium. Here, the requisite equations of p-T equilibration are explored for Mie-Grüneisen Equations of State (EOS). Examples with Davis reactants and Davis products EOSs will be given. The Davis forms are used in both the WSD and AWSD reactive flow models.

42 ENGINEERING↗
DOE OSTI2021

Quantum Fokker-Planck modeling of degenerate electrons

In this work, an implicit and conservative numerical scheme is proposed for the isotropic quantum Fokker-Planck equation describing the evolution of degenerate electrons subject to elastic collisions with other electrons and ions. The electron-ion and electron-electron collision operators are discretized using a discontinuous Galerkin method, and the electron energy distribution is updated by an implicit time integration method. The numerical scheme is designed to satisfy all conservation laws exactly. Numerical tests and comparisons with other modeling approaches are shown to demonstrate the accuracy and conservation properties of the proposed method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

GenASiS Mathematics: Object-oriented manifolds, operations, and solvers for large-scale physics simulations (version 2)

We report GenASiS Mathematics provides modern Fortran classes furnishing extensible object-oriented functionality for the solution of fields governed by selected partial differential equations. The initial release included extensible object-oriented implementations of simple meshes and the evolution of generic conserved currents thereon. This revision - Version 2 of Mathematics - includes significant reorganization and streamlining of these classes, higher-order reconstruction by a different method, a Poisson solver, coarsening to avoid Courant time step limitations near coordinate singularities, and the offloading of computational kernels to GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The collisional particle-in-cell method for the Vlasov–Maxwell–Landau equations

We introduce an extension of the particle-in-cell method that captures the Landau collisional effects in the Vlasov–Maxwell–Landau equations. The method arises from a regularisation of the variational formulation of the Landau equation, leading to a discretisation of the collision operator that conserves mass, charge, momentum and energy, while increasing the (regularised) entropy. The collisional effects appear as a fully deterministic effective force, thus the method does not require any transport–collision splitting. The scheme can be used in arbitrary dimension, and for a general interaction, including the Coulomb case. We validate the scheme on scenarios such as the Landau damping, the two-stream instability and the Weibel instability, demonstrating its effectiveness in the numerical simulation of plasma.

Bailo, Rafael (ORCID:0000000180183799)↗
DOE OSTI2020

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $k^2_⊥ρ^2_s$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mutually magnetically insulated two-species Brillouin flow

In this work, we use the Brillouin flow analytic framework to examine the physics of Magnetically Insulated Transmission Lines (MITL). Here, we derive a model applicable to any particle species, including both positive and negative ions, in planar and cylindrical configurations. We then show how to self-consistently solve for two-species simultaneously, using magnetically insulated electrons and positive ions as an example. We require both layers to be spatially separated and magnetically insulated (mutually magnetically insulated); for a 7.5 cm gap with a 2 MV bias voltage, this condition requires magnetic fields in excess of 2.73 T. We see a close match between mutually insulated MITL performance and “superinsulated” (high degree of magnetic insulation) electron-only theory, as may be expected for these high magnetic fields. However, the presence of ions leads to several novel effects: (1) Opposite to electron-only theory, total electron currents increase rather than decrease as the degree of magnetic insulation becomes stronger. The common assumption of neglecting electrons for superinsulated MITL operation must be revisited when ions are present—we calculate up to 20× current enhancement. (2) The electron flow layer thickness increases up to double, due to ion space-charge enhancement. (3) The contributions from both ions and electrons to the MITL flow impedance are calculated. The flow impedance drops by over 50% when ions fill the gap, which can cause significant reflections at the load if not anticipated and degrade performance. Additional effects and results from the inclusion of the ion layer are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Novel Relaxation Time Approximation to the Relativistic Boltzmann Equation

In this work, we show that the widely used relaxation time approximation to the relativistic Boltzmann equation contains basic flaws, being incompatible with micro- and macroscopic conservation laws if the relaxation time depends on energy or general matching conditions are applied. We propose a new approximation that fixes such fundamental issues and maintains the basic properties of the linearized Boltzmann collision operator. We show how this correction affects transport coefficients, such as the bulk viscosity and particle diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗