Higher-order compact finite difference scheme

From HandWiki
The nine-point HOC stencil

High-order compact finite difference schemes are used for solving third-order differential equations created during the study of obstacle boundary value problems. They have been shown to be highly accurate and efficient. They are constructed by modifying the second-order scheme that was developed by Noor and Al-Said in 2002. The convergence rate of the high-order compact scheme is third order, the second-order scheme is fourth order.[1]

Differential equations are essential tools in mathematical modelling. Most physical systems are described in terms of mathematical models that include convective and diffusive transport of some variables. Finite difference methods are amongst the most popular methods that have been applied most frequently in solving such differential equations. A finite difference scheme is compact in the sense that the discretised formula comprises at most nine point stencils which includes a node in the middle about which differences are taken. In addition, greater order of accuracy (more than two) justifies the terminology 'higher-order compact finite difference scheme' (HOC). This can be achieved in several ways. The higher-order compact scheme considered here [2] is by using the original differential equation to substitute for the leading truncation error terms in the finite difference equation. Overall, the scheme is found to be robust, efficient and accurate for most computational fluid dynamics (CFD) applications discussed here further.

The simplest problem for the validation of the numerical algorithms is the Lid Driven cavity problem. Computed results in form of tables, graphs and figures for a fluid with Prandtl number = 0.71 with Rayleigh number (Ra) ranging from 103 to 107 are available in the literature.[2] The efficacy of the scheme is proved when it very clearly captures the secondary and tertiary vortices at the sides of the cavity at high values of Ra.

Another milestone was the development of these schemes for solving two dimensional steady/unsteady convection diffusion equations. A comprehensive study of flow past an impulsively started circular cylinder was made.[3] The problem of flow past a circular cylinder has continued to generate tremendous interest[clarification needed] amongst researchers working in CFD mainly because it displays almost all the fluid mechanical phenomena for incompressible, viscous flows in the simplest of geometrical settings. It was able to analyze and visualize the flow patterns more accurately for Reynold's number (Re) ranging from 10 to 9500 compared to the existing numerical results. This was followed by its extension to rotating counterpart of the cylinder surface for Re ranging from 200 to 1000.[4] More complex phenomenon that involves a circular cylinder undergoing rotational oscillations while translating in a fluid is studied for Re as high as 500.[5] [6]

Another benchmark in the history is its extension to multiphase flow phenomena. Natural processes such as gas bubble in oil, ice melting, wet steam are observed everywhere in nature. Such processes also play an important role with the practical applications in the area of biology, medicine, environmental remediation. The scheme has been successively implemented to solve one and two dimensional elliptic and parabolic equations with discontinuous coefficients and singular source terms.[7] These type of problems hold importance numerically because they usually lead to non-smooth or discontinuous solutions across the interfaces. Expansion of this idea from fixed to moving interfaces with both regular and irregular geometries is currently going on.[8][9]

References

  1. Xie, S.; Li, P.; Gao, Z.; Wang, H. (2012). "High order compact finite difference schemes for a system of third order boundary value problem". Applied Mathematics and Computation 219 (5): 2564. doi:10.1016/j.amc.2012.08.091. 
  2. 2.0 2.1 Kalita JC, Dalal DC and Dass AK., A class of higher-order compact schemes for the unsteady two-dimensional convection-diffusion equations with variable convection coefficients., Int. J. Numer. Meth. Fluids, Vol. 101,(2002), pp. 1111–1131
  3. J. C. Kalita, and R. K. Ray., A transformation-free HOC scheme for incompressible viscous flows past an impulsively started circular cylinder, Int. J. Numer. Meth. Fluids, Vol. 228,(2009), pp. 5207–5236
  4. R. K. Ray., A transformation free HOC scheme for incompressible viscous flow past a rotating and translating circular cylinder, J. Sci. Comput., Vol. 46,(2011), pp. 265–293
  5. H. V. R. Mittal, Rajendra K. Ray and Qasem M. Al-Mdallal, A numerical study of initial flow past an impulsively started rotationally oscillating circular cylinder using a transformation-free HOC scheme, Physics of Fluids, vol. 29, no. 9 (2017), pp. 093603
  6. H. V. R. Mittal, Qasem M. Al-Mdallal and Rajendra K. Ray, Locked-on vortex shedding modes from a rotationally oscillating circular cylinder, Ocean Engineering, vol. 146 (2017), pp. 324-338
  7. Rajendra K. Ray, J. C. Kalita, and A. K. Dass, An efficient HOC scheme for transient convection–diffusion reaction equations with discontinuous coefficients and singular source terms, Proc. Appl. Math. Mech., vol. 7, no. 1(2007), pp. 1025603–1025604
  8. H. V. R. Mittal, Jiten C. Kalita and Rajendra K. Ray, A class of finite difference schemes for interface problems with an HOC approach, International Journal for Numerical Methods in Fluids, vol. 82, no. 9 (2016), pp. 567-606
  9. H. V. R. Mittal, Ray, Rajendra K. Ray, Solving Immersed Interface Problems Using a New Interfacial Points-Based Finite Difference Approach, SIAM Journal on Scientific Computing, vol. 40, no. 3 (2018), pp. A1860-A1883