Numerical simulation of unsteady incompressible viscous flows in generalized coordinate systems

Several numerical solutions of the three dimensional unsteady incompressible Navier-Stokes equations in generalized coordinate systems are presented. The governing equations are discretized by finite volumes with special care to the accurate approximation of the geometric quantities. The unknowns are the pressure and the volume fluxes over the computational cell faces. This formulation results in a robust fractional step solution method for solving discrete equations. Although this method is formulated for the three dimensional case, only two dimensional unsteady results are given. Results are presented for a lid driven two dimensional cavity flow at Reynolds number of 10,000, and for the flow over a circular cylinder with vortex shedding for several Reynolds numbers in the range 100 less than Re less than 1000

CFD - Mature Technology?

Over the past 30 years, numerical methods and simulation tools for fluid dynamic problems have advanced as a new discipline, namely, computational fluid dynamics (CFD). Although a wide spectrum of flow regimes are encountered in many areas of science and engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to a large demand for predicting the aerodynamic performance characteristics of flight vehicles, such as commercial, military, and space vehicles. As flow analysis is required to be more accurate and computationally efficient for both commercial and mission-oriented applications (such as those encountered in meteorology, aerospace vehicle development, general fluid engineering and biofluid analysis) CFD tools for engineering become increasingly important for predicting safety, performance and cost. This paper presents the author's perspective on the maturity of CFD, especially from an aerospace engineering point of view

An upwind-differencing scheme for the incompressible Navier-Stokes equations

The steady state incompressible Navier-Stokes equations in 2-D are solved numerically using the artificial compressibility formulation. The convective terms are upwind-differenced using a flux difference split approach that has uniformly high accuracy throughout the interior grid points. The viscous fluxes are differenced using second order accurate central differences. The numerical system of equations is solved using an implicit line relaxation scheme. Although the current study is limited to steady state problems, it is shown that this entire formulation can be used for solving unsteady problems. Characteristic boundary conditions are formulated and used in the solution procedure. The overall scheme is capable of being run at extremely large pseudotime steps, leading to fast convergence. Three test cases are presented to demonstrate the accuracy and robustness of the code. These are the flow in a square-driven cavity, flow over a backward facing step, and flow around a 2-D circular cylinder

Computation of viscous incompressible flows

Incompressible Navier-Stokes solution methods and their applications to three-dimensional flows are discussed. A brief review of existing methods is given followed by a detailed description of recent progress on development of three-dimensional generalized flow solvers. Emphasis is placed on primitive variable formulations which are most promising and flexible for general three-dimensional computations of viscous incompressible flows. Both steady- and unsteady-solution algorithms and their salient features are discussed. Finally, examples of real world applications of these flow solvers are given

Computation of incompressible viscous flows through turbopump components

A finite-difference, three-dimensional, incompressible Navier-Stokes formulation for calculating the flow through turbopump components is presented. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line-relaxation method. Both steady and unsteady flow calculations can be performed using the presented algorithm. In this paper, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using the Baldwin-Lomax model. Numerical results are compared with experimental measurements and good agreement is found between the two. Time-accurate calculations will be reported in future publications

Development of a fractional-step method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A fractional step method is developed for solving the time-dependent three-dimensional incompressible Navier-Stokes equations in generalized coordinate systems. The primitive variable formulation uses the pressure, defined at the center of the computational cell, and the volume fluxes across the faces of the cells as the dependent variables, instead of the Cartesian components of the velocity. This choice is equivalent to using the contravariant velocity components in a staggered grid multiplied by the volume of the computational cell. The governing equations are discretized by finite volumes using a staggered mesh system. The solution of the continuity equation is decoupled from the momentum equations by a fractional step method which enforces mass conservation by solving a Poisson equation. This procedure, combined with the consistent approximations of the geometric quantities, is done to satisfy the discretized mass conservation equation to machine accuracy, as well as to gain the favorable convergence properties of the Poisson solver. The momentum equations are solved by an approximate factorization method, and a novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two- and three-dimensional laminar test cases are computed and compared with other numerical and experimental results to validate the solution method. Good agreement is obtained in all cases

Current CFD Practices in Launch Vehicle Applications

The quest for sustained space exploration will require the development of advanced launch vehicles, and efficient and reliable operating systems. Development of launch vehicles via test-fail-fix approach is very expensive and time consuming. For decision making, modeling and simulation (M&S) has played increasingly important roles in many aspects of launch vehicle development. It is therefore essential to develop and maintain most advanced M&S capability. More specifically computational fluid dynamics (CFD) has been providing critical data for developing launch vehicles complementing expensive testing. During the past three decades CFD capability has increased remarkably along with advances in computer hardware and computing technology. However, most of the fundamental CFD capability in launch vehicle applications is derived from the past advances. Specific gaps in the solution procedures are being filled primarily through "piggy backed" efforts.on various projects while solving today's problems. Therefore, some of the advanced capabilities are not readily available for various new tasks, and mission-support problems are often analyzed using ad hoc approaches. The current report is intended to present our view on state-of-the-art (SOA) in CFD and its shortcomings in support of space transport vehicle development. Best practices in solving current issues will be discussed using examples from ascending launch vehicles. Some of the pacing will be discussed in conjunction with these examples

Computation of Lifting Wing-Flap Configurations

Research has been carried out on the computation of lifting wing-flap configurations. The long term goal of the research is to develop improved computational tools for the analysis and design of high lift systems. Results show that state-of-the-art computational methods are sufficient to predict time-averaged lift and overall flow field characteristics on simple high-lift configurations. Recently there has been an increased interest in the problem of airframe generated noise and experiments carried out in the 7 x 10 wind tunnel at NASA Ames have identified the flap edge as an important source of noise. A follow-on set of experiments will be conducted toward the end of 1995. The computations being carried out under this project are coordinated with these experiments. In particular, the model geometry being used in the computations is the same as that in the experiments. The geometry consists of a NACA 63-215 Mod B airfoil section which spans the 7 x lO tunnel. The wing is unswept and has an aspect ratio of two. A 30% chord Fowler flap is deployed modifications of the flap edge geometry have been shown to be effective in reducing noise and the existing code is currently being used to compute the effect of a modified geometry on the edge flow

Efficient simulation of incompressible viscous flow over multi-element airfoils

The incompressible, viscous, turbulent flow over single and multi-element airfoils is numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The computer code uses the method of pseudo-compressibility with an upwind-differencing scheme for the convective fluxes and an implicit line-relaxation solution algorithm. The motivation for this work includes interest in studying the high-lift take-off and landing configurations of various aircraft. In particular, accurate computation of lift and drag at various angles of attack, up to stall, is desired. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil; an accurate prediction of stall is obtained. The approach used for multi-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared: a patched system of grids, and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time (on a CRAY YMP) per element in the airfoil configuration

High-Fidelity Simulation in Biomedical and Aerospace Engineering

Contents include the following: Introduction / Background. Modeling and Simulation Challenges in Aerospace Engineering. Modeling and Simulation Challenges in Biomedical Engineering. Digital Astronaut. Project Columbia. Summary and Discussion