Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Verification of iterative matrix solutions for multipoint kinetics equations

Multipoint kinetic equation systems have been solved numerically using matrix algebra software and 1st to 4th order implicit schemes based on single-step matrix propagators. These matrix schemes have been validated successfully on demanding point kinetics benchmarks with various prescribed reactivity insertions. Verification tests have also been performed on a simple 3-region fast reactor core model with asymmetric step or ramp reactivity insertions and cross-checked with the multipoint kinetics SACRE code developed at INL and based on existing and validated stiff ODE solver packages. Accurate results can be obtained at a limited time expense. The reason for intriguing results obtained in some transients has been elucidated and linked to the multipoint matrix coefficient interpolation method during the transient.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2021

Lie group analysis of the point-reactor neutron kinetics equations for various reactivity models

This work applies Lie Group Theory to the nonlinear point-reactor neutron kinetics equations for several reactivity feedback and reactivity insertion models. The goal of this work is to connect known analytical solutions to their corresponding Lie groups as well as to obtain new analytical solutions and the circumstances from which they arise. The reactivity models we consider are the Nordheim-Fuchs Model, an arbitrary-order polynomial reactivity insertion model, and the Fuchs Ramp-Insertion Model. In all three models, we generalize the analysis when possible and obtain analytical solutions in every case. For the first two models, we succeed in finding the Lie groups corresponding to classical solutions and find additional solutions that were otherwise unknown. Finally, for the final model, we obtain a reduced-order nonlinear kinetics equation which is not analytically solvable, but by applying a second-order truncation, we find a novel closed-form solution and demonstrate its validity with a numerical example.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A Numerical Scheme for Wave Turbulence: 3-Wave Kinetic Equations

Here, we introduce a finite volume scheme to solve a special case of isotropic 3-wave kinetic equations. We test our numerical solution against theoretical results concerning the long time behavior of the energy and observe that our solutions verify the energy cascade phenomenon. To our knowledge, this is the first numerical scheme that can capture the long time asymptotic behavior of solutions to those isotropic 3-wave kinetic equations, where the energy cascade can be observed. Our numerical energy cascade rates are in good agreement with previously obtained theoretical results. The finite volume scheme given here relies on a new identity, allowing one to reduce the number of terms needed in the collision operators.

3-wave equation↗
DOE OSTI2024

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalised kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, B $\vartheta$ /B, it generalises the drift kinetic theory of [Hazeltine Phys. Plasmas 15 (1973) 77] to the limit of B $\vartheta$ /B ~ 1 (e.g. to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalised gyrokinetic formalism of [Dudkovskaia et al. Plasma Phys. Control. Fusion 65 (2023) 045010], derived specifically to capture neoclassical equilibrium currents in gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of [Dudkovskaia et al. Nucl. Fusion 63 (2023) 016020] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. In conclusion, resolving ρ ~ ρ $\vartheta$ (or B $\vartheta$ ~ B), where ρ $\vartheta$ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A Regularized Entropy-Based Moment Method for Kinetic Equations

Here, we present a new entropy-based moment method for the velocity discretization of kinetic equations. This method is based on a regularization of the optimization problem defining the original entropy-based moment method, and this gives the new method the advantage that the moment vectors of the solution do not have to take on realizable values. We show that this equation still retains many of the properties of the original equations, including hyperbolicity, an entropy-dissipation law, and rotational invariance. The cost of the regularization is mismatch between the moment vector of the solution and that of the ansatz returned by the regularized optimization problem. However, we show how to control this error using the parameter defining the regularization. This suggests that with proper choice of the regularization parameter, the new method can be used to generate accurate solutions of the original entropy-based moment method, and we confirm this with numerical simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Multiscale convergence properties for spectral approximations of a model kinetic equation

Here in this work, we prove rigorous convergence properties for a semi-discrete, moment-based approximation of a model kinetic equation in one dimension. This approximation is equivalent to a standard spectral method in the velocity variable of the kinetic distribution and, as such, is accompanied by standard algebraic estimates of the form N -q , where N is the number of modes and q > 0 depends on the regularity of the solution. However, in the multiscale setting, the error estimate can be expressed in terms of the scaling parameter ϵ, which measures the ratio of the mean-free-path to the characteristic domain length. We show that, for isotropic initial conditions, the error in the spectral approximation is $\mathcal{O}$($ϵ^{N+1}$). More surprisingly, the coefficients of the expansion satisfy super convergence properties. In particular, the error of the ℓ th coefficient of the expansion scales like $\mathcal{O}$($ϵ^{2N}$ ) when ℓ = 0 and $\mathcal{O}$($ϵ^{2N+2-ℓ}$ ) for all 1 ≤ ℓ ≤ N. This result is significant, because the low-order coefficients correspond to physically relevant quantities of the underlying system. All the above estimates involve constants depending on N, the time t, and the initial condition. We investigate specifically the dependence on N, in order to assess whether increasing N actually yields an additional factor of ϵ in the error. Numerical tests will also be presented to support the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

High-order accurate solutions of the point kinetics equations with the spectral deferred correction method

Solving initial value problems with higher-order methods can improve the accuracy of the simulation results or the efficiency of the calculation. In this paper, we apply the spectral deferred correction (SDC) method to solve the initial value problem of the point kinetics equations (PKE). SDC is a stable, robust, and efficient high-order time-integration method capable of an arbitrary order of accuracy. For our implementation we show that it is A-stable for orders up to 8 and the order of accuracy is verified for PKE problems with a range of different reactivities. A 5.-order SDC method was then implemented to solve the exact PKE (EPKE) in the Transient Multilevel (TML) method of MPACT. The error from solutions of the EPKE is shown to be negligible. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Physics-Informed Neural Network Solution of Point Kinetics Equations for a Nuclear Reactor Digital Twin

A digital twin (DT) for nuclear reactor monitoring can be implemented using either a differential equations-based physics model or a data-driven machine learning model. The challenge of a physics-model-based DT consists of achieving sufficient model fidelity to represent a complex experimental system, whereas the challenge of a data-driven DT consists of extensive training requirements and a potential lack of predictive ability. We investigate the performance of a hybrid approach, which is based on physics-informed neural networks (PINNs) that encode fundamental physical laws into the loss function of the neural network. We develop a PINN model to solve the point kinetic equations (PKEs), which are time-dependent, stiff, nonlinear, ordinary differential equations that constitute a nuclear reactor reduced-order model under the approximation of ignoring spatial dependence of the neutron flux. The PINN model solution of PKEs is developed to monitor the start-up transient of Purdue University Reactor Number One (PUR-1) using experimental parameters for the reactivity feedback schedule and the neutron source. The results demonstrate strong agreement between the PINN solution and finite difference numerical solution of PKEs. We investigate PINNs performance in both data interpolation and extrapolation. For the test cases considered, the extrapolation errors are comparable to those of interpolation predictions. Extrapolation accuracy decreases with increasing time interval.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Quarkonium Polarization Kinetic Equation from Open Quantum Systems and Effective Field Theories

Recent measurements of polarization phenomena in relativistic heavy ion collisions have aroused a great interest in understanding dynamical spin evolution of the QCD matter. In particular, the spin alignment signature of J/ψ has been recently observed in Pb-Pb collisions at LHC, which may infer nontrivial spin transport of quarkonia in quark gluon plasmas. Motivated by this, we study the spin-dependent in-medium dynamics of quarkonia by using the potential nonrelativistic QCD (pNRQCD) and the open quantum system framework. By applying the Markovian approximation and Wigner transformation, we systematically derive the Boltzmann transport equation for vector quarkonia with polarization dependence in the quantum optical limit. As opposed to the previous study for the spin-independent case where the collision terms depend on chromoelectric correlators, the new kinetic equation incorporates gauge invariant correlators of chromomagnetic fields that determine the recombination and dissociation terms with polarization dependence at the order we are working in the multipole expansion. In the quantum Brownian motion limit, the Lindblad equation with new transport coefficients defined in terms of the chromomagnetic field correlators have also been derived. Our formalism is generic and valid for both weakly-coupled and strongly-coupled quark gluon plasmas. It may be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Wave kinetic equation for inhomogeneous drift-wave turbulence beyond the quasilinear approximation

Here, the formation of zonal flows from inhomogeneous drift-wave (DW) turbulence is often described using statistical theories derived within the quasilinear approximation. However, this approximation neglects wave–wave collisions. Hence, some important effects such as the Batchelor–Kraichnan inverse-energy cascade are not captured within this approach. Here we derive a wave kinetic equation that includes a DW collision operator in the presence of zonal flows. Our derivation makes use of the Weyl calculus, the quasinormal statistical closure and the geometrical-optics approximation. The obtained model conserves both the total enstrophy and energy of the system. The derived DW collision operator breaks down at the Rayleigh–Kuo threshold. This threshold is missed by homogeneous-turbulence theory but expected from a full-wave quasilinear analysis. In the future, this theory might help better understand the interactions between drift waves and zonal flows, including the validity domain of the quasilinear approximation that is commonly used in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗