Stable coupling between vector and scalar variables for the IDO scheme on collocated grids

The Interpolated Differential Operator (IDO) scheme on collocated grids provides fourth-order discretizations for all the terms of the fluid flow equations. However, computations of fluid flows on collocated grids are not guaranteed to produce accurate solutions because of the poor coupling between velocity vector and scalar variables. A stable coupling method for the IDO scheme on collocated grids is proposed, where a new representation of first-order derivatives is adopted. It is important in deriving the representation to refer to the variables at neighboring grid points, keeping fourth-order truncation error. It is clear that accuracy and stability are drastically improved for shallow water equations in comparison with the conventional IDO scheme. The effects of the stable coupling are confirmed in incompressible flow calculations for DNS of turbulence and a driven cavity problem. The introduction of a rational function into the proposed method makes it possible to calculate shock waves with the initial conditions of extreme density and pressure jumps.

[1]  D. Gottlieb,et al.  Numerical analysis of spectral methods , 1977 .

[2]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

[3]  Tayfan E. Tezduyar,et al.  Stabilized Finite Element Formulations for Incompressible Flow Computations , 1991 .

[4]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[5]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[6]  A. A. Amsden,et al.  A simplified MAC technique for incompressible fluid flow calculations , 1970 .

[7]  Clive Temperton,et al.  Semi-Implicit Semi-Lagrangian Integration Schemes for a Barotropic Finite-Element Regional Model , 1986 .

[8]  Clive Temperton,et al.  An Efficient Two‐Time‐Level Semi‐Lagrangian Semi‐Implicit Integration Scheme , 1987 .

[9]  Takayuki Aoki.  3D simulation for falling papers , 2001 .

[10]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[11]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[12]  Takayuki Aoki,et al.  Poisson equation solver with fourth-order accuracy by using interpolated differential operator scheme☆☆☆ , 2002 .

[13]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[14]  Mark L. Wilkins,et al.  Use of artificial viscosity in multidimensional fluid dynamic calculations , 1980 .

[15]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[16]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

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

[18]  Feng Xiao,et al.  Constructing oscillation preventing scheme for advection equation by rational function , 1996 .

[19]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[20]  T. Aoki.  MULTI-DIMENSIONAL ADVECTION OF CIP (CUBIC-INTERPOLATED PROPAGATION) SCHEME , 1995 .

[21]  Takayuki Aoki,et al.  Interpolated differential operator (IDO) scheme for solving partial differential equations , 1997 .

[22]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation. II, Two- and three-dimensional solvers , 1991 .

[23]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[24]  Takayuki Aoki,et al.  A numerical scheme for strong blast wave driven by explosion , 2006 .

[25]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[26]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[27]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[28]  J. Quirk.  An alternative to unstructured grids for computing gas dynamic flows around arbitrarily complex two-dimensional bodies , 1994 .