The aim of this study is to develop a code that is capable of solving three-dimensional compressible flows which are viscous and turbulent, and parallelization of this code. Purpose of parallelization is to obtain a computational efficiency in time respect which enables the solution of complex flow problems in reasonable computational times. In the first part of the study, which is the development of a three-dimensional Navier-Stokes solver for turbulent flows, first step is to develop a two-dimensional Euler code using Roe flux difference splitting method. This is followed by addition of sub programs involving calculation of viscous fluxes. Third step involves implementation of Baldwin-Lomax turbulence model to the code. Finally, the Euler code is generalized to three-dimensions. At every step, code validation is done by comparing numerical results with theoretical, experimental or other numerical results, and adequate consistency between these results is obtained. In the second part, which is the parallelization of the developed code, two-dimensional code is parallelized by using Message Passing Interface (MPI), and important improvements in computational times are obtained.

An axisymmetric, Navier-Stokes finite volume flow solver, which uses Harten, Lax and van Leer (HLL) and Harten, Lax and van Leer&ndash / Contact (HLLC) upwind flux differencing scheme for spatial and uses Runge-Kutta explicit multi-stage time stepping scheme for temporal discretization on unstructured meshe is developed. Developed solver can solve the compressible axisymmetric flow. The spatial accuracy of the solver can be first or second order accurate. Second order accuracy is achieved by piecewise linear reconstruction. Gradients of flow variables required for piecewise linear reconstruction are calculated by Green-Gauss theorem. Baldwin-Lomax turbulent model is used to compute the turbulent viscosity. Approximate Riemann solver of HLL and HLLC implemented in solver are validated by solving a cylindrical explosion case. Also the solver&rsquo / s capability of solving unstructured, multi-zone domain is investigated by this problem. First and second order results of solver are compared by solving the flow over a circular bump. Axisymmetric flow in solid propellant rocket motor is solved in order to validate the axisymmetric feature of solver. Laminar flow over flat plate is solved for viscous terms validation. Turbulent model is studied in the flow over flat plate and flow with mass injection test cases.

In this thesis we describe the development of Reynolds- averaged Navier Stokes code for the flow past two- dimensional configuration. Particularly, emphasis has been laid on the study of low Reynolds number airfoil aerodynamics. The thesis consists of five chapters covering the back ground history, problem formulation, method of solution and discussion of the results and conclusion. Chapter I deals with a detailed background history of low Reynolds number aerodynamics, problem associated with it, state of the art, its importance in practical applications in aircraft industries. Chapter II describes the mathematical model of the flow physics and various levels of approximations. Also it gives an account of complexity of the equations at low Reynolds number regarding flow separation, transition and reattachment. Chapter III describes method of solution, numerical algorithm developed, description of various upwind schemes, grid system, finite volume discrieti-zation of the governing equations described in Chapter II. Chapter IV describes the application of the newly developed Navier Stokes code for the test cases from GAMM Workshop proceedings. Also it describes validation of the code for Euler solutions, Blasius solution for the flow past flat plate and compressible Navier Stokes solution for the flow past NACA 0012 Airfoil at low Reynolds number. Chapter V describes the application of the Navier Stokes code for the more test cases of current practical interest . In this chapter laminar separation bubble characteristics are investigated in detail regarding formation, growth and shedding in an unsteady environment. Finally the conclusion is drawn regarding the robustness of the newly developed code in predicting the airfoil aerodynamic characteristics at low Reynolds number both in steady and unsteady environment. Lastly, suggestion for future work has been highlighted.

For turbulent flow calculations, some of the well-known turbulence models in the literature are applied on a previously developed Navier-Stokes solver designed to handle laminar flows. A finite volume formulation, which is cell-based for inviscid terms and cell-vertex for viscous terms, is used for numerical discretization of the Navier-Stokes equations in conservative form. This formulation is combined with one-step, explicit time marching Lax-Wendroff numerical scheme that is second order accurate in space. To minimize non-physical oscillations resulting from the numerical scheme, second and fourth order artificial smoothing terms are added. To increase the convergence rate of the solver, local time stepping technique is applied. Before applying turbulence models, Navier-Stokes solver is tested for a case of subsonic, laminar flow over a flat plate. The results are in close agreement with Blasius similarity solutions. To calculate turbulent flows, Boussinesq eddy-viscosity approach is utilized. The eddy viscosity (also called turbulent viscosity), which arises as a consequence of this approach, is calculated using Cebeci-Smith, Michel et. al., Baldwin-Lomax, Chien&rsquo / s k-epsilon and Wilcox&rsquo / s k-omega turbulence models. To evaluate the performances of these turbulence models and to compare them with each other, the solver has been tested for a case of subsonic, laminar - transition fixed - turbulent flow over a flat plate. The results are verified by analytical solutions and empirical correlations.