FVM Interpolation and Gradient Schemes literature review

             FVM:: Interpolation and Gradient Schemes literature review

FVM-Flux limiters and interpolation schemes

What is FVM?

FVM is a finite volume method for representing and evaluating PDE's in the form of algebraic equations which are similar to the FDM or FEM. Values are calculate at discrete places on a meshed geometry. Finite volume refers to the small volume surrounding each node point on a mesh. In this method, volume integrals in a partial differential equation that contains a divergence term are converted to surface integrals,using the divergence theorem. Those terms are then evaluated as fluxes at the surface of each finite volume, these methods are conservative.

This method is helpful if the mesh is not even or unstructured meshes, like a mesh we give in openFoam using simple grading which is not possible in matlab. In short FVM computes the solution at each finite volume where FDM computes the solution at each finite points. Also FVM seems to conserve the properties in a better way than the FDM because of the fact that the governing equations solved in an integral form.

So while writing the governing equation in integral form we come up with suitable methods to calculate the flux terms. In this process we use interpolation and gradient scheme. We use the flux limiters to bound the values.


1. Interpolation scheme:-

The approximation of surface and volume integrals may require values of the variable at locations other than the computational nodes of CV. Values at these nodes are obtained using the interpolation formaule.

For example in 1D linear heat conduction equation in steady state, the source and the alpha  values are functions of temperature, We need to calculate them seperately of temperature to that particular faces. To calcualte the temperature at that face, we use the interpolation scheme which takes the temperature fro the adjoining cells.

We mention few of the possibilities:-

a. Upwind interpolation:-

It is used for calculating approximate value of the variable at the east face of the control volume, it is first order accurate and is numerically diffusive with coefficient of numerical diffusion. 

`f_e={(f_P if (v*n)_e>0 @ f_E if (v*n)_e<0)`

It is equivalant to using a forward or backward finite difference method . The Upwind relation is given below:-


b. Linear interpolation:-

To approximate the value of variable at control volume face centre by two nearest computational nodes location e.

It is equivalant to central differencing scheme of first order derivative. It is second order accurate and produce oscillatory solutions. The rekation for linear interpolation is given below:-


c. Quadratic upwind interpolation:-

It is third order accurate but prone to oscillations it is derived from polinomial fittings. The QUICK interpolation on the uniform cartesian grid is given by


Where D,U,UU denotes the downstream, first upstream and second upstream node respectively (E,P and W or P,E and EE depending on the flow direction)

A quadratic curve is fitted with two upstream nodes and one downstram node but in regions with strong gradients, overshoots and undershoots can occure. This can lead to stability problem in the calculation.

The other interpolation schemes are:-

d. Hybrid interpolation schemes:-(CDS and UDS)

e. Total variable diminishing (TVD) scheme


2. Flux limiters:-

Flux limiters are used in high resolution schemes-Numerical schemes used to solve the problems in science and engineering, particularly fluid dyanamics, described by PDE's. They are used in high resolution schemes, such as the MUSCL scheme, to avoid the spurious oscillations that would otherwise occure with high order spatial discritization schemes due to shocks, discountinuities or sharp changes in the solution domain. Use of flux limiters, together with an appropriate high resolution schemes, make the solutions total variable diminishing(TVD).

The main idea behind the flux limiter schemes is to limit the spatial derivatives to realistic values- for scientific and engineering problems this usually means physically realisable and meaningfull values. They are used in hifh resolution schemes for solving problems described ny PDE's and only come into operation when sharp wave fronts are present. For smoothly changing waves, the flux limiters do not operate and and the spatial derivatives can be represented by higher order approximations without introducing spurious oscillations. Consider the 1D semidiscrete scheme below,


Where `F(u_(i+1/2)) and F(u_(i-1/2))` represents egde fluxes for the ith cell. These edge fluxes can be represented by low and high resolution schemes, the a flux limiter can switch between these schemes depending upon the gradients close to the particular cell, as follows:-



Where `f^(low)`= low precision, high resoluation flux 

`f^(high)`=High precision low resolution flux

`phi(r)`=flux limiter function

and r represents the ratio of successive gradients on the solution mesh.


The limiter function is constrained to be grater than or equal to zero, `phi(r)>=0`. Therefore when the limiter is equal to zero(sharp gradient, opposite slopes and zero gradient), the flux is represented by a low resolulation scheme. similarly when the limiter is equal to 1(smooth solution), it is represented by a high resolutiton schemes. The various limiters have differing swithing  characterstics and are selected according to the particular proble and solution scheme. No particular limiter has been found to work well foe all problems, and a particular choice is usually made on a trial and error basis.


Projects by Saurabh kumar sharma

Explore Tutorial No 1- SI final Introduction:- A spark ignition IC engine model take for my case study is single cylinder model. It is port injected and fuel and air enters in to the combustion chamber as homogenious mixture. This model of our engine runs at 1800 rpm Read more

  Assignment 1: Understanding GT Power   1. Introduction:- GT Suite is multiphysics CAE system simulation software. The foundation of GT suite is a versatile multi-physics plateform for constructing models of general systems based on many underlying Read more

Gearbox Sloshing effect The objective is to perform gear sloshing simulation on ansys fluent for two different gearboxes with different clearance. The User defined function is used which provides the angular velocity of 200rpm. Four simulations were performed:- Case Read more

  Challenge on Combustion Objective:- Perform a combustion simulation on the combuster model and plot the variation of the mass fraction of the different species in the simulation using line probes at different locations of the combuster. You need to plot Read more


Cyclone Separator Challenge



We have to perform an analysis on the cyclone seperator model for  different number of inectioned particles.

1. The four imperical models use Read more

  Gate Valve Parametric Study A gate valve, also known as a sluice valve, is a valve that opens by lifting a barrier(gate) out of the path of the fluid.Gate valves are used to shut off the flow of liquids rather than for flow regulation. Gate valve require Read more

Conjugate Heat Transfer Analysis on a graphics card. Description of a conjugate heat transfer analysis:- Conjugate heat transfer represents the combination of heat transfer in liquids and solids simulteneously. The mode of heat transfer in solids is conduction and con Read more

    Rayleigh Taylor Instability Challenge SOLUTION-1:- Models used for Rayleigh-Taylor instability:- Early-time growth Potential flow Buoyancy-drag balance Energy budget   1. Early-time growth:- It was noted in the taylor (1950) that t Read more

              Exhaust Port Challenge SOLUTION1:- Why and Where a CHT analysis is used? Conjugate Heat Transfer analysis:- Conjugate heat transfer analysis is generally used when there is a teperature variation during heating and co Read more

              Ahmed Body Challenge  SOLUTION1: Ahmed body:- The Ahmed body was originally described By S.R. Ahmed in 1984. Ahmed body is a simplified car body. The shape provides a  model to study geometric effects Read more


The End