Chapter Two Literature Survey
Chapter Two Literature Survey
Chapter Two Literature Survey
Literature survey
In this chapter general review which deal with numerical method of flow over a body are presented, the review is classified into: 1- Grid generation and body fitted coordinates. 2- Numerical method in solving flow fields problems for compressible viscous and non-viscous fluid flow. 3- Boundary condition. 4- Turbulence modeling. 5- Artificial dissipation. 6-The scope of the present work.
2.1
simple uniform boundary, a numerical procedure is needed with regular rectangular grids . Since engineers very rarely come across such boundaries in reality, it is then necessary to find out a procedure which can map the irregular physical reality into the rectangular regions that can easily be used to realize the computational scheme. Generating a regular grid over an irregular computation region makes the specification of boundary conditions easier, which influences the accuracy of the solution significantly regardless of the approximation accuracy inside the computation region. The use of curvilinear coordinate system with coordinate lines coincident with the boundary of physical domain is called body fitted coordinate, which is a coordinate system that coincides with body
xx + yy = P ( , )
xx + yy = Q( , )
(2-1) (2-2)
P = a( x . y y . x ) 2 Q = b( x . y y . x ) 2
(2-3)
Where a and b are constants. The results of generation potential flow about different airfoils showed excellent comparison with the analytic solutions. Barfield et al [8] achieved control of the interior coordinate line spacing by solving the following system of elliptic PDEs.
X + aX = 0 X + bX = 0
constants a and b.
(2-4)
Then the interior coordinate line spacing is changed by varying the Many researchers proposed a method for controlling grid using different grid clustering[9][10][11]. Hyperbolic schemes required only the inner and outer boundaries for closed domain. Marching outward from the inner boundary until the outer boundary is reached generates the grid. The procedure is non iterative but provides no control over the position and shape of outer boundary. Steger and Chausse [12] presented a hyperbolic grid generation scheme for airfoil. The parabolic grid generation which combines the attributes of both elliptic and hyperbolic grid generation, was derived by Nakamura [13]. It was found that this procedure is suitable for coordinate line spacing and orientation and has been applied successfully to the simple wing-body configuration.
2.2
Numerical Method
Number of investigators proposed different numerical approaches
for the prediction of viscous and inviscid flow over different configuration hence computation of fluid dynamics schemes (CFD) contains a wide range of applications. However the survey here is limited to the compressible supersonic viscous and inviscid flows. The use of numerical methods in solving the flow problems started early 1950. It was first used by Courant et al. [2]. In recent years there is a considerable number of solutions to some equations of fluid flow in regions with curved or complex shaped boundaries. Contributions have come from different branches of physics and engineering. Several method have been presented to treat irregular boundaries fluid flow problems [1726]. Lax and Wendroff [27], [28] and [29] solved time Euler and compressible Navier-Stocks equations . The second-order Lax-Wendroff method has led to a whole family of variants when applied to non-linear systems, the most propular of these variants is due to MacCormack [30]. In the present work MacCormack explicit technique is used . Courant et al [2] has introduced for the first time the explicit up-wind scheme . Later several extensions to second order accuracy and implicit time integration have been developed . D.Migdal, et al. [31] found that where conservative form of the Euler equations required fine mesh size, large computer storage and long computational time are also needed. J.D, Anderson Jr., [32] used explicit finite difference MacCormak predictor-Corrector method for unsteady quasi-one-dimensional equations in conservative form . R.A Serra [33] developed a time-dependant numerical method for the solution of inviscid mixed subsonic-supersonic flows with shocks.
10
11
2.3
Boundary conditions
In order to obtain a unique solution of a partial differential equation
a set of supplementing conditions must be provided to determines the arbitrary function. The supplementary conditions are classified as boundary and initial conditions. Initial conditions are required for which the dependent variable is specified at some initial state. Boundary conditions are required for which the dependent or its derivative must be satisfied on the boundaries of the domain. The types of boundary condition are [44]. 1234The Dirichlet boundary conditions if the dependent variable along the boundary is prescribed. The Neumann boundary conditions if the normal gradient of the dependent variable along the boundary is specified. The Robain boundary condition, if the imposed boundary condition is a linear combination of the Dirichlet and Neumann types The mixed boundary condition if a certain portion of domain is the DirIchlet B.C. and another portion is a Neumann B.C.This is known as a mixed B.C. The B.C. treatment are perhaps the most important factor for accuracy of the computation flow results. In general the boundary condition can be classified according to the type of the flow, viscous on
12
13
14
15
16