Search NASA⌕ Search

SEARCH · Search NASA

Results for “source term”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗
NASA NTRS1990

Theoretical models of gas dynamics and star formation in interacting ring galaxies

A series of one and two dimensional hydrodynamic simulations of a ring wave in interstellar gas disks was completed. These calculations included nonlinear source terms to model the effects of interstellar interactions and star formation, as well as the spatial-temporal gas flow. Toomre's kinematical model was merged with the Arnold, Shandarin, and Zeldovich 'pancake' theory of caustics in galaxy formation. The resulting theory can describe almost all the structure in restricted three-body simulations of single-pass collisions, even with multi-component potentials. Off-center galactic collisions were studied to understand the dynamics involved. Multi-color optical and near-infrared observations of faint tidal features were performed in about two dozen interacting galaxies selected from the Arp atlas. This sample provided evidence for ongoing star formation in tidal structures, and even enhancements of star formation in some cases. The task of assembling the data for gas-rich, late-type galaxies, was undertaken to see if a more coherent picture of the gas distribution would emerge from the more complete data. Analytic solutions of the equations with subsonic flows to balance gas consumption for expulsion form a galactic fountain were also derived.

Struck-Marcell, Curtis↗
NASA NTRS1990

Near-field noise of a single rotation propfan at an angle of attack

The near-field noise characteristics of a propfan operating at an angle of attack are examined utilizing the unsteady pressure field obtained from a three-dimensional Euler simulation of the propfan flowfield. The near-field noise is calculated employing three different procedures: a direct computation method in which the noise field is extracted directly from the Euler solution, and two acoustic-analogy-based frequency domain methods which utilize the computed unsteady pressure distribution on the propfan blades as the source term. The inflow angles considered are -0.4, 1.6, and 4.6 degrees. The results of the direct computation method and one of the frequency domain methods show qualitative agreement with measurements. They show that an increase in the inflow angle is accompanied by an increase in the sound pressure level at the outboard wing boom locations and a decrease in the sound pressure level at the (inboard) fuselage locations. The trends in the computed azimuthal directivities of the noise field also conform to the measured and expected results.

Nallasamy, M.↗
NASA NTRS1991

Radiative interactions in a hydrogen-fueled supersonic combustor

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with finite-rate chemistry and radiation, for hydrogen-air systems. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed flow in a channel with a 10-deg compression ramp. Three different chemistry models are used, accounting for increasing number of reactions and participating species. Two chemistry models assume nitrogen as inert, while the third model accounts for nitrogen reactions and NO(x) formation. The tangent slab approximation is used in the radiative flux formulation. A pseudogray model is used to represent the absorption-emission characteristics of the participating species. Results otained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species, and that this can have a significant influence on the flowfield.

Chandrasekhar, R.↗
NASA NTRS1991

Computational aspects of chemically reacting flows

An implicit, time-dependent approach is taken to solve the inviscid conservative equations in a vertex-based finite-volume formulation. The convective fluxes are estimated by a modified Riemann solver, and higher-order accuracy is obtained by means of the MUSCL. Implicit techniques are adopted to deal with the equation stiffness as a result of the clustered grid and chemical source term. The real-gas Roe matrix and the ADI and LU schemes have been tested on the sample problems: hypersonic flow over blunt configurations, such as the forebody of a double ellipse, and the near wake of a conic aerobrake. Temperature distribution and standoff distance are validated against the shock-fitting, central-differencing results. Equilibrium temperature may be recovered by artificially increasing the chemical rates at the wall. The present method satisfactorily resolves the flow structure, including the strong expansion behind the flat base. The relative merits and performance of the real-gas TVD and implicit schemes are discussed briefly.

Li, C. P.↗
NASA NTRS1991

Application of a Reynolds Stress turbulence model to a supersonic hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows undergoing hydrogen-air chemical reactions. An assumed Beta Probability Density Function is applied to account for the chemical source terms in the conservation equations. An algebraic Reynolds Flux model is used for the fluctuating density-velocity as well as the species mass fraction-velocity correlations. The variances of temperature and species fluctuations are also modelled using an algebraic flux technique. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion processes. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant.

Chandrasekhar, R.↗
NASA NTRS1991

Nonlinear analysis of longitudinal-mode liquid propellant rocket combustion instability

A computational technique for the nonlinear analysis of longitudinal-mode liquid propellant rocket combustion instability is examined based on the unsteady, quasi-one-dimensional Euler equations with appropriate source terms introduced to account for interphase transport coupling with the spray. The method is first assessed for unsteady, nonreacting, isentropic duct flow with specified admittances at the outflow and inflow boundaries. For small amplitude disturbances, numerical results based on a sufficient number of grid points compare favorably with predictions from a small-disturbance linearized analysis.

Jeng, San-Mou↗
NASA NTRS1991

Simulation of cryogenic liquid flows with vapor bubbles

Liquid flows in rocket engine components (such as bearings, seals, and pumps) often involve the formation of vapor bubbles due to local superheating of the fluid (either boiling or cavitation). Under the present effort, an analysis has been developed for liquid flows with vapor bubbles, based on a combined Eulerian-Lagrangian technique, in which the continuous (liquid) phase is treated by solving a system of Eulerian conservation equations, while the discrete (vapor bubble) phase is dealt with by integrating Lagrangian equations of motion in computational coordinates. Vapor bubbles of changing size can be accommodated easily by this analysis, and models for the simulation of bubble formation, growth, and motion have been included. The effect of bubble motion and other bubble processes on the continuous (liquid) phase has been accounted for by appropriate bubble mass, momentum, and energy interchange source terms in the Eulerian conservation equations. To demonstrate the viability of the resulting procedure, the cavitating flow of liquid oxygen through a simplified model of a labyrinth seal has been successfully calculated.

De Jong, Frederik J.↗
NASA NTRS1991

Mass-loading and parallel magnetized shocks

Recent observations at comets Giacobini-Zinner and Halley suggest that simple nonreacting gas dynamics or MHD is an inappropriate description for the bow shock. The thickness of the observed (sub)shock implies that mass-loading is an important dynamical process within the shock itself, thereby requiring that the Rankine-Hugoniot conditions possess source terms. This leads to shocks with properties similar to those of combustion shocks. The paper considers parallel magnetized shocks subjected to mass-loading, describes some properties which distinguish them from classical MHD parallel shocks, and establishes the existence of a new kind of MHD compound shock. These results will be of importance both to observations and numerical simulations of the comet-solar wind interaction.

Zank, G. P.↗
NASA NTRS1992

Spurious Solutions Of Nonlinear Differential Equations

Report utilizes nonlinear-dynamics approach to investigate possible sources of errors and slow convergence and non-convergence of steady-state numerical solutions when using time-dependent approach for problems containing nonlinear source terms. Emphasizes implications for development of algorithms in CFD and computational sciences in general. Main fundamental conclusion of study is that qualitative features of nonlinear differential equations cannot be adequately represented by finite-difference method and vice versa.

Yee, H. C.↗
NASA NTRS1991

Direct computation of turbulence and noise

Jet exhaust turbulence noise is computed using a time dependent solution of the three dimensional Navier-Stokes equations to supply the source terms for an acoustic computation based on the Phillips convected wave equation. An extrapolation procedure is then used to determine the far field noise spectrum in terms of the near field sound. This will lay the groundwork for studies of more complex flows typical of noise suppression nozzles.

Berman, C.↗
NASA NTRS1991

Turbulence modeling in supersonic combusting flows

To support the National Aerospace Plane project, the RPLUS3D CFD code has been developed at NASA Lewis. The code has the ability to solve three-dimensional flowfields with finite rate combustion of hydrogen and air. The combustion processes of the hydrogen-air system are simulated by an 18-reaction path, 8-species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though it uses a fully implicit numerical algorithm. A k-epsilon (two equation) turbulence model is incorporated into the RPLUS3D code.

Chitsomboon, Tawit↗
NASA NTRS1992

A compressible Navier-Stokes solver with two-equation and Reynolds stress turbulence closure models

This report outlines the development of a general purpose aerodynamic solver for compressible turbulent flows. Turbulent closure is achieved using either two equation or Reynolds stress transportation equations. The applicable equation set consists of Favre-averaged conservation equations for the mass, momentum and total energy, and transport equations for the turbulent stresses and turbulent dissipation rate. In order to develop a scheme with good shock capturing capabilities, good accuracy and general geometric capabilities, a multi-block cell centered finite volume approach is used. Viscous fluxes are discretized using a finite volume representation of a central difference operator and the source terms are treated as an integral over the control volume. The methodology is validated by testing the algorithm on both two and three dimensional flows. Both the two equation and Reynolds stress models are used on a two dimensional 10 degree compression ramp at Mach 3, and the two equation model is used on the three dimensional flow over a cone at angle of attack at Mach 3.5. With the development of this algorithm, it is now possible to compute complex, compressible high speed flow fields using both two equation and Reynolds stress turbulent closure models, with the capability of eventually evaluating their predictive performance.

Navier-Stoke solver↗
NASA NTRS1992

Upwind schemes and bifurcating solutions in real gas computations

The area of high speed flow is seeing a renewed interest due to advanced propulsion concepts such as the National Aerospace Plane (NASP), Space Shuttle, and future civil transport concepts. Upwind schemes to solve such flows have become increasingly popular in the last decade due to their excellent shock capturing properties. In the first part of this paper the authors present the extension of the Osher scheme to equilibrium and non-equilibrium gases. For simplicity, the source terms are treated explicitly. Computations based on the above scheme are presented to demonstrate the feasibility, accuracy and efficiency of the proposed scheme. One of the test problems is a Chapman-Jouguet detonation problem for which numerical solutions have been known to bifurcate into spurious weak detonation solutions on coarse grids. Results indicate that the numerical solution obtained depends both on the upwinding scheme used and the limiter employed to obtain second order accuracy. For example, the Osher scheme gives the correct CJ solution when the super-bee limiter is used, but gives the spurious solution when the Van Leer limiter is used. With the Roe scheme the spurious solution is obtained for all limiters.

Suresh, Ambady↗
NASA NTRS1992

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulent models. The probability density function (PDF) method offers an attractive alternative: in a PDF model, the chemical source terms are closed and do not require additional models. Because the number of computational operations grows only linearly in the Monte Carlo scheme, it is chosen over finite differencing schemes. A grid dependent Monte Carlo scheme following J.Y. Chen and W. Kollmann has been studied in the present work. It was found that in order to conserve the mass fractions absolutely, one needs to add further restrictions to the scheme, namely alpha(sub j) + gamma(sub j) = alpha(sub j - 1) + gamma(sub j + 1). A new algorithm was devised that satisfied this restriction in the case of pure diffusion or uniform flow problems. Using examples, it is shown that absolute conservation can be achieved. Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗
NASA NTRS1992

Implementation of a kappa-epsilon turbulence model to RPLUS3D code

The RPLUS3D code has been developed at the NASA Lewis Research Center to support the National Aerospace Plane (NASP) project. The code has the ability to solve three dimensional flowfields with finite rate combustion of hydrogen and air. The combustion process of the hydrogen-air system are simulated by an 18 reaction path, 8 species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though a fully implicit numerical algorithm is used. A k-epsilon turbulence model has recently been incorporated into the code. Initial validations have been conducted for a flow over a flat plate. Results of the validation studies are shown. Some difficulties in implementing the k-epsilon equations to the code are also discussed.

Chitsomboon, Tawit↗
NASA NTRS1991

Transient cooling of a square region of radiating medium

A time-accurate numerical solution was carried out for transient radiative cooling of a gray emitting and absorbing medium in a square two-dimensional region. The integro-differential energy equation for transient temperature distributions was solved in two stages. At each time increment, the local radiative source term was obtained by numerical integration of the temperature field using two-dimensional Gaussian integration over rectangular subregions. Then the differential portion of the equation was integrated forward in time by use of the local first and second time derivatives. The results were compared with available limiting case, and excellent agreement was obtained. Transient results are given for a wide range of optical thicknesses of the region. Optimum transient cooling is obtained when the optical side length is about 4.

Siegel, Robert↗
NASA NTRS1992

Finite difference solution for transient cooling of a radiating-conducting semitransparent layer

Transient solutions were obtained for cooling a semitransparent material by radiation and conduction. The layer is in a vacuum environment so the only means for heat dissipation is by radiation from within the medium leaving through the boundaries. Heat conduction serves only to partially equalize temperatures across the layer. As the optical thickness is increased, steep temperature gradients exist near the boundaries when conduction is relatively small. A solution procedure is required that will provide accurate temperature distributions adjacent to the boundaries, or radiative heat losses will be in error. The approach utilized numerical Gaussian integration to obtain the local radiative source term, and a finite difference procedure with variable space and time increments to solve the transient energy equation.

Siegel, Robert↗