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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Full potential solution of a transonic quasi-3-D flow through a cascade using artificial compressibility

A reliable method is presented for calculating the flowfield about a cascade of arbitrary 2-D airfoils. The method approximates the three-dimensional flow in a turbomachinery blade row by correcting for streamtube convergence and radius change in the throughflow direction. The method is a fully conservative solution of the full potential equation incorporating the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to ensure stability and the capture of shock waves. Comparison of results for several supercritical blades shows good agreement with their hodograph solutions. Other calculations for these profiles as well as standard NACA blade sections indicate that this is a useful scheme for analyzing both the design and off-design performance of turbomachinery blading.

Farrell, C.↗
NASA NTRS1981

Numerical method for gas dynamics combining characteristic and conservation concepts

An efficient implicit numerical method that solves the compressible Navier-Stokes equations in arbitrary curvilinear coordinates by the finite-volume technique is presented. An intrinsically dissipative difference scheme and a fully implicit treatment of boundary conditions, based on characteristic and conservation concepts, are used to improve stability and accuracy. Efficiency is achieved by using a diagonal form of the implicit algorithm and spatially varying time-steps. Comparisons of various schemes and methods are presented for one- and two-dimensional flows, including transonic separated flow past a thick circular-arc airfoil in a channel. The new method is equal to or better than a version of MacCormack's hybrid method in accuracy and it converges to a steady state up to an order of magnitude faster.

Coakley, T. J.↗
NASA NTRS1981

Calculations of transonic potential flow over cascades

Transonic flow through a cascade was studied by using the full potential equation and the finite volume method of Jameson and Caughey. The C-type computational grid is generated by an electrostatic analogy and simple shearing transformation. The solution algorithm includes an option of using either an artificial density or an artificial viscosity formulation of the dissipative term. Using the developed code, flows through a cascade of NACA 0012 airfoils and flows through a cascade of shockless blades were computed. It is found that the designed flow through the shockless blade is accurately predicted, the artificial density formulation shows more tolerance to the mesh irregularity, and the C-type mesh does not extend very far upstream for a small pitch-cord ratio.

Jou, W. H.↗
NASA NTRS1982

Computer program for calculating full potential transonic, quasi-three-dimensional flow through a rotating turbomachinery blade row

A fast, reliable computer code is described for calculating the flow field about a cascade of arbitrary two dimensional airfoils. The method approximates the three dimensional flow in a turbomachinery blade row by correcting for stream tube convergence and radius change in the throughflow direction. A fully conservative solution of the full potential equation is combined with the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. The instructions required to set up and use the code are included. The name of the code is QSONIC. A numerical example is also given to illustrate the output of the program.

Farrell, C. A.↗
NASA NTRS1982

Singularity embedding method in potential flow calculations

The so-called H-type mesh is used in a finite-element (or finite-volume) calculation of the potential flow past an airfoil. Due to coordinate singularity at the leading edge, a special singular trial function is used for the elements neighboring the leading edge. The results using the special singular elements are compared to those using the regular elements. It is found that the unreasonable pressure distribution obtained by the latter is removed by the embedding of the singular element. Suggestions to extend the present method to transonic cases are given.

Jou, W. H.↗
NASA NTRS1982

Transonic potential flow and coordinate generation for bodies in a wind tunnel

An accurate method has been developed for computation of transonic potential flow about a 2-D lifting airfoil or an axisymmetric body in a wind tunnel. The computational mesh is nearly orthogonal everywhere and is generated by a sequence of Schwarz-Christoffel transformations and shearings to obtain an '0' grid near the body. A conservative finite-volume scheme for the full-potential equation and exact boundary conditions is used together with 'retarded' density to solve efficiently transonic flow with embedded shocks and large regions of supersonic flow. Supersonic free-stream flows are also solved with captured bow shocks and embedded subsonic regions.

Doria, M. L.↗
NASA NTRS1982

A full potential inverse method based on a density linearization scheme for wing design

A mixed analysis inverse procedure based on the full potential equation in conservation form was developed to recontour a given base wing to produce density linearization scheme in applying the pressure boundary condition in terms of the velocity potential. The FL030 finite volume analysis code was modified to include the inverse option. The new surface shape information, associated with the modified pressure boundary condition, is calculated at a constant span station based on a mass flux integration. The inverse method is shown to recover the original shape when the analysis pressure is not altered. Inverse calculations for weakening of a strong shock system and for a laminar flow control (LFC) pressure distribution are presented. Two methods for a trailing edge closure model are proposed for further study.

Shankar, V.↗
NASA NTRS1983

TAWFIVE: A user's guide

The Transonic Analysis of a Wing and Fuselage with Interacted Viscous Effects (TAWFIVE) was developed. A finite volume full potential method is used to model the outer inviscid flow field. First-order viscous effects are modeled by a three dimensional integral boundary layer method. Both turbulent and laminar boundary layers are treated. Wake thickness and curvature effects are modeled using a two dimensional strip method. A very brief discussion of the engineering aspects of the program is given. The input and use of the program are covered in great detail.

Melson, N. D.↗
NASA NTRS1983

Multi-grid calculation of three-dimensional transonic potential flows

A multi-grid algorithm has been developed to speed the iterative convergence of calculations for the transonic potential flow past swept wings and wing-fuselage combinations. The method is based upon a fully-conservative, finite-volume approximation to the steady potential equation which is second-order accurate everywhere in the flow field except near shock waves. The multi-grid scheme is incorporated within the framework of an alternating successive-line-overrelaxation (SLOR) solver of the difference equations. Computed results confirm the second-order accuracy of the scheme, and demonstrate the effectiveness of the multi-grid procedure.

Caughey, D. A.↗
NASA NTRS1983

Transonic flow simulation of prop-fan configurations

This paper presents the results of a study that used a three-dimensional transonic analysis method to analyze flow around a wing-mounted prop-fan configuration. The transonic analysis method employs a cell-oriented finite volume approach to solve the full potential equation in conservative form over a numerically generated surface-fitted grid. The propeller slipstream effects are simulated using a linearized transpiration boundary condition. The results of analyses of flow over clean wing/body, wing/body/nacelle, and wing/body nacelle with slipstream are compared to test data.

Samant, S. S.↗
NASA NTRS1982

An experience in mesh generation for three-dimensional calculation of potential flow around a rotating propeller

An attempt is made to develop a three-dimensional, finite volume computational code for highly swept, twisted, small aspect ratio propeller blades with supersonic tip speeds, in a way that accounts for cascade effects, hub-induced flow, and nonlinear transonic effects. Attention is presently given to the generation of a computational mesh for such a complex propeller configuration, with the aim of sharing developmental process experience. The problem treated is unique, in that blade chord, blade length, hub length and blade-to-blade distance represent several characteristic length scales among which there is considerable disparity. An ad hoc mesh-generation scheme is accordingly developed.

Jou, W.-H.↗
NASA NTRS1983

Computational technique for three-dimensional compressible flow past wings at high angles of attack

A computational technique, based on the integral solution of the full potential equation, has been developed for the solution of three-dimensional high-subsonic flows past wings at high angles of attack. The problem includes two sources of nonlinearities: a boundary oriented nonlinearity (separated flow roll up) and a region oriented nonlinearity (flow compressibility). The former is represented by a nonlinear vortex lattice while the latter is represented by a source distribution inside a finite volume. The solution is obtained by using double iteration cycles; a separated flow (wake) iteration cycle and a compressibility iteration cycle. The computational technique is applied to a delta wing. The results show that the technique is very promising and efficient.

Kandil, O. A.↗
NASA NTRS1983

Computational methods for ideal compressible flow

Conservative dissipative difference schemes for computing one dimensional flow are introduced, and the recognition and representation of flow discontinuities are discussed. Multidimensional methods are outlined. Second order finite volume schemes are introduced. Conversion of difference schemes for a single linear convection equation into schemes for the hyperbolic system of the nonlinear conservation laws of ideal compressible flow is explained. Approximate Riemann solvers are presented. Monotone initial value interpolation; and limiters, switches, and artificial dissipation are considered.

Vanleer, B.↗
NASA NTRS1983

Three-dimensional unsteady Euler equation solutions using flux vector splitting

A method for numerically solving the three dimensional unsteady Euler equations using flux vector splitting is developed. The equations are cast in curvilinear coordinates and a finite volume discretization is used. An explicit upwind second-order predictor-corrector scheme is used to solve the discretized equations. The scheme is stable for a CFL number of two and local time stepping is used to accelerate convergence for steady-state problems. Characteristic variable boundary conditions are developed and used in the far field and at surfaces. No additional dissipation terms are included in the scheme. Numerical results are compared with results from an existing three dimensional Euler code and experimental data.

Whitfield, D. L.↗
NASA NTRS1984

Error reduction program: A progress report

Five finite differences schemes were evaluated for minimum numerical diffusion in an effort to identify and incorporate the best error reduction scheme into a 3D combustor performance code. Based on this evaluated, two finite volume method schemes were selected for further study. Both the quadratic upstream differencing scheme (QUDS) and the bounded skew upstream differencing scheme two (BSUDS2) were coded into a two dimensional computer code and their accuracy and stability determined by running several test cases. It was found that BSUDS2 was more stable than QUDS. It was also found that the accuracy of both schemes is dependent on the angle that the streamline make with the mesh with QUDS being more accurate at smaller angles and BSUDS2 more accurate at larger angles. The BSUDS2 scheme was selected for extension into three dimensions.

Syed, S. A.↗
NASA NTRS1984

A compressible Navier-Stokes code for turbulent flow modeling

An implicit, finite volume code for solving two dimensional, compressible turbulent flows is described. Second order upwind differencing of the inviscid terms of the equations is used to enhance stability and accuracy. A diagonal form of the implicit algorithm is used to improve efficiency. Several zero and two equation turbulence models are incorporated to study their impact on overall flow modeling accuracy. Applications to external and internal flows are discussed.

Coakley, T. J.↗
NASA NTRS1984

Investigation of the three-dimensional flow field within a transonic fan rotor: Experiment and analysis

An experimental investigation of the three dimensional flow field through a low aspect ratio, transonic, axial flow fan rotor has been conducted using an advanced laser anemometer (LA) system. Laser velocimeter measurements of the rotor flow field at the design operating speed and over a range of through flow conditions are compared to analytical solutions. The numerical technique used herein yields the solution to the full, three dimensional, unsteady Euler equations using an explicit time marching, finite volume approach. The numerical analysis, when coupled with a simplified boundary layer calculation, generally yields good agreement with the experimental data. The test rotor has an aspect ratio of 1.56, a design total pressure ratio of 1.629 and a tip relative Mach number of 1.38. The high spatial resolution of the LA data matrix (9 radial by 30 axial by 50 blade to blade) permits details of the transonic flow field such as shock location, turning distribution and blade loading levels to be investigated and compared to analytical results.

Pierzga, M. J.↗
NASA NTRS1983

Numerical approximations to nonlinear conservation laws with locally varying time and space grids

Numerical approximations to the initial value problem for nonlinear systems of conservation laws are considered. The considered system is said to be hyperbolic when all eigenvalues of every real linear combination of the Jacobian matrices are real. Solutions may develop discontinuities in finite time, even when the initial data are smooth. In the investigation, explicit finite difference methods which use locally varying time grids are considered. The global CFL restriction is replaced by a local restriction. The numerical flux function is studied from a finite volume viewpoint, and a differencing technique is developed at interface points between regions of distinct time increments.

Osher, S.↗