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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Full-Scale Direct Numerical Simulation of Two- and Three-Dimensional Instabilities and Rivulet Formulation in Heated Falling Films

A thin film draining on an inclined plate has been studied numerically using finite element method. Three-dimensional governing equations of continuity, momentum and energy with a moving boundary are integrated in an arbitrary Lagrangian Eulerian frame of reference. Kinematic equation is solved to precisely update interface location. Rivulet formation based on instability mechanism has been simulated using full-scale computation. Comparisons with long-wave theory are made to validate the numerical scheme. Detailed analysis of two- and three-dimensional nonlinear wave formation and spontaneous rupture forming rivulets under the influence of combined thermocapillary and surface-wave instabilities is performed.

Krishnamoorthy, S.↗
NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗
NASA NTRS1985

The numerical simulation of steady transonic rotational flow using a dual potential formulation

A finite-difference method is presented that simulates steady transonic rotational flow of an inviscid fluid by representing the velocity field as the sum of scalar and vector potentials. This dual potential velocity decomposition extends the validity of the scalar (full) velocity potential to include vorticity. The inclusion of a vector potential also permits an alternate treatment of lift that does not require a circulation wake cut. This is accomplished by specifying the vector potential as a constant on the airfoil surface in order to satisfy a Kutta condition. The governing equations are solved as iteratively decoupled scalar equations using approximate factorization techniques, and the overall efficiency approaches that of the full potential equation. The governing equations are able to convect entropy and vorticity throughout the flow field and are equivalent to the Euler equations in continuous flow domains, however at shocks the Rankine-Hugoniot entropy jump must be supplied. An entropy correction method is presented and verified with transonic airfoil solutions of the Euler equations.

Chaderjian, N. M.↗
NASA NTRS1952

Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System

The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.

Mager, Arthur↗
NASA NTRS1994

Auroral plasma transport processes in the presence of kV potential structures

We have simulated plasma transport processes in the presence of a quasi-two-dimensional current filament, that generated kV potential structure in the auroral region. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The model uses the 16 moment system of equations and simultaneously solves coupled continuity and momentum equations and equations describing the transport along the magnetic field lines of parallel and perpendicular thermal energy and heat flows for each species. The lower end of the simulation is at an altitude of 800 km, in the collisional topside ionosphere, while the upper end is at 10 R(sub E) in the magnetosphere. The plasma consists of hot electrons and protons of magnetospheric origin and low-energy electrons, protons, and oxygen ions of ionospheric origin. The dynamical interaction of the individual current filaments with ionospheric and magnetospheric plasma generates a potential structure in the horizontal direction and kilovolt field-aligned potential drops along the field lines. The side-by-side display exhibits the evolution of the implied potential structure in the horizontial direction. In the presence of this potential structure and parallel electric field ionospheric plasma density is depleted and velocity is reduced, while density enhancement and increased velocity is observed in magnetospheric plasma. The ionospheric and magnetospheric electron temperatures increase below 2 R(sub E) due to magnetic mirror force on converging geomagnetic field lines. The primary cross-field motion produced by the horizontal E field (E x B drift) is perpendicular to both of the significant spatial directions and is thus ignorable in this geometry. The effects of other cross-field drift processes are discussed. The simulation thus provides insight into the dynamical evolution of two-dimensional potential structures driven by an imposed finite width, field-aligned current profile.

Ganguli, Supriya B.↗
NASA NTRS1992

Towards the development of micromechanics equations for ceramic matrix composites via fiber substructuring

A generic unit cell model which includes a unique fiber substructuring concept is proposed for the development of micromechanics equations for continuous fiber reinforcement ceramic composites. The unit cell consists of three constituents: fiber, matrix, and an interphase. In the present approach, the unit cell is further subdivided into several slices and the equations of micromechanics are derived for each slice. These are subsequently integrated to obtain ply level properties. A stand alone computer code containing the micromechanics model as a module is currently being developed specifically for the analysis of ceramic matrix composites. Towards this development, equivalent ply property results for a SiC/Ti-15-3 composite with 0.5 fiber volume ratio are presented and compared with those obtained from customary micromechanics models to illustrate the concept. Also, comparisons with limited experimental data for the ceramic matrix composite, SiC/RBSN (Reaction Bonded Silicon Nitride) with a 0.3 fiber volume ratio are given to validate the concepts.

Murthy, P. L. N.↗
NASA NTRS1989

Estimates of equatorial upwelling between 140 deg and 110 deg W during 1984

The equation of continuity is used to estimate profiles of vertical velocity between 25 and 120 m from moored current measurements in arrays nested within the triangle with vertices located at 1 deg 30 arcmins S, 140 deg W and along the equator at 140 deg and 110 deg W during December 1983 through March 1984 and May-September 1984. All directions of the 4- to 5-month mean values were upward. The ensemble averaged mean vertical velocity was 0.000022 m/s. Upwelling speeds decreased eastward.

Halpern, David↗
NASA NTRS1967

Theoretical Study of Natural Convection Flows in Closed-End Cylindrical Vessels Final Report

This study is concerned with the analytical solutions of natural convection flows in closed-end cylindrical vessels to obtain exact solutions of temperature and velocity distribution in the laminar flow region under steady state condition. The temperature and velocity distributions, in general power series of displacements, are substituted into the three basic equations of continuity, momentum and energy. The relationships between coefficients of all powers are obtained through initial and boundary conditions and recurrence formulae. For constant wall temperature conditions the general solution of temperature and velocity may be expressed as a function of displacement, geometry ratio and Rayleigh number. Prandtl number does not enter as an independent parameter.

Convective flow↗
NASA NTRS1972

6300 A quantum efficiency of the recombination mechanism in the night-time F layer.

Simultaneous airglow and electron content measurements made at Hawaii are used to infer the number of 6300 and 6364 A quanta produced per electron lost in the nighttime F layer of the ionosphere. The equation of continuity of electrons is then solved numerically to obtain the electron density profile, and the amount of quenching is estimated. This leads to the number of excitations of O(super-1 D) per O2(+) recombination (epsilon). We find, for an exospheric temperature of 1100 K, epsilon is equal to 1.1 plus or minus 0.6, in good agreement with Zipf's laboratory measurement at 300 K.

Brown, W. E.↗
NASA NTRS1973

The role of eddy turbulence in the development of self-consistent models of the lower and upper thermosphere

Numerical solutions of mutually coupled time dependent equations of continuity, momentum and energy balance are presented to illustrate the effect of eddy turbulence on the neutral composition and temperature of the lower and upper atmosphere. The procedure adopted eliminates the necessity of making ad-hoc assumptions about the conditions at the turbopause level and allows the study of the two regions in a self-consistent manner. From the illustrative examples comprising parametric changes in the eddy diffusion coefficient, the specific roles of eddy turbulence in the development of theoretical models of the thermosphere are discussed.

Chandra, S.↗
NASA NTRS1974

Pressure Drop Across Woven Screens Under Uniform and Nonuniform Flow Conditions

Tests were conducted to determine the experimental pressure drop and velocity data for water flowing through woven screens. The types of materials used are dutch twill and square weave fabrics. Pressure drop measures were made at four locations in a rectangular channel. The data are presented as change in pressure compared with the average entry velocity and the numerical relationship is determined by dividing the volumetric flow rate by the screen area open to flow. The equations of continuity and momentum are presented. A computer program listing an extension of a theoretical model and data from that computer program are included.

Ludewig, M.↗
NASA NTRS1978

The Martian ionosphere in light of the Viking observations

A theoretical model has been constructed in which the ion density and the ion and electron temperature distributions are calculated by solving the coupled continuity-momentum equations and the coupled energy equations. The latest experimental results from the Viking 1 and 2 landers are used to vary some of the parameters in the model in order to obtain agreement between the theoretical and experimental results. It is found that solar EUV radiation alone is not able to maintain the observed high ion temperatures. It was also established that the energy coupling between the electron and ion gas is insufficient to account for the measured ion temperatures even in the presence of very large electron temperatures. Direct heat input to the ion gas, probably due to solar wind-ionosphere interactions, can result in ion temperature values in reasonable agreement with the observations. The ion densities calculated with the present model agree well with the Viking observations in the chemically controlled region, but at higher altitudes, dynamic transport processes need to be invoked to achieve consistency among the observed and calculated temperature and density values.

Chen, R. H.↗
NASA NTRS1980

Energetics and thermal structure of the middle atmosphere

A study of a large number of temperature measurements in the middle atmosphere shows a much more complex thermal structure of this region than described in the U.S. Standard Atmosphere, 1976. The mesopause height which is generally assumed to be at 80 km varies between 70 and 100 km, often with two minima in temperature at about 70 and 100 km and a maximum between 80-85 km. By solving the energy balance equation and the equations of continuity, the physical significance of the observed thermal structure is discussed in terms of the energetics of the various regions of the middle atmosphere. It is show that the solar UV radiation plays a major role only in the energy budget of the stratosphere and the lower thermosphere. The energetics of the mesosphere is primarily influenced by the dissipation of eddy energy. The temperature in the region is a good indicator of the eddy diffusivity and can be used in deriving the eddy diffusion coefficient.

Chandra, S.↗
NASA NTRS1981

An algorithm to generate input data from meteorological and space shuttle observations to validate a CH4-CO model

Objective procedures to analyze data from meteorological and space shuttle observations to validate a three dimensional model were investigated. The transport and chemistry of carbon monoxide and methane in the troposphere were studied. Four aspects were examined: (1) detailed evaluation of the variational calculus procedure, with the equation of continuity as a strong constraint, for adjustment of global tropospheric wind fields; (2) reduction of the National Meteorological Center (NMC) data tapes for data input to the OSTA-1/MAPS Experiment; (3) interpolation of the NMC Data for input to the CH4-CO model; and (4) temporal and spatial interpolation procedures of the CO measurements from the OSTA-1/MAPS Experiment to generate usable contours of the data.

Peters, L. K.↗
NASA NTRS1982

Modelling of dimensional stability of fiber reinforced composite materials

Various methods of predicting the expansion and diffusion properties of composite laminates are reviewed. The prediction equations for continuous fiber composites can be applied to SMC composites as the effective fiber aspect ratio in the latter is large enough. The effect of hygrothermal expansion on the dimensional stability of composite laminates was demonstrated through the warping of unsymmetric graphite/epoxy laminates. The warping is very sensitive to the size of the panel, and to the moisture content which is in turn sensitive to the relative humidity in the environment. Thus, any long term creep test must be carried out in a humidity-controlled environment. Environmental effects in SMC composites and bulk polyester were studied under seven different environments. The SMC composites chosen are SMC-R25, SMC-R40, and SMC-R65.

Hahn, H. T.↗