Interpolated differential operator (IDO) scheme for solving partial differential equations

Abstract We present a numerical scheme applicable to a wide variety of partial differential equations (PDEs) in space and time. The scheme is based on a high accurate interpolation of the profile for the independent variables over a local area and repetitive differential operations regarding PDEs as differential operators. We demonstrate that the scheme is uniformly applicable to hyperbolic, ellipsoidal and parabolic equations. The equations are solved in terms of the primitive independent variables, so that the scheme has flexibility for various types of equations including source terms. We find out that the conservation holds accurate when a Hermite interpolation is used. For compressible fluid problems, the shock interface is found to be sharply described by adding an artificial viscosity term.

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

[2]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[3]  Takashi Yabe,et al.  A New Higher-Order Godunov Method for General Hyperbolic Equations , 1988 .

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

[5]  H. Takewaki,et al.  The cubic-interpolated Pseudo particle (cip) method: application to nonlinear and multi-dimensional hyperbolic equations , 1987 .

[6]  Feng Xiao,et al.  Description of Complex and Sharp Interface during Shock Wave Interaction with Liquid Drop , 1993 .

[7]  N. Zabusky,et al.  Interaction of "Solitons" in a Collisionless Plasma and the Recurrence of Initial States , 1965 .

[8]  Takashi Yabe,et al.  Implicit CIP (Cubic-Interpolated Propagation) method in one dimension , 1995 .

[9]  R. D. Richtmyer,et al.  A Method for the Numerical Calculation of Hydrodynamic Shocks , 1950 .

[10]  E. Tadmor,et al.  Non-oscillatory central differencing for hyperbolic conservation laws , 1990 .

[11]  Smadar Karni,et al.  Multicomponent Flow Calculations by a Consistent Primitive Algorithm , 1994 .

[12]  Huanan Yang,et al.  An artificial compression method for ENO schemes - The slope modification method. [essentially nonoscillatory , 1990 .

[13]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[14]  A. Harten.  ENO schemes with subcell resolution , 1989 .

[15]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[16]  Takashi Yabe,et al.  Cubic interpolated pseudo-particle method (CIP) for solving hyperbolic-type equations , 1985 .

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

[18]  Y. Kondoh,et al.  Kernel optimum nearly-analytical discretization (KOND) algorithm applied to parabolic and hyperbolic equations , 1992 .