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 · Search NASA
Results for “MOMENT EQUATION”
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.
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.↗
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.↗
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.↗
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↗
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↗
Electromagnetic Field Penetration Studies
A numerical method is presented to determine electromagnetic shielding effectiveness of rectangular enclosure with apertures on its wall used for input and output connections, control panels, visual-access windows, ventilation panels, etc. Expressing EM fields in terms of cavity Green's function inside the enclosure and the free space Green's function outside the enclosure, integral equations with aperture tangential electric fields as unknown variables are obtained by enforcing the continuity of tangential electric and magnetic fields across the apertures. Using the Method of Moments, the integral equations are solved for unknown aperture fields. From these aperture fields, the EM field inside a rectangular enclosure due to external electromagnetic sources are determined. Numerical results on electric field shielding of a rectangular cavity with a thin rectangular slot obtained using the present method are compared with the results obtained using simple transmission line technique for code validation. The present technique is applied to determine field penetration inside a Boeing-757 by approximating its passenger cabin as a rectangular cavity filled with a homogeneous medium and its passenger windows by rectangular apertures. Preliminary results for, two windows, one on each side of fuselage were considered. Numerical results for Boeing-757 at frequencies 26 MHz, 171-175 MHz, and 428-432 MHz are presented.
Deshpande, M.D.↗
Moment method with isoparametric elements for three-dimensional anisotropic scatterers
A novel method for computing the frequency-domain electromagnetic fields scattered from, and penetrating into, arbitrarily shaped, three-dimensional, lossy, inhomogeneous anisotropic scatters is presented. The method is based on a general volume integrodifferential formulation of the scattering problem and consists of the numerical solution of the coupled integral equations by the moment method and point matching. The numerical model of the scatterer is obtained by parametric volume elements, and the basis functions used to represent the field within each element are the same used in the finite-element method. Element integration problems due to the singular kernel of the integral equations are treated in some detail. Numerical results for both the isotropic and the anisotropic spherical scatterer are presented, including comparisons with results obtained by different numerical methods for the isotropic cases considered. The capability of the numerical code presented to deal with cases where the material parameters of the scatterer are given by singular matrices is discussed for two particular examples.
Graglia, Roberto D.↗
A moment method solution of a volume-surface integral equation using isoparametric elements and point matching
Three integral equations for TE (transverse electric) scattering by dielectric cylinders having large values of permittivity are examined. A moment-method solution of the volume-surface integral equation is developed employing isoparametric elements and point matching. The solution is shown to be more accurate and stable than the traditional solutions using pulse basis, particularly in the case of scatterers having large refractive indices.
Jin, Jian-Ming↗
Applications of two-dimensional integral equation theory to reflector antenna analysis
Applying method of moments to solution of integral equations for current induced on conducting cylinders in two-dimensional field
Rusch, W. V. T.↗
Solution of Euler's Equations of Motion and Eulerian Angles for near symmetric rigid bodies subject to constant moments
Analytic expressions are found for Euler's Equations of Motion and for the Eulerian Angles for both symmetric and near symmetric rigid bodies under the influence of arbitrary constant body-fixed torques. These solutions provide the body-fixed angular velocities and the attitude of the body, respectively, as functions of time. They are of special interest in applications to spinning spacecraft (such as the Galileo Spacecraft to be launched in 1984) because they include the effect of time-varying spin rate. Thus they can be applied to spin-up and spin-down maneuvers as well as to error analysis for thruster misalignments. The solutions are given for arbitrary initial conditions in terms of Fresnel, Sine and Cosine Integrals. Numerical integration of the governing differential equations has verified that the approximate analytic solutions are very accurate in many physical situations of interest.
Longuski, J. M.↗
On the joint bimodality of temperature and moisture near stratocumulus cloud tops
The observed distributions of the thermodynamic variables near stratocumulus top are highly bimodal. Two simple models of sub-grid fractional cloudiness motivated by this observed bimodality are examined. In both models, certain low order moments of two independent, moist-conservative thermodynamic variables are assumed to be known. The first model is based on the assumption of two discrete populations of parcels: a warm-day population and a cool-moist population. If only the first and second moments are assumed to be known, the number of unknowns exceeds the number of independent equations. If the third moments are assumed to be known as well, the number of independent equations exceeds the number of unknowns. The second model is based on the assumption of a continuous joint bimodal distribution of parcels, obtained as the weighted sum of two binormal distributions. For this model, the third moments are used to obtain 9 independent nonlinear algebraic equations in 11 unknowns. Two additional equations are needed to determine the covariance within the two subpopulations. In case these two internal covariance vanish, the system of equations can be solved analytically.
Randall, D. A.↗
Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes
Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.
Groth, Clinton P. T.↗