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

Application of program LAURA to thermochemical nonequilibrium flow through a nozzle

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.

Gnoffo, Peter A.↗
NASA NTRS1993

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS1993

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially nonoscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.

Casper, J.↗
NASA NTRS1993

Effects of bleed-hole geometry and plenum pressure on three-dimensional shock-wave/boundary-layer/bleed interactions

A numerical study was performed to investigate 3D shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. This study was focused on how bleed-hole geometry and pressure ratio across bleed holes affect the bleed rate and the physics of the flow in the vicinity of the holes. The aspects of the bleed-hole geometry investigated include angle of bleed hole and the number of bleed holes. The plenum/freestream pressure ratios investigated range from 0.3 to 1.7. This study is based on the ensemble-averaged, 'full compressible' Navier-Stokes (N-S) equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the ensemble-averaged N-S equations were obtained by an implicit finite-volume method using the partially-split, two-factored algorithm of Steger on an overlapping Chimera grid.

Chyu, Wei J.↗
NASA NTRS1993

Supersonic vortex breakdown over a delta wing in transonic flow

The effects of freestream Mach number and angle of attack on the leading-edge vortex breakdown due to the terminating shock on a 65-degree, sharp-edged, cropped delta wing are investigated computationally, using the time-accurate solution of the laminar unsteady compressible full Navier-Stokes equations with the implicit upwind flux-difference splitting, finite-volume scheme. A fine O-H grid consisting of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively, is used for all the flow cases. Keeping the Reynolds number fixed at 3.23 x 10 exp 6, the Mach number is varied from 0.85 to 0.9 and the angle of attack is varied from 20 to 24 deg. The results show that, at 20-deg angle of attack, the increase of the Mach number from 0.85 to 0.9 results in moving the location of the terminating shock downstream. The results also show that, at 0.85 Mach number, the increase of the angle of attack from 20 to 24 deg results in moving the location of the terminating shock upstream. The results are in good agreement with the experimental data.

Kandil, Hamdy A.↗
NASA NTRS1993

Reduction of numerical diffusion in three-dimensional vortical flows using a coupled Eulerian/Lagrangian solution procedure

A new approach based on the coupling of an Eulerian and a Lagrangian solver, aimed at reducing the numerical diffusion errors of standard Eulerian time-marching finite-volume solvers, is presented. The approach is applied to the computation of the secondary flow in two bent pipes and the flow around a 3D wing. Using convective point markers the Lagrangian approach provides a correction of the basic Eulerian solution. The Eulerian flow in turn integrates in time the Lagrangian state-vector. A comparison of coarse and fine grid Eulerian solutions makes it possible to identify numerical diffusion. It is shown that the Eulerian/Lagrangian approach is an effective method for reducing numerical diffusion errors.

Felici, Helene M.↗
NASA NTRS1993

Active control of asymmetric conical flow using spinning and rotatory oscillations

The effectiveness of active control on asymmetric flows around circular cones is investigated computationally using cone spinning and rotatory oscillation around its axis. The investigation uses the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, flux-difference splitting, finite-volume scheme. The present solutions are obtained under the locally-conical-flow assumption in order to understand the flow physics using very fine grids for reasonable flow resolution at low computational cost. For all the computational solutions, a grid of 241 x 81 x 2 points in the wrap-around, normal and axial directions, respectively, is used. The grid is spinning or oscillating rigidly with the cone according to its motion and the kinematical and dynamical boundary conditions are modified accordingly. The computational applications include the effects of uniform spinning rates and periodic rotatory oscillations at different amplitudes and frequencies on the flow asymmetry.

Kandil, Osama A.↗
NASA NTRS1993

Shock-vortex interaction over a 65-degree delta wing in transonic flow

Transonic flow over a 65-deg swept-back, sharp-edged, cropped delta wing is investigated computationally using the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with an implicit, upwind, flux-difference splitting, finite-volume scheme. Coarse and fine O-H grids are used to obtain the solution. The grid consists of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively. The results are presented for an angle of attack of 20 deg Mach number of 0.85 and Reynolds number of 3.23 x 10 exp 6. With the fine grid, the results show that a system of shocks has been captured over the upper wing surface and that the leading-edge vortex core experiences an unsteady supersonic vortex breakdown after passing through a spanwise shock near the wing trailing edge. The computed results at a certain time are in good agreement with the experimental data. Topological aspects of the vortex breakdown flowfield are also presented and discussed.

Kandil, Osama A.↗
NASA NTRS1993

A detailed study of mean-flow solutions for stability analysis of transitional flows

A finite-volume upwind-difference parabolized Navier-Stokes code is utilized to obtain laminar mean-flow solutions at Mach 3.5 on a half-angle cone of 5 deg at an angle-of-attack of 2 deg. A detailed study is conducted on this configuration; the main focus is the velocity profiles in the leeward and windward symmetry planes at various axial locations. Comparisons of the solution profiles are made with both a central-difference code that incorporates scalar and matrix dissipation models and another state-of-the-art upwind-difference finitevolume code. The results obtained emphasize the importance of using matrix dissipation models for schemes that require explicit artificial dissipation. These results also illustrate the accuracy and efficiency of the planeby-plane marching procedure for computing mean-flow solutions for predicting the onset of transition with linear instability.

Ramakrishnan, R.↗
NASA NTRS1993

Simulation of tail buffet using delta wing-vertical tail configuration

Computational simulation of the vertical tail buffet problem is accomplished using a delta wing-vertical tail configuration. Flow conditions are selected such that the wing primary-vortex cores experience vortex breakdown and the resulting flow interacts with the vertical tail. This multidisciplinary problem is solved successively using three sets of equations for the fluid flow, aeroelastic deflections and grid displacements. For the fluid dynamics part, the unsteady, compressible, full Navier-Stokes equations are solved accurately in time using an implicit, upwind, flux-difference splitting, finite-volume scheme. For the aeroelastic part, the aeroelastic equation for bending vibrations is solved accurately in time using the Galerkin method and the four-stage Runge-Kutta scheme. The grid for the fluid dynamics computations is updated every few time steps using a third set of interpolation equations. The computational application includes a delta wing of aspect ratio 1 and a rectangular vertical tail of aspect ratio 2, which is placed at 0.5 root-chord length downstream of the wing trailing edge. The wing angle of attack is 35 deg and the flow Mach number and Reynolds number are 0.4 and 10,000, respectively.

Kandil, Osama A.↗
NASA NTRS1992

Point-implicit relaxation strategies for viscous, hypersonic flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's symmetric total variation diminishing scheme. Newton relaxation of the fully coupled equation set is employed on a cell-to-cell basis. Under-relaxation of the inviscid and over-relaxation of the viscous contributions to the residual are implemented. Computational work is easily partitioned among many processors in an asynchronous, dynamic mode for convergence acceleration. An overview of the physical models employed herein for thermochemical nonequilibrium is included. Several test cases and comparisons with experimental data are presented involving hypersonic flow over blunt bodies which illustrate the qualitative and quantitative capabilities of this approach.

Gnoffo, P. A.↗
NASA NTRS1993

Numerical study of the transient flow in the driven tube and the nozzle section of a shock tunnel

The initial flow in a shock tunnel was examined numerically using computational fluid dynamics (CFD). A finite-volume total variation diminishing (TVD) scheme was used to calculate the transient flow in a shock tunnel. Both viscous and inviscid, chemically nonreacting flows were studied. The study consisted of two parts, the first dealt with the transient flow in the driven-tube/nozzle interface region (inviscid calculations). The effects of varying the geometry in this region was examined. The second part of the study examined the transient flow in the nozzle (viscous calculations). The results were compared to experimental data.

Tokarcik-Polsky, Susan↗
NASA NTRS1993

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Dasgupta, A.↗
NASA NTRS1989

Upwind-biased, point-implicit relaxation strategies for hypersonic flowfield simulations on supercomputers

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗
NASA NTRS1993

Energy considerations in computational aeroacoustics

A finite-volume multistage time-stepping Euler code is used to investigate the use of CFD algorithms for the direct calculation of acoustics. The 2D compressible inviscid flow about an accelerating or decelerating circular cylinder is used as a model problem. The time evolution of the energy transfer from the cylinder to the fluid, as the cylinder is moved from rest to some nonnegligible velocity, is clearly seen. By examining the temporal and spatial characteristics of the numerical solution, a distinction can be made between the propagating acoustic energy, the convecting energy associated with the entropy change in the fluid, and the energy contained in the local aerodynamic field. Systematic variation of the cylinder acceleration shows that the radiated acoustic energy depends strongly upon the rate of acceleration or deceleration. The computational grid has a large effect on the ratio of acoustic energy to nonphysical entropy associated energy, while the role of the explicit artificial viscosity seems to be of second order. The entropy term was nearly negligible in all cases the cylinder was started slowly.

Brentner, Kenneth S.↗
NASA NTRS1993

Large time-step stability of explicit one-dimensional advection schemes

There is a wide-spread belief that most explicit one-dimensional advection schemes need to satisfy the so-called 'CFL condition' - that the Courant number, c = udelta(t)/delta(x), must be less than or equal to one, for stability in the von Neumann sense. This puts severe limitations on the time-step in high-speed, fine-grid calculations and is an impetus for the development of implicit schemes, which often require less restrictive time-step conditions for stability, but are more expensive per time-step. However, it turns out that, at least in one dimension, if explicit schemes are formulated in a consistent flux-based conservative finite-volume form, von Neumann stability analysis does not place any restriction on the allowable Courant number. Any explicit scheme that is stable for c is less than 1, with a complex amplitude ratio, G(c), can be easily extended to arbitrarily large c. The complex amplitude ratio is then given by exp(- (Iota)(Nu)(Theta)) G(delta(c)), where N is the integer part of c, and delta(c) = c - N (less than 1); this is clearly stable. The CFL condition is, in fact, not a stability condition at all, but, rather, a 'range restriction' on the 'pieces' in a piece-wise polynomial interpolation. When a global view is taken of the interpolation, the need for a CFL condition evaporates. A number of well-known explicit advection schemes are considered and thus extended to large delta(t). The analysis also includes a simple interpretation of (large delta(t)) total-variation-diminishing (TVD) constraints.

Leonard, B. P.↗
NASA NTRS1993

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially non-oscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted.

Casper, Jay↗
NASA NTRS1993

Computation of confined coflow jets with three turbulence models

A numerical study of confined jets in a cylindrical duct is carried out to examine the performance of two recently proposed turbulence models: an RNG-based K-epsilon model and a realizable Reynolds stress algebraic equation model. The former is of the same form as the standard K-epsilon model but has different model coefficients. The latter uses an explicit quadratic stress-strain relationship to model the turbulent stresses and is capable of ensuring the positivity of each turbulent normal stress. The flow considered involves recirculation with unfixed separation and reattachment points and severe adverse pressure gradients, thereby providing a valuable test of the predictive capability of the models for complex flows. Calculations are performed with a finite-volume procedure. Numerical credibility of the solutions is ensured by using second-order accurate differencing schemes and sufficiently fine grids. Calculations with the standard K-epsilon model are also made for comparison. Detailed comparisons with experiments show that the realizable Reynolds stress algebraic equation model consistently works better than does the standard K-epsilon model in capturing the essential flow features, while the RNG-based K-epsilon model does not seem to give improvements over the standard K-epsilon model under the flow conditions considered.

Zhu, J.↗