Week 9 - FVM Literature Review

Study of various interpolation schemes and flux limiters in FVM.

Finite volume method :

The finite volume method (FVM) is a method for representing and evaluating partial diffrential 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 theorm. 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.

Interpolation Schemes in Finite Volume Method (FVM):

The interpolation schemes used in FVM method are as follows

1. Upwind Interpolation Scheme. (UDS)
2. Central Differencing Approximation. (CDS)
3. Quadratic Upwind Differencing Scheme. (QUICK)
4. Hybrid Interpolation Scheme. (CDS & UDS)
5. Total variation Diminishing Scheme. (TVD)
   
   

1. Upwind Interpolation Scheme(UDS):

In this scheme we use the node in the upwind or upstream of the given node. It is generally used for convection dominated problems (i.e Pe>1). ,

fe={fp,if(v.n)e>0fE,if(v.n)e<0fe={fp,if(v.n)e>0fE,if(v.n)e<0

The values for this at the interface is given by

For v>0;ui−12≈ui−1;ui+12≈uiv>0;ui-12≈ui-1;ui+12≈ui
Ic≈vui−ui−1ΔxIc≈vui-ui-1Δx - (Backward differencing scheme)

For v>0;ui−12≈ui−1;ui+12≈uiv>0;ui-12≈ui-1;ui+12≈ui
Ic≈vui+1−uiΔxIc≈vui+1-uiΔx - (Forward differencing scheme)

Taylor series expansion is

v.ui+12=v.ui−(vΔx2).(∂(u)∂(x))i+12v.ui+12=v.ui-(vΔx2).(∂(u)∂(x))i+12 where the second term has the truncation error which resembles a diffusive Flux term.

Upwind schemes are the approximations that satisfy the boundedness unconditionally and is First order accurate in space.

2.Central Differencing Approximation (CDS):

Approximate values of the variable at the control volume (CV) face centre by the linear interpolation of the values at the two nearest computational nodes.

fe=fEλe+fp(1−λe)fe=fEλe+fp(1-λe) where λe=xe−xpxE−xPλe=xe-xpxE-xP

The linear interpolation is equivalent to that of using central difference formula of the first order derivative and hence this scheme is also termed as Central Differencing Scheme(CDS).

3. Interpolation Polynomial:

p1(x)=uL.xR−xxr+x+uR.x−xLxR−xLp1(x)=uL.xR-xxr+x+uR.x-xLxR-xL

The interface values are given by:

ui−12≈ui−1+ui2ui-12≈ui-1+ui2 and ui+12≈ui+ui+12ui+12≈ui+ui+12

Ic≈v.ui+1−ui−12ΔxIc≈v.ui+1-ui-12Δx - (Central differencing scheme)

Taylor Series Approximation Polynomial:

ui+12=ui+ui+12−((Δx)28).(∂2(u)∂(x)2)i+12+...ui+12=ui+ui+12-((Δx)28).(∂2(u)∂(x)2)i+12+...

The approximation is second order precise. It is the simplest and most widely used interpolation technique for the evaluation of gradients required to measure diffusive fluxes.

FLUX LIMITERS:

Flux limiters are used in the numerical schemes to solve problems of fluid dynamics, described by highly coupled non-linear Partial Differential Equations. If we tend to use low order numerical schemes for the solving of the governing equations then we might get a highly oscillatory and unstable solution near the discontinuity if present. Whereas if we tend to use higher-order numerical schemes to capture the phenomenon more accurately we tend to exploit the computational power unnecessarily. The main idea behind the use of the flux limiters is to limit the spatial derivative terms near the discontinuities to physical real values so as to maintain physical consistency. There are several types of flux/slope limiter functions. They come in handy when sharp wave-fronts/peaks are present. For smoothly changing waves/functions, the flux limiters do not operate and the spatial derivatives can be represented normally.

It is imperative to note that flux limiters are also referred to as slope limiters because they both have the same mathematical form, and both have the effect of limiting the solution gradient near shocks or discontinuities.

F(ui+12)=flow(i+12)−ϕ(ri)(flow(i+12)−fhigh(i+12))F(ui+12)=flow(i+12)-ϕ(ri)(flow(i+12)-fhigh(i+12))

F(ui−12)=flow(i−12)−ϕ(ri−1)(flow(i−12)−fhigh(i−12))F(ui-12)=flow(i-12)-ϕ(ri-1)(flow(i-12)-fhigh(i-12))

Hereflowflow = Low precision flux (1st order precise)

fhighfhigh= High precision flux (Highest order precision)

ϕ(r)ϕ(r) = Flux limiter function and r=ui−ui−1ui+1−uir=ui-ui-1ui+1-ui.