Search NASA⌕ Search

SEARCH · Search NASA

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

At least 379 records · Page 21

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, where droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Computational Fluid Dynamics↗

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, wherein droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Self-Pressurization↗

Test-particle method in kinetic theory of a plasma.

The introduction of a test particle into a system is considered. The system may be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy. The field particles form a cloud which surrounds the test particle. The cloud is described by a conditional probability function which satisfies a certain equation. A generalization of the superposition principle reported by Rostoker (1964) to higher order correlation functions is discussed. Kinetic equations with the generalized Lenard-Balescu term are obtained, taking into account also diffusion by waves. The characteristics regarding the absorption or emission of waves by particles can be calculated.

Matsuda, K.↗

Kinetic effects on Alfven wave nonlinearity. II - The modified nonlinear wave equation

A previously developed Vlasov theory is used here to study the role of resonant particle and other kinetic effects on Alfven wave nonlinearity. A hybrid fluid-Vlasov equation approach is used to obtain a modified version of the derivative nonlinear Schroedinger equation. The differences between a scalar model for the plasma pressure and a tensor model are discussed. The susceptibilty of the modified nonlinear wave equation to modulational instability is studied. The modulational instability normally associated with the derivative nonlinear Schroedinger equation will, under most circumstances, be restricted to left circularly polarized waves. The nonlocal term in the modified nonlinear wave equation engenders a new modulational instability that is independent of beta and the sense of circular polarization. This new instability may explain the occurrence of wave packet steepening for all values of the plasma beta in the vicinity of the earth's bow shock.

Spangler, Steven R.↗

A second-order accurate kinetic-theory-based method for inviscid compressible flows

An upwind method for the numerical solution of the Euler equations is presented. This method, called the kinetic numerical method (KNM), is based on the fact that the Euler equations are moments of the Boltzmann equation of the kinetic theory of gases when the distribution function is Maxwellian. The KNM consists of two phases, the convection phase and the collision phase. The method is unconditionally stable and explicit. It is highly vectorizable and can be easily made total variation diminishing for the distribution function by a suitable choice of the interpolation strategy. The method is applied to a one-dimensional shock-propagation problem and to a two-dimensional shock-reflection problem.

Deshpande, Suresh M.↗

Stably stratified shear turbulence: A new model for the energy dissipation length scale

A model is presented to compute the turbulent kinetic energy dissipation length scale l(sub epsilon) in a stably stratified shear flow. The expression for l(sub epsilon) is derived from solving the spectral balance equation for the turbulent kinetic energy. The buoyancy spectrum entering such equation is constructed using a Lagrangian timescale with modifications due to stratification. The final result for l(sub epsilon) is given in algebraic form as a function of the Froude number Fr and the flux Richardson number R(sub f), l(sub epsilon) = l(sub epsilon)(Fr, R(sub f). The model predicts that for R(sub f) less than R(sub fc), l(sub epsilon) decreases with stratification. An attractive feature of the present model is that it encompasses, as special cases, some seemingly different models for l(sub epsilon) that have been proposed in the past by Deardorff, Hunt et al., Weinstock, and Canuto and Minotti. An alternative form for the dissipation rate epsilon is also discussed that may be useful when one uses a prognostic equation for the heat flux. The present model is applicable to subgrid-scale models, which are needed in large eddy simulations (LES), as well as to ensemble average models. The model is applied to predict the variation of l(sub epsilon) with height z in the planetary boundary layer. The resulting l(sub epsilon) versus z profile reproduces very closely the nonmonotonic profile of l(sub epsilon) exhibited by many LES calculations, beginning with the one by Deardorff in 1974.

Cheng, Y.↗

Theoretical studies of solar lasers and converters

The geometry and setup for the n-C3F7I iodine laser are illustrated. The mathematical modeling of this system is described. The chemical kinetics are summarized. A sensitivity analysis was performed on the parameters occurring in the differential equations describing the chemical kinetics.

Heinbockel, John↗

Cosmic ray propagation in interplanetary space

The validity of the test-particle picture, the approximation of static fields, and the spatial-diffusion approximation are discussed in a general way before specific technical assumptions are introduced. It is argued that the spatial-diffusion equation for the intensity per unit energy has a much wider range of applicability than the kinetic (Fokker-Planck) equation it is derived from. This gives strong weight to the phenomenological propagation theory. The general success (and possible failure at small energies) of the phenomenological theory for the modulation of galactic cosmic rays and solar events is described. Apparent effects such as the 'free boundary' are given disproportionate weight since they establish the connection with the detailed plasma physics of the solar wind. Greatest attention is paid to the pitch-angle diffusion theory. A general theory is presented which removes the well-known secularities of the quasi-linear approximation. The possible breakdown of any pitch-angle diffusion theory at very small energies is perhaps connected with the observed 'turn up' of the spectrum at low energies. A first attempt to derive the spatial dependence of the diffusion coefficient in the solar cavity, using such a divergence free scattering theory, is described and compared with recent observations out to 5 AU.

Voelk, H. J.↗

Anisotropic interface kinetics and tilted cells in unidirectional solidification

A nonlinear evolution equation is derived which governs the cellular structure of a binary alloy with a small segregation coefficient, including the effects of anisotropic interface kinetics. This equation, applicable to long-wave instabilities of a planar interface, describes the spatial pattern of the growing disturbances. The presence of anisotropy causes the cells to grow at an angle to the normal of the planar front. This transition to a cellular morphology is shown to be a subcritical instability.

Young, G. W.↗

Missile impact craters (White Sands Missile Range, New Mexico) and applications to lunar research: Contributions to astrogeology

Craters in natural materials at White Sands Missile Range, N. Mex., were produced by the impact of high-velocity to hypervelocity missiles traveling along oblique trajectories with kinetic energies between 2.1 and 81 × 1014 ergs. The oblique impacts produce craters 2 to 10 m across with morphologies and ejecta that are bilaterally symmetrical with respect to the plane of the missile trajectory. Rims are high and the amount of ejecta large in down-trajectory and lateral directions, whereas rims are low to nonexistent and ejecta thin to absent up-trajectory. Symmetry development and modifications of the symmetry are a function of target material, local topography, and angle of impact. Seven mappable units can be recognized in and around the craters. Three of these are ejecta: thick ejecta near the crater, thin to discontinuous ejecta at greater distances, and scattered ejecta at the greatest distances to the limit of throwout. These ejecta units may be absent on the up-trajectory side; if present, they are rarely as thick or continuous as on other sides of the crater. Three units are target materials: undeformed target material exposed in local patches through thin to discontinuous ejecta and everywhere between the fragments of scattered ejecta, tilted and broken target material exposed in upper crater walls, and shattered and fractured target material exposed on the up-trajectory crater wall. The seventh unit is slope material composed of talus and fallback within the crater. Development, character, and exposure of these units varies chiefly with the target material. Ejecta from the craters is chiefly broken but relatively undeformed target material that may range in size from very fine grained debris to large blocks. Where the target is porous, significant amounts of the ejecta are composed of sheared and compressed fragments, some coated with dark layers of mixed projectile pieces, powder, and fused metal mixed with crushed target material. For layered targets, the original stratigraphic sequence is crudely preserved and in inverted order in thick ejecta. Secondary impact craters are produced by the impact of ejected fragments when the surrounding surface materials are sufficiently weak. A wide variety of secondary impact crater relations may result. Secondary craters nearest the primary crater have blocks in them that are larger than or the same size as the crater they produced. Farther from the primary crater, the fragments are generally smaller than the secondary crater and are ejected from it. Excavation of four craters revealed a mixed breccia beneath the crater floor composed of missile pieces, sheared and compressed target material, and crushed debris. Banded, disaggregated target material and nonmixed breccia surrounded the mixed breccia, and these breccias were surrounded by a zone of conjugate fractures. Beneath the ejecta on the lateral and down-trajectory crater flanks, the target materials were tilted upward and broken. Up-trajectory, open fractures and downward displacement occurred in two of the craters. No displacement was observed for the other two. Beneath the down-trajectory rims of craters with distinct layering, overturned synclines were observed. Missile breakup and behavior during cratering are a function of target and missile properties. Missile breakup depends on missile velocity and is most extensive at high velocities, where the missile is fragmented, powdered, and partly fused. Burial of missile or its fragmented, powdered, and fused remains is greatest for porous targets and least for dense cohesive targets. For very porous targets, camouflet structures containing the fragmented missile may form. Least squares fit to the data on craters in dry to moist targets indicate V(a) = 10^(-11.433)E(p)^(1.205) where V(a) is the volume of the apparent crater and E(p) is the kinetic energy of the missile. This equation is consistent with expectations of the equations relating apparent depth and radius to kinetic energy. Extrapolation of displaced masses and kinetic energies for laboratory impacts with sand and rock converge near 10^(15) to 10^(16) ergs, where the extrapolations are near the data on missile impact craters, corrected for impact angle. Displaced masses of craters produced by missile impacts and by chemical explosives with small scaled depths of burial are about the same when the kinetic energies of the missiles (corrected for angle of impact) are equal to the TNT equivalent energy of the explosive. The problem of equivalent scaled depth of burst for an impact crater is complicated and not entirely resolved, however. Both missile impact craters and chemical explosive craters in water-saturated targets are larger than their counterparts in dry to moist materials. Data collected during the study of missile impact craters have helped resolve a number of problems in lunar research: (1) the soillike nature of lunar surface materials was predicted, (2) sizes of craters produced by artificial impacts were correctly predicted, (3) certain features imaged by Surveyor were found to be analogous to features associated with missile impact craters, {4) missile impacts were used in support of the Apollo passive seismic experiment, (5) craters seen in Apollo orbital photographs were found to be similar to some missile impact craters, (6) missile impact craters supplied data on sample collection and crater phenomenology used in training astronauts, and (7) some returned lunar samples are similar to coated, sheared, and compressed fragments ejected from missile impact craters.

H. J. Moore↗

Equations of motion and two-equation turbulence model for plane or axisymmetric turbulent flows in body-oriented orthogonal curvilinear coordinates and mass-averaged dependent variables

The full Navier-Stokes time-dependent, compressible, turbulent, mean-flow equations in mass-averaged variables for plane or axisymmetric flow are presented. The equations are derived in a body-oriented, orthogonal, curvilinear coordinate system. Turbulence is modelled by a system of two equations for mass-averaged turbulent kinetic energy and dissipation rate proposed. These equations are rederived and some new features are discussed. A system of second order boundary layer equations is then derived which includes the effects of longitudinal curvature and the normal pressure gradient. The Wilcox and Chambers approach is used in considering effects of streamline curvature on turbulence phenomena in turbulent boundary layer type flows. Their two-equation turbulence model with curvature terms are rederived for the cases considered in the present report. The derived system equations serves as a basis for an investigation of problems where streamline curvature is of the order of the characteristic length in the longitudinal direction.

Sislian, J. P.↗

Second-Order Modeling of Low-Reynolds-Number Turbulence Near Walls

This paper presents a set of second-order closure models for low-Reynolds-number turbulence near the wall. Existing closure models for the Reynolds-stress equations were modified to show proper near-wall behavior. A dissipation-rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The equations are solved for the mean velocity, the Reynolds stresses, and the dissipation rate of the turbulent kinetic energy. The calculations are compared with both direct numerical simulations and with measurements. It is shown that the present models perform well in predicting the behavior of the turbulence near a wall. Significant improvements over previous models in predicting the components of the Reynolds stress tensor are obtained in the present models.

Shih, T.-H.↗

Budget of Turbulent Kinetic Energy in a Shock Wave Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Quantities present in the exact equation of the turbulent kinetic energy transport were accumulated and used to calculate terms like production, dissipation, molecular diffusion, and turbulent transport. The present results for a turbulent boundary layer were validated by comparison with direct numerical simulation data. It was found that a longer development domain was necessary for the boundary layer to reach an equilibrium state and a finer mesh resolution would improve the predictions. In spite of these findings, trends of the present budget match closely with that of the direct numerical simulation. Budgets for the SBLI region are presented at key axial stations. These budgets showed interesting dynamics as the incoming boundary layer transforms and the terms of the turbulent kinetic energy budget change behavior within the interaction region.

shock wave/boundary-layer interaction↗

Stellar Mixing: I. Formalism

In this paper we use the Reynolds stress models (RSM) to derive algebraic expressions for the following variables: a) heat fluxes; b) J.l fluxes; and c) momentum fluxes. These relations, which are fully 3D, include: 1) stable and unstable stratification, represented by the Brunt-Vaislila frequency, N(exp 2) =-g/H(sub p_(del - del(sub ad))(1 - RI(sub mu)); 2) double diffusion, salt-fingers, and semi-convection, represented by the density ratio R(sub mu) = del(sub mu)/(del - del(sub ad)); 3) shear (differential rotation), represented by the mean squared shear Sigma(exp 2) or by the Richardson number, Ri =N(exp 2)Sigma(exp -2); 4) radiative losses represented by a Peclet number, Pe; 5) a complete analytical solution of the ID version of the model. In general, the model requires the solution of two differential equations for the eddy kinetic energy K and its rate of dissipation, epsilon. In the local and stationary cases, when production equals dissipation, the model equations are all algebraic.

Canuto, V .M.↗

NCC: A Physics-Based Design and Analysis Tool for Combustion Systems

The National Combustion Code (NCC) is an integrated system of computer codes for physics-based design and analysis of combustion systems. It uses unstructured meshes and runs on parallel computing platforms. The NCC is composed of a set of distinct yet closely related modules. They are: (1) a gaseous flow module solving 3-D Navier-Stokes equations; (2) a turbulence module containing the non-linear k-epsilon models; (3) a chemistry module using either the conventional reduced kinetics approach of solving species equations or the Intrinsic Low Dimensional Manifold (ILDM) kinetics approach of table looking up in conjunction with solving the equations of the progressive variables; (4) a turbulence-chemistry interaction module including the option of solving the joint probability density function (PDF) for species and enthalpy; and (5) a spray module for solving the liquid phase equations. In early 1995, an industry-government team was formed to develop the NCC. In July 1998, the baseline beta version was completed and presented in two NCC sessions at the 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, July 1998. An overview of this baseline beta version was presented at the NASA HPCCP/CAS Workshop 98, August 1998. Since then, the effort has been focused on the streamlining, validation, and enhancement of the th baseline beta version. The progress is presented in two NCC sessions at the AIAA 38 Aerospace Sciences Meeting & Exhibit, January 2000. At this NASA HPCCP/CAS Workshop 2000, an overview of the NCC papers presented at the AIAA 38 th Aerospace Sciences Meeting & Exhibit is presented, with emphasis on the reduction of analysis time of simulating the (gaseous) reacting flows in full combustors. In addition, results of NCC simulation of a modern turbofan combustor will also be reported.

Liu, Nan-Suey↗

Waveguide CO2 laser gain - Dependence on gas kinetic and discharge properties

A simple rate-equation approach is used to examine the gas-kinetic and discharge properties of waveguide CO2 lasers. The dependence of the population inversion and laser small-signal gain on gas pressure, gas mixture, pumping rate (discharge current), tube bore diameter, and wall temperature is calculated along with the dependence of laser tunability on the gas-kinetic properties and cavity losses. It is found that for low-loss cavities, the laser tunability may substantially exceed the molecular full width at half-maximum. Furthermore, the more helium-rich gas mixtures give greater tunability when cavity losses are small, and less tunability when cavity losses are large. By contrast with conventional lasers, the waveguide-laser transition is homogeneously broadened. Thus, the dependence of gain on gas pressure and other kinetic properties differs substantially from that predicted by scaling results from conventional low-pressure lasers.

Cohen, S. C.↗