Conservative form of interpolated differential operator scheme for compressible and incompressible fluid dynamics

The proposed scheme, which is a conservative form of the interpolated differential operator scheme (IDO-CF), can provide high accurate solutions for both compressible and incompressible fluid equations. Spatial discretizations with fourth-order accuracy are derived from interpolation functions locally constructed by both cell-integrated values and point values. These values are coupled and time-integrated by solving fluid equations in the flux forms for the cell-integrated values and in the derivative forms for the point values. The IDO-CF scheme exactly conserves mass, momentum, and energy, retaining the high resolution more than the non-conservative form of the IDO scheme. A direct numerical simulation of turbulence is carried out with comparable accuracy to that of spectral methods. Benchmark tests of Riemann problems and lid-driven cavity flows show that the IDO-CF scheme is immensely promising in compressible and incompressible fluid dynamics studies.

[1]  Datta V. Gaitonde,et al.  Optimized Compact-Difference-Based Finite-Volume Schemes for Linear Wave Phenomena , 1997 .

[2]  Yohsuke Imai,et al.  Stable coupling between vector and scalar variables for the IDO scheme on collocated grids , 2006, J. Comput. Phys..

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

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

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

[6]  Xu-Dong Liu,et al.  Solution of Two-Dimensional Riemann Problems of Gas Dynamics by Positive Schemes , 1998, SIAM J. Sci. Comput..

[7]  H. Joachim Schroll.  Relaxed High Resolution Schemes for Hyperbolic Conservation Laws , 2004, J. Sci. Comput..

[8]  Eitan Tadmor,et al.  Solution of two‐dimensional Riemann problems for gas dynamics without Riemann problem solvers , 2002 .

[9]  Tomoaki Kunugi,et al.  Stability and accuracy of the cubic interpolated propagation scheme , 1997 .

[10]  Takashi Yabe,et al.  Constructing exactly conservative scheme in a non-conservative form , 2000 .

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

[12]  Feng Xiao,et al.  A Simple CIP Finite Volume Method for Incompressible Flows , 2004 .

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

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

[15]  Feng Xiao,et al.  An efficient method for capturing free boundaries in multi‐fluid simulations , 2003 .

[16]  T. Yabe,et al.  Exactly conservative semi-Lagrangian scheme for multi-dimensional hyperbolic equations with directional splitting technique , 2001 .

[17]  T. Yabe,et al.  An Exactly Conservative Semi-Lagrangian Scheme (CIP–CSL) in One Dimension , 2001 .

[18]  T. Yabe,et al.  Conservative and oscillation-less atmospheric transport schemes based on rational functions , 2002 .

[19]  Yohsuke Imai,et al.  Accuracy study of the IDO scheme by Fourier analysis , 2006, J. Comput. Phys..

[20]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[21]  Marcelo H. Kobayashi,et al.  Regular Article: On a Class of Padé Finite Volume Methods , 1999 .

[22]  Atsushi Yashima,et al.  Numerical simulation of flow failure of geomaterials based on fluid dynamics , 2005 .

[23]  Takehiro Himeno,et al.  Numerical Analysis for Propellant Management in Rocket Tanks , 2005 .

[24]  James P. Collins,et al.  Numerical Solution of the Riemann Problem for Two-Dimensional Gas Dynamics , 1993, SIAM J. Sci. Comput..

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

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

[27]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[28]  Feng Xiao,et al.  A New Paradigm of Computer Graphics by Universal Solver for Solid, Liquid and Gas , 2004 .

[29]  Masatoshi Sekine,et al.  Accurate basis set by the CIP method for the solutions of the Schrödinger equation , 2004 .

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

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

[32]  Kenji Takizawa,et al.  Simulation and experiment on swimming fish and skimmer by CIP method , 2005 .

[33]  Marcelo H. Kobayashi,et al.  A fourth-order-accurate finite volume compact method for the incompressible Navier-Stokes solutions , 2001 .

[34]  P. Chu,et al.  A Three-Point Combined Compact Difference Scheme , 1998 .

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

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