# Diffusion term

### From CFD-Wiki

Line 45: | Line 45: | ||

</math> <br> | </math> <br> | ||

where <math> \bar \nabla \phi _f </math> and <math> \Gamma _f </math> are suitable face averages. <br> | where <math> \bar \nabla \phi _f </math> and <math> \Gamma _f </math> are suitable face averages. <br> | ||

+ | |||

+ | === Orthogonal Correction Approaches === | ||

+ | In non-orthogonal grids, the gradient direction that will yield an expression involving the values at the neighboring control volumes will have to be along the line joining the centroids of the two control volumes. If this direction has a unit vector denoted by <math>\vec e</math> then, by definition <br> | ||

+ | <math> | ||

+ | \vec e {\rm{ = }} \frac{{{\rm{\overrightarrow{d_{PN}}}}}} {\left| {\overrightarrow{d_{PN}}} \right|} | ||

+ | </math><br> | ||

+ | then the gradient in the direction of <math>\vec e</math> can be written as <br> | ||

+ | <math> \nabla \phi _f \cdot \vec e = \frac {\partial \phi_f} {\partial e} = \frac { \phi_N - \phi_P} {\left| {\overrightarrow{d_{PN}}} \right|} </math> <br> | ||

+ | |||

+ | If the surface vector <math>\vec {S_f}</math> is written as the summation of two vectors <math>\vec {E}</math> and <math>\vec {T}</math> <br> | ||

+ | <math>\vec {S_f} = \vec {E} + \vec {T}</math> <br> | ||

+ | where <math>\vec {E}</math> is in the direction joining the centroids of the two control volumes, we will then be able to express the diffusive flux in terms of the neighboring control volumes plus an additional correction. This is done as follows <br> | ||

+ | |||

+ | <math> \nabla \phi_f \cdot \vec {S_f} = \nabla \phi_f \cdot \vec {E} + \nabla \phi_f \cdot \vec {T} </math> <br> | ||

+ | |||

+ | <math> \nabla \phi_f \cdot \vec {S_f} = E \nabla \phi_f \cdot \vec {e} + \nabla \phi_f \cdot \vec {T} </math> .... (where E is the magnitude of <math> \vec E </math><br> | ||

+ | At the outset, one obtains <br> | ||

+ | <math> \nabla \phi_f \cdot \vec {S_f} = E \frac { \phi_N - \phi_P} {\left| {\overrightarrow{d_{PN}}} \right|} + \nabla \phi_f \cdot \vec {T} </math> <br> | ||

+ | <br> | ||

+ | The first term in the above equation can be thought of as the orthogonal contribution to the diffusive flux, while the second term represents the non-orthogonal effects. At this point, the vector <math>\vec {T}</math> has not been defined yet. There are three main methods to define this vector. | ||

+ | |||

+ | ==== Minimum Correction ==== | ||

+ | ==== Orthogonal Correction ==== | ||

+ | ==== Over Relaxed Correction ==== | ||

== References == | == References == |

## Revision as of 05:08, 5 December 2005

## Contents |

## Discretisation of the Diffusion Term

### Description

For a general control volume (orthogonal, non-orthogonal), the discretization of the diffusion term can be written in the following form

where

- S denotes the surface area of the control volume
- denotes the area of a face for the control volume

As usual, the subscript f refers to a given face. The figure below describes the terminology used in the framework of a general non-orthogonal control volume

**A general non-orthogonal control volume**

Note: The approaches those are discussed here are applicable to non-orthoganal meshes as well as orthogonal meshes.

A control volume in mesh is made up of set of faces enclosing it. Where represents the magnitude of area of the face. And **n** represents the normal unit vector of the face under consideration.

If and are position vector of centroids of cells P and N respectively. Then, we define

We wish to approaximate the diffusive flux at the face.

### Approach 1

A first approach is to use a simple expression for estimating the gradient of a scalar normal to the face.

where is a suitable face average.

This approach is not very good when the non-orthogonality of the faces increases. If this is the case, it is advisable to use one of the following approaches.

### Approach 2

We define the vector

giving us the expression:

where and are suitable face averages.

### Orthogonal Correction Approaches

In non-orthogonal grids, the gradient direction that will yield an expression involving the values at the neighboring control volumes will have to be along the line joining the centroids of the two control volumes. If this direction has a unit vector denoted by then, by definition

then the gradient in the direction of can be written as

If the surface vector is written as the summation of two vectors and

where is in the direction joining the centroids of the two control volumes, we will then be able to express the diffusive flux in terms of the neighboring control volumes plus an additional correction. This is done as follows

.... (where E is the magnitude of

At the outset, one obtains

The first term in the above equation can be thought of as the orthogonal contribution to the diffusive flux, while the second term represents the non-orthogonal effects. At this point, the vector has not been defined yet. There are three main methods to define this vector.

#### Minimum Correction

#### Orthogonal Correction

#### Over Relaxed Correction

## References

**Ferziger, J.H. and Peric, M. (2001)**,*Computational Methods for Fluid Dynamics*, ISBN 3540420746, 3rd Rev. Ed., Springer-Verlag, Berlin..**Hrvoje, Jasak (1996)**, "Error Analysis and Estimation for the Finite Volume Method with Applications to Fluid Flows", PhD Thesis, Imperial College, University of London (download).

* Return to Numerical Methods *