SEARCH · Engineering Papers

Results for “Source Terms”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Extension of a semi-implicit shock-capturing algorithm for 3-D fully coupled, chemically reacting flows in generalized coordinates

A semiimplicit high-resolution shock-capturing method for multidimensional systems of hyperbolic conservation laws with stiff source terms has been developed by Yee and Shinn (1987). The goal of this work is to extend this method to solve the three-dimensional fully coupled Navier-Stokes equations for a hypersonic chemically reacting flow in generalized coordinates. In this formulation, the global continuity equation was replaced by all the species continuity equations. The shock-capturing technique is a second-order-accurate, symmetric total-variation-diminishing method which accounts fully and directly for the coupling among the fluid and all the species. To verify the current approach, it was implemented into an existing computer code which contained the MacCormack method. Test results for a five-species reacting flow are shown to be oscillation-free around the shock, and the time spent per iteration only doubles when compared to the result using classical way of supplying numerical dissipation. The extra computation is more than justified by the elimination of spurious oscillation and nonlinear instability associated with the classical shock-capturing schemes in computing hypersonic reacting flows.

Shinn, J. L.↗

NASA NTRS1985

A hybrid numerical technique for predicting the aerodynamic and acoustic fields of advanced turboprops

A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.

Homicz, G. F.↗

NASA NTRS1988

Investigation of supersonic chemically reacting and radiating channel flow

The 2-D time-dependent Navier-Stokes equations are used to investigate supersonic flows undergoing finite rate chemical reaction and radiation interaction for a hydrogen-air system. The explicit multistage finite volume technique of Jameson is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The multidimensional radiative transfer equations for a nongray model are provided for a general configuration and then reduced for a planar geometry. Both pseudo-gray and nongray models are used to represent the absorption-emission characteristics of the participating species. The supersonic inviscid and viscous, nonreacting flows are solved by employing the finite volume technique of Jameson and the unsplit finite difference scheme of MacCormack. The specified problem considered is of the flow in a channel with a 10 deg compression-expansion ramp. The calculated results are compared with those of an upwind scheme. The problem of chemically reacting and radiating flows are solved for the flow of premixed hydrogen-air through a channel with parallel boundaries, and a channel with a compression corner. Results obtained for specific conditions indicate that the radiative interaction can have a significant influence on the entire flow field.

Mani, Mortaza↗

NASA NTRS1987

Ring vortex/cylinder sound production revisited

The acoustic radiation produced by the passage of a vortex ring near a circular cylinder is calculated utilizing the Coriolis acceleration of the vortex as the source term and the low frequency Green's function approach. Necessary integrals are evaluated exactly via complex integration techniques. The results exhibit good agreement (on a Pascal basis) with experimental data. Of particular interest is the rotation of the directivity pattern during this transient encounter and the extreme sensitivity of the results to the distance at which the vortex passes by the cylinder.

Hardin, Jay C.↗

NASA NTRS1987

Blade-vortex interaction noise predictions using measured blade surface pressures

The generation of helicopter noise by blade-vortex interactions during descent under impulsive conditions is investigated analytically. A noise-prediction technique is developed on the basis of the dipole source term of the Ffowcs-Williams/Hawkings equation and applied to data from simultaneous blade-pressure and acoustic measurements obtained by Cowan et al. (1986) on a 10-ft-diameter 4-blade rotor model in a wind tunnel. Preliminary results show that input-blade-airload azimuth resolution of 1 deg or better and computational azimuth step size of 2 deg or less are required to achieve good agreement between predicted and recorded acoustic time histories. The need for more sophisticated methods to model chordwise input data and for a more extensive experimental data base is indicated.

Ziegenbein, Perry R.↗

NASA NTRS1988

A space-marching method for the computation of viscous internal flows

A space-marching method has been developed to compute 3-D viscous flows in internal geometries. The Navier-Stokes equations have been posed as an initial-value problem by neglecting the effects of streamwise diffusion and treating the streamwise pressure gradient as a known source term. The fully coupled system of equations has been solved by a noniterative algorithm at each streamwise step of the computation. A low Mach number formulation of the equations has been used to compute incompressible flow fields. A computer program has been written to implement all aspects of the space-marching algorithm. The program is modular and is easily adapted to the widely varying geometries of internal flows. The space-marching algorithm has been tested by computing simple flows with known analytical solutions. The method has been used to predict complex 3-D turbulent flows. The algorithm is stable and very economical. A single sweep of the flow field by the space-marching method is approximately equivalent to one time-step of the time-marching method.

Govindan, T. R.↗

NASA NTRS1988

Navier-Stokes solutions of flowfield characteristics produced by ice accretion

The flowfield and resultant heat transfer rates over a series of ice accretion shapes have been obtained through numerical solutions of the Navier-Stokes equations. The influence of roughness is modeled by including blockage, form drag and stagnation heating effects as source terms in the governing equations. Using the flowfield information obtained from the Navier-Stokes equations the droplet impingement efficiencies are computed using a PNS-type solving scheme. Good agreement is achieved between the numerical results and experimental data.

Scott, J. N.↗

NASA NTRS1988

Spatio-temporal adaptation algorithm for two-dimensional reacting flows

A spatio-temporal adaptive algorithm for solving the unsteady Euler equations with chemical source terms is presented. Quadrilateral cells are used in two spatial dimensions which allow for embedded meshes tracking moving flow features with spatially varying time-steps which are multiples of global minimum time-steps. Blast wave interactions corresponding to a perfect gas (frozen) and a Lighthill dissociating gas (nonequilibrium) are considered for circular arc cascade and 90 degree bend duct geometries.

Pervaiz, Mehtab M.↗

NASA NTRS1988

Supersonic combusting flow - A comparison of numerical and experimental results

The flow field of an unconfined coaxial supersonic burner (SSB) was obtained numerically, using a CFD code. The explicit MacCormack algorithm, with point implicit treatment of chemistry source terms, was employed in the CFD code. Combustion of hydrogen and air was simulated by a two-step finite-rate combustion model whereas turbulence was accounted for by a Prandtl mixing length scheme. Boundary conditions and some important features of the numerical scheme were discussed. The qualitative features of the numerical solution were compared with those obtained by regular as well as schlieren photographs. General agreements of the CFD solution and the experimental photographs were observed.

Chitsomboon, T.↗

NASA NTRS1988

CFD prediction of the reacting flow field inside a subscale scramjet combustor

A three-dimensional, Reynolds-averaged Navier-Stokes CFD code has been used to calculate the reacting flowfield inside a hydrogen-fueled, subscale scramjet combustor. Pilot fuel was injected transversely upstream of the combustor and the primary fuel was injected transversely downstream of a backward facing step. A finite rate combustion model with two-step kinetics was used. The CFD code used the explicit MacCormack algorithm with point-implicit treatment of the chemistry source terms. Turbulent mixing of the jets with the airstream was simulated by a simple mixing length scheme, whereas near wall turbulence was accounted for by the Baldwin-Lomax model. Computed results were compared with experimental wall pressure measurements.

Chitsomboon, T.↗

NASA NTRS1988

Transient heat transfer from shrinking LOX-drop

In order to achieve prescribed experimental conditions in wind tunnels, the nitrogen-oxygen mixture is enriched by injection of liquid oxygen (LOX) upstream of methane burners. The objective of the study is to determine the LOX drop evaporation rate, which is dominated by heat transfer from the air stream (mass transfer mechanism can be later coupled with the solution from the present work). Since the concentration of drops in the surrounding gas is quite high, the concept of infinite medium cannot be utilized. Drop evaporation, as part of system's mass balance, is an important source-term in the finite-difference 2/3-dimensional network terms of laboratory-system (Eulerian) coordinates, while individual drop behavior, including its thermal history, is analyzed in Lagrangian coordinates.

Kondic, Nenad N.↗

NASA NTRS1989

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗

NASA NTRS1988

Navier-Stokes calculation of solid-propellant rocket motor internal flowfields

A comprehensive numerical analysis has been carried out to study the detailed physical and chemical processes involved in the combustion of homogeneous propellant in a rocket motor. The formulation is based on the time-dependent full Navier-Stokes equations, with special attention devoted to the chemical reactions in both gas and condensed phases. The turbulence closure is achieved using both the Baldwin-Lomax algebraic model and a modified k-epsilon two-equation scheme with a low Reynolds number and near-wall treatment. The effects of variable thermodynamic and transport properties are also included. The system of governing equations are solved using a multi-stage Runge-Kutta shceme with the source terms treated implicitly. Preliminary results clearly demonstrate the presence of various combustion regimes in the vicinity of propellant surface. The effects of propellant combustion on the motor internal flowfields are investigated in detail.

Hsieh, Kwang-Chung↗

NASA NTRS1989

A set of strongly coupled, upwind algorithms for computing flows in chemical nonequilibrium

Two new algorithms have been developed to predict the flow of viscous, hypersonic, chemically reacting gases over three-dimensional bodies. Both take advantage of the benefits of upwind differencing, Total Variation Diminishing (TVD) techniques and of a finite-volume framework, but obtain their solution in two separate manners. The first algorithym is a time-marching scheme, and is generally used to obtain solutions in the subsonic portions of the flow field. The second algorithm is a much less expensive, space-marching scheme and can be used for the computation of the larger, supersonic portion of the flow field. Both codes compute their interface fluxes with a new temporal Riemann solver and the resulting schemes are made fully implicit including the chemical source terms.

Molvik, Gregory A.↗

NASA NTRS1989

Parabolized Navier-Stokes algorithm for chemically reacting flows

A second-order parabolized Navier-Stokes algorithm based on the MacCormack (1969) explicit scheme is used to study three-dimensional chemically reacting flows with finite-rate chemistry. The method can treat the chemical source term implicitly, and it accounts for the multicomponent diffusion and convection of the chemical species. The method is demonstrated with the nonreacting case of a Mach-3 flow over a double-wedge compression corner and the case of streamwise hydrogen injection at sonic velocity in a Mach-2.44 vitiated air stream.

Kamath, H.↗

NASA NTRS1989

A numerical investigation of the influence of surface roughness on heat transfer in ice accretion

The flowfield and resulting heat transfer rate over a series of ice accretion shapes is obtained by solving the Navier-Stokes equations. The influence of surface roughness on surface heat transfer is examined by including blockage, form drag, and stagnation heating effects as source terms in the governing equations. The results indicate increases of a factor of three in cooling rates due to distributed roughness compared to smotth surfaces. In addition, droplet impingement efficiencies are studied for the same series of ice accretion shapes using a time-dependent solution procedure.

Scott, J. N.↗

NASA NTRS1989

Dynamic loading on impact surfaces of a high subsonic elliptic jet

The feedback resonance in a high subsonic velocity range was investigated for the case where a small aspect ratio (2:1) elliptic jet impinged perpendicularly on a flat plate. This resonance was observed for Mach numbers greater than 0.7 for plate locations between 4 and 9 major radii downstream of the exhaust plane. At the predominant resonant conditions, the loading region and sound source terms were found to be concentrated along the major axis direction on the impact surface.

Wat, J. K.↗

NASA NTRS1988

Effects of furnace temperature profile on the interface shape during Bridgman crystal growth

The effects of asymmetry in furnace temperature profile and pulling velocity on the crystal interface shape are demonstrated while neglecting the latent heat of solidification. It is concluded that the furnace temperature profile may be varied in order to influence the interface shape of the crystal. An exact thermal analysis is performed on the Bridgman technique by including the latent heat of solidification as a source term. The exact temperature field is obtained for the case of a flat melt-crystal interface. The earlier observation regarding the influence of furnace temperature profile on the crystal interface shape is confirmed and a criterion for a flat crystal interface is obtained. Various furnace temperature profiles are selected and their corresponding results are presented.

Taghavi, K.↗