Finite volume method for three-dimensional diffusion problem
Finite volume method (FVM) is a numerical method. FVM in computational fluid dynamics is used to solve the partial differential equation which arises from the physical conservation law by using discretisation. Convection is always followed by diffusion and hence where convection is considered we have to consider combine effect of convection and diffusion. But in places where fluid flow plays a non-considerable role we can neglect the convective effect of the flow. In this case we have to consider more simplistic case of only diffusion. The general equation for steady convection-diffusion can be easily derived from the general transport equation for property  by deleting transient.

General transport equation is defined as:

 …………………………………………….1

Where,

 is a conservative form of all fluid flow,	

 is density,

 is a net rate of flow of  out of fluid element represents convective term,

 is a transient term,

 is a rate of change   of due to diffusion,

 is a rate of increase of  due  to source.

Due to steady state condition transient term becomes zero and due to absence of convection convective term becomes zero, therefore steady state three- dimensional convection and diffusion equation becomes:

………………………………………………………….2

Therefore,

…………………………………………………………………….3

Flow should also satisfy continuity equation therefore,

 ………………………………………………………………………………………………………4

To solve the problem we will follow following general steps

Grid formation:

1. Divide the domain into discrete control volume.

2. Place the nodal point between end points defining the physical boundaries. Boundaries/ faces of the control volume are created midway between adjacent nodes.

3. Set up the control volume near the edge of domain such that physical as well as control volume boundaries will coincide with each other.

4. Considering a general nodal point P accompanied by six neighboring nodal point ‘E’ (which represent east), ‘W’ (which represent west), ‘N’ (which represent north), ‘S’ (which represent south), ‘T’ (which represent Top), ‘B’ (which represent bottom). In the considered control volume east side face is referred by ‘e’, west side face is referred by ‘w’, north side face is referred by ‘n’, south side face is referred by ‘s’, top side face is referred by ‘t’, bottom side face is referred by ‘b’.

5. Now the distance between nodes W and P, between nodes P and E, between nodes P and N, between nodes S and P, between nodes P and T, between nodes B and P are denoted as respectively.

Discretisation:

On integration of equation 3 in one dimension over the general control volume gives:

 [

Now using central differencing method we can rewrite above equation as

[

This can be rearranged to give the discretised equation for interior nodes:

Where

Solution of equation:

1.	For solving the one- dimensional convection- diffusion problem we have to express equation (8) at all the grid nodes.

2.	Now obtained set of algebraic equations is then solved to obtain the distribution of the transported property .

See also
Finite volume method
Computational fluid dynamics
Finite volume method for one-dimensional steady state diffusion
Convection
Control volume
Central differencing scheme

External links
http://mathworld.wolfram.com/FiniteVolumeMethod.html
 The finite volume method by R. Eymard, T Gallouët and R. Herbin, update of the article published in Handbook of Numerical Analysis, 2000
https://web.archive.org/web/20140210101323/http://s6.aeromech.usyd.edu.au/aero/cvanalysis/integral_approach.pdf
http://www.phy.davidson.edu/fachome/dmb/py200/centraldiff.htm
http://opencourses.emu.edu.tr/course/view.php?id=27&lang=en

References

Computational fluid dynamics
Mathematical problems