SEARCH · Engineering Papers

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

32 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Theory of the Lattice Boltzmann Equation: Symmetry properties of Discrete Velocity Sets

In the lattice Boltzmann equation, continuous particle velocity space is replaced by a finite dimensional discrete set. The number of linearly independent velocity moments in a lattice Boltzmann model cannot exceed the number of discrete velocities. Thus, finite dimensionality introduces linear dependencies among the moments that do not exist in the exact continuous theory. Given a discrete velocity set, it is important to know to exactly what order moments are free of these dependencies. Elementary group theory is applied to the solution of this problem. It is found that by decomposing the velocity set into subsets that transform among themselves under an appropriate symmetry group, it becomes relatively straightforward to assess the behavior of moments in the theory. The construction of some standard two- and three-dimensional models is reviewed from this viewpoint, and procedures for constructing some new higher dimensional models are suggested.

Rubinstein, Robert↗

NASA NTRS1998

Direct Acceleration of Pickup Ions at The Solar Wind Termination Shock: The Production of Anomalous Cosmic Rays

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

Ellison, Donald C.↗

NASA NTRS2011

Gas Flux and Density Surrounding a Cylindrical Aperture in the Free Molecular Flow Regime

The equations for rigorously calculating the particle flux and density surrounding a cylindrical aperture in the free molecular flow regime are developed and presented. The fundamental equations for particle flux and density from a reservoir and a diffusely reflecting surface will initially be developed. Assumptions will include a Maxwell-Boltzmann speed distribution, equal particle and wall temperatures, and a linear flux distribution along the cylindrical aperture walls. With this information, the equations for axial flux and density surrounding a cylindrical aperture will be developed. The cylindrical aperture will be divided into multiple volumes and regions to rigorously determine the surrounding axial flux and density, and appropriate limits of integration will be determined. The results of these equations will then be evaluated. The linear wall flux distribution assumption will be assessed. The axial flux and density surrounding a cylindrical aperture with a thickness-to-radius ratio of 1.25 will be presented. Finally, the equations determined in this study will be verified using multiple methods.

Soulas, George C.↗

NASA NTRS2002

Lattice Boltzmann Equation On a 2D Rectangular Grid

We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitras~y inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.

Bouzidi, MHamed↗

NASA NTRS1961

Consistent Pl Analysis of Aqueous Uranium-235 Critical Assemblies

The lethargy-dependent equations of the consistent Pl approximation to the Boltzmann transport equation for slowing down neutrons have been used as the basis of an IBM 704 computer program. Some of the effects included are (1) linearly anisotropic center of mass elastic scattering, (2) heavy element inelastic scattering based on the evaporation model of the nucleus, and (3) optional variation of the buckling with lethargy. The microscopic cross-section data developed for this program covered 473 lethargy points from lethargy u = 0 (10 Mev) to u = 19.8 (0.025 ev). The value of the fission neutron age in water calculated here is 26.5 square centimeters; this value is to be compared with the recent experimental value given as 27.86 square centimeters. The Fourier transform of the slowing-down kernel for water to indium resonance energy calculated here compared well with the Fourier transform of the kernel for water as measured by Hill, Roberts, and Fitch. This method of calculation has been applied to uranyl fluoride - water solution critical assemblies. Theoretical results established for both unreflected and fully reflected critical assemblies have been compared with available experimental data. The theoretical buckling curve derived as a function of the hydrogen to uranium-235 atom concentration for an energy-independent extrapolation distance was successful in predicting the critical heights of various unreflected cylindrical assemblies. The critical dimensions of fully water-reflected cylindrical assemblies were reasonably well predicted using the theoretical buckling curve and reflector savings for equivalent spherical assemblies.

Fieno, Daniel↗

NASA NTRS2000

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗

NASA NTRS1974

Acoustical theory of turbulence

Fluid functions are expanded in a series of functions constituting a complete and orthogonal system of wave solutions of the linearized Navier-Stokes system. The motion of characteristic waves is described by Hamiltonian equations of ray acoustics. Amplitude of the waves is described with the aid of nonlinear interaction terms. The distribution functions (squares of amplitudes) satisfy Boltzmann-type equations and completely describe the mean turbulence properties.

Kentzer, C. P.↗

NASA NTRS1996

Cross-Effects in Microgravity Flows

Film growth by chemical/physical vapor deposition is a process of considerable interest in microgravity experiments. The absence of natural convection should allow better control of film growth processes but, in highly non-isothermal ampoules, thermal slip (creep) can become a matter of significant concern. The reported research is a theoretical and experimental investigation of the flow of gas/vapor mixtures under non-continuum conditions. The Boltzmann equation has been solved for a monatomic gas under non-condensing conditions and the various phenomenological coefficients have been computed. Computations for realistic potentials as well as for velocity and creep slip have been completed and the creep slip has been found to be dependent on the type of gas confirming the accuracy of previous variational results. The variational technique has been extended and planar flows calculated via the Burnett solutions. Velocity, diffusion and creep slips have been computed for gas mixtures and previously unknown dependencies of the creep slip on the mixture properties have been observed. Also for gas mixtures, an integral representation of the linearized Boltzmann operator has been developed for use in numerical and variational calculations for all intermolecular force laws. Two, two-bulb capillary systems have been designed, built and tested for the measurements of cross-flows; one of glass for isothermal measurements and one of stainless steel for non-isothermal measurements. Extensive data have been collected for Ar-He and N2-He mixtures at a variety of pressures and mole ratios. Viscosity, velocity slip coefficients and tangential momentum accommodation coefficients have been obtained from measurements with a spinning rotor gauge via a new theory that has been formulated for the spinning rotor gauge in the slip regime. The FIDAP fluid dynamics code has been applied to condensing flows in ampoules in the continuum regime and agreement obtained with the earlier work of Duval.

Loyalka, Sudarshan K.↗

NASA NTRS1991

Lattice Boltzmann model for simulation of magnetohydrodynamics

A numerical method, based on a discrete Boltzmann equation, is presented for solving the equations of magnetohydrodynamics (MHD). The algorithm provides advantages similar to the cellular automaton method in that it is local and easily adapted to parallel computing environments. Because of much lower noise levels and less stringent requirements on lattice size, the method appears to be more competitive with traditional solution methods. Examples show that the model accurately reproduces both linear and nonlinear MHD phenomena.

Chen, Shiyi↗

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 NTRS1983

Ion-acoustic solitary waves in a magnetized plasma with arbitrary electron equation of state

The oblique propagation of fully nonlinear, slow ion-acoustic solitary waves in a collisionless, low-beta, magnetized plasma is examined. The analysis includes the effects of a finite ion pressure, electron trapping, and multicomponent particle populations. The existence of both compressive and rarefactive modes propagating obliquely to the magnetic field in a plasma with two distinct Boltzmann electron populations and cold ions is demonstrated. It is shown that paired electrostatic shocks observed over the earth's auroral zone may be closely related to the rarefactive modes. As a measure of the collisionless dissipation encountered by the solitary waves, the linear response of the plasma to slow ion-acoustic waves is also examined.

Witt, E.↗

NASA NTRS1996

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗

NASA NTRS1985

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗

NASA NTRS2004

A space radiation transport method development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest-order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard finite element method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 ms and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re-configurable computing and could be utilized in the final design as verification of the deterministic method optimized design. Published by Elsevier Ltd on behalf of COSPAR.

NASA Center JSC↗