Search NASA⌕ Search

SEARCH · Search NASA

Results for “energy 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 217 records · Page 12

Two-Flux Method for Transient Radiative Transfer in a Semitransparent Layer

The two-flux method was used to obtain transient solutions for a plane layer including internal reflections and scattering. The layer was initially at uniform temperature, and was heated or cooled by external radiation and convection. The two-flux equations were examined as a means for evaluating the radiative flux gradient in the transient energy equation. Comparisons of transient temperature distributions using the two-flux method were made with results where the radiative flux gradient was evaluated from the exact radiative transfer equations. Good agreement was obtained for optical thicknesses from 0.5 to 5 and for refractive indices of 1 and 2. Illustrative results obtained with the two-flux method demonstrate the effect of isotropic scattering coupled with changing the refractive index. For small absorption with large scattering the maximum layer temperature is increased when the refractive index is increased. For larger absorption the effect is opposite, and the maximum temperature decreases with increased refractive index .

Siegel, Robert↗

Validation of a Node-Centered Wall Function Model for the Unstructured Flow Code FUN3D

In this paper, the implementation of two wall function models in the Reynolds averaged Navier-Stokes (RANS) computational uid dynamics (CFD) code FUN3D is described. FUN3D is a node centered method for solving the three-dimensional Navier-Stokes equations on unstructured computational grids. The first wall function model, based on the work of Knopp et al., is used in conjunction with the one-equation turbulence model of Spalart-Allmaras. The second wall function model, also based on the work of Knopp, is used in conjunction with the two-equation k-! turbulence model of Menter. The wall function models compute the wall momentum and energy flux, which are used to weakly enforce the wall velocity and pressure flux boundary conditions in the mean flow momentum and energy equations. These wall conditions are implemented in an implicit form where the contribution of the wall function model to the Jacobian are also included. The boundary conditions of the turbulence transport equations are enforced explicitly (strongly) on all solid boundaries. The use of the wall function models is demonstrated on four test cases: a at plate boundary layer, a subsonic di user, a 2D airfoil, and a 3D semi-span wing. Where possible, different near-wall viscous spacing tactics are examined. Iterative residual convergence was obtained in most cases. Solution results are compared with theoretical and experimental data for several variations of grid spacing. In general, very good comparisons with data were achieved.

Carlson, Jan-Renee↗

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, it is assumed that the total temperature is constant. Therefore, the energy equation can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed. It is shown that the system is strictly hyperbolic and well posed for the initial value problems. Boundary conditions are described such that the linearized system is well posed. The Hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗

On the Navier-Stokes equations with constant total temperature

For various applications in fluid dynamics, one can assume that the total temperature is constant. Therefore, the energy equations can be replaced by an algebraic relation. The resulting set of equations in the inviscid case is analyzed in this paper. It is shown that the system is strictly hyperbolic and well posed for the initial-value problem. Boundary conditions are described such that the linearized system is well posed. The hopscotch method is investigated and numerical results are presented.

Gottlieb, D.↗

Modeling scalar flux and the energy and dissipation equations

Closure models derived from the Two-Scale Direct-Interaction Approximation were compared with data from direct simulations of turbulence. Attention was restricted to anisotropic scalar diffusion models, models for the energy dissipation equation, and models for energy diffusion.

Yoshizawa, A.↗

The role of latent heat release in baroclinic waves - Without beta-effect

An attempt is made to understand the role of latent heat release in the development of baroclinic waves by developing the analytical theory of two-level, quasi-geostrophic baroclinic waves proposed by Phillips (1954) without beta effect, for the case of an atmosphere initially saturated with water vapor and subject to pseudo-adiabatic lifting and dry adiabatic subsistence. The saturated atmosphere condition precludes consideration of conditional convective instability. A consideration of disturbances at each level which depends only on the zonal coordinates and time results in a set of differential equations which is linear in all respects except for the latent heating term in the thermodynamic energy equation. A fundamental nonlinearity is present in this model which is due to heating via the release of latent heat when vertical motion is upward.

Tang, C.-M.↗

Modeling Candle Flame Behavior In Variable Gravity

The burning of a candle, as typical non-propagating diffusion flame, has been used by a number of researchers to study the effects of electric fields on flame, spontaneous flame oscillation and flickering phenomena, and flame extinction. In normal gravity, the heat released from combustion creates buoyant convection that draws oxygen into the flame. The strength of the buoyant flow depends on the gravitational level and it is expected that the flame shape, size and candle burning rate will vary with gravity. Experimentally, there exist studies of candle burning in enhanced gravity (i.e. higher than normal earth gravity, g(sub e)), and in microgravity in drop towers and space-based facilities. There are, however, no reported experimental data on candle burning in partial gravity (g < g(sub e)). In a previous numerical model of the candle flame, buoyant forces were neglected. The treatment of momentum equation was simplified using a potential flow approximation. Although the predicted flame characteristics agreed well with the experimental results, the model cannot be extended to cases with buoyant flows. In addition, because of the use of potential flow, no-slip boundary condition is not satisfied on the wick surface. So there is some uncertainty on the accuracy of the predicted flow field. In the present modeling effort, the full Navier-Stokes momentum equations with body force term is included. This enables us to study the effect of gravity on candle flames (with zero gravity as the limiting case). In addition, we consider radiation effects in more detail by solving the radiation transfer equation. In the previous study, flame radiation is treated as a simple loss term in the energy equation. Emphasis of the present model is on the gas-phase processes. Therefore, the detailed heat and mass transfer phenomena inside the porous wick are not treated. Instead, it is assumed that a thin layer of liquid fuel coated the entire wick surface during the burning process. This is the limiting case that the mass transfer process in the wick is much faster than the evaporation process at the wick surface.

Alsairafi, A.↗

Studies of earth simulation experiments

The low gravity environment of earth orbit offers the potential for performing experiments involving baroclinic Geophysical Fluid Dynamics (GFD) on spherical surfaces. These experiments in turn have the potential for providing deeper understanding of large scale planetary and solar circulations. However, to perform these experiments, one requires an experimental technique whereby a radially directed body force can be generated to simulate a radial gravitational force field. One viable technique is the use of dielectric fluids with temperature dependent dielectric permittivity in a radially directed electric field. Application of the Boussinesq approximation to the equations of motion for this system and restrictions on the size of certain electrodynamic terms in the energy equations yields a set of equations which are analogous to the equations of motions of geophysical systems like the earth's atmosphere on term by term basis. The theoretical design of GFD experiments for performance in earth orbit are described along with results of preliminary tests of a prototype.

Hart, J. E.↗

A numerical study of the unsteady leading edge separation bubble on an oscillating airfoil

A preliminary study of the unsteady viscous flow in the region of an airfoil leading edge was performed, in which the interaction between the viscous and inviscid flow fields is neglected. The solution method uses the finite difference form of the governing equations throughout the separated flow field and incorporates a transition model based on the integral turbulence kinetic energy equation. The validity of the numerical procedure is verified by making comparisons with analytical solutions to several test problems, including unsteady flow over a plate oscillating in its own plane. The method was then applied to the problem of unsteady viscous flow over a NACA 0012 airfoil oscillating sinusoidally in pitch. The flow field characteristics were in qualitative agreement with experimental results. The bubble moved forward on the airfoil and decreased in size as incidence was increased. Viscous flow in the leading edge region was found to be quasi-steady, while bubble height varied inversely with Reynolds number.

Kreskovsky, J. P.↗

Turbulent heat flux measurements in a transitional boundary layer

During an experimental investigation of the transitional boundary layer over a heated flat plate, an unexpected result was encountered for the turbulent heat flux (bar-v't'). This quantity, representing the correlation between the fluctuating normal velocity and the temperature, was measured to be negative near the wall under certain conditions. The result was unexpected as it implied a counter-gradient heat transfer by the turbulent fluctuations. Possible reasons for this anomalous result were further investigated. The possible causes considered for this negative bar-v't' were: (1) plausible measurement error and peculiarity of the flow facility, (2) large probe size effect, (3) 'streaky structure' in the near wall boundary layer, and (4) contributions from other terms usually assumed negligible in the energy equation including the Reynolds heat flux in the streamwise direction (bar-u't'). Even though the energy balance has remained inconclusive, none of the items (1) to (3) appear to be contributing directly to the anomaly.

Sohn, K. H.↗

Heat transfer to the transpired turbulent boundary layer.

This paper contains a summarization of five years work on an investigation on heat transfer to the transpired turbulent boundary layer. Experimental results are presented for friction coefficient and Stanton number over a wide range of blowing and suction for the case of constant free-stream velocity, holding certain blowing parameters constant. The problem of the accelerated turbulent boundary layer with transpiration is considered, experimental data are presented and discussed, and theoretical models for solution of the momentum equation under these conditions are presented. Data on turbulent Prandtl number are presented so that solutions to the energy equation may be obtained. Some examples of boundary layer heat transfer and friction coefficient predictions are presented using one of the models discussed, employing a finite difference solution method.

Kays, W. M.↗

The large Reynolds number - Asymptotic theory of turbulent boundary layers.

A self-consistent, asymptotic expansion of the one-point, mean turbulent equations of motion is obtained. Results such as the velocity defect law and the law of the wall evolve in a relatively rigorous manner, and a systematic ordering of the mean velocity boundary layer equations and their interaction with the main stream flow are obtained. The analysis is extended to the turbulent energy equation and to a treatment of the small scale equilibrium range of Kolmogoroff; in velocity correlation space the two-thirds power law is obtained. Thus, the two well-known 'laws' of turbulent flow are imbedded in an analysis which provides a great deal of other information.

Mellor, G. L.↗

A method for developing K/S boundary conditions

A method for developing the missing general K/S (Kustaanheimo/Stiefel) boundary conditions is presented, with use of the formalism of optimal control theory. As illustrative examples, the method is applied to the transfer between two position and velocity vectors and to the K/S Lambert problem to derive the missing terminal conditions. The necessary equations for a solution are then developed to the K/S Lambert problem with both the fictitious time, s, and the generalized eccentric anomaly, E, as the independent variables. The latter formulation, requiring the solution of only one nonlinear, well-behaved equation in one unknown, E, results in considerable simplification of the problem. This simplification is possible because the energy equation, in the E-formulation, is separable.

Jezewski, D. J.↗

Comparison between the South Pacific Convergence Zone and South Atlantic Convergence Zone in January 1979

Using one of the FGGE data sets, the eddy potential and eddy kinetic energy budgets during the period of January 10-27, 1979 were derived for the regions containing the South Pacific Convergence Zone (SPCZ) and the South Atlantic Convergence Zone. A complete set of energy equations developed by Brennan and Vincent (1980) was partitioned into zonal and eddy components; the eddy forms of energy were derived from a grid point's departure from the zonal average over a 5-deg latitude belt. The results of the analysis indicate that the eddy kinetic energy content (KE) is greater than the eddy available potential energy content (AE) throughout the period. Both AE and KE show large losses beginning when the SPCZ enters a decaying state: AE shows a gradual but slow decrease until after January 22, when it decreases rapidly; KE oscillates, reaching a maximum on January 22 before a rapid decrease.

Vincent, Dayton G.↗

Computer program for solving compressible nonsimilar-boundary-layer equations for laminar, transitional, or turbulent flows of a perfect gas

A computer program is described which solves the compressible laminar, transitional, or turbulent boundary-layer equations for planar or axisymmetric flows. Three-point implicit difference relations are used to reduce the momentum and energy equations to finite-difference form. These equations are solved simultaneously without iteration. Turbulent flow is treated by the inclusion of either a two-layer eddy-viscosity model or a mixing-length formulation. The eddy conductivity is related to the eddy viscosity through a static turbulent Prandtl number which may be an arbitrary function of the distance from the wall boundary. The transitional boundary layer is treated by the inclusion of an intermittency function which modifies the fully turbulent model. The laminar-boundary-layer equations are recovered when the intermittency is zero, and the fully turbulent equations are solved when the intermittency is unity.

Price, J. M.↗

Population inversion calculations using near resonant charge exchange as a pumping mechanism

Near resonance charge exchange between ions of a large ionization potential gas such as helium or neon and vapors of metals such as zinc, cadmium, selenium, or tellurium has produced laser action in the metal ion gas. The possibility of obtaining population inversions in near resonant charge exchange systems (Xe-Ca, Xe-Mg, Xe-Sr, Xe-Ba, Ar-Mg, N-Ca) was investigated. The analysis is an initial value problem that utilizes rate equations for the densities of relevant levels of the laser gas (Ca, Ba, Mg, or Sr) and an electron energy equation. Electron excitation rates are calculated using the Bohr-Thomson approximation for the cross section. Approximations to experimental values of the electron ionization cross section and the ion-atom charge exchange cross section are used. Preliminary results have been obtained for the Ca-Xe system and show that it is possible to obtain gains greater than 10 to the 14th power/m with inversion times up to 8x10 to the minus 7th power second. A possible charge exchange laser system using a MPD arc plasma accelerator is also described.

Chubb, D. L.↗

Wake Statistics of Different-Scale Wind Turbines under Turbulent Boundary Layer Inflow

Subscale wind turbines can be installed in the field for the development of wind technologies, for which the blade aerodynamics can be designed in a way similar to that of a full-scale wind turbine. However, it is not clear whether the wake of a subscale turbine, which is located closer to the ground and faces different incoming turbulence, is also similar to that of a full-scale wind turbine. In this work we investigate the wakes from a full-scale wind turbine of rotor diameter 80 m and a subscale wind turbine of rotor diameter of 27 m using large-eddy simulation with the turbine blades and nacelle modeled using actuator surface models. The blade aerodynamics of the two turbines are the same. In the simulations, the two turbines also face the same turbulent boundary inflows. The computed results show differences between the two turbines for both velocity deficits and turbine-added turbulence kinetic energy. Such differences are further analyzed by examining the mean kinetic energy equation.

17 WIND ENERGY↗

Supercritical time-dependent accretion onto compact objects. I - Neutron stars

The time-dependent flow resulting from spherically symmetric, supercritical accretion onto a nonmagnetic neutron star is computed for accretion rates between 3 and 30. In addition, the fully time-dependent self-consistent set of coupled radiation-hydrodynamic equations governing such flow has been solved. Effects that have not previously been treated are taken into account, including separate energy equations for ions and electrons allowing for nonequilibrium, variable Eddington factors in the radiation transfer to account for anisotropy of the radiation field due to highly spherically extended envelopes, and convection of ions, electrons and photons in the accretion envelope. It is found that supercritical flow, for material falling from a radius of less than 10 to the 10th cm, is characterized by stable accretion with no evidence of stellar wind outflow at any phase of the evolution.

Klein, R. I.↗