Analysis of non-conservative interpolation techniques in overset grid finite-volume methods

Abstract This paper reports on the challenges of using an overset grid approach when simulating incompressible flows with a cell-centred finite-volume solver. The focal point is on the grid coupling approach. Usually, overset grid methods for unstructured three-dimensional grids apply a local interpolation of field values onto the partner grid as a coupling procedure. Such local interpolation is obviously insufficient to guarantee that the global sum of the mass fluxes across the overlapping interfaces vanishes. Hence, the inherent mass conservation of finite volume methods is violated by the non-conservative property of the overset grid coupling. Since incompressible finite-volume solvers directly use the mass defect when solving for the pressure, severe pressure fluctuations may be provoked. A sensitivity study of the residual mass defect on overlapping grids and the related pressure fluctuations is performed for a variety of coupling strategies. Emphasis is given to five different aspects, i.e. relative grid motion, the order of accuracy of the employed baseline interpolation, the formulation of global mass conservation enforcing practices, the resolution quality of the overlapping grids and multiphase-flow issues. Examples included refer to a simple 2D cylinder flow, a 2D/3D lid-driven cavity flow, a 2D body force disturbed channel flow and 2D oscillating hydrofoil flows in fully wetted single-phase and immiscible two-phase conditions. Results reveal that transient flows and simulations with relative grid motion are subject to significant mass and pressure fluctuations. High-order interpolation and grid refinement prove beneficial, but can still be afflicted with severe disturbances, particularly when the resolution properties of the overlapping grids disagree. Two-phase flows with a high density ratio show a distinct compensating behaviour, while the introduction of flux or wetted volume correction practices leads to notable improvements for single-phase flows in both homogeneous and heterogeneous resolution conditions.

[1]  Kazuhiro Nakahashi An Intergrid-Boundary Definition Method for Overset Unstructured Grid Approach , 1999 .

[2]  James H. Duncan,et al.  The breaking and non-breaking wave resistance of a two-dimensional hydrofoil , 1983, Journal of Fluid Mechanics.

[3]  William D. Henshaw,et al.  A Scheme for Conservative Interpolation on Overlapping Grids , 1994, SIAM J. Sci. Comput..

[4]  Frederick Stern,et al.  Turn and zigzag maneuvers of a surface combatant using a URANS approach with dynamic overset grids , 2013 .

[5]  M. Berger ON CONSERVATION AT GRID INTERFACES. , 1987 .

[6]  Jun Shao,et al.  Quantitative V&V of CFD simulations and certification of CFD codes , 2006 .

[7]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[8]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[9]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[10]  J. Azevedo,et al.  Three-dimensional viscous flow simulations over the VLS using overset grids , 2004 .

[11]  Rainald Lhner Applied Computational Fluid Dynamics Techniques , 2008 .

[12]  Sergey Yakubov,et al.  Experience using pressure-based CFD methods for Euler–Euler simulations of cavitating flows , 2015 .

[13]  J. Benek,et al.  A flexible grid embedding technique with application to the Euler equations , 1983 .

[14]  Z. J. Wang,et al.  A Fully Conservative Interface Algorithm for Overlapped Grids , 1995 .

[15]  A. Popel,et al.  Application of Chimera grid to modelling cell motion and aggregation in a narrow tube , 2007, International journal for numerical methods in fluids.

[16]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[17]  C. J. Freitas,et al.  SIMULATION OF FLUID–STRUCTURE INTERACTION USING PATCHED-OVERSET GRIDS , 1999 .

[18]  E. Atta,et al.  Component-adaptive grid interfacing , 1981 .

[19]  Björn Sjögreen,et al.  Conservative and Non-Conservative Interpolation between Overlapping Grids for Finite Volume Solutions of Hyperbolic Problems , 1994 .

[20]  Fotis Sotiropoulos,et al.  An overset-grid method for 3D unsteady incompressible flows , 2003 .

[21]  I. Demirdzic,et al.  Space conservation law in finite volume calculations of fluid flow , 1988 .

[22]  R. Rannacher,et al.  Benchmark Computations of Laminar Flow Around a Cylinder , 1996 .

[23]  O. Ubbink Numerical prediction of two fluid systems with sharp interfaces , 1997 .

[24]  T. Rung,et al.  Hybrid MPI/OpenMP parallelization of an Euler–Lagrange approach to cavitation modelling , 2013 .