Transport equations for gases and plasmas obtained by the 13-moment method - A summary
Momentum transfer solutions to Boltzmann transport equations for gas mixture and plasmas
Harstad, K. G.↗
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.
Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection
A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.
79 ASTRONOMY AND ASTROPHYSICS↗
Non-relativistic equations of bulk motion of a relativistic gas. II.
Nonrelativistic bulk motion equations of relativistic gas deduced from moments of collisionless Boltzmann equation, discussing applications to galactic cosmic ray gas
Lerche, I.↗
Closures and Simulation for Thermal Radiation Transport in Stochastic Media with Nonlinear Temperature Dependence
Because of the practical challenge of rendering very complex realistic spatial structures for numerical work, it is common practice to resort to characterizing such media as stochastic mixtures of materials, ideally parametrized with low order statistics such as the mean, variance, and correlation functions of the now random material properties. This enables realizations of the medium to be repeatedly generated and radiation transport computations to, in principle, be performed for a large ensemble of these realizations to obtain a statistically well-characterized radiation field. Statistical post-processing yields desired quantities such as conditional and unconditional mean radiation flux and probability distributions of transmitted radiation. However, such computations prove expensive for all but the simplest stochastic geometries and are most suited for benchmarking approximate models. The most common approximations lead to homogenized media so that transport computations are required only on a single medium realization but by construct provide only limited statistical information on the radiation field. Almost all approximate approaches to this problem attempt to develop equations for low order moments of the radiation intensity (mean, second moment, correlation function) but inevitably encounter a closure problem: the equation for any statistical moment will contain terms depending on unknown higher-order moments. Thus, the challenge shifts to one of developing closure relations that relate the unknown moments to the lower order moments. Under very special conditions, an exact closure can be derived but in general closures are heuristically stated constitutive relations. Also, closure approaches depend on whether the mixing statistics are spatially and/or temporally continuous as in fluctuating turbulent fields, or discontinuous as in randomly mixed solid chunks of material. Thus, unconditional averaging is generally applied in the former case but conditional averaging is more appropriate when the mixing is discontinuous. In this work, the emphasis is on binary statistical mixtures of immiscible fluids as as such the quantities of interest are averages (flux, temperature) conditioned on the material type.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Scatter-free collimated convection and cosmic-ray transport at 1 AU
Observations of solar cosmic rays between 0.1 MeV and 1 GeV are more consistent with scatter-free than diffusive propagation near 1 AU. Since these particles also are constrained to interplanetary magnetic field lines, the parallel and transverse flow are decoupled by collimated convection. Local propagation is therefore more appropriately described by Liouville's equation, and transport properties can be examined from the exact zeroeth and first moments of the equation averaged over all directions in momentum space. Observational results which may be explained directly without invoking local scattering include the decay of solar flare events as well as the radial gradient, diurnal and semidiurnal variation of galactic cosmic rays.
Roelof, E. C.↗
Stochastic Modeling of the Joint Neutron Number-Cumulative Fission Fragment Kinetic Energy Deposition Distribution and its Statistical Moments [Slides]
We investigate the joint distribution of the neutron number and cumulative fission-fragment kinetic energy (FKE) deposition, with a specific focus on low-order statistical moments: the mean, variance, and correlation. Starting from a point-kinetic framework, we derive a forward Master equation (FME) for the joint distribution and develop the corresponding moment equations.
42 ENGINEERING↗
A critical analysis of the Grad approximation for closing out the magnetohydrodynamic equations for plasmas
Error analysis on moment distribution functions of magnetohydrodynamic plasma equations
Meador, W. E.↗
Quantum theory of an optical maser. VI - Transient behavior.
The transient behavior of a laser is discussed using the quantum theory as did Scully and Lamb. The formal solution of the density-matrix equation is expressed in terms of exponentially decaying eigenmodes. Some of the lower decay constants are obtained numerically. The equations for the moments of the density matrix are then derived and solved by a truncation method. The equations of motion are integrated numerically for the case where the average number of photons in a laser cavity has the realistically large value 1.3 x 100,000. An alternative Fokker-Planck-equation approach is discussed.
Wang, Y. K.↗
Third-moment closure of turbulence for predictions of separating and reattaching shear flows: A study of Reynolds-stress closure model
A numerical study of computations in backward-facing steps with flow separation and reattachment, using the Reynolds stress closure is presented. The highlight of this study is the improvement of the Reynold-stress model (RSM) by modifying the diffusive transport of the Reynolds stresses through the formulation, solution and subsequent incorporation of the transport equations of the third moments, bar-u(i)u(j)u(k), into the turbulence model. The diffusive transport of the Reynolds stresses, represented by the gradients of the third moments, attains greater significance in recirculating flows. The third moments evaluated by the development and solution of the complete transport equations are superior to those obtained by existing algebraic correlations. A low-Reynolds number model for the transport equations of the third moments is developed and considerable improvement in the near-wall profiles of the third moments is observed. The values of the empirical constants utilized in the development of the model are recommended. The Reynolds-stress closure is consolidated by incorporating the equations of k and e, containing the modified diffusion coefficients, and the transport equations of the third moments into the Reynolds stress equations. Computational results obtained by the original k-e model, the original RSM and the consolidated and modified RSM are compared with experimental data. Overall improvement in the predictions is seen by consolidation of the RMS and a marked improvement in the profiles of bar-u(i)u(j)u(k) is obtained around the reattachment region.
Amano, R. S.↗
Second-order closure PBL model with new third-order moments: Comparison with LES data
This paper contains two parts. In the first part, a new set of diagnostic equations is derived for the third-order moments for a buoyancy-driven flow, by exact inversion of the prognostic equations for the third-order moment equations in the stationary case. The third-order moments exhibit a universal structure: they all are a linear combination of the derivatives of all the second-order moments, bar-w(exp 2), bar-w theta, bar-theta(exp 2), and bar-q(exp 2). Each term of the sum contains a turbulent diffusivity D(sub t), which also exhibits a universal structure of the form D(sub t) = a nu(sub t) + b bar-w theta. Since the sign of the convective flux changes depending on stable or unstable stratification, D(sub t) varies according to the type of stratification. Here nu(sub t) approximately equal to wl (l is a mixing length and w is an rms velocity) represents the 'mechanical' part, while the 'buoyancy' part is represented by the convective flux bar-w theta. The quantities a and b are functions of the variable N(sub tau)(exp 2), where N(exp 2) = g alpha derivative of Theta with respect to z and tau is the turbulence time scale. The new expressions for the third-order moments generalize those of Zeman and Lumley, which were subsequently adopted by Sun and Ogura, Chen and Cotton, and Finger and Schmidt in their treatments of the convective boundary layer. In the second part, the new expressions for the third-order moments are used to solve the ensemble average equations describing a purely convective boundary laye r heated from below at a constant rate. The computed second- and third-order moments are then compared with the corresponding Large Eddy Simulation (LES) results, most of which are obtained by running a new LES code, and part of which are taken from published results. The ensemble average results compare favorably with the LES data.
Canuto, V. M.↗
Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model
The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Code-verification techniques for the method-of-moments implementation of the combined-field integral equation
Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the combined-field integral equation. Here, we demonstrate the effectiveness of these approaches for cases with and without coding errors.
97 MATHEMATICS AND COMPUTING↗
Discontinuous Galerkin Methods for NonLinear Differential Systems
This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.
Barth, Timothy↗
Stability of matrix equations arising in electromagnetics
Many problems in electromagnetics can be formulated as a Fredholm equation of the first kind. In these problems it is desired to determine some function f(x) from its integral transform g(y). Except for special cases for which an explicit inverse operator is known, the equation is usually discretized by the method of moments into a matrix equation Af = g, which is then solved on a computer. A method of calculating the sensitivity of the solution of Af = g due to errors in either the vector g or the matrix A is developed. The usefulness of this method is demonstrated with problems from antenna theory, scattering theory, and remote sensing.
Klein, C.↗
On the 2018 Outburst of the Accreting Millisecond X-Ray Pulsar Swift J1756.9-2508 As Seen with NICER
We report on the coherent timing analysis of the 182 Hz accreting millisecond X-ray pulsar SwiftJ1756.92508during its 2018 outburst as observed with the Neutron Star Interior Composition Explorer (NICER). Combiningour NICER observations with Rossi X-ray Timing Explorer observations of the 2007 and 2009 outbursts, we alsostudied the long-term spin and orbital evolution of this source. We find that the binary system is well describedby a constant orbital period model, with an upper limit on the orbital period derivative of Pb < 7.4 ´ 10-13 ss1.Additionally, we improve upon the source coordinates through astrometric analysis of the pulse arrival times,finding R.A.=17h56m57 18±0 08 and decl.=25°0627 8±3 5, while simultaneously measuring thelong-term spin frequency derivative as n = -7.3 ´ 10-16 Hzs1. We briefly discuss the implications of thesemeasurements in the context of the wider population of accreting millisecond pulsars. We reported on the coherent timing analysis of the 2018 outburst of Swift J1756 as observed with NICER. Consistent with analyses of the previous outbursts (Krimm et al. 2007b; Patruno et al. 2010), we find that the X-ray pulsations have energy dependent amplitudes; the fractional amplitude of the fundamental increases with energy, whereas the fractional amplitude of the harmonic shows a slight decline with energy. This energy dependent behavior is not unusual in AMXPs (Patruno & Watts 2012) and can be interpreted in terms of the thermal emission from the stellar hotspot and reprocessing in the accretion column (e.g., Gierliński et al. 2002; Ibragimov & Poutanen 2009). The pulse arrival times of the 2018 outburst are well described by a timing model consisting of a circular orbit with a constant spin frequency. The pulse phases with respect to this model do not show spurious residuals with time or orbital phase, and no evidence is found that the pulse arrival times exhibit an additional delay associated with passing through the gravitational well of the companion star (Shapiro delay). We note, however, that the expected Shapiro delay is given as (Shapiro et al. 1971) Equation (5) where Φ is the orbital phase, G is the gravitational constant, c is the speed of light, and i is the inclination. Even for the maximum allowed companion mass, ${M}_{C}=0.030\,{M}_{\odot }$ (Krimm et al. 2007b, but see Section 4.2 for more details) and an inclination of 90°, the largest delay we can expect is only 4 μs. As this time-delay is smaller than the uncertainty on our phase residuals by nearly two orders of magnitude (see Figure 1), we are not sensitive to Shapiro delays in Swift J1756. Comparing our measurements for the 2018 outburst with those of the 2007 and 2009 outbursts as observed with RXTE, we analyzed the long-term evolution of this source. We found that the binary system is consistent with having a constant orbital period and that the pulsar shows a spin frequency derivative of $\dot{\nu }=-7.3\times {10}^{-16}\,\mathrm{Hz}\,{{\rm{s}}}^{-1}$. 4.1. Spin-down Evolution The long-term spin frequency derivative measured in Swift J1756 is of the same order as the spin frequency derivatives measured in other AMXPs (Hartman et al. 2008; Patruno 2010; Riggio et al. 2011). This frequency change is most likely driven by the neutron star's loss of rotational energy. If so, then the spin-down luminosity is given as Equation (6) where I represents the neutron star moment of inertia. The long-term spin-down of a neutron star is usually assumed to be dominated by the braking torque associated with a spinning magnetic field. Assuming this mechanism is responsible for the observed spin-down in Swift J1756, we can compute the magnetic dipole moment as (Spitkovsky 2006) Equation (7) where α is the misalignment angle between the rotational and magnetic poles. Considering α = 0°–90°, we then find a magnetic field strength of $B\simeq (4\mbox{--}6)\times {10}^{8}$ G at the stellar magnetic poles. This magnetic field strength estimate is in line with those obtained for other accreting millisecond pulsars (see Mukherjee et al. 2015 and references therein). 4.2. Orbit Evolution The observed long-term binary evolution of Swift J1756 is consistent with this source having a constant orbital period and a lower limit on the evolutionary timescale of Equation (8) Binary evolution theory predicts that systems of this type evolve due to angular momentum loss through gravitational radiation (Kraft et al. 1962; Rappaport et al. 1982; Verbunt 1993). For conservative mass transfer, the binary period derivative is given by di Salvo et al. (2008), Equation (9) where MNS is the neutron star mass, $q={M}_{C}/{M}_{\mathrm{NS}}$ is the binary mass ratio, and −1/3 < n < 1 is the mass–radius index of the companion star. Depending on the source inclination, Krimm et al. (2007b) derived a companion mass of ${M}_{C}\,=0.007\mbox{--}0.022\,{M}_{\odot }$ for a neutron star mass of 1.4 ${M}_{\odot }$. For a neutron star mass of 2.2 ${M}_{\odot }$, the allowed range increased to ${M}_{C}=0.009\mbox{--}0.030\,{M}_{\odot }$. In both cases, they assumed an upper limit on the inclination of i < 85°, motivated by the fact that Swift J1756 does not show eclipses in its light curve. Accounting for the extreme cases of stellar masses and n, the binary may either be contracting or expanding. In either case, however, the rate of change is limited to $| {\dot{P}}_{b}| \lesssim 7\times {10}^{-14}$ s s−1, which is well below the upper limit obtained in this work. Although the binary evolution timescale we obtain for Swift J1756 is consistent with theory, it is worth noting that this is not generally true for low-mass X-ray binaries (see Patruno et al. 2017, for a comprehensive discussion). The AMXP SAX J1808.4–3658, in particular, has been found to evolve on a much shorter timescale, with a first derivative on the orbital period of $3.5\times {10}^{-12}$ s s−1 (Hartman et al. 2008; Patruno et al. 2012; Sanna et al. 2017a). Two models have been proposed to explain this discrepancy: highly nonconservative mass transfer due to irradiation of the companion star by the pulsar (di Salvo et al. 2008; Burderi et al. 2009), and spin–orbit coupling in the companion star (Hartman et al. 2008, 2009). While the latter depends on the companion star, and may vary from source to source, the former should operate in all AMXPs (see also Patruno 2017; Sanna et al. 2017c), including Swift J1756. The spin-down luminosity impinging on the companion star can be estimated as Equation (10) where ${\dot{E}}_{\mathrm{abl}}$ is the ablation luminosity, RL2 is the Roche lobe radius of the companion (Eggleton 1983), and a the binary separation. The irradiation fraction is $f={\dot{E}}_{\mathrm{abl}}/{\dot{E}}_{\mathrm{sd}}$, which, accounting for the range of allowed neutron star and companion masses, evaluates to f = 0.15%–0.35%. The associated mass loss for the companion is given by Equation (11) such that, assuming an efficiency of η = 100%, ${\dot{M}}_{C}\,\sim -3\times {10}^{-10}\,{M}_{\odot }$ yr−1. The effect of this mass loss on the orbital period follows through the relation (Frank et al. 2002) Equation (12) giving a period derivative due to mass loss of ${\dot{P}}_{b,\mathrm{ML}}\,=5\times {10}^{-12}$ s s−1. This value is well above our limit on the period derivative. Hence, in order for this mechanism to be consistent with our observations of Swift J1756, the efficiency at which the companion star converts the incident luminosity into mass loss must be η < 15%. This value is very different from the 40% required in SAX J1808.4–3658 (Patruno et al. 2016) and is instead in line with the <5% efficiency determined for IGR J00291+5934 (Patruno 2017). This work was supported by NASA through the NICER mission and the Astrophysics Explorers Program, and made use of data and software provided by the High Energy Astrophysics Science Archive Research Center (HEASARC). P.B. was supported by an NPP fellowship at NASA Goddard Space Flight Center. D.A. acknowledges support from the Royal Society.
Bult, Peter↗
Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma
A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.
neoclassical transport↗
Equation of State Constraints from Nuclear Physics, Neutron Star Masses, and Future Moment of Inertia Measurements
We explore constraints on the equation of state (EOS) of neutron-rich matter based on microscopic calculations up to nuclear densities and observations of neutron stars. In a previous work we showed that predictions based on modern nuclear interactions derived within chiral effective field theory and the observation of two-solar-mass neutron stars result in a robust uncertainty range for neutron star radii and the EOS over a wide range of densities. In this work we extend this study, employing both the piecewise polytrope extension from Hebeler et al. as well as the speed of sound model of Greif et al., and show that moment of inertia measurements of neutron stars can significantly improve the constraints on the EOS and neutron star radii.
79 ASTRONOMY AND ASTROPHYSICS↗
Numerical investigations of low-density nozzle flow by solving the Boltzmann equation
A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.
Deng, Zheng-Tao↗