Development of an incompressible navier-stokes solver with alternating cell direction implicit method on structured and unstructured quadrilateral grids

Download
2007
Baş, Onur
In this research, the Alternating Cell Direction Implicit method is used in temporal discretisation of the incompressible Navier-Stokes equations and compared with the well known and widely used Point Gauss Seidel scheme on structured and quadrilateral unstructured meshes. A two dimensional, laminar and incompressible Navier-Stokes solver is developed for this purpose using the artificial compressibility formulation. The developed solver is used to obtain steady-state solutions with implicit time stepping methods and a third order data reconstruction scheme (U-MUSCL) is added to obtain high order spatial accuracy. The Alternating Cell Directions Implicit method and Point Gauss Seidel scheme is compared in terms of convergence iteration number and total computation time using test cases with growing complexity, including laminar flat plate, single and multi-element airfoil calculations. Both structured and quadrilateral unstructured grids are used in single element airfoil calculations. In these test cases, it is seen that a reduction between 13% and 20% is obtained in total computation time by usage of Alternating Cell Directions Implicit method when compared with the Point Gauss Seidel method.

Suggestions

Implementation of different flux evaluation schemes into a two-dimensional Euler solver
Eraslan, Elvan; Aksel, Mehmet Haluk; Department of Mechanical Engineering (2006)
This study investigates the accuracy and efficiency of several flux splitting methods for the compressible, two-dimensional Euler equations. Steger-Warming flux vector splitting method, Van Leer flux vector splitting method, The Advection Upstream Splitting Method (AUSM), Artificially Upstream Flux Vector Splitting Scheme (AUFS) and Roe’s flux difference splitting schemes were implemented using the first- and second-order reconstruction methods. Limiter functions were embedded to the second-order reconstruc...
A discontinuous subgrid eddy viscosity method for the time-dependent Navier-Stokes equations
Kaya Merdan, Songül (Society for Industrial & Applied Mathematics (SIAM), 2005-01-01)
In this paper we provide an error analysis of a subgrid scale eddy viscosity method using discontinuous polynomial approximations for the numerical solution of the incompressible Navier-Stokes equations. Optimal continuous in time error estimates of the velocity are derived. The analysis is completed with some error estimates for two fully discrete schemes, which are first and second order in time, respectively.
Error analysis for the numerical evaluation of the diagonal forms of the scalar spherical addition theorem
Koc, S; Song, JM; Chew, WC (Society for Industrial & Applied Mathematics (SIAM), 1999-04-29)
The numerical solution of wave scattering from large objects or from a large cluster of scatterers requires excessive computational resources and it becomes necessary to use approximate-but fast-methods such as the fast multipole method; however, since these methods are only approximate, it is important to have an estimate for the error introduced in such calculations. An analysis of the error for the fast multipole method is presented and estimates for truncation and numerical integration errors are obtain...
Two dimensional finite volume weighted essentially non-oscillatory euler schemes with uniform and non-uniform grid coefficients
Elfarra, Monier Ali; Akmandor, İbrahim Sinan; Department of Aerospace Engineering (2005)
In this thesis, Finite Volume Weighted Essentially Non-Oscillatory (FV-WENO) codes for one and two-dimensional discretised Euler equations are developed. The construction and application of the FV-WENO scheme and codes will be described. Also the effects of the grid coefficients as well as the effect of the Gaussian Quadrature on the solution have been tested and discussed. WENO schemes are high order accurate schemes designed for problems with piecewise smooth solutions containing discontinuities. The key ...
Inverse Sturm-Liouville Systems over the whole Real Line
Altundağ, Hüseyin; Taşeli, Hasan; Department of Mathematics (2010)
In this thesis we present a numerical algorithm to solve the singular Inverse Sturm-Liouville problems with symmetric potential functions. The singularity, which comes from the unbounded domain of the problem, is treated by considering the limiting case of the associated problem on the symmetric finite interval. In contrast to regular problems which are considered on a finite interval the singular inverse problem has an ill-conditioned structure despite of the limiting treatment. We use the regularization t...
Citation Formats
O. Baş, “Development of an incompressible navier-stokes solver with alternating cell direction implicit method on structured and unstructured quadrilateral grids,” M.S. - Master of Science, Middle East Technical University, 2007.