Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Volume”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Prediction of unsteady loads on maneuvering delta wings using time-accurate Euler schemes

Three-dimensional steady and unsteady vortex-dominated flows around sharp-edged delta wings are considered in this paper. The problem is formulated by using the unsteady conservative Euler equations for the flow relative motion with respect to a moving frame of reference. An implicit approximately-factored finite volume scheme is used to solve the resulting equations on a three-dimensional computational grid which is generated by using a modified Joukowski transformation in cross-flow planes at the grid chord stations. The scheme is applied to a delta wing undergoing pitching oscillation around a large angle of attack. The initial conditions correspond to a steady flow around a delta wing of aspect ratio of one, freestream Mach number of 0.3 and mean angle of attack of 20.5. The steady flow results are compared with those of an explicit computational scheme and the experimental data, and they are in good agreement.

Kandil, Osama A.↗
NASA NTRS1987

Flow solution on a dual-block grid around an airplane

The compressible flow around a complex fighter-aircraft configuration (fuselage, cranked delta wing, canard, and inlet) is simulated numerically using a novel grid scheme and a finite-volume Euler solver. The patched dual-block grid is generated by an algebraic procedure based on transfinite interpolation, and the explicit Runge-Kutta time-stepping Euler solver is implemented with a high degree of vectorization on a Cyber 205 processor. Results are presented in extensive graphs and diagrams and characterized in detail. The concentration of grid points near the wing apex in the present scheme is shown to facilitate capture of the vortex generated by the leading edge at high angles of attack and modeling of its interaction with the canard wake.

Eriksson, Lars-Erik↗
NASA NTRS1987

Adaptive grid generation

The fundamental principles of adaptive grid generation for the numerical analysis of physical phenomena described by systems of partial differential equations are examined in an analytical review. Topics addressed include weight functions, equidistribution in one dimension, the specification of coefficients in the linear weight, the attraction to a given grid on a curve, evolutionary forces, and metric notation. Consideration is given to curve-by-curve methods, finite-volume methods, variational methods, and temporal aspects.

Eiseman, Peter R.↗
NASA NTRS1988

Calculations of three-dimensional flows using the isenthalpic Euler equations with implicit flux-vector splitting

A numerical method for solving the isenthalpic form of the Euler equations is developed. The method is based on the concept of flux vector splitting in its implicit form applied to a cell centered finite volume scheme. Approximate factorization is implemented in solving the implicit part of the governing equations. Time marching to a steady state solution requires short computational times due to the relative efficiency of the basic method. Computational times are further reduced by the implementation of multigrid. Results for several basic cases are shown.

Cannizzaro, Frank E.↗
NASA NTRS1988

Computational simulation of turbulent flow in Space Shuttle Main Engine turnaround ducts

Two axisymmetric 180-deg turnaround ducts are used in the Space Shuttle Main Engine (SSME) to connnect the preburners with the main thrust chamber. The prediction of the flow field and heat transfer within the turnaround ducts is important to minimize the pressure drops, size of recirculation and stagnation zones, local overheating, etc. This paper presents the computational approach and selected results for the turbulent flow in the turnaround ducts of the fuel and oxidizer sides of the SSME. The time-averaged, Navier-Stokes equations for the viscous, compressible, turbulent flow are solved in body-fitted-coordinates by using a finite-volume approach. Two turbulence models, viz: the k-epsilon model and a multiple scale turbulence model, are used to examine the sensitivity of calculated flows. Both models produce almost identical solutions for the fuelside turnaround duct (which has no recirculation region). However, for the oxidizer-side duct which has a large recirculation region, the two models show quite different results.

Przekwas, Andrzej J.↗
NASA NTRS1988

The multi-zone calculation of turbomachinery flows. II - The multi-zone calculation of the turbulent, two-specie flow through the SSME HPFTP first and second stage cavities

A multi-zone Navier-Stokes methodology to calculate the two-specie flow through the first and second stage cooling cavities of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump (HPFTP) is developed. A simplified two-component fluid formulation is used to model the interaction of coolant and hot gas. Johnston's secant approximation is used to define an appropriate near wall velocity for use in a three-dimensional law of the wall. The basic Navier-Stokes algorithm used is a finite-volume, predictor-corrector algorithm which uses a pressure correction technique. A multi-zone method is used to partition each cavity into easily handled subdomains. The results show that coolant flow is pumped up the turbine wheel for both cavities, creating a region of large temperature gradients on the turbine shank.

Williams, M.↗
NASA NTRS1988

Unsteady delta-wing flow computation using an implicit factored Euler scheme

The conservative unsteady Euler equations for the flow relative motion in the moving frame of reference are used to solve for the steady and unsteady flows around sharp-edged delta wings. The resulting equations are solved by using an implicit approximately-factored finite-volume scheme. Implicit second-order and explicit second- and fourth-order dissipations are added to the scheme. The boundary conditions are explicitly satisfied. The grid is generated by locally using a modified Joukowski transformation in cross-flow planes at the grid chord stations. The computational applications cover a steady flow around a delta wing whose results serve as the initial conditions for the unsteady flow around a pitching delta wing about a large angle of attack. The steady results are compared with the experimental data and the periodic solution is achieved within the third cycle of oscillation.

Kandil, Osama A.↗
NASA NTRS1988

Control of supersonic intersection flowfields through filleting and sweep

The problem of supersonic flow control using fillets and sweep for a wing/body junction has been investigated numerically using a three-dimensional Navier-Stokes code which employs the MacCormack time-split finite-volume technique. An elliptic grid generation technique with direct control over spacing has been incorporated for constructing the grid at a filleted wing/body junction. The computed results for pressure distribution, particle paths, and limiting streamlines on the flat plate and fin surface for a swept fin show a decrease in the peak pressure on the fin leading edge and in the extent of the separated flow region. Moreover, the results for the filleted juncture clearly show that the flow streamline patterns lose much of their vortical character with proper filleting. It is demonstrated that fillets with a radius of 3.5 times the fin leading edge diameter are required to weaken the usual necklace vortex interaction.

Lakshmanan, B.↗
NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗
NASA NTRS1988

Computation of steady and unsteady vortex-dominated flows with shock waves

The unsteady Euler equations have been derived in the conservation form for the flow relative motion with respect to a rotating frame of reference. The resulting equations are solved by using a central-difference finite-volume scheme with four-state Runge-Kutta time stepping. For steady flow applications local time stepping is used, and for unsteady applications the minimum global time stepping is used. A three-dimensional fully vectorized computer program has been developed and applied to steady and unsteady maneuvering delta wings. The capability of the three-dimensional program has been demonstrated for a rigid sharp-edged delta wing undergoing uniform rolling in a conical flow and rolling oscillations in a locally conical flow.

Kandil, Osama A.↗
NASA NTRS1988

Numerical simulation of underexpanded plumes using upwind algorithms

The solution of the flow characteristics of highly underexpanded plumes is presently sought by means of a high-resolution total variation-diminishing upwind scheme featuring finite-volume metrics, Roe averaging with entropy correction, and monotonic upstream schemes for conservation laws. The existing approximate factorization schemes are used both as preconditioners and in order to apply the conjugate residual algorithm. The plumes treated are axisymmetric and laminar, and for an ideal gas that is allowed to be thermally incomplete. Simulations of the SSME and SRB plumes have been successfully conducted.

Obayashi, Shigeru↗
NASA NTRS1988

Interactive grid generation for turbomachinery flow field simulations

The control point form of algebraic grid generation presented provides the means that are needed to generate well structured grids for turbomachinery flow simulations. It uses a sparse collection of control points distributed over the flow domain. The shape and position of coordinate curves can be adjusted from these control points while the grid conforms precisely to all boundaries. An interactive program called TURBO, which uses the control point form, is being developed. Basic features of the code are discussed and sample grids are presented. A finite volume LU implicit scheme is used to simulate flow in a turbine cascade on the grid generated by the program.

Choo, Yung K.↗
NASA NTRS1988

Higher order turbulence closure models

Theoretical models are developed and numerical studies conducted on various types of flows including both elliptic and parabolic. The purpose of this study is to find better higher order closure models for the computations of complex flows. This report summarizes three new achievements: (1) completion of the Reynolds-stress closure by developing a new pressure-strain correlation; (2) development of a parabolic code to compute jets and wakes; and, (3) application to a flow through a 180 deg turnaround duct by adopting a boundary fitted coordinate system. In the above mentioned models near-wall models are developed for pressure-strain correlation and third-moment, and incorporated into the transport equations. This addition improved the results considerably and is recommended for future computations. A new parabolic code to solve shear flows without coordinate tranformations is developed and incorporated in this study. This code uses the structure of the finite volume method to solve the governing equations implicitly. The code was validated with the experimental results available in the literature.

Amano, Ryoichi S.↗
NASA NTRS1988

Three-dimensional multigrid algorithms for the flux-split Euler equations

The Full Approximation Scheme (FAS) multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each of the splitting algorithms uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computation required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined.

Anderson, W. Kyle↗
NASA NTRS1989

Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a tubine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of independence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Choo, Yung K.↗
NASA NTRS1989

Unsteady transonic flow using Euler equations

An implicit, two factor, split flux, finite volume Euler equations solution algorithms is applied to the time accurate solution of transonic flow about an NACA 0012 airfoil and a rectangular planform supercritical wing undergoing pitch oscillations. Accuracy for Courant numbers greater than one is analyzed. Freezing the flux Jacobians can result in significant savings for steady state solutions; the accuracy of freezing flux Jacobians for unsteady results is investigated. The Euler algorithm results are compared with experimental results for an NACA 0012 and a rectangular planform supercritical wing.

Belk, Dave M.↗
NASA NTRS1989

Unsteady Navier-Stokes computations over airfoils using both fixed and dynamic meshes

A finite volume implicit approximate factorization method which solves the thin layer Navier-Stokes equations was used to predict unsteady turbulent flow airfoil behavior. At a constant angle of attack of 16 deg, the NACA 0012 airfoil exhibits an unsteady periodic flow field with the lift coefficient oscillating between 0.89 and 1.60. The Strouhal number is 0.028. Results are similar at 18 deg, with a Strouhal number of 0.033. A leading edge vortex is shed periodically near maximum lift. Dynamic mesh solutions for unstalled airfoil flows show general agreement with experimental pressure coefficients. However, moment coefficients and the maximum lift value are underpredicted. The deep stall case shows some agreement with experiment for increasing angle of attack, but is only qualitatively comparable past stall and for decreasing angle of attack.

Rumsey, Christopher L.↗
NASA NTRS1989

Transonic flutter calculations using the Euler equations

In transonic flutter problems where shock motion plays an important part, it is believed that accurate predictions of the flutter boundaries will require the use of codes based on the Euler equations. Only Euler codes can obtain the correct shock location and shock strength, and the crucially important shock excursion amplitude and phase lag. The present study is based on the finite volume scheme developed by Jameson and Venkatakrishnan for the 2-D unsteady Euler equations. The equations are solved in integral form on a moving grid. The variable are pressure, density, Cartesian velocity components, and total energy.

Bendiksen, Oddvar O.↗