Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Volume Scheme”

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

Computations of vortical supersonic flows over a blunt-nose-cylinder at high angles of attack

Three-dimensional, viscous, and separated flows over a blunt-nose-cylinder at 20, 32, and 44-deg attack angles were computed. The approaching freestream was at a Mach number of 1.6 and a unit Reynolds number of 2 x 10 to the 6th/ft with a total temperature of 585 R. The cylinder used for the computations had a length-to-diameter ratio of 6.67 with a base diameter of 3 in. The flowfield was dominated by large-scale and multiple vortices generated by crossflow separation. The effect of turbulence on the flow structure of one case was modeled algebraically with modifications to correct the length and velocity scales in the regions of separation. The mass averaged Navier-Stokes equations were solved by an approximately factored, upwind-biased, implicit, finite volume scheme. The initialization of the flows was enhanced by a mesh sequencing strategy applied to the diagonalized form of the discretized equations. The convergence to steady-state was accelerated by a multigrid algorithm and using the block inversions for the discretized equations. Calculations were compared with experimental results.

Baysal, O.↗
NASA NTRS1988

A longitudinally-patched grid approach with applications to high speed flows

The combined use of an implicit, upwind finite-volume scheme with an 'equivalent' gamma formulation (for real gas effects) and with a three-dimensional conservative patched grid formulation is discussed. Results are presented over a wide Mach number range on both single and patched grids.

Walters, R. W.↗
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

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 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.↗
NASA NTRS1988

Zonal techniques for flowfield simulation about aircraft

A technique for performing conservative flowfield calculations on zonal meshes is described. The underlying flow solver is an implicit, upwind finite volume scheme which can incorporate either a perfect gas or an equilibrium air equation of state. Two different approaches which yield identical results, in terms of performing a conservative flux calculation on a zonal interface, are described and compared in terms of numerical efficiency. The capability of the method to handle relatively complex geometries is demonstrated by considering the flowfield about a model SR71 aircraft.

Walters, Robert W.↗
NASA NTRS1989

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗
NASA NTRS1989

Thickening oscillation of a delta wing using Navier-Stokes and Navier-displacement equations

The problem of unsteady, supersonic, locally conical, vortical flow around a delta wing undergoing thickening oscillation is solved using the unsteady, thin-layer Navier-Stokes equations and the unsteady linearized Navier-displacement equations. The unsteady, thin-layer Navier-Stokes equations are solved using the implicit approximate-factorization finite-volume scheme to compute the conservative components of the flow vector field. With the conservative components known at any time step, the linearized, Navier-displacement equations are solved using the alternating, direction-implicit scheme to obtain the grid points displacements due to known displacement boundary conditions. A grid-displacements limiter, in the form of a low mesh Reynolds number, is used to limit grid-folding in regions of highly reversed flow.

Chuang, Hsin-Kung A.↗
NASA NTRS1989

Navier-Stokes computations of symmetric and asymmetric vortex shedding around slender bodies

The unsteady, compressible Reynolds-averaged thin-layer Navier-Stokes equations are used to solve for symmetric and asymmetric vortical flows around slender pointed bodies of revolution. The modified Baldwin and Lomax algebraic two-layer turbulent model is used for the eddy viscosity calculation in prescribed turbulent regions. The implicit, upwind flux-difference splitting finite-volume scheme is used to obtain the solutions. Computational results are presented for a low-speed symmetrical vortical flow around a 3.5-caliber tangent-ogive cylinder for fully laminar flow. Computed results are presented for a high-speed asymmetric vortex-shedding flow around a cone at 3.286 relative incidence, using a fine computational grid. The asymmetric vortex-shedding flow has been produced through a small asymmetric transitional perturbation in the cross-flow plane.

Kandil, Osama A.↗
NASA NTRS1990

Two-dimensional Euler and Navier-Stokes Time accurate simulations of fan rotor flows

Two numerical methods are presented which describe the unsteady flow field in the blade-to-blade plane of an axial fan rotor. These methods solve the compressible, time-dependent, Euler and the compressible, turbulent, time-dependent, Navier-Stokes conservation equations for mass, momentum, and energy. The Navier-Stokes equations are written in Favre-averaged form and are closed with an approximate two-equation turbulence model with low Reynolds number and compressibility effects included. The unsteady aerodynamic component is obtained by superposing inflow or outflow unsteadiness to the steady conditions through time-dependent boundary conditions. The integration in space is performed by using a finite volume scheme, and the integration in time is performed by using k-stage Runge-Kutta schemes, k = 2,5. The numerical integration algorithm allows the reduction of the computational cost of an unsteady simulation involving high frequency disturbances in both CPU time and memory requirements. Less than 200 sec of CPU time are required to advance the Euler equations in a computational grid made up of about 2000 grid during 10,000 time steps on a CRAY Y-MP computer, with a required memory of less than 0.3 megawords.

Boretti, A. A.↗
NASA NTRS1990

Analysis and mitigation of numerical dissipation in inviscid and viscid computation of vortex-dominated flows

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 NTRS1990

Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction

High order accurate finite-volume schemes for solving the Euler equations of gasdynamics are developed. Central to the development of these methods are the construction of a k-exact reconstruction operator given cell-averaged quantities and the use of high order flux quadrature formulas. General polygonal control volumes (with curved boundary edges) are considered. The formulations presented make no explicit assumption as to complexity or convexity of control volumes. Numerical examples are presented for Ringleb flow to validate the methodology.

Barth, Timothy J.↗
NASA NTRS1990

Unsteady inviscid and viscous computations for vortex-dominated flows

The unsteady supersonic flow around a rigid sharp-edged delta wing is solved using the unsteady Euler and thin-layer Navier-Stokes equations. The problem is formulated relative to a moving frame of reference so that computation of the grid motion is not required. A three-dimensional, implicit, aproximately factored, central-differencing, finite-volume scheme is used to obtain the time-accurate, unsteady, locally conical flow. The inviscid and viscous results are compared during the periodic response of a wing undergoing forced rolling-oscillation motion. The results show substantial differences in the surface pressure results and other distributed flow characteristics. These differences are attributed to the inviscid and viscous flow predictions in the vortex-dominated regions. The results of the total loads do not show any significant differences.

Kandil, Osama A.↗
NASA NTRS1990

Computation of vortex-dominated flow for a delta wing undergoing pitching oscillation

The conservative, unsteady Euler equations for the flow relative to a moving frame of reference are used to solve for the three-dimensional steady and unsteady flows around a sharp-edged delta wing. 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 crossflow 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 at a large mean angle of attack. The steady results are compared with the experimental data, and the unsteady results are compared with results of a flux-difference splitting scheme.

Kandil, Osama A.↗
NASA NTRS1990

Numerical computation of viscous flow about unconventional airfoil shapes

A new two-dimensional computer code was developed to analyze the viscous flow around unconventional airfoils at various Mach numbers and angles of attack. The Navier-Stokes equations are solved using an implicit, upwind, finite-volume scheme. Both laminar and turbulent flows can be computed. A new nonequilibrium turbulence closure model was developed for computing turbulent flows. This two-layer eddy viscosity model was motivated by the success of the Johnson-King model in separated flow regions. The influence of history effects are described by an ordinary differential equation developed from the turbulent kinetic energy equation. The performance of the present code was evaluated by solving the flow around three airfoils using the Reynolds time-averaged Navier-Stokes equations. Excellent results were obtained for both attached and separated flows about the NACA 0012 airfoil, the RAE 2822 airfoil, and the Integrated Technology A 153W airfoil. Based on the comparison of the numerical solutions with the available experimental data, it is concluded that the present code in conjunction with the new nonequilibrium turbulence model gives excellent results.

Ahmed, S.↗
NASA NTRS1990

Asymmetric supersonic flow around cones with noncircular sections

The unsteady, compressible, single and double thin-layer, Navier-Stokes equations are used to solve for steady and unsteady, asymmetric, supersonic flow around pointed bodies with noncircular sections at high incidence and zero side slip. The equations are solved by using an implicit, upwind, flux-difference splitting finite-volume scheme. Since the flow is locally conical, the solutions are presented on a cross-flow plane at the axial station of unit. The grid is generated by using a modified Joukowski transformation. The computational applications cover noncircular sections with elliptic and diamond shapes. Unsteady asymmetric vortex shedding has been captured at large angles of attack. It is shown that for the same flow conditions and same cross-section fineness ratio, the diamond-section cones develop less flow asymmetry than the elliptic-section cones. Passive control of flow asymmetry was also demonstrated computationally.

Kandil, Osama A.↗
NASA NTRS1991

Computation of compressible quasi-axisymmetric slender vortex flow and breakdown

The unsteady, compressible Navier-Stokes equations are used to compute and analyze compressible quasi-axisymmetric isolated vortices. The Navier-Stokes equations are solved using an implicit, upwind, flux difference splitting finite volume scheme. The developed three dimensional solver was verified by comparing its solution profiles with those of a slender, quasi-axisymmetric vortex solver for a subsonic, quasi-axisymmetric vortex in an unbounded domain. The Navier-Stokes solver is then used to solve for a supersonic, quasi-axisymmetric vortex flow in a configured circular duct. Steady and unsteady vortex-shock interactions and breakdown were captured. The problem was also calculated using the Euler solver of the same code; the results were compared with those of the Navier-Stokes solver. The effect of the initial swirl was investigated.

Kandil, Osama A.↗
NASA NTRS1991

Supersonic quasi-axisymmetric vortex breakdown

An extensive computational study of supersonic quasi-axisymmetric vortex breakdown in a configured circular duct is presented. The unsteady, compressible, full Navier-Stokes (NS) equations are used. The NS equations are solved for the quasi-axisymmetric flows using an implicit, upwind, flux difference splitting, finite volume scheme. The quasi-axisymmetric solutions are time accurate and are obtained by forcing the components of the flowfield vector to be equal on two axial planes, which are in close proximity of each other. The effect of Reynolds number, for laminar flows, on the evolution and persistence of vortex breakdown, is studied. Finally, the effect of swirl ration at the duct inlet is investigated.

Kandil, Osama A.↗