# Finite Volume Method (FVM)

This section describes the formulation and methodology of finite volume method to solve the governing equations on a computational domain. It also describes the cell centered and face centered approaches used for finite volume formulation.

The finite volume formulation is based on the approximate solution of the integral form of the conservation equations. The problem domain is divided into a set of non overlapping control volumes referred to as finite volumes, where the variable of interest is usually taken at the centroid of the finite volume. The finite volume is also referred to as a cell or element.

The governing equations are integrated over each finite volume and interpolation profiles are assumed in order to describe the variation of the concerned variable of interest between the cell centroids. The resulting discretization equation expresses the conservation principle for the variable inside the finite volume.

The governing differential equations discussed in this manual each have a dependent variable that obeys a generalized conservation principle. If the dependent variable (scalar or vector) is denoted by $\phi $ , the generic differential equation is:

- $\u03f5$ is the diffusion coefficient.
- The transient term $\frac{\partial \rho \phi}{\partial t}$ accounts for the rate of change of $\phi $ inside the control volume.
- The convection term $\nabla \cdot \left(\rho \overrightarrow{u}\phi \right)$ accounts for the transport of $\phi $ due to the velocity field $\overrightarrow{u}$ .
- The diffusion term $\nabla \cdot \left(\u03f5\nabla \phi \right)$ accounts for the transport of $\phi $ due to its gradients.
- The source term ${S}_{\phi}$ accounts for any sources or sinks that affect the quantity $\phi $ .

When the above equation is integrated over a three dimensional control volume it yields

The divergence terms (convection and diffusion) can be converted into surface integrals using Gauss divergence theorem and the above equation can be rewritten as

- $\widehat{n}$ is the unit vector normal to the surface $dA$
- $S$ is the boundary of the control volume
- $\widehat{n}\cdot \left(\rho \overrightarrow{u}\phi \right)$ is the convective flux of $\phi $ across the boundary $dS$
- $\widehat{n}\cdot \left(\u03f5\nabla \phi \right)$ is the diffusive flux of $\phi $ across the boundary $dS$

The integral conservation in the scalar transport equation applies to each control volume as well as the complete solution domain thus satisfying the global conservation of quantities such as mass, momentum and energy. These quantities can be evaluated as fluxes at the surfaces of each control volume.

In order to obtain an algebraic (discretised) equation for each control volume the surface and volume integrals are approximated using quadrature formulae in terms of function values at the storage location. This may require values of variable at points other than the computational nodes of the control volume. Values at these locations are obtained using interpolation schemes.

- Cell centered approach: The domain is defined by a suitable grid of finite volumes and the computational nodes are assigned at the centroid of the control volumes.
- Face/Node centered approach: The domain is first defined by a set of nodes and control volumes are constructed around these nodes as cell centers such that the faces of the control volume lie between these nodes.

The cell centered approach and node centered approach have nearly the same accuracy and efficiency for most of the cases which use a structured grid.

A cell centered approach is employed with the points P, W, E, N, S representing the cell centers and the notations n, e, s, w representing the faces of cell P. The notation N (n), E(e), S(s), W(w) represent the north, east, south and west directions, respectively.

The velocity $u$ is stored at the nodes N, E, S, W and it is represented as ${u}_{N}$ , ${u}_{E}$ , ${u}_{S}$ , ${u}_{W}$ , respectively. The finite volume approximation starts by integrating the x-momentum equation over the control volume P.

- Evaluation of Source Terms
- The source terms can be approximated by evaluating them at cell centers and multiplying them by the volume of the cell.
- Evaluation of Diffusive Fluxes
- The diffusive fluxes can be approximated as
- Evaluation of Convective Fluxes
- The convective fluxes can be approximated as