SEARCH · Search NASA

Results for “KROOK 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1990

A phase-space fluid simulation of a two-component narrow planetary ring - Particle size segregation, edge formation, and spreading rates

The Krook kinetic equation for identical planetary ring particles is presently generalized for the case of two-component systems, and the equations are numerically solved on the basis of Brophy and Esposito's (1989) phase-space CFD method. Attention is given to the simulation results obtained for a two-component narrow ring, in which the large particles are eight times as massive as the small particles. This ring's unconstrained edge dynamics are resolved by the simulation, and are found to exhibit a sharpening that would not have been expected in single-component rings.

Brophy, Thomas G.↗

NASA NTRS1992

Numerical simulation of satellite-ring interactions - Resonances and satellite-ring torques

The Krook kinetic equation for planetary rings is numerically solved in two spatial dimensions and in time, with (1) interparticle collisions and (2) satellite-forcing, but (3) without self-gravity, for the case of a flattened planetary ring that undergoes gravitational perturbation by a nearby satellite. It is noted that the amplitude of wakes is limited by purely kinematic effects, even in the absence of collisions. Attention is given to the results of a simulation of an inner Lindblad-resonance location, as the distribution approaches steady state; these simulations do not show an increase in velocity dispersion in the resonance zone, obviating a net torque.

Brophy, Thomas G.↗

NASA NTRS1973

Shock-wave structure in a fully ionized gas.

The structure of a steady plane shock in a fully ionized gas has been investigated using the Boltzmann equation with the Gross-Krook type model as the governing equation and the discrete-ordinate method as a tool. The present results agree well with the results obtained by the continuum approach for weak shocks. For strong shocks the present approach gives results that are considerably different from those yielded by the continuum approach, particularly in the high pressure region. Also, the present method gives smooth and continuous shock profiles for high Mach numbers. On the other hand, the results from the continuum approach contain discontinuities in slope. The ion and electron distribution functions across the shock are found to be singly peaked for all cases considered and the effects of the induced electric field is small.

Lu, C. S.↗

NASA NTRS1970

High frequency Hall current instability

The Hall current electrojet instability was studied to determine if it generates VLF hiss in the auroral zone of magnetic substorms. Using the kinetic equations with particle conserving Krook collision operator for electrons and ions, and Poisson's equation, it is shown that in the low frequency region, the range of frequency of unstable modes increases as the electron density increases, and that the high frequency component must be taken into account in calculations of turbulent saturation of the Hall current instability driven beyond threshold. It is concluded that Hall currents can generate waves in the VLF hiss frequency band which are related to local auroral electrojet activity.

Lee, K.↗

DOE OSTI2021

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1973

Test of statistical models for gases with and without internal energy states.

The problem of nonlinear rarefied Couette flow with heat transfer has been studied for both monatomic and diatomic gases using the Boltzmann equation with the Bhatnagar-Gross-Krook type models as the governing equation and the method of discrete ordinates as a tool. The calculated results have been compared with the existing experimental data in order to test the accuracy and the applicability of the statistical models for this one-dimensional problem. The calculated density results are found to be in good agreement with available experimental data; the calculated heat flux solution for the linear case is found to always be lower than the experimental data of Teagan and Springer. The comparisons made here indicate that the statistical models are indeed reasonably accurate so that their use is justified in the type of problems investigated.

Huang, A. B.↗

NASA NTRS1973

High-speed leading edge problem.

The sharp leading edge problem has been studied for both monatomic and diatomic gases using the Boltzmann equation with the Bhatnagar-Gross-Krook type models as the governing equation and the discrete ordinate method with a closed-boundary value approach as a tool. Plate length relative to the freestream mean free path is taken to be 52. The gas-surface interaction law is assumed to be diffuse reflection. The local distribution functions of molecular velocities and internal energies (for the diatomic gas) for the entire flowfield have been calculated for a freestream Mach number of 6.1. Comparisons are made between the calculated results and experimental data.

Huang, A. B.↗

NASA NTRS1979

A theory of local and global processes which affect solar wind electrons. I - The origin of typical 1 AU velocity distribution functions - Steady state theory

A kinetic theory for the velocity distribution of solar wind electrons which illustrates the global and local properties of the solar wind expansion is proposed. By means of the Boltzmann equation with the Krook collision operator accounting for Coulomb collisions, it is found that Coulomb collisions determine the population and shape of the electron distribution function in both the thermal and suprathermal energy regimes. For suprathermal electrons, the cumulative effects of Coulomb interactions are shown to take place on the scale of the heliosphere itself, whereas the Coulomb interactions of thermal electrons occur on a local scale near the point of observation (1 AU). The bifurcation of the electron distribution between thermal and suprathermal electrons is localized to the deep solar corona (1 to 10 solar radii).

Scudder, J. D.↗

NASA NTRS1979

High-altitude gas releases - Transition from collisionless flow to diffusive flow in a nonuniform atmosphere

The paper is motivated by a need for describing high-altitude gas releases from rockets or the Space Shuttle. In the tenuous upper atmosphere, the injected gases expand from a collisionless to a collision-dominated state. The analysis presented extends the efforts of Baum (1973-1974) to include the effects of a nonuniform background atmosphere and a velocity-dependent collision frequency. The unsteady expansion of gas releases is analyzed using gas kinetic theory. The Boltzmann equation with the Krook collision integral is solved numerically. At late times (after many collisions with the background atmosphere) the solution is identical to the one given by diffusion in a nonuniform atmosphere (Bernhardt, 1979). The theoretical model predicts elongation and heating of the vapor trails due to collisions.

Bernhardt, P. A.↗

NASA NTRS1968

Small perturbation method in kinetic theory.

Oseen equations derived for compressible flow by applying small perturbation method and Chapman- Enskog procedure to Bhatnagar-Gross-Krook model of Boltzmann equation

Chi, L. K.↗

NASA NTRS1993

Assessment of three numerical methods for the computation of a low-density plume flow

Results from three numerical methods including one based on the Navier-Stokes equations, one based on kinetic theory using the DSMC method, and one based on the Boltzmann equation with a Krook-type collision term are compared to each other and to experimental data for a model problem of heated nitrogen flow in a conical nozzle expanding into a vacuum. The problem simulates flow in a resistojet, a low-thrust, electrothermal rocket. The continuum method is applied to both the internal flow and near-field plume. The DSMC and Boltzmann methods are applied primarily to the plume. Experimental measurements of Pitot pressure and flow angle, taken with an apparatus that duplicates the model nozzle flow, are used in the comparisons.

Penko, Paul F.↗

DOE OSTI2024

Kinetic closures for unmagnetized and magnetized plasmas

Parallel and perpendicular closures with cyclotron resonance effects retained for the five-moment (density, temperature, and flow velocity) fluid equations are derived by solving the kinetic equation with the Bhatnagar–Gross–Krook operator in Fourier space. For parallel propagation, the parallel closures are reduced to those of Ji et al. [Phys. Plasmas 20, 082121 (2013)]. The closures when combined to the fluid equations reproduce the fully kinetic dispersion relation that can be directly derived from the kinetic equation. The closures for the five-moment fluid system can be utilized to derive closures for the extended fluid system, which is demonstrated by deriving closures for the ten-moment system consisting of density, flow velocity, temperature, and viscosity tensor equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2000

A gas-kinetic BGK scheme for the compressible Navier-Stokes equations

This paper presents an improved gas-kinetic scheme based on the Bhatnagar-Gross-Krook (BGK) model for the compressible Navier-Stokes equations. The current method extends the previous gas-kinetic Navier-Stokes solver developed by Xu and Prendergast by implementing a general nonequilibrium state to represent the gas distribution function at the beginning of each time step. As a result, the requirement in the previous scheme, such as the particle collision time being less than the time step for the validity of the BGK Navier-Stokes solution, is removed. Therefore, the applicable regime of the current method is much enlarged and the Navier-Stokes solution can be obtained accurately regardless of the ratio between the collision time and the time step. The gas-kinetic Navier-Stokes solver developed by Chou and Baganoff is the limiting case of the current method, and it is valid only under such a limiting condition. Also, in this paper, the appropriate implementation of boundary condition for the kinetic scheme, different kinetic limiting cases, and the Prandtl number fix are presented. The connection among artificial dissipative central schemes, Godunov-type schemes, and the gas-kinetic BGK method is discussed. Many numerical tests are included to validate the current method.

Xu, Kun↗

DOE OSTI2024

Semi-Lagrangian nodal discontinuous Galerkin method for the BGK model

In this work, we propose a semi-Lagrangian (SL) nodal discontinuous Galerkin (DG) solver for the BGK equation. The BGK model was introduced by Bhatnagar, Gross, and Krook [1] as a relaxation model for the fundamental Boltzmann equation [5], which describes the kinetic dynamic of rarefied gases with a probability distribution function. The challenges of designing efficient numerical schemes for the Boltzmann equation mainly come from its high dimensionality and complicated nonlinear collision operator. The BGK model gains interests since it has much lower computational cost, due to the relatively simple structure of the relaxation operator in replacement of the collision operator, while simultaneously preserving several important physical properties, such as macroscopic quantities and dissipation of entropy.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1992

Numerical solution of a three-dimensional cubic cavity flow by using the Boltzmann equation

A three-dimensional cubic cavity flow has been analyzed for diatomic gases by using the Boltzmann equation with the Bhatnagar-Gross-Krook (B-G-K) model. The method of discrete ordinate was applied, and the diffuse reflection boundary condition was assumed. The results, which show a consistent trend toward the Navier-Stokes solution as the Knudson number is reduced, give us confidence to apply the method to a three-dimensional geometry for practical predictions of rarefied-flow characteristics. The CPU time and the main memory required for a three-dimensional geometry using this method seem reasonable.

Hwang, Danny P.↗

DOE OSTI2024

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗

NASA NTRS1971

Evaluation of two statistical models using the shock structure problem.

The accuracy of two statistical models for the collision integral of the Boltzmann equation has been evaluated by applying the models to the solution of the problem of shock structure in a monatomic gas and then comparing the theoretical results with available ex perimental data. The two models considered here are the Bhatnagar-Gross-Krook and ellipsoidal statistical models. The Mach number range covered is 1.59-10.7 and profiles for density and, where available, temperature are compared. The method of numerical solution is the discrete ordinate technique which looks quite promising for application to more complicated models. The results indicate that the ellipsoidal statistical model, which gives a correct value for the Prandtl number, gives accurate results for a low Mach number shock. However, the accuracy degenerates as the Mach number increases. The Bhatnagar-Gross-Krook model gives poorer agreement with experimental data in all cases examined.

Giddens, D. P.↗