Numerical Solution of Compressible Steady Flows around the NACA 0012

The article presents results of a numerical solution of subsonic and transonic flows described by the system of Navier-Stokes equations in 2D laminar compressible flows around the NACA 0012 airfoil. Authors used Runge-Kutta method to numerically solve the flows around the NACA 0012 airfoil.


Introduction
A numerical code has been developed for simulating transonic flow field around the NACA 0012 airfoil.In these simulations the meshes type C has been used, which was created by mesh generator and has been described in [1][2].In this case mesh has been created for numerical solution over profile NACA 0012.

Mathematical models
The 2D flow of a viscid compressible fluid (the Newtonian fluid) is described by the system of Navier-Stokes equations where and [3]     In the above equations, W is conservative q is heat flux (given by Fourier's Law assuming Mayer's formula),  represents dynamical viscosity (calculated by Sutherland formula [3]) and Pr is laminar Prandtl number.
All variables were considered dimensionless and inflow variables (with the suffix ∞) were used as the reference variable for the transformation of the equations to dimensionless form.Thus after the transformation the system of Navier-Stokes equations is written where Re  is a Reynolds number.

Specification of test case
We selected for numerical solution a structured mesh formed by quadrilateral finite volumes.The grid around profiles (wing) usually consists of a C-grid in the flow direction.In the case of the C-topology the aerodynamics body is enclosed by one family grid lines, which also form the wake region.The situation is sketched in figure 1.The initial grid is generated algebraically by using the linear TFI method [3].Afterwards, elliptic PDE's are employed to produce boundary-orthogonal grid with specific wall spacing.The NACA0012 airfoil contour is approximated by a Bezier spline [4].
The grid with 640x64 elements created by ours program for this simulation were presented in figure 2. The detail of grid near wall was presented in Fig. 3.We used also another grid that had 192x64 elements.

Numerical methods
For modeling of the mentioned flow case, numerical scheme of the finite volume method Runge-Kutta (RK) is used on non-orthogonal structured grids of quadrilateral cells ij D . , Re 0, 25; 0,1667; 0,375; 0,5; 1 where , i j D represents the surface or volume of the cell.
Numerical approximations F  , G  of convective terms F , G were considered in the forward form of the first order of accuracy.Numerical approximations R  , S  of dissipative terms R , S were approximated in the central form of the second order of accuracy and by using dual cells and applying Green's formula.

01121-p.2 EFM 2012
The scheme was extended to include Jameson's artificial dissipation because of the stability of the method The Jameson's artificial dissipation damps undesirable oscillations and improves the stability of the method.The convergence to the steady state is followed by log L2 residual defined by 1 Rez where N is a number of all elements in the computational domain.

Boundary conditions
For supersonic inflow the conservative variables on the boundary are determined by freestream values only.The values are specified based on the given Mach number Ma  and two flow angles (angle of attack and side-slip angle).
For supersonic outflow all conservative variables at the boundary must be determined from the solution inside domain.
Three characteristic variables are prescribed based on the free stream values at far field subsonic inflow.One characteristic variable is extrapolated from interior of the physical domain.
In the case of subsonic outflow, three flow variables (density and the two velocity components) have to be extrapolated from interior of the physical domain.The remaining four variables (pressure) must be specified externally (as a multiple of the input pressure).
On wall zero derivatives of velocity vector along normal is considered.This condition is implemented by using virtual cells.Such cells adjoin from outside on the boundary cells and we prescribe values of unknowns inside of them to obtain the desired effect.

Numerical results and conclusion
As mentioned above, the authors have dealt with the numerical simulation of 2D laminar transonic flows around the NACA 0012 airfoil and applied the multistage Runge-Kutta method and orthogonal structured grids C type with 640x64 (192x64) cells for its created by the linear TFI method [3].