View source for Finite volume method for two dimensional diffusion problem

{{Context|date=October 2024}}
The methods used for solving two dimensional [[Diffusion]] problems are similar to those used for one dimensional problems. The general equation for steady diffusion can be easily derived from the general transport equation for property ''Φ'' by deleting transient and convective terms<ref>{{cite web|author= |url=http://www.efunda.com/formulae/fluids/navier_stokes.cfm |title=Navier-Stokes Equations in Fluid Mechanics |publisher=Efunda.com |date= |accessdate=2013-10-29}}</ref>
<math>\frac{\partial{}}{\partial{}x}\left(\Gamma{}\frac{\partial{}\phi{}}{\partial{}x}\right)+\frac{\partial{}}{\partial{}y}\left(\Gamma{}\frac{\partial{}\phi{}}{\partial{}y}\right)+S=0</math>
where,<br />
<math>\Gamma</math> is the Diffusion coefficient<ref>{{cite web|url=http://www.life.illinois.edu/crofts/bioph354/diffusion1.html |title=Diffusion – useful equations |publisher=Life.illinois.edu |date= |accessdate=2013-10-29}}</ref> and <math>S</math> is the Source term.<ref>{{cite web|url=http://www.physics.drexel.edu/~valliere/PHYS405/DiffusionEquation/diffusion.html |title=SSCP: Programming Strategies |publisher=Physics.drexel.edu |date= |accessdate=2013-10-29}}</ref>
A portion of the two dimensional [[Polygon mesh|grid]] used for [[Discretization]] is shown below:
[[File:Cfd graph.jpg|thumb|Graph of 2 dimensional plot]]
In addition to the east (E) and west (W) neighbors, a general grid node P, now also has north (N) and south (S) neighbors. The same notation is used
here for all faces and cell dimensions as in one dimensional analysis. When the above equation is formally integrated over the [[Control volume]], we obtain
<math>\int_{\Delta{v}}\frac{\partial{}}{\partial{}x}\left(\Gamma{}\frac{\partial{}\phi{}}{\partial{}x}\right)dxdy+\int_{\Delta{v}}\frac{\partial{}}{\partial{}y}\left(\Gamma{}\frac{\partial{}\phi{}}{\partial{}y}\right)dxdy+\int_{\Delta{v}}S_{\phi{}}dV=0</math>
Using the divergence theorem, the equation can be rewritten as :
<math>
\left[{\Gamma{}}_eA_e\left(\frac{\partial{}\phi{}}{\partial{}x}\right)_{e}-{\Gamma{}}_wA_w\left(\frac{\partial{}\phi{}}{\partial{}x}\right)_{w}\right]+\left[{\Gamma{}}_nA_n\left(\frac{\partial{}\phi{}}{\partial{}y}\right)_{n}-{\Gamma{}}_sA_s\left(\frac{\partial{}\phi{}}{\partial{}y}\right)_{s}\right]+\bar{S}\Delta{}V=0
</math>
This equation represents the balance of generation of the property φ in a [[Control volume]] and the [[flux]]es through its cell faces. The derivatives can by represented as follows by using [[Taylor series]] approximation:
<math>
{{\Gamma{}}_wA_w\left(\frac{\partial{}\phi{}}{\partial{x}}\right)}_w=
{\Gamma{}}_wA_w\frac{({\phi{}}_p-{\phi{}}_w)}{{\delta{}x}_{WP}}
</math>
Flux across the east face =
<math>
{{\Gamma{}}_eA_e\left(\frac{\partial{}\phi{}}{\partial{x}}\right)}_e=
{\Gamma{}}_eA_e\frac{({\phi{}}_e-{\phi{}}_p)}{{\delta{}x}_{PE}}
</math>
Flux across the south face =
<math>
{{\Gamma{}}_sA_s\left(\frac{\partial{}\phi{}}{\partial{y}}\right)}_s=
{\Gamma{}}_sA_s\frac{({\phi{}}_p-{\phi{}}_s)}{{\delta{}y}_{SP}}
</math>
Flux across the north face =
<math>
{{\Gamma{}}_nA_n\left(\frac{\partial{}\phi{}}{\partial{y}}\right)}_n=
{\Gamma{}}_nA_n\frac{({\phi{}}_n-{\phi{}}_p)}{{\delta{}y}_{PN}}
</math>
Substituting these expressions in equation (2) we obtain
<math>
{\Gamma{}}_eA_e\frac{({\phi{}}_e-{\phi{}}_p)}{{\delta{}x}_{PE}} - {\Gamma{}}_wA_w\frac{({\phi{}}_p-{\phi{}}_w)}{{\delta{}x}_{WP}} + {\Gamma{}}_nA_n\frac{({\phi{}}_n-{\phi{}}_p)}{{\delta{}y}_{PN}} - {\Gamma{}}_sA_s\frac{({\phi{}}_p-{\phi{}}_s)}{{\delta{}y}_{SP}}+ \bar{S}\Delta{}V=0
</math>
When the source term is represented in linearized form <math> \bar{S}\Delta{}V = S_u + S_p*S_{\varphi}</math> ,
this equation can be rearranged as,
<math>\left[
\frac{{\Gamma{}}_eA_e}{{\delta{}x}_{PE}}
+ \frac{{\Gamma{}}_wA_w}{{\delta{}x}_{WP}}
+ \frac{{\Gamma{}}_nA_n}{{\delta{}y}_{PN}}
+ \frac{{\Gamma{}}_sA_s}{{\delta{}y}_{SP}}
- S_p\right] \phi{}_P </math> =
<math>
\frac{{\Gamma{}}_eA_e}{{\delta{}x}_{PE}}\phi{}_E
+ \frac{{\Gamma{}}_wA_w}{{\delta{}x}_{WP}}\phi{}_W
+ \frac{{\Gamma{}}_nA_n}{{\delta{}x}_{PN}}\phi{}_N
+ \frac{{\Gamma{}}_sA_s}{{\delta{}x}_{SP}}\phi{}_S
+ S_u  </math>
This equation can now be expressed in a general [[Discretization|discretized]] equation form for internal nodes, i.e.,
<math>a_P \phi{}_P= a_W \phi{}_W + a_E \phi{}_E + a_S \phi{}_S + a_N \phi{}_N + S_u</math>
Where,
{| class="wikitable"
|-
! <math>a_W</math> !! <math>a_E</math> !! <math>a_S</math> !! <math>a_N</math> !! <math>a_P</math>
|-
|<math>\frac{{\Gamma{}}_wA_w}{{\delta{}x}_{WP}}
</math>  || <math>\frac{{\Gamma{}}_eA_e}{{\delta{}x}_{PE}}
</math> || <math>\frac{{\Gamma{}}_sA_s}{{\delta{}x}_{SP}}
</math> || <math>\frac{{\Gamma{}}_nA_n}{{\delta{}x}_{PN}}
</math> || <math>a_W+a_E+a_S+a_N-S_p</math>
|}
The face areas in y two dimensional case are :
<math> A_{w} = A_{e} = \Delta{}y </math>
and
<math> A_{n} = A_{s} = \Delta{}x</math>.
We obtain the distribution of the property <math>\varphi{}</math> i.e. a given two dimensional situation by writing [[Discretization|discretized]] equations of the form of equation (3) at each grid node of the subdivided domain. At the boundaries where the temperature or fluxes are known the discretized equation are modified to incorporate the [[Boundary value problem|boundary conditions]]. The boundary side coefficient is set to zero (cutting the link with the boundary) and the flux crossing this boundary is introduced as a source which is appended to any existing <math>s_{u} </math> and <math> S_{p}</math> terms. Subsequently the resulting set of equations is solved to obtain the two dimensional distribution of the property <math>\varphi{}</math>
==References==
* Patankar, Suhas V. (1980), Numerical Heat Transfer and Fluid Flow, Hemisphere.
* Hirsch, C. (1990), Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows, Wiley.
* Laney, Culbert B.(1998), Computational Gas Dynamics, Cambridge University Press.
* LeVeque, Randall(1990), Numerical Methods for Conservation Laws, ETH Lectures in Mathematics Series, Birkhauser-Verlag.
* Tannehill, John C., et al., (1997), Computational Fluid mechanics and Heat Transfer, 2nd Ed., Taylor and Francis.
* Wesseling, Pieter(2001), Principles of Computational Fluid Dynamics, Springer-Verlag.
* Carslaw, H. S. and Jager, J. C. (1959). Conduction of Heat in Solids. Oxford: Clarendon Press
* Crank, J. (1956). The Mathematics of Diffusion. Oxford: Clarendon Press
* Thambynayagam, R. K. M (2011). The Diffusion Handbook: Applied Solutions for Engineers: McGraw-Hill
{{reflist}}
== External links ==
* http://opencourses.emu.edu.tr/course/view.php?id=27&lang=en
* https://web.archive.org/web/20120303230200/http://nptel.iitm.ac.in/courses/112105045/
* http://ingforum.haninge.kth.se/armin/CFD/dirCFD.htm {{Webarchive|url=https://web.archive.org/web/20120713003851/http://ingforum.haninge.kth.se/armin/CFD/dirCFD.htm |date=2012-07-13 }}
* [http://www.math.uci.edu/~chenlong/226/FVM.pdf Finite volume method, Cheng Long]
* [http://www.scholarpedia.org/article/Finite_volume_method Finite volume method, Robert Eymard et al. (2010), Scholarpedia,5(6):9835]
== See also ==
*[[Computational fluid dynamics]]
* [[Finite difference]]
* [[Heat equation]]
* [[Fokker–Planck equation]]
* [[Fick's laws of diffusion]]
* [[Maxwell–Stefan equation]]
*[[Diffusion equation]]
* [[Convection–diffusion equation]]
[[Category:Computational fluid dynamics]]