Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moment 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Formulation of 8-moment plasma transport with application to the Nernst effect

8-moment plasma models using two different distribution functions are used to study the Nernst effect and heat transport in dense plasma. These models are presented in hyperbolic form in contrast to traditional parabolic systems derived from perturbing the distribution function, as in Braginskii [Rev. Plasma Phys. 1, 205 (1965)]. The hyperbolic moment formulation can be solved implicitly in time with straightforward and fast local solvers. The numerical implementation of 8-moment models with the relaxation method in the PERSEUS code is also presented. To test 8-moment PERSEUS compared to Braginskii’s transport equations, a verification test for the Nernst thermomagnetic wave by Velikovich et al. [Phys. Plasmas 26, 112702 (2019)] is performed that confirms the presence of the same physics, but with slight differences in the transport coefficients, which are tabulated in the limits of high and low magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Level 2 Milestone: Develop and Incorporate Novel Code Verification of Fundamental Equations in Gemma and Set Up Appropriate Tests

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementations of the electric-, magnetic-, and combinedfield integral equations, are frequently used to solve Maxwell’s equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this report, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches in Gemma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Orbital angular momentum at small $x$ revisited

We revisit the problem of the small Bjorken-x asymptotics of the quark and gluon orbital angular momentum (OAM) distributions in the proton utilizing the revised small-x helicity evolution derived recently in [1]. We relate the quark and gluon OAM distributions at small x to the polarized dipole amplitudes and their (first) impact-parameter moments. To obtain the OAM distributions, we derive novel small-x evolution equations for the impact-parameter moments of the polarized dipole amplitudes in the double-logarithmic approximation (summing powers of α s ln 2 (1/x) with α s the strong coupling constant). We solve these evolution equations numerically and extract the leading large-N c , small-x asymptotics of the quark and gluon OAM distributions, which we determine to be ${L}_{q+\overline{q}}\left(x,{Q}^2\right)\sim {L}_G\left(x,{Q}^2\right)\sim \Delta \Sigma \left(x,{Q}^2\right)\sim \Delta G\left(x,{Q}^2\right)\sim {\left(\frac{1}{2}\right)}^{3.66\sqrt{\frac{\alpha_s{N}_c}{2_π}}}$in agreement with [2] within the precision of our numerical evaluation. (Here N c is the number of quark colors.) We also investigate the ratios of the quark and gluon OAM distributions to their helicity distribution counterparts in the small-x region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Gravothermal evolution of dark matter halos with differential elastic scattering

Here, we study gravothermal evolution of dark matter halos in the presence of differential self-scattering that has strong velocity and angular dependencies. We design controlled N-body simulations to model Rutherford and Møller scatterings in the halo, and follow its evolution in both core-expansion and -collapse phases. The simulations show the commonly-used transfer cross section underestimates the effects of dark matter self-interactions, but the viscosity cross section provides an accurate approximation for modeling angular-dependent dark matter scattering. We investigate thermodynamic properties of the halo, and find that the three moments of the Boltzmann equation under the fluid approximation are satisfied. We further propose a constant effective cross section, which integrates over the halo's characteristic velocity dispersion with weighting kernels motivated by kinetic theory of heat conduction. The effective cross section provides a good approximation to differential self-scattering for most of the halo evolution. It indicates that we can map astrophysical constraints on a constant self-interacting cross section to an SIDM model with velocity- and angular-dependent scatterings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Implementation of Monte Carlo Transport in the General Relativistic SpEC Code

Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the last is a subdominant source of error at the accuracy reached by current simulations and for the interactions currently included in our code. We also provide tests of the most important features of this code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗
DOE OSTI2024

Compton rocket effect due to the action of radiation reaction force in degenerate plasma

A closed set of fluid equations with radiation reaction force (RRF) are constructed from the moments of the appropriate single particle kinetic equation describing a relativistic degenerate (high density) electron plasma. The closure, in analogy with the Maxwellian closure for non-degenerate plasmas, is affected via a parametrized Fermi-Dirac distribution. It is shown that the degeneracy increases RRF just as will be predicted from the so-called “Compton Rocket” effect.

Physics↗
DOE OSTI2020

Analytical model for quasi-linear flow response to resonant magnetic perturbation in resistive-inertial and viscous-resistive regimes

In this work, a quasi-linear model for plasma flow response to the resonant magnetic perturbation (RMP) in a tokamak has been rigorously developed in the resistive-inertial (RI) and viscousresistive (VR) regimes purely from the two-field reduced MHD model. Models for plasma response to RMP are commonly composed of equations for the resonant magnetic field response (i.e. the magnetic island) and the torque balance of plasma flow. However, in previous plasma response models, the magnetic island and the torque balance equations are often derived separately from reduced MHD and full MHD equations, respectively. By contrast, in this work we derive both the magnetic island response and the torque balance equations in a quasi-linear model for plasma flow response entirely from a set of two-field reduced MHD equations. Such a quasi-linear model can recover previous plasma flow response models within certain limits and approximations. Furthermore, the physical origins of quasi-linear forces and moments in the flow response equation are also accurately calculated and clarified self-consistently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Plasma Modeling Hierarchy and Verification Approach

This report reviews a hierarchy of formal mathematical models for describing plasma phenomena. Starting with the Boltzmann equation, a sequence of approximations and modeling assumptions can be made that progressively reduce to the equations for magnetohydrodynamics. Understanding the assumptions behind each of these models and their mathematical form is essential to appropriate use of each level of the hierarchy. A sequence of moment models of the Boltzmann equation are presented, then focused into a generalized three-fluid model for neutral species, electrons, and ions. This model is then further reduced to a two-fluid model, for which Braginskii described a useful closure. Further reduction of the two-fluid model yields a Generalized Ohm's Law model, which provides a connection to magnetohydrodynamic approaches. A verification approach based on linear plasma waves is presented alongside the model hierarchy, which is intended as an initial and necessary but not sufficient step for verification of plasma models within this hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Data Science Enabled Enabled Discovery of Superconductors (Final Progress Report)

This Final Technical Report describes efforts by 4 PIs at the University of Florida (Peter Hirschfeld, Richard Hennig, Greg Stewart and James Hamlin), over the period September 2019-August 2023, to use data science and machine learning techniques to discover new conventional superconductors. The PIs constructed a discovery loop with two theorists and two experimentalists to: develop algorithms to machine learn descriptors correlating strongly with the critical temperature Tc (PI's Peter Hirschfeld, UF Physics and Richard Hennig, UF Materials Science and En), synthesize and measure properties of promising materials, and feed back the knowledge gained into the prediction algorithm. This work was motivated by the theoretical prediction and experimental discovery of high-pressure, high-pressure hydride superconductors, and to find ways to recreate the high critical temperatures in these systems at ambient pressure. Highlights from the grant include: 1) a new equation for Tc in terms of moments of the electron-phonon spectral function, improving on the so-called Allen-Dynes equation (1975); 2) study of the metastable A15 superconductor Nb3Si, formed under explosive compression at ~1000GPa to determine the kinetic barrier to the ground state structure; 3) the development of ultra-fast machine-learned atomic potentials for molecular dynamics, and 4) the discovery of superconductivity at 19K in WB2 arising from metastable defect structures in the crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗
DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fluid structure interaction of a subaqueous pendulum: Analyzing the effect of wake correction via large eddy simulations

The dynamic behavior of a subaqueous cylindrical pendulum and corresponding flow dynamics are investigated. The objectives were twofold: (i) to examine whether the two-dimensional model equations sufficiently capture the three-dimensional dynamics and (ii) to investigate the emerging three-dimensional vortical flow structures. Large eddy simulations with two-way coupling fluid structure interaction were carried out using the immersed boundary method to simulate the motion of the pendulum and its interactions with the initially stagnant water. The resulting pendulum motion is compared against measured data obtained in a series of experimental tests to validate the simulation results and the model equations with and without wake corrections. An analysis of the flow vorticity revealed the development of a vortex ring during the first swing and the formation of tip vortices. The evolution of the vortex rings emerging from the motion of the subaqueous cylindrical pendulum was visualized using Q-criteria showing a reasonable agreement with vortical structures observed in the experiment using particle imaging velocimetry. The hydrodynamic moments acting on the simulated pendulum and the moments calculated from the model equations are analyzed. Using the insights from these numerical simulations, a modification of the wake correction is proposed to enhance the accuracy of the rate of decay and period. The transient effect of coherent flow on pendulum dynamics, especially the added mass effect, is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Geometric integration of classical spin dynamics via a mean-field Schrödinger equation

The Landau-Lifshitz equation describes the time evolution of magnetic dipoles and can be derived by taking the classical limit of a quantum mechanical spin Hamiltonian. To take this limit, one constrains the many-body quantum state to a tensor product of coherent states, thereby neglecting entanglement between sites. Expectation values of the quantum spin operators produce the usual classical spin dipoles. One may also consider expectation values of polynomials of the spin operators, leading to quadrupole and higher-order spin moments, which satisfy a dynamical equation of motion that generalizes the Landau-Lifshitz dynamics [Zhang and Batista, Phys. Rev. B 104, 104409 (2021)]. Here we reformulate the dynamics of these N 2 –1 generalized spin components as a mean-field Schrödinger equation on the N-dimensional coherent state. Furthermore, this viewpoint suggests efficient integration methods that respect the local symplectic structure of the classical spin dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Modelling the ejection of primary aerosols during the fast pyrolysis of biomass anisotropic particles

A model for the fast pyrolysis of anisotropic biomass particles is presented which considers bubbling dynamics within the liquid intermediate phase (metaplast) and aerosol ejection from this phase. The model employs the population balance equation and the method of moments to estimate the production rate and resultant size distribution of aerosol ejections, incorporating a detailed CRECK reaction mechanism, and considers the effect of anisotropic biomass microstructure on the intraparticle transport of mass and energy. Here, this study investigates the impact of particle size, heating rate (heat transfer coefficient), and lignocellulosic composition on aerosol ejection. The model predicts that, at high heating rates (convective heat transfer coefficient of 359 W/m 2 .K), aerosols can contribute over 20% to the heavy fraction yield in bio-oil for small particles (1 mm diameter, 4 mm length). The model can predict aerosol size distribution and surface area, indicating an average size of 20 μm for bubbles and 5 μm for aerosols during increased bubble production and aerosol ejection rates. These findings are consistent with prior experimental results and provide essential information for future modeling of extra-particle reactions of the aerosols as they progress through the reactor.

09 BIOMASS FUELS↗