Week 9 - FVM Literature Review

Finite volume Method:

• The finite volume method (FVM) is a method for representing and evaluating partial differential equations in the form of algebraic equations. Similar to the finite difference method (FDM) or finite element method (FEM), values are calculated at discrete places on a meshed geometry.
• One important feature of finite volume schemes is their conservation properties. Since they are based on applying conservation principles over each small control volume, global conservation is also ensured.
• \"Finite volume\" refers to the small volume surrounding each node point on a mesh means the grid points where variables are stored are typically defined as being at the centre of each control volume.
• The method starts by dividing the flow domain into a number of small control volumes.
• 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.
• The advantage of the finite volume method is that it is easily formulated to allow for unstructured meshes, and this method is used in many computational fluid dynamic packages.

Interpolation schemes in Finite Volume Method:

Interpolation is a process to estimate the values at unknown points by using the points with known value and sample points on either sides of the unknown point. In FVM, interpolation schemes are uThe relation obtained as:

ϕe=68ϕc+38ϕD−18ϕUU

4. Hybrid Interpolation Scheme: It is the blend of CDS and UDS scheme, based on Peclet number

fe=γfCDS+(1−γ)fUDS

5. Total variation Diminishing scheme:

fe=fp+12ψ(r)(fE−fp)

where ψis the flux limiter.

FLUX LIMITERS:

For convective fluid flow, it is observed that the low order schemes are usually stable but quite dissipative in nature around the points of discontinuity or shocks while the higher order schemes are unstable in nature and show oscillations in the vicinity of discontinuity or shocks. Highly accurate and oscillation free schemes are known as High Resolution Schemes.

Flux Limiter functions are used to fine tune the higher order and lower order schemes in such a way that the resulting scheme gives a higher order accuracy in the smooth region of the flow and maintains first order accuracy in the vicinity of shocks and discontinuities. For such a scheme TVD scheme is employed.

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

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

where flow = Low precision flux (1st order accurate)

fhig = High precision flux (Highest order accurate)

ϕ(r) = Flux limiter function where, r = ui−ui−1ui+1−ui