Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity 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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗
NASA NTRS1992

FDDO and DSMC analyses of rarefied gas flow through 2D nozzles

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas expanding through a two-dimensional nozzle and into a surrounding low-density environment. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, the variable hard sphere model is used as a molecular model and the no time counter method is employed as a collision sampling technique. The results of both the FDDO and the DSMC methods show good agreement. The FDDO method requires less computational effort than the DSMC method by factors of 10 to 40 in CPU time, depending on the degree of rarefaction.

Chung, Chan-Hong↗
NASA NTRS1983

A 4-cylinder Stirling engine computer program with dynamic energy equations

A computer program for simulating the steady state and transient performance of a four cylinder Stirling engine is presented. The thermodynamic model includes both continuity and energy equations and linear momentum terms (flow resistance). Each working space between the pistons is broken into seven control volumes. Drive dynamics and vehicle load effects are included. The model contains 70 state variables. Also included in the model are piston rod seal leakage effects. The computer program includes a model of a hydrogen supply system, from which hydrogen may be added to the system to accelerate the engine. Flow charts are provided.

Daniele, C. J.↗
NASA NTRS1990

Direct simulation for the instability and breakup of laminar liquid jets

A direct numerical simulation method is described for predicting the deformation of laminar liquid jets. In the present nonlinear direct simulation, the convective term, which has been discarded in past linear analyses by Rayleigh and others, is included in the hydrodynamic equations. It is shown that only by maintaining full complexity of the nonlinear surface tension term accurate drop formation can be predicted. The continuity and momentum equations in the transient form are integrated on an adaptive grid, conforming the jet and surface wave shape. The equations, which are parabolic in time and elliptic in space, are solved by a TVD scheme with characteristic flux splitting. The results of the present work are discussed and compared with available measurements and other analyses. The comparison shows that among the predictions, the current 1-D direct simulation results agree best with the experimental data. Furthermore, the computer time requirements are much (an order of magnitude) smaller than those of previously reported multidimensional analyses.

Chuech, S. G.↗
NASA NTRS1988

Analytical study of the liquid phase transient behavior of a high temperature heat pipe

The transient operation of the liquid phase of a high temperature heat pipe is studied. The study was conducted in support of advanced heat pipe applications that require reliable transport of high temperature drops and significant distances under a broad spectrum of operating conditions. The heat pipe configuration studied consists of a sealed cylindrical enclosure containing a capillary wick structure and sodium working fluid. The wick is an annular flow channel configuration formed between the enclosure interior wall and a concentric cylindrical tube of fine pore screen. The study approach is analytical through the solution of the governing equations. The energy equation is solved over the pipe wall and liquid region using the finite difference Peaceman-Rachford alternating direction implicit numerical method. The continuity and momentum equations are solved over the liquid region by the integral method. The energy equation and liquid dynamics equation are tightly coupled due to the phase change process at the liquid-vapor interface. A kinetic theory model is used to define the phase change process in terms of the temperature jump between the liquid-vapor surface and the bulk vapor. Extensive auxiliary relations, including sodium properties as functions of temperature, are used to close the analytical system. The solution procedure is implemented in a FORTRAN algorithm with some optimization features to take advantage of the IBM System/370 Model 3090 vectorization facility. The code was intended for coupling to a vapor phase algorithm so that the entire heat pipe problem could be solved. As a test of code capabilities, the vapor phase was approximated in a simple manner.

Roche, Gregory Lawrence↗
NASA NTRS1990

Properties of large scale plasma flow during the early stage of the plasmaspheric refilling

The objective is to better characterize the macroscopic properties of the interhemisphere plasma flow by solving a more complete set of hydrodynamic equations than that solved previously. Specifically, the ion continuity, momentum and energy equations were solved for the plasma flow along the closed magnetic field lines. During the initial stage of the supersonic outflow in the equatorial region, the ions cool substantially. Using the hydrodynamic model for the large-scale plasma flow, the dynamics of shocks was examined which form in the geomagnetic flux tubes during the early stages of refilling. These shocks are more like those forming in neutral gases than the electrostatic shocks driven by microinstabilities involving ion-ion interaction. Therefore, the shocks seen in the hydrodynamic model are termed as hydrodynamic shocks. Such shocks are generally unsteady and therefore the usual shock jump conditions given by Rankine-Hugoniot relations are not strictly applicable to them. The density, flow velocity and temperature structures associated with the shocks are examined for both asymmetrical and symmetrical flows. In the asymmetrical flow the outflow from one of two conjugate ionospheres is dominant. On the other hand, in the symmetrical case outflows from the two ionospheric sources are identical. Both cases are treated by a two-stream model. In the late type of flow, the early-time refilling shows a relaxation type of oscillation, which is driven by the large-scale interactions between the two identical streams. After this early stage, the resulting temperature structure shows some interesting features. In the equatorial region the streams are isothermal, but in the off-equatorial regions the streams have quite different temperatures, and also densities and flow velocities. The dense and slow stream is found to be warmer than the low-density fast stream. In the late stage of refilling, the temperature is found to steadily increase from the conjugate ionospheres towards the equator; the equatorial temperature is found to be as high as about 8000 K compared to the ionospheric temperature of 3600 K.

Singh, Nagendra↗
NASA NTRS1988

Naturally occurring continued fractions in the variation of Kepler's equation

A family of functions involving integrals of universal functions is introduced. These functions have some interesting mathematical properties including the fact that they may be expressed as Gaussian continued fractions. A unique method of performing the integration is demonstrated which indicates why these functions may be important in the variation of Kepler's equation.

Shepperd, Stanley W.↗
NASA NTRS1992

Error analysis of finite element results on plates with nonuniform grids

Discrete Dirac delta functions in two dimensions and the governing plate differential equations are used to produce continuous approximations from discrete finite-element data on nonuniform grids. The resulting continuous approximation when compared to the discrete data provides an error analysis. Moreover, the continuous approximation can be differentiated to compute continuous stresses needed for component failure prediction. Stresses can then be calculated accurately at boundaries and at interfaces where finite-element analysis alone may not provide reliable results. The continuous solution along with the differential equations also provides a residual error which is used in the error analysis. This is a unique feature of this approach. Two application examples of the procedure are provided. One example involves the linear response of a clamped square plate under uniform transverse load. Results from finite-element analysis and the continuous solution from the present analysis are compared to a closed-form elasticity solution. The continuous solution is used in a Zienkiewicz-Zhu error estimator to assess the performance of the present analysis. The continuous solution is observed to dramatically correct the finite-element results. In a second example, results for the nonlinear response of a stiffened panel are compared to the finite-element solution.

Sistla, Rajaram↗
NASA NTRS1993

Structure of confined laminar spray diffusion flames: Numerical investigation

The structure of confined laminar spray diffusion flames is investigated numerically by solving the gas-phase conservation equations for mass species, continuity, momentum, and energy and the liquid-phase equations for droplet position, velocity, size, and temperature. A one-step global reaction scheme along with six equilibrium reactions are employed to model the flame chemistry. Monodisperse as well as polydisperse sprays are considered. The numerical results demonstrate that liquid spray flames substantially differ from gaseous flames in their structure, i.e., temperature, concentration, and velocity fields, shape, and dimensions under the same conditions. Spray flames are predicted to be taller and narrower than their counterpart gaseous ones and their shapes are almost cylindrical. This is in agreement with experimental observations. The numerical computations also show that the use of the equilibrium reactions with the one-step reaction scheme decreases the flame temperature compared to the one-step reaction scheme without the equilibrium reactions and more importantly increases the surface area of the flame zone due to a phenomenon termed 'equilibrium broadening.' The spray flames also possess a finite thickness with minimal overlap of the fuel and oxygen species. A case for which a fuel-mixture consisting of 20 to 80 percent gas-liquid by mass is introduced into the combustor is also investigated and compared with predictions using only gaseous or liquid fuel.

Mawid, M. A.↗
NASA NTRS1987

Influence of large deflection and transverse shear on random response of rectangular symmetric composite laminates to acoustic loads

Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived accounting for von Karman strains. The effect of transverse shear is included in the formulation and the rotatory inertia effect is ignored. Using a single-mode Galerkin procedure the nonlinear modal equation is obtained. Direct equivalent linearization is employed. The response of acoustic excitation on moderately thick composite panels is studied. Further, the effects of transverse shear on large deflection vibration of laminates under random excitation are studied. Mean-square deflection and mean-square inplane stresses are obtained for some symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for simply supported plate. The analogous equations for a clamped case are given in the appendix.

Mei, Chuh↗
NASA NTRS1987

Effects of large deflection and transverse shear on response of rectangular symmetric composite laminates subjected to acoustic excitation

Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived accounting for von Karman strains. The effect of transverse shear is included in the formulation. Using a single-mode Galerkin procedure, the nonlinear modal equation is obtained. The direct equivalent linearization method is employed for solution of this equation. The response to acoustic excitation of moderately thick composite panels is studied. Further, the effects of transverse shear on large deflection vibration of laminates under random excitation are studied. Mean-square deflections and mean-square inplane stresses are obtained for symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for simply supported plates. The analogous equations for a clamped case are also obtained.

Mei, Chuh↗
NASA NTRS1998

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗
NASA NTRS1983

Thermodynamic analysis of journal bearings operating under steady state loading in laminar regime

Thermohydrodynamic effects in journal bearings operating under steady load in laminar regime are investigated. An analytical model for the finite and infinitely long journal bearings is formulated. The model includes correction factors for the cavitation effects in the unloaded region of the bearing and the mixing of the recirculating oil and supply oil at the oil inlet. A finite difference computer program is developed to numerically solve the governing equations of the continuity, Reynolds, energy, Laplace heat conduction, and a viscosity-temperature relation simultaneously. The program includes a numerical technique for obtaining an isothermal shaft temperature. The numerical results of temperature distribution and the heat effects on the bearing load carrying capacity agree closely with those of experimental findings. Several different sets of simpler boundary conditions for the energy equation are studied.

Khonsari, M. M.↗
NASA NTRS1989

Hypervelocity atmospheric flight: Real gas flow fields

Flight in the atmosphere is examined from the viewpoint of including real gas phenomena in the flow field about a vehicle flying at hypervelocity. That is to say, the flow field is subject not only to compressible phenomena, but is dominated by energetic phenomena. There are several significant features of such a flow field. Spatially, its composition can vary by both chemical and elemental species. The equations which describe the flow field include equations of state and mass, species, elemental, and electric charge continuity; momentum; and energy equations. These are nonlinear, coupled, partial differential equations that have been reduced to a relatively compact set of equations in a self-consistent manner (which allows mass addition at the surface at a rate comparable to the free-stream mass flux). The equations and their inputs allow for transport of these quantities relative to the mass-average behavior of the flow field. Thus transport of mass by chemical, thermal, pressure, and forced diffusion; transport of momentum by viscosity; and transport of energy by conduction, chemical considerations, viscosity, and radiative transfer are included. The last of these complicate the set of equations by making the energy equations a partial integrodifferential equation. Each phenomenon is considered and represented mathematically by one or more developments. The coefficients which pertain are both thermodynamically and chemically dependent. Solutions of the equations are presented and discussed in considerable detail, with emphasis on severe energetic flow fields. Hypervelocity flight in low-density environments where gaseous reactions proceed at finite rates chemical nonequilibrium is considered, and some illustrations are presented. Finally, flight where the flow field may be out of equilibrium, both chemically and thermodynamically, is presented briefly.

Howe, John T.↗
NASA NTRS1990

Hypervelocity atmospheric flight: Real gas flow fields

Flight in the atmosphere is examined from the viewpoint of including real gas phenomena in the flow field about a vehicle flying at hypervelocity. That is to say, the flow field is subject not only to compressible phenomena, but is dominated by energetic phenomena. There are several significant features of such a flow field. Spatially, its composition can vary by both chemical and elemental species. The equations which describe the flow field include equations of state and mass, species, elemental, and electric charge continuity; momentum; and energy equations. These are nonlinear, coupled, partial differential equations that were reduced to a relatively compact set of equations of a self-consistent manner (which allows mass addition at the surface at a rate comparable to the free-stream mass flux). The equations and their inputs allow for transport of these quantities relative to the mass-averaged behavior of the flow field. Thus transport of mass by chemical, thermal, pressure, and forced diffusion; transport of momentum by viscosity; and transport of energy by conduction, chemical considerations, viscosity, and radiative transfer are included. The last of these complicate the set of equations by making the energy equation a partial integrodifferential equation. Each phenomenon is considered and represented mathematically by one or more developments. The coefficients which pertain are both thermodynamically and chemically dependent. Solutions of the equations are presented and discussed in considerable detail, with emphasis on severe energetic flow fields. For hypervelocity flight in low-density environments where gaseous reactions proceed at finite rates, chemical nonequilibrium is considered and some illustrations are presented. Finally, flight where the flow field may be out of equilibrium, both chemically and thermodynamically, is presented briefly.

Howe, John T.↗