Preconditioned characteristic boundary conditions based on artificial compressibility method for solution of incompressible flows

Abstract The preconditioned characteristic boundary conditions based on the artificial compressibility (AC) method are implemented at artificial boundaries for the solution of two- and three-dimensional incompressible viscous flows in the generalized curvilinear coordinates. The compatibility equations and the corresponding characteristic variables (or the Riemann invariants) are mathematically derived and then applied as suitable boundary conditions in a high-order accurate incompressible flow solver. The spatial discretization of the resulting system of equations is carried out by the fourth-order compact finite-difference (FD) scheme. In the preconditioning applied here, the value of AC parameter in the flow field and also at the far-field boundary is automatically calculated based on the local flow conditions to enhance the robustness and performance of the solution algorithm. The code is fully parallelized using the Concurrency Runtime standard and Parallel Patterns Library (PPL) and its performance on a multi-core CPU is analyzed. The incompressible viscous flows around a 2-D circular cylinder, a 2-D NACA0012 airfoil and also a 3-D wavy cylinder are simulated and the accuracy and performance of the preconditioned characteristic boundary conditions applied at the far-field boundaries are evaluated in comparison to the simplified boundary conditions and the non-preconditioned characteristic boundary conditions. It is indicated that the preconditioned characteristic boundary conditions considerably improve the convergence rate of the solution of incompressible flows compared to the other boundary conditions and the computational costs are significantly decreased.

[1]  M. Y. Hashemi,et al.  Efficient and non-reflecting far-field boundary conditions for incompressible flow calculations , 2014, Appl. Math. Comput..

[2]  Y. Adam,et al.  A Hermitian finite difference method for the solution of parabolic equations , 1975 .

[3]  Andrew J. Majda,et al.  Absorbing Boundary Conditions for Numerical Simulation of Waves , 1977 .

[4]  Philip Beran,et al.  Steady and unsteady solutions of the Navier-Stokes equations for flows about airfoils at low speeds , 1991 .

[5]  Xi-Yun Lu,et al.  Large-eddy simulation of the compressible flow past a wavy cylinder , 2010, Journal of Fluid Mechanics.

[6]  Arnaud G. Malan,et al.  An improved unsteady, unstructured, artificial compressibility, finite volume scheme for viscous incompressible flows: Part II. Application , 2002 .

[7]  Kazem Hejranfar,et al.  Preconditioned characteristic boundary conditions for solution of the preconditioned Euler equations at low Mach number flows , 2012, J. Comput. Phys..

[8]  S. Chakravarthy.  Euler equations - Implicit schemes and boundary conditions , 1983 .

[9]  Tim Colonius,et al.  MODELING ARTIFICIAL BOUNDARY CONDITIONS FOR COMPRESSIBLE FLOW , 2004 .

[10]  Eli Turkel,et al.  Local preconditioners for steady and unsteady flow applications , 2005 .

[11]  E. Turkel,et al.  Preconditioned methods for solving the incompressible low speed compressible equations , 1987 .

[12]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[13]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[14]  Jeffrey L. Young,et al.  Practical aspects of higher-order numerical schemes for wave propagation phenomena , 1999 .

[15]  M. Hafez,et al.  Improved Numerical Simulations of Incompressible Flows Based on Viscous/Inviscid Interaction Procedures , 2007 .

[16]  D. Tritton.  Experiments on the flow past a circular cylinder at low Reynolds numbers , 1959, Journal of Fluid Mechanics.

[17]  Michael B. Giles,et al.  Nonreflecting boundary conditions for Euler equation calculations , 1990 .

[18]  Pseudo-Compressibility Methods for the Incompressible Flow Equations , 1993 .

[19]  S. Dennis,et al.  Numerical solutions for steady flow past a circular cylinder at Reynolds numbers up to 100 , 1970, Journal of Fluid Mechanics.

[20]  K. Lam,et al.  Three-dimensional nature of vortices in the near wake of a wavy cylinder , 2004 .

[21]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[22]  Guowei Yang,et al.  Computational-Fluid-Dynamics Solver with Preconditioned Method for Aerodynamic Simulation of High-Lift Configuration , 2016 .

[23]  K. Hejranfar,et al.  Application of a preconditioned high‐order accurate artificial compressibility‐based incompressible flow solver in wide range of Reynolds numbers , 2018 .

[24]  V. Vatsa,et al.  Pressure updating methods for the steady-state fluid equations , 1995 .

[25]  Antony Jameson,et al.  A new implicit algorithm with multigrid for unsteady incompressible flow calculations , 1995 .

[26]  A. Bayliss,et al.  Radiation boundary conditions for wave-like equations , 1980 .

[27]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[28]  Quasi-three dimensional analysis of global instabilities: onset of vortex shedding behind a wavy cylinder , 2011, Journal of Fluid Mechanics.

[29]  E. Turkel,et al.  PRECONDITIONING TECHNIQUES IN COMPUTATIONAL FLUID DYNAMICS , 1999 .

[30]  Eli Turkel,et al.  Review of preconditioning methods for fluid dynamics , 1993 .

[31]  Stuart E. Rogers,et al.  Upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations , 1990 .

[32]  Kazem Hejranfar,et al.  Assessment of three preconditioning schemes for solution of the two‐dimensional Euler equations at low Mach number flows , 2012 .

[33]  Krishnan Mahesh,et al.  High order finite difference schemes with good spectral resolution , 1997 .

[34]  K. Thompson.  Time-dependent boundary conditions for hyperbolic systems, II , 1990 .

[35]  Herman Deconinck,et al.  Application of time preconditioning and high‐order compact discretization method for low Mach number flows , 2013 .

[36]  G. W. Hedstrom,et al.  Nonreflecting Boundary Conditions for Nonlinear Hyperbolic Systems , 1979 .

[37]  A. S. Grove,et al.  An experimental investigation of the steady separated flow past a circular cylinder , 1964, Journal of Fluid Mechanics.

[38]  Kwok-cheung Lam,et al.  Effects of large spanwise wavelength on the wake of a sinusoidal wavy cylinder , 2016 .

[39]  Kaveh Parseh,et al.  Assessment of Characteristic Boundary Conditions Based on the Artificial Compressibility Method in Generalized Curvilinear Coordinates for Solution of the Euler Equations , 2017, Comput. Methods Appl. Math..

[40]  R. Hirsh,et al.  Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique , 1975 .

[41]  Wayne A. Smith,et al.  Preconditioning Applied to Variable and Constant Density Flows , 1995 .

[42]  Miguel R. Visbal,et al.  High-Order-Accurate Methods for Complex Unsteady Subsonic Flows , 1999 .

[43]  L. K. Taylor,et al.  Numerical solution of the two-dimensional time-dependent incompressible Euler equations , 1994 .

[44]  Anwar Ahmed,et al.  Transverse flow over a wavy cylinder , 1992 .

[45]  F. Hu.  Absorbing Boundary Conditions , 2004 .

[46]  Dimitris E. Papantonis,et al.  A characteristic-based method for incompressible flows , 1994 .

[47]  C. L. Merkle,et al.  The application of preconditioning in viscous flows , 1993 .

[48]  Jonathan G. Tyler.  Analysis and Implementation of High-Order Compact Finite Difference Schemes , 2007 .

[49]  Kazem Hejranfar,et al.  Implementation of a high-order compact finite-difference lattice Boltzmann method in generalized curvilinear coordinates , 2014, J. Comput. Phys..

[50]  P. Nithiarasu.  An efficient artificial compressibility (AC) scheme based on the characteristic based split (CBS) method for incompressible flows , 2003 .

[51]  Arnaud G. Malan,et al.  An improved unsteady, unstructured, artificial compressibility, finite volume scheme for viscous incompressible flows: Part I. Theory and implementation , 2002 .

[52]  A. Majda,et al.  Radiation boundary conditions for acoustic and elastic wave calculations , 1979 .

[53]  R. Bouard,et al.  Experimental determination of the main features of the viscous flow in the wake of a circular cylinder in uniform translation. Part 1. Steady flow , 1977, Journal of Fluid Mechanics.

[54]  D. Gottlieb,et al.  The Stability of Numerical Boundary Treatments for Compact High-Order Finite-Difference Schemes , 1993 .