The discrete ordinate method for unsteady linearized Boltzmann Bhatnagar-Gross-Krook equation.
Discrete ordinate method for unsteady linearized Boltzmann Bhatnagar-Gross-Krook equation of rarefied gas flow
SEARCH · Search NASA
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.
Discrete ordinate method for unsteady linearized Boltzmann Bhatnagar-Gross-Krook equation of rarefied gas flow
Iterative solution to Krook-Boltzmann kinetic equation noting Navier-Stokes numerical solution, computation and analysis of distribution function within shock wave
Numerical convergence of iterative solution to complete Bhatnagar-Gross-Krook equation for strong shock waves
This article provides a concise exposition of the multiple-relaxation-time lattice Boltzmann equation, with examples of fifteen-velocity and nineteen-velocity models in three dimensions. Simulation of a diagonally lid-driven cavity flow in three dimensions at Re=500 and 2000 is performed. The results clearly demonstrate the superior numerical stability of the multiple-relaxation-time lattice Boltzmann equation over the popular lattice Bhatnagar-Gross-Krook equation.
Slip coefficient of gas calculated from linearized Boltzmann-Bhatnagar-Gross-Krook equation for slip velocity problem
It is shown that the use of the Krook equation greatly simplifies the discussion of the collisional dynamics of particulate disks. Collisions are assumed to be inelastic from the outset. A simple heuristic argument is used to compute the effective collision rate in a disk of spherical particles with a power-law distribution of sizes. For Saturn's rings, the effective collision rate for momentum transport is substantially lower than that conventionally estimated on the basis of an observed optical depth at visual wavelengths. The closure of the vertically integrated set of moment equations without discarding the third-order moments at the outset is discussed, allowing for the possibility of a bent disk. In the limit that the collision frequency is much larger than the orbit frequency, the usual Navier-Stokes equations of viscous hydrodynamics for a thin disk are recovered with an explicit expression for the shear viscosity. The method is generalized to include the effects of gravitational scattering.
A quantitative formalism for Keplerian particulate disk dynamics having its basis in a heuristic description of viscous transport is presently applied to the case of a ring with optical thickness radial gradient. A steady-state velocity distribution solution directly yields both the radial mass transport and the viscosity. The analytical method employed involves the solution of a novel form of the Krook equation through the separation of the phase-space collisional-products distribution into a symmetrical component and a delta function-approachable remainder. While this model is simplified for the case treated, its general approach may be extended to less artificially restricted cases.
Gas dynamics - structure of strong shock waves in krook collision equation
Steady and unsteady state problems and shock wave structure using Krook model of Boltzmann equation
A system of nonlinear integral equations equivalent to the Krook kinetic equation for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust two-dimensional nozzles. The method of characteristics was used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it exhausts into a vacuum. Results predict backscatter and show the effect of the inside wall boundary layer on the flowfields external to the nozzle.
A system of nonlinear integral equations equivalent to the Krook kinetic equations for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust three-dimensional axisymmetric nozzles. The method of characteristics is used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it expands into a vacuum. Results predict backscatter and show the effect of the nozzle wall boundary layer on the external flowfields.
The back-scattering of outgas contamination near an orbiting spacecraft due to intermolecular collisions was analyzed. Analytical tools were developed for making reasonably accurate quantitative estimates of the outgas contamination return flux, given a knowledge of the pertinent spacecraft and orbit conditions. Two basic collision mechanisms were considered: (1) collisions involving only outgas molecules (self-scattering) and (2) collisions between outgas molecules and molecules in the ambient atmosphere (ambient-scattering). For simplicity, the geometry was idealized to a uniformly outgassing sphere and to a disk oriented normal to the freestream. The method of solution involved an integration of an approximation of the Boltzmann kinetic equation known as the BGK (or Krook) model equation. Results were obtained in the form of simple equations relating outgas return flux to spacecraft and orbit parameters. Results were compared with previous analyses based on more simplistic models of the collision processes.
A system of nonlinear integral equations equivalent to the steady-state Krook kinetic equation was used to model the flow from a low-thrust axisymmetric nozzle. The mathematical model was used to numerically calculate the number density, temperature, and velocity of a simple gas as it expands into a near vacuum. With these quantities the gas pressure and flow directions of the gas near the exit plane were calculated and compared with experimental values for a low-thrust nozzle of the same geometry and mass flow rate.
The thermal and suprathermal processes involving solar wind electrons are discussed from a theoretical point of view. A model for the electron distribution function, f(e), based on the solutions of the Boltzmann equation in Krook's approximation is outlined: the angular and energy dependences of for various distances from the Sun between the coronal base and the Earth are presented with the express purpose of obtaining the radial profile of the heat flux generated by the suprathermals. The basic fluid-dynamical equations for the solar wind plasma as a whole along any given magnetic field line are solved. Results are in good agreement with observations. The predicted density and temperature profiles for positive ions exhibit the characteristics of recent measurements, both in the corona (above the coronal holes) and between 0.3 and 1 A.U. (the Helios region). It is concluded that the suprathermal electrons play an important role in the formation and the dynamics of the solar wind.
Shock wave profile computation using krook collision model equation in iteration scheme starting from navier-stokes solution
In this article we present a class of high order unconditionally strong stability preserving (SSP) implicit two-derivative Runge--Kutta schemes and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$ and IMEX Runge--Kutta schemes of order up to $p=3$. For the multiderivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar--Gross--Krook kinetic equation. We present numerical results to support the theoretical results on a variety of problems.
In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.
Numerical solution of Bhatnagar-Gross-Krook model for shock wave structure problem, noting Chapman-Enskog approximation