Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservation form”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Low-dissipation centred schemes for hyperbolic equations in conservative and non-conservative form

We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗
DOE OSTI2026

Entropy Stable Conservative Flux Form Neural Networks

We propose an entropy-stable conservative flux form neural network (CFN) to predict the dynamics of unknown governing conservation laws. The design of the network is based on the entropy-stable, second-order, and non-oscillatory Kurganov-Tadmor (KT) scheme. The proposed entropy-stable CFN, hereafter referred to as ESCFN, uses slope limiting as a denoising mechanism, ensuring accurate predictions in both noisy and sparse observation environments, as well as in both smooth and discontinuous regions. Importantly, our method is designed to predict long term dynamics of the unknown conservation law exclusively from a short temporal window of observed data, that is, without oracle knowledge of the PDE or later-time solution profiles. Numerical experiments demonstrate that the ESCFN achieves both stability and conservation while maintaining accuracy over extended time domains, and successfully predicts shock propagation speeds in long-term simulations. Furthermore, it is also robust to both noisy and sparse data environments.

Hyperbolic conservation laws↗
DOE OSTI2021

Fluid Dynamics in Curvilinear Coordinates without Fictitious Forces

The use of curvilinear coordinates is sometimes indicated by the inherent geometry of a fluid dynamics problem, but this introduces fictitious forces into the momentum equations that spoil the strict conservative form. If one is willing to work in three dimensions, these fictitious forces can be eliminated by solving for rectangular (Cartesian) momentum components on a curvilinear mesh. A thoroughly geometric approach to fluid dynamics on spacetime demonstrates this transparently, while also giving insight into a greater unity of the relativistic and nonrelativistic cases than is usually appreciated.

pseudo-forces↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2021

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗
DOE OSTI2025

Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator

In our previous work Walton and Tran (2023), numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the long time dynamics of the equation given the multiple blow-up time phenomenon. In this second work in the series, we propose numerical treatments for the complete 3-wave kinetic equations, in which the complete, much more complicated collision operators are fully considered based on a novel conservative form of the equation. Here we then derive an implicit finite volume scheme to solve the equation. The new discretization uses an adaptive time-stepping method which allows for the simulations to be carried to very long times. Our computed solutions are compared with previously derived long-time asymptotic estimates for the decay rate of total energy of time-dependent solutions of 3-wave kinetic equations and found to be in excellent agreement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Generalized Kompaneets formalism for inelastic neutrino-nucleon scattering in supernova simulations

Based on the Kompaneets approximation, we develop a robust methodology to calculate spectral redistribution via inelastic neutrino-nucleon scattering in the context of core-collapse supernova simulations. The resulting equations conserve lepton number to machine precision and scale linearly, not quadratically, with number of energy groups. The formalism also provides an elegant means to derive the rate of energy transfer to matter which, as it must, automatically goes to zero when the neutrino radiation field is in thermal equilibrium. Furthermore, we derive the next higher order in $ϵ/mc^2$ correction to the neutrino Kompaneets equation. Unlike other Kompaneets schema, ours also generalizes to the case of anisotropic angular distributions, while retaining the conservative form that is a hallmark of the classical Kompaneets equation. In this work, our formalism enables immediate incorporation into supernova codes that follow the spectral angular moments of the neutrino radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the n e = d ln T e /d ln n e = 0 case, but substantially larger than the Sugama model for the n e = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Current-conserving relativistic linear response for collisional plasmas

Here we investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi–Dirac equilibrium distribution in terms of the average relaxation rate k. We study the resulting covariant, gauge invariant, and current conserving form of the polarization tensor in the ultrarelativistic and non-relativistic limits. We evaluate the susceptibility in the ultrarelativistic limit and explore their dependence on k. Finally, we study the dispersion relations for the longitudinal and transverse poles of the propagator. We show that for k >2ω p , where ω p is the plasma frequency, the plasma wave modes are overdamped. In the opposite case, k$\ll$ ω p , the propagating plasma modes are weakly damped.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The TPR domain of PgaA is a multifunctional scaffold that binds PNAG and modulates PgaB-dependent polymer processing

The synthesis of exopolysaccharides as biofilm matrix components by pathogens is a crucial factor for chronic infections and antibiotic resistance. Many periplasmic proteins involved in polymer processing and secretion in Gram-negative synthase dependent exopolysaccharide biosynthetic systems have been individually characterized. The operons responsible for the production of PNAG, alginate, cellulose and the Pel polysaccharide each contain a gene that encodes an outer membrane associated tetratricopeptide repeat (TPR) domain containing protein. While the TPR domain has been shown to bind other periplasmic proteins, the functional consequences of these interactions for the polymer remain poorly understood. Herein, we show that the C-terminal TPR region of PgaA interacts with the de-N-acetylase domain of PgaB, and increases its deacetylase activity. Additionally, we found that when the two proteins form a complex, the glycoside hydrolase activity of PgaB is also increased. To better understand structure-function relationships we determined the crystal structure of a stable TPR module, which has a conserved groove formed by three repeat motifs. Tryptophan quenching, mass spectrometry analysis and molecular dynamics simulation studies suggest that the crystallized TPR module can bind PNAG/dPNAG via its electronegative groove on the concave surface, and potentially guide the polymer through the periplasm towards the porin for export. Our results suggest a scaffolding role for the TPR domain that combines PNAG/dPNAG translocation with the modulation of its chemical structure by PgaB.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Opportunities for agrivoltaic systems to achieve synergistic food-energy-environmental needs and address sustainability goals

Achieving decarbonization goals to address global climate change and increasing energy needs requires significant continued investments in solar energy. The expansion of utility-scale solar development across the globe has increased the pressure on land resources for energy generation and other land uses (e.g., agriculture, biodiversity conservation). To address this growing issue, greater emphasis has been placed on solar development strategies that maximize the benefits of solar energy generation and multiple ecosystem services, such as the development of agrivoltaics systems that co-locate solar energy production and various forms of conservation and agricultural land uses. The purpose of this paper is to systematically synthesize the potential ecosystem services of agrivoltaics and summarize how these development strategies could address several United Nations Sustainable Development Goals (SDGs). Our review will focus on four broad potential ecosystem services of agrivoltaics: (1) energy and economic benefits; (2) agricultural provisioning services of food production and animal husbandry; (3) biodiversity conservation; and (4) regulating ecosystem services such as carbon sequestration and water and soil conservation. In particular, we will highlight the state of the science, challenges, and knowledge gaps that represent opportunities for further study to better understand how solar energy deployment can facilitate sustainable development.

14 SOLAR ENERGY↗
DOE OSTI2020

Deeply virtual Compton scattering off the neutron

The three-dimensional structure of nucleons (protons and neutrons) is embedded in so-called generalized parton distributions, which are accessible from deeply virtual Compton scattering. In this process, a high-energy electron is scattered off a nucleon by exchanging a virtual photon. Then, a highly energetic real photon is emitted from one of the quarks inside the nucleon, which carries information on the quark's transverse position and longitudinal momentum. By measuring the cross-section of deeply virtual Compton scattering, Compton form factors related to the generalized parton distributions can be extracted. Here, we report the observation of unpolarized deeply virtual Compton scattering off a deuterium target. From the measured photon-electroproduction cross-sections, we have extracted the cross-section of a quasifree neutron and a coherent deuteron. Due to the approximate isospin symmetry of quantum chromodynamics, we can determine the contributions from the different quark flavours to the helicity-conserved Compton form factors by combining our measurements with previous ones probing the proton's internal structure. These results advance our understanding of the description of the nucleon structure, which is important to solve the proton spin puzzle. The internal structure of the neutron has now been probed by highly energetic photons scattering off it. Combined with previous results for protons, these measurements reveal the contributions of quark flavours to the nucleon structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗