Engineering Papers⌕ Search

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The dependence of the dissipation on rotation

A form of the evolution equation for the dissipation showing the effects of rotation is derived directly from the Navier-Stokes equations. Rotation affects the dissipation if the flow is inhomogeneous or compressible. This is due to the fact that pressure scales with inverse Rossby number and certain pressure correlations usually neglected have to be retained when Re Ro is approximately equal to 1. The term reflecting the effects of rotation on the dissipation is a volume integral of the two-point velocity correlation similar in structure to the 'rapid-pressure' correlation appearing in the second-moment equations. Some of its properties are indicated.

Ristorcelli, J. R.↗
NASA NTRS1996

Non-Parabolic Hydrodynamic Formulations for the Simulation of Inhomogeneous Semiconductor Devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models cannot fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations or the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship ((hk)(exp 2)/2m = W(1 + alphaW). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(exp y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships; parabolic, Kane dispersion and power law dispersion.

Smith, A. W.↗
NASA NTRS1995

Non-parabolic hydrodynamic formulations for the simulation of inhomogeneous semiconductor devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models can not fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations of the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship (hk)(exp 2)/2m = W(1 + alpha(W)). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(sup y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships: parabolic, Kane dispersion, and power low dispersion.

Smith, Arlynn W.↗
NASA NTRS1991

Near-critical spherical accretion by neutron stars - General relativistic treatment

Time-independent spherical accretion by a neutron star is studied using general relativistic radiation hydrodynamics. Numerical integrations of the flow equations are presented. These show that when the luminosity is sufficiently close to (but below) the Eddington limit, the flow velocity increases with decreasing radius far from the neutron star, reaches a maximum at an intermediate radius, and decreases at small radii. A large fraction of the binding energy of the flow is transferred to the radiation through scattering before the flow strikes the surface of the neutron star. Following Miller's treatment of accretion at luminosites near the Eddington limit (which neglected general relativistic effects), analytic approximations for the decelerating phase of the flow's velocity profile are derived. The dependence of the solutions on the variable Eddington factor prescription chosen to close the radiation moment equations is also examined.

Park, Myeong-Gu↗
DOE OSTI2020

High‐order cubature rules for tetrahedra

Summary In this paper, we construct new high‐order numerical integration schemes for tetrahedra, with positive weights and integration points that are in the interior of the domain. The construction of cubature rules is a challenging problem, which requires the solution of strongly nonlinear algebraic (moment) equations with side conditions given by affine inequality constraints. We present a robust algorithm based on a sequence of three modified Newton procedures to solve the constrained minimization problem. In the literature, numerical integration rules for the tetrahedron are available up to order p =15. We obtain integration rules for the tetrahedron from p =2 to p =20, which are computed using multiprecision arithmetic. For p ≤15, our approach provides integration rules that have the same or fewer number of integration points than existing rules; for p =16 to p =20, our rules are new. Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, rational, weakly singular and trigonometric test functions over tetrahedra with flat and curved faces. In all tests, improvements in accuracy is realized as p is increased, though in some cases nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2023

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1971

Theory of the unmagnetized plasma.

The Vlasov mathematical model of a plasma, which has come to be thought more useful than any other in describing the dynamical behavior of the majority of plasmas of interest, is first examined. Macroscopic variables and moment equations; linear electrostatics solutions; plasma oscillations, ion acoustic waves, and linear instabilities are treated, as well as external fields, 'test' charges, and nonlinear Vlasov phenomena. Plasmas are statistically described, and attention is given to the kinetic theory of the stable, uniform plasma and the Balescu-Lenard equation; two-time ensemble averages and fluctuation spectra in stable plasmas; the kinetic theory of the unstable plasma; and ensembles of Vlasov plasmas. Some illustrative experiments are described. Four appendixes deal with the electrostatic approximation and transverse waves; solution of the linearized Vlasov equation in a magnetic field; estimates of correlation functions from thermal equilibrium; and equivalence of spatially uniform BBGKY and Klimontovich correlations.

Montgomery, D. C.↗
NASA NTRS1972

A solar-wind model including proton thermal anisotropy.

A model of the solar wind is presented which divides the interplanetary space into two regions. It assumes that the solar wind is one-fluid in an 'inner region' of radius less than 0.4 A, and two-fluid and collisionless in an 'outer region' of radius greater than 0.4 A. The second and third moment equations of the Vlasov equation together with conservation equations are used to govern the solar wind flow in the outer region. The model produces solutions for all macroscopic quantities as well as for the microscopic proton distribution function, both as functions of heliocentric distances up to 50 A. Calculated results at 1 A include the solar wind velocity, electron temperature, proton temperature, proton thermal anisotropy and magnetic field magnitude. They are in good agreement with observations.

Whang, Y. C.↗
NASA NTRS1972

Low-frequency macroscopic instabilities of fully ionized magnetoplasma

Studies are described of low-frequency quasi-static instabilities in a fully ionized plasma. The plasma is assumed to be immersed in a uniform magnetic field, and is either uniform or has a number density gradient perpendicular to the magnetic field. A moment equation description of the ion and electron dynamics is used; collisions are assumed to have a strong effect on electron motion along the magnetic field. Before considering specific modes, a stability analysis is developed which allows a classification of wave growth characteristics to be made for a bounded system from solutions to the dispersion relation for an infinite system. Also, a method is given for calculating the normal mode frequencies and wave profiles by using the reflection coefficients at the boundaries. For wave propagation perpendicular to the magnetic field, the flute wave is studied in cylindrical geometry. The destabilizing effect of a radial electric field is considered by solving a differential equation.

Rognlien, T. D.↗
NASA NTRS1974

The stability of motion of satellites with cavities partially filled with liquid

The purpose of the present paper is to determine rigorously the effect of viscous dissipation on the stability and time-dependent motion of a spinning rigid body with a spherical cavity partially filled with liquid. The cavity is such that its center does not lie on the body axis of rotation. The work includes the problem formulation, consisting of the boundary-value problem for the liquid and moment equations for the entire system. Because of large Reynold's numbers involved, viscosity effects are negligible everywhere except for a thin boundary layer near the wetted surface. Using a boundary-layer analysis, the liquid problem is reduced to the solution of the inviscid equations subject to modified boundary conditions.

Nayfeh, A. H.↗
NASA NTRS1975

The stability of motion of satellites with cavities partially filled with liquid

The stability and time dependent motion of a spinning satellite, simulated by a rigid body with a cavity partially filled with liquid is examined. The problem formulation, consisting of the boundary-value problem for the liquid and moment equations for the entire system is presented. Because of large Reynold's numbers involved, viscosity effects are negligible everywhere except for a thin boundary layer near the wetted surface. Using a boundary-layer analysis, the effect of the boundary layer is replaced by modified boundary conditions for the liquid. The solution of the differential equations for the inviscid problem is solved in closed form. A semi-analytical numerical solution of the inviscid equations subject to the viscous boundary condition has proved unsucessful.

Nayfeh, A. H.↗
NASA NTRS1975

Strong scintillations in astrophysics. III - The fluctuations in intensity

Fluctuations in intensity due to wave propagation in a wide range of turbulent media, are discussed on the basis of a previously developed theory of strong scintillations. An asymptotic analytic solution to the fourth-moment equation is presented, and its range of validity is discussed. Various parameters characterizing the intensity fluctuations, such as scintillation index and correlation scale, are expressed in terms of turbulence parameters for both Gaussian and power-law irregularity spectra as well as for thick and thin turbulent media.

Lee, L. C.↗
NASA NTRS1976

Analysis of bilinear noise models in circuits and devices

There are a number of applications in which linear noise models are inappropriate. In the paper, the use of bilinear noise models in circuits and devices is considered. Several physical problems are studied in this framework. These include circuits involving varying parameters (such as variable resistance circuits constructed using field-effect transistors), the effect of switching jitter on sampled data system performance and communication systems involving voltage-controlled oscillators and phase-lock loops. In addition, several types of analytical techniques for stochastic bilinear systems are considered. Specifically, the moment equations of Brockett for bilinear systems driven by white noise are discussed, and closed-form expressions for certain bilinear systems (those that evolve an Abelian or solvable Lie groups) driven by white or colored noise are derived. In addition, an approximate statistical technique involving the use of harmonic expansions is described.

Willsky, A. S.↗
NASA NTRS1977

Prediction of Skylab attitude motion for space shuttle revisit

Computer simulations were employed to generate results, based on a fairly complicated analytical model. Environmental effects considered include atmospheric drag and gravity gradients. Mathematical models were based on Euler's moment equations and transformations from the body frame to inertial frame. A range of altitudes was considered, from 278 kilometers to 417 kilometers. Skylab is predicted to be in this orbital range during the 1979-1980 period. Results indicated that there will be a cyclic attitude motion due to the combined effect of gravity gradient and atmosphere drag. This motion can be characterized as being a slow tumble.

Alexander, S. G.↗
NASA NTRS1977

Kinetic theory analysis of rarefied gas flow through finite length slots

A kinetic-theory analysis is made of the flow of a rarefied monatomic gas through a two-dimensional slot connecting two reservoirs. Numerical solutions are obtained by the moment and discrete-ordinate methods. The former method portrays the transition-regime characteristics well but has limitations in the free-molecule regime. The latter method gives accurate results in the free-molecule and slip regimes and bolsters confidence in the accuracy of the transition-regime results. The numerical solution for the mass flux through the slot agrees well with an approximate analytical solution of the moment equations for length-to-width ratios from 6 to 0.5, pressure ratios from 0.8 to 0.1, and Knudsen numbers from 5 to 0.5.

Raghuraman, P.↗
NASA NTRS1978

Analysis of dynamic system response to product random processes

The response of dynamic systems to the product of two independent Gaussian random processes is developed by use of the Fokker-Planck and associated moment equations. The development is applied to the amplitude modulated process which is used to model atmospheric turbulence in aeronautical applications. The exact solution for the system response is compared with the solution obtained by the quasi-steady approximation which omits the dynamic properties of the random amplitude modulation. The quasi-steady approximation is valid as a limiting case of the exact solution for the dynamic response of linear systems to amplitude modulated processes. In the nonlimiting case the quasi-steady approximation can be invalid for dynamic systems with low damping.

Sidwell, K.↗