Higher-order compact finite difference scheme

Last updated
The nine-point HOC stencil HOC stencil.jpg
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]

Related Research Articles

<span class="mw-page-title-main">Fluid dynamics</span> Aspects of fluid mechanics involving flow

In physics, physical chemistry and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluids—liquids and gases. It has several subdisciplines, including aerodynamics and hydrodynamics. Fluid dynamics has a wide range of applications, including calculating forces and moments on aircraft, determining the mass flow rate of petroleum through pipelines, predicting weather patterns, understanding nebulae in interstellar space and modelling fission weapon detonation.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

In physics, the Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances, named after French engineer and physicist Claude-Louis Navier and Anglo-Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842-1850 (Stokes).

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

<span class="mw-page-title-main">Finite volume method</span>

The finite volume method (FVM) is a method for representing and evaluating partial differential equations in the form of algebraic equations. In the finite volume method, volume integrals in a partial differential equation that contain a divergence term are converted to surface integrals, using the divergence theorem. These terms are then evaluated as fluxes at the surfaces of each finite volume. Because the flux entering a given volume is identical to that leaving the adjacent volume, these methods are conservative. Another advantage of the finite volume method is that it is easily formulated to allow for unstructured meshes. The method is used in many computational fluid dynamics packages. "Finite volume" refers to the small volume surrounding each node point on a mesh.

<span class="mw-page-title-main">D'Alembert's paradox</span>

In fluid dynamics, d'Alembert's paradox is a contradiction reached in 1752 by French mathematician Jean le Rond d'Alembert. D'Alembert proved that – for incompressible and inviscid potential flow – the drag force is zero on a body moving with constant velocity relative to the fluid. Zero drag is in direct contradiction to the observation of substantial drag on bodies moving relative to fluids, such as air and water; especially at high velocities corresponding with high Reynolds numbers. It is a particular example of the reversibility paradox.

<span class="mw-page-title-main">Level-set method</span>

Level-set methods (LSM) are a conceptual framework for using level sets as a tool for numerical analysis of surfaces and shapes. The advantage of the level-set model is that one can perform numerical computations involving curves and surfaces on a fixed Cartesian grid without having to parameterize these objects. Also, the level-set method makes it very easy to follow shapes that change topology, for example, when a shape splits in two, develops holes, or the reverse of these operations. All these make the level-set method a great tool for modeling time-varying objects, like inflation of an airbag, or a drop of oil floating in water.

Fluid mechanics is the branch of physics concerned with the mechanics of fluids and the forces on them. It has applications in a wide range of disciplines, including mechanical, aerospace, civil, chemical and biomedical engineering, geophysics, oceanography, meteorology, astrophysics, and biology.

<span class="mw-page-title-main">Sergei K. Godunov</span>

Sergei Konstantinovich Godunov is a Soviet and Russian professor at the Sobolev Institute of Mathematics of the Russian Academy of Sciences in Novosibirsk, Russia.

<span class="mw-page-title-main">High-resolution scheme</span> Scheme used in the numerical solution of partial differential equations

High-resolution schemes are used in the numerical solution of partial differential equations where high accuracy is required in the presence of shocks or discontinuities. They have the following properties:

<span class="mw-page-title-main">Bram van Leer</span>

Bram van Leer is Arthur B. Modine Emeritus Professor of aerospace engineering at the University of Michigan, in Ann Arbor. He specializes in Computational fluid dynamics (CFD), fluid dynamics, and numerical analysis. His most influential work lies in CFD, a field he helped modernize from 1970 onwards. An appraisal of his early work has been given by C. Hirsch (1979)

Pressure-correction method is a class of methods used in computational fluid dynamics for numerically solving the Navier-Stokes equations normally for incompressible flows.

In computational physics, the term upwind scheme typically refers to a class of numerical discretization methods for solving hyperbolic partial differential equations, in which so-called upstream variables are used to calculate the derivatives in a flow field. That is, derivatives are estimated using a set of data points biased to be more "upwind" of the query point, with respect to the direction of the flow. Historically, the origin of upwind methods can be traced back to the work of Courant, Isaacson, and Rees who proposed the CIR method.

In computational fluid dynamics, shock-capturing methods are a class of techniques for computing inviscid flows with shock waves. The computation of flow containing shock waves is an extremely difficult task because such flows result in sharp, discontinuous changes in flow variables such as pressure, temperature, density, and velocity across the shock.

<span class="mw-page-title-main">Volume of fluid method</span> Free-surface modelling technique

In computational fluid dynamics, the volume of fluid (VOF) method is a free-surface modelling technique, i.e. a numerical technique for tracking and locating the free surface. It belongs to the class of Eulerian methods which are characterized by a mesh that is either stationary or is moving in a certain prescribed manner to accommodate the evolving shape of the interface. As such, VOF is an advection scheme—a numerical recipe that allows the programmer to track the shape and position of the interface, but it is not a standalone flow solving algorithm. The Navier–Stokes equations describing the motion of the flow have to be solved separately. The same applies for all other advection algorithms.

<span class="mw-page-title-main">Coolfluid</span>

COOLFluiD is a component based scientific computing environment that handles high-performance computing problems with focus on complex computational fluid dynamics (CFD) involving multiphysics phenomena.

In applied mathematics, the name finite pointset method is a general approach for the numerical solution of problems in continuum mechanics, such as the simulation of fluid flows. In this approach the medium is represented by a finite set of points, each endowed with the relevant local properties of the medium such as density, velocity, pressure, and temperature.

False diffusion is a type of error observed when the upwind scheme is used to approximate the convection term in convection–diffusion equations. The more accurate central difference scheme can be used for the convection term, but for grids with cell Peclet number more than 2, the central difference scheme is unstable and the simpler upwind scheme is often used. The resulting error from the upwind differencing scheme has a diffusion-like appearance in two- or three-dimensional co-ordinate systems and is referred as "false diffusion". False-diffusion errors in numerical solutions of convection-diffusion problems, in two- and three-dimensions, arise from the numerical approximations of the convection term in the conservation equations. Over the past 20 years many numerical techniques have been developed to solve convection-diffusion equations and none are problem-free, but false diffusion is one of the most serious problems and a major topic of controversy and confusion among numerical analysts.

Fluid motion is governed by the Navier–Stokes equations, a set of coupled and nonlinear partial differential equations derived from the basic laws of conservation of mass, momentum and energy. The unknowns are usually the flow velocity, the pressure and density and temperature. The analytical solution of this equation is impossible hence scientists resort to laboratory experiments in such situations. The answers delivered are, however, usually qualitatively different since dynamical and geometric similitude are difficult to enforce simultaneously between the lab experiment and the prototype. Furthermore, the design and construction of these experiments can be difficult, particularly for stratified rotating flows. Computational fluid dynamics (CFD) is an additional tool in the arsenal of scientists. In its early days CFD was often controversial, as it involved additional approximation to the governing equations and raised additional (legitimate) issues. Nowadays CFD is an established discipline alongside theoretical and experimental methods. This position is in large part due to the exponential growth of computer power which has allowed us to tackle ever larger and more complex problems.

The finite point method (FPM) is a meshfree method for solving partial differential equations (PDEs) on scattered distributions of points. The FPM was proposed in the mid-nineties in, and with the purpose to facilitate the solution of problems involving complex geometries, free surfaces, moving boundaries and adaptive refinement. Since then, the FPM has evolved considerably, showing satisfactory accuracy and capabilities to deal with different fluid and solid mechanics problems.

<span class="mw-page-title-main">Numerical modeling (geology)</span> Technique to solve geological problems by computational simulation

In geology, numerical modeling is a widely applied technique to tackle complex geological problems by computational simulation of geological scenarios.

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. 1 2 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