Computational Fluid Dynamics
Computational Fluid Dynamics
Computational Fluid Dynamics
ON
COMPUTATIONAL
FLUID
DYNAMICS
PRESENTED BY
CHINDU DAS
AJEESH KUMAR.G
III MECHANICAL
LORD JEGANNATH
COLLEGE OF ENGG & TECH
PSN NAGAR
INTRODUCTION
Computational fluid dynamics (CFD) is one of the branches of fluid mechanics that
uses numerical methods and algorithms to solve and analyze problems that involve fluid flows.
Computers are used to perform the millions of calculations required to simulate the interaction of
liquids and gases with surfaces defined by boundary conditions. Even with high-speed
supercomputers only approximate solutions can be achieved in many cases. Ongoing research,
however, may yield software that improves the accuracy and speed of complex simulation
scenarios such as transonic or turbulent flows. Initial validation of such software is often
performed using a wind tunnel with the final validation coming in flight tests.
In the two-dimensional realm, quite a number of Panel Codes have been developed for
airfoil analysis and design. These codes typically have a boundary layer analysis included, so
that viscous effects can be modeled. Professor Richard Eppler of the University of Stuttgart
developed the PROFIL code, partly with NASA funding, which became available in the early
1980s. This was soon followed by MIT Professor Mark Drela's XFOIL code. Both PROFIL and
XFOIL incorporate two-dimensional panel codes, with coupled boundary layer codes for airfoil
analysis work. PROFIL uses a conformal transformation method for inverse airfoil design, while
XFOIL has both a conformal transformation and an inverse panel method for airfoil design. Both
codes are widely used.
An intermediate step between Panel Codes and Full Potential codes were codes that used
the Transonic Small Disturbance equations. In particular, the three-dimensional WIBCO code,
developed by Charlie Boppe of Grumman Aircraft in the early 1980s has seen heavy use.
Developers next turned to Full Potential codes, as panel methods could not calculate the
non-linear flow present at transonic speeds. The first description of a means of using the Full
Potential equations was published by Earll Murman and Julian Cole of Boeing in 1970. Frances
Bauer, Paul Garabedian and David Korn of the Courant Institute at New York University (NYU)
wrote a series of two-dimensional Full Potential airfoil codes that were widely used, the most
important being named Program H. A further growth of Program H was developed by Bob
Melnik and his group at Grumman Aerospace as Grum foil. Antony Jameson, originally at
Grumman Aircraft and the Courant Institute of NYU, worked with David Caughey to develop
the important three-dimensional Full Potential code FLO22 in 1975. Many Full Potential codes
emerged after this, culminating in Boeing's Tranair (A633) code, which still sees heavy use.
The next step was the Euler equations, which promised to provide more accurate
solutions of transonic flows. The methodology used by Jameson in his three-dimensional FLO57
code (1981) was used by others to produce such programs as Lockheed's TEAM program and
IAI/Analytical Methods' MGAERO program. MGAERO is unique in being a structured cartesian
mesh code, while most other such codes use structured body-fitted grids (with the exception of
NASA's highly successful CART3D code, Lockheed's SPLITFLOW code and Georgia Tech's
NASCART-GT). Antony Jameson also developed the three-dimensional AIRPLANE code
(1985) which made use of unstructured tetrahedral grids.
In the two-dimensional realm, Mark Drela and Michael Giles, then graduate students at
MIT, developed the ISES Euler program (actually a suite of programs) for airfoil design and
analysis. This code first became available in 1986 and has been further developed to design,
analyze and optimize single or multi-element airfoils, as the MSES program. MSES sees wide
use throughout the world. A derivative of MSES, for the design and analysis of airfoils in a
cascade, is MISES, developed by Harold "Guppy" Youngren while he was a graduate student at
MIT.
The Navier–Stokes equations were the ultimate target of developers. Two-dimensional
codes, such as NASA Ames' ARC2D code first emerged. A number of three-dimensional codes
were developed (OVERFLOW, CFL3D are two successful NASA contributions), leading to
numerous commercial packages.
Technicalities
The most fundamental consideration in CFD is how one treats a continuous fluid in a
discretized fashion on a computer. One method is to discretize the spatial domain into small cells
to form a volume mesh or grid, and then apply a suitable algorithm to solve the equations of
motion (Euler equations for inviscid, and Navier–Stokes equations for viscous flow). In addition,
such a mesh can be either irregular (for instance consisting of triangles in 2D, or pyramidal
solids in 3D) or regular; the distinguishing characteristic of the former is that each cell must be
stored separately in memory. Where shocks or discontinuities are present, high resolution
schemes such as Total Variation Diminishing (TVD), Flux Corrected Transport (FCT),
Essentially NonOscillatory (ENO), or MUSCL schemes are needed to avoid spurious oscillations
(Gibbs phenomenon) in the solution.
If one chooses not to proceed with a mesh-based method, a number of alternatives exist,
notably :
• Spectral methods, a technique where the equations are projected onto basis functions like
the spherical harmonics and Chebyshev polynomials,
Methodology
In all of these approaches the same basic procedure is followed.
• During preprocessing
o The geometry (physical bounds) of the problem is defined.
o The volume occupied by the fluid is divided into discrete cells (the mesh). The
mesh may be uniform or non uniform.
o The physical modeling is defined – for example, the equations of motions +
enthalpy + radiation + species conservation
o Boundary conditions are defined. This involves specifying the fluid behaviour and
properties at the boundaries of the problem. For transient problems, the initial
conditions are also defined.
• The simulation is started and the equations are solved iteratively as a steady-state or
transient.
• Finally a postprocessor is used for the analysis and visualization of the resulting solution.
Discretization methods
The stability of the chosen discretization is generally established numerically rather than
analytically as with simple linear problems. Special care must also be taken to ensure that the
discretization handles discontinuous solutions gracefully. The Euler equations and Navier–
Stokes equations both admit shocks, and contact surfaces.
Some of the discretization methods being used are:
• Finite volume method (FVM). This is the "classical" or standard approach used most
often in commercial software and research codes. The governing equations are solved on
discrete control volumes. FVM recasts the PDE's (Partial Differential Equations) of the
N-S equation in the conservative form and then discretize this equation. This guarantees
the conservation of fluxes through a particular control volume. Though the overall
solution will be conservative in nature there is no guarantee that it is the actual solution.
Moreover this method is sensitive to distorted elements which can prevent convergence if
such elements are in critical flow regions. This integration approach yields a method that
is inherently conservative (i.e. quantities such as density remain physically meaningful)
[citation needed]:
where Q is the vector of conserved variables, F is the vector of fluxes (see Euler
equations or Navier–Stokes equations), V is the cell volume, and is the cell surface area.
• Finite element method (FEM). This method is popular for structural analysis of solids,
but is also applicable to fluids. The FEM formulation requires, however, special care to
ensure a conservative solution. The FEM formulation has been adapted for use with the
Navier–Stokes equations. Although in FEM conservation has to be taken care of, it is
much more stable than the FVM. approach[4] Consequently it is the new direction in
which CFD is moving[citation needed]. Generally stability/robustness of the solution is better in
FEM though for some cases it might take more memory than FVM methods.[5]
In this method, a weighted residual equation is formed:
where Ri is the equation residual at an element vertex i , Q is the conservation equation
expressed on an element basis, Wi is the weight factor and Ve is the volume of the
element.
• Finite difference method. This method has historical importance and is simple to
program. It is currently only used in few specialized codes. Modern finite difference
codes make use of an embedded boundary for handling complex geometries making these
codes highly efficient and accurate. Other ways to handle geometries are using
overlapping-grids, where the solution is interpolated across each grid.
Where Q is the vector of conserved variables, and F, G, and H are the fluxes in the x, y,
and z directions respectively.
• Boundary element method. The boundary occupied by the fluid is divided into surface
mesh.
• High-resolution schemes are used where shocks or discontinuities are present. To capture
sharp changes in the solution requires the use of second or higher order numerical
schemes that do not introduce spurious oscillations. This usually necessitates the
application of flux limiters to ensure that the solution is total variation diminishing.
Turbulence models
Turbulent flow produces fluid interaction at a large range of length scales. This problem
means that it is required that for a turbulent flow regime calculations must attempt to take this
into account by modifying the Navier–Stokes equations. Failure to do so may result in an
unsteady simulation. When solving the turbulence model there exists a trade-off between
accuracy and speed of computation.
In fact, statistically unsteady (or non-stationary) flows can equally be treated. This is
sometimes referred to as URANS. There is nothing inherent in Reynolds averaging to preclude
this, but the turbulence models used to close the equations are valid only as long as the time over
which these changes in the mean occur is large compared to the time scales of the turbulent
motion containing most of the energy.
Vortex method
The Vortex method is a grid-free technique for the simulation of turbulent flows. It uses
vortices as the computational elements, mimicking the physical structures in turbulence. Vortex
methods were developed as a grid-free methodology that would not be limited by the
fundamental smoothing effects associated with grid-based methods.
To be practical, however, vortex methods require means for rapidly computing velocities
from the vortex elements – in other words they require the solution to a particular form of the N-
body problem (in which the motion of N objects is tied to their mutual influences). A long-
sought breakthrough came in the late 1980’s with the development of the Fast Multi pole Method
(FMM), an algorithm by V. Rokhlin (Yale) and L. Greengard (Courant Institute) that has been
heralded as one of the top ten advances in numerical science of the 20th century. This
breakthrough paved the way to practical computation of the velocities from the vortex elements
and is the basis of successful algorithms.
Software based on the Vortex method offer the engineer a new means for solving tough
fluid dynamics problems with minimal user intervention. All that is required is specification of
problem geometry and setting of boundary and initial conditions. Among the significant
advantages of this modern technology;
• It is practically grid-free, thus eliminating numerous iterations associated with RANS and
LES.
• All problems are treated identically. No modeling or calibration inputs are required.
• Time-series simulations, which are crucial for correct analysis of acoustics, are possible.
• The small scale and large scale are accurately simulated at the same time.
Solution algorithms
Discretization in space produces a system of ordinary differential equations for unsteady
problems and algebraic equations for steady problems. Implicit or semi-implicit methods are
generally used to integrate the ordinary differential equations, producing a system of (usually)
nonlinear algebraic equations.
Applying a Newton or Picard iteration produces a system of linear equations which is
non-symmetric in the presence of advection and indefinite in the presence of incompressibility.
Such systems, particularly in 3D, are frequently too large for direct solvers, so iterative methods
are used, either stationary methods such as successive over relaxation or Krylov subspace
methods. Krylov methods such as GMRES, typically used with preconditioning, operate by
minimizing the residual over successive subspaces generated by the preconditioned operator.
Multi grid is especially popular, both as a solver and as a pre conditioner, due to its
asymptotically optimal performance on many problems. Traditional solvers and pre conditioners
are effective at reducing high-frequency components of the residual, but low-frequency
components typically require many iterations to reduce. By operating on multiple scales, multi
grid reduces all components of the residual by similar factors, leading to a mesh-independent
number of iterations.
For indefinite systems, pre-conditioners such as incomplete LU factorization, additive
Schwarz, and multi grid perform poorly or fail entirely, so the problem structure must be used for
effective pre-conditioning. The traditional methods commonly used in CFD are the SIMPLE and
Uzawa algorithms which exhibit mesh-dependent convergence rates, but recent advances based
on block LU factorization combined with multi grid for the resulting definite systems, have led
to pre conditioners which deliver mesh-independent convergence rates.