Accuracy study of the IDO scheme by Fourier analysis

The numerical accuracy of the Interpolated Differential Operator (IDO) scheme is studied with Fourier analysis for the solutions of Partial Differential Equations (PDEs): advection, diffusion, and Poisson equations. The IDO scheme solves governing equations not only for physical variable but also for first-order spatial derivative. Spatial discretizations are based on Hermite interpolation functions with both of them. In the Fourier analysis for the IDO scheme, the Fourier coefficients of the physical variable and the first-order derivative are coupled by the equations derived from the governing equations. The analysis shows the IDO scheme resolves all the wavenumbers with higher accuracy than the fourth-order Finite Difference (FD) and Compact Difference (CD) schemes for advection equation. In particular, for high wavenumbers, the accuracy is superior to that of the sixth-order Combined Compact Difference (CCD) scheme. The diffusion and Poisson equations are also more accurately solved in comparison with the FD and CD schemes. These results show that the IDO scheme guarantees highly resolved solutions for all the terms of fluid flow equations.

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

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

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

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

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

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

[7]  P. Moin,et al.  Turbulence statistics in fully developed channel flow at low Reynolds number , 1987, Journal of Fluid Mechanics.

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

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

[10]  P. Spalart.  Direct simulation of a turbulent boundary layer up to Rθ = 1410 , 1988, Journal of Fluid Mechanics.

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

[12]  J. Ferziger,et al.  An adaptive multigrid technique for the incompressible Navier-Stokes equations , 1989 .

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

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

[15]  Kenji Takizawa,et al.  Multi-dimensional semi-Lagrangian scheme that guarantees exact conservation , 2002 .

[16]  Solutions of Partial Differential Equations with the CIP-BS Method , 2004 .

[17]  YOICHI OGATA,et al.  Multi-dimensional semi-lagrangian characteristic approach to the shallow water equations by the CIP method , 2004, Int. J. Comput. Eng. Sci..

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

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