A higher-order implicit IDO scheme and its CFD application to local mesh refinement method

The Interpolated Differential Operator (IDO) scheme has been developed for the numerical solution of the fluid motion equations, and allows to produce highly accurate results by introducing the spatial derivative of the physical value as an additional dependent variable. For incompressible flows, semi-implicit time integration is strongly affected by the Courant and diffusion number limitation. A high-order fully-implicit IDO scheme is presented, and the two-stage implicit Runge-Kutta time integration keeps over third-order accuracy. The application of the method to the direct numerical simulation of turbulence demonstrates that the proposed scheme retains a resolution comparable to that of spectral methods even for relatively large Courant numbers. The scheme is further applied to the Local Mesh Refinement (LMR) method, where the size of the time step is often restricted by the dimension of the smallest meshes. In the computation of the Karman vortex street problem, the implicit IDO scheme with LMR is shown to allow a conspicuous saving of computational resources.

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

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

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

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

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

[6]  Yohsuke Imai,et al.  A Numerical Method on Cartesian Grid for Image-based Blood Flow Simulation , 2004 .

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

[8]  Tayfun E. Tezduyar,et al.  The multi-domain method for computation of the aerodynamics of a parachute crossing the far wake of an aircraft , 2001 .

[9]  J. Butcher.  Implicit Runge-Kutta processes , 1964 .

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

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

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

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

[14]  Marek Behr,et al.  Enhanced-Discretization Interface-Capturing Technique (EDICT) for computation of unsteady flows with interfaces , 1998 .

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

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

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