# Finite volume

### From CFD-Wiki

(→Convection Schemes) |
|||

(48 intermediate revisions not shown) | |||

Line 1: | Line 1: | ||

- | + | The Finite Volume Method (FVM) is one of the most versatile discretization techniques used in CFD. Based on the control volume formulation of analytical fluid dynamics, the first step in the FVM is to divide the domain into a number of control volumes (aka cells, elements) where the variable of interest is located at the centroid of the control volume. The next step is to integrate the differential form of the governing equations (very similar to the control volume approach) over each control volume. Interpolation profiles are then assumed in order to describe the variation of the concerned variable between cell centroids. The resulting equation is called the discretized or discretization equation. In this manner, the discretization equation expresses the conservation principle for the variable inside the control volume. | |

+ | The most compelling feature of the FVM is that the resulting solution satisfies the conservation of quantities such as mass, momentum, energy, and species. This is '''exactly''' satisfied for any control volume as well as for the whole computational domain and '''for any number of control volumes'''. Even a coarse grid solution exhibits exact integral balances.<br> | ||

- | + | FVM is the ideal method for computing discontinuous solutions arising in compressible flows. Any discontinuity must satisfy the [[Rankine-Hugoniot jump condition]] which is a consequence of conservation. Since finite volume methods are conservative they automatically satisfy the jump conditions and hence give physically correct weak solutions. | |

- | + | FVM is also preferred while solving partial differential equations containing discontinuous coefficients. | |

- | == | + | ==Integral form of conservation law== |

- | + | Consider a single partial differential equation in conservation form | |

- | + | :<math> | |

- | :<math> | + | \frac{\partial u}{\partial t} + \nabla \cdot f(u) = 0 |

- | + | ||

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

- | |||

- | + | where <math>u</math> is a conserved quantity and | |

- | + | :<math> | |

- | :<math> | + | f(u) = [ f_1(u), \ldots, f_d(u)] |

- | + | ||

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

- | |||

- | + | is the flux vector. Integrating this over any volume <math>V</math> in <math>R^d</math> we get the '''integral form''' of the conservation law | |

- | + | ||

- | + | :<math> | |

- | + | \frac{\partial}{\partial t} \int_V u dx + \oint_{\partial V} f_i n_i ds = 0 | |

- | + | ||

- | :<math> | + | |

- | + | ||

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

- | |||

- | |||

- | |||

- | |||

- | |||

- | + | where the [[Greens theorem | divergence theorem]] has been used to convert the divergence term to a surface integral over the surface <math>\partial V</math> and <math>(n_1, \ldots, n_d)</math> is the unit outward normal to <math>\partial V</math>. The above equation is the starting point for the finite volume method. | |

- | + | ||

- | + | ||

- | </math> | + | |

- | < | + | |

- | + | ||

- | + | ||

- | + | ||

- | + | ==Finite volume formulation== | |

- | + | ||

- | + | ||

- | + | ||

- | + | The computational domain <math>V_h</math> is divided into non-overlapping cells or finite volumes | |

- | + | ||

- | + | :<math>V_r, r=1,\ldots, N, \quad V_h = \cup_r V_r | |

- | + | ||

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

- | |||

- | + | Usually these cells are polygons (triangles, quadrilaterals) in 2-D and polyhedra (tetrahedron, hexahedron, prims, etc) in 3-D. | |

- | + | Introduce the cell average value over the cell <math>V_r</math> | |

- | + | :<math> | |

+ | u_r = \frac{1}{|V_r|} \int_{V_r} u(x) dx | ||

+ | </math> | ||

- | + | which is the basic unknown quantity in the finite volume method. Let | |

- | + | ||

- | + | ||

- | + | ||

- | + | ||

- | + | ||

- | + | ||

- | + | ||

- | + | :<math> | |

+ | N(r) = \{ \mbox{Set of cells which share a common face with } V_r \} | ||

+ | </math> | ||

- | + | Then the integral conservation law for cell <math>V_r</math> is | |

- | + | :<math> | |

- | + | |V_r| \frac{du_r}{dt} + \sum_{s \in N(r)} \int_{V_r \cap V_s} f_i n_i ds = 0 | |

- | + | ||

- | + | ||

- | + | ||

- | :<math> | + | |

- | + | ||

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

- | |||

- | |||

- | |||

- | + | It now remains to approximate the flux integral. This can be achieved using gaussian quadrature. Taking <math>p</math> gaussian points we get | |

- | + | :<math> | |

- | + | \int_{V_r \cap V_s} f_i n_i ds = \Delta s_{rs} \sum_{m=1}^p \omega_m F^m_{rs} | |

- | + | ||

- | :<math> | + | |

- | + | ||

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

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | |||

- | + | where <math>\omega_m</math> are the gaussian weights and <math>F</math> is an approximation to <math>f_i n_i</math>. | |

- | + | == External links == | |

+ | * [http://www.imtek.uni-freiburg.de/simulation/mathematica/imsReferencePointers/FVM_introDocu.html The Finite Volume Method (FVM) - An introduction] by Oliver Rübenkönig of Albert Ludwigs University of Freiburg, available under the GNU Free Document License|GFDL. | ||

- | |||

- | + | <i> Return to [[Numerical methods | Numerical Methods]] </i> | |

- | + | {{stub}} |

## Latest revision as of 09:15, 3 January 2012

The Finite Volume Method (FVM) is one of the most versatile discretization techniques used in CFD. Based on the control volume formulation of analytical fluid dynamics, the first step in the FVM is to divide the domain into a number of control volumes (aka cells, elements) where the variable of interest is located at the centroid of the control volume. The next step is to integrate the differential form of the governing equations (very similar to the control volume approach) over each control volume. Interpolation profiles are then assumed in order to describe the variation of the concerned variable between cell centroids. The resulting equation is called the discretized or discretization equation. In this manner, the discretization equation expresses the conservation principle for the variable inside the control volume.

The most compelling feature of the FVM is that the resulting solution satisfies the conservation of quantities such as mass, momentum, energy, and species. This is **exactly** satisfied for any control volume as well as for the whole computational domain and **for any number of control volumes**. Even a coarse grid solution exhibits exact integral balances.

FVM is the ideal method for computing discontinuous solutions arising in compressible flows. Any discontinuity must satisfy the Rankine-Hugoniot jump condition which is a consequence of conservation. Since finite volume methods are conservative they automatically satisfy the jump conditions and hence give physically correct weak solutions.

FVM is also preferred while solving partial differential equations containing discontinuous coefficients.

## Integral form of conservation law

Consider a single partial differential equation in conservation form

where is a conserved quantity and

is the flux vector. Integrating this over any volume in we get the **integral form** of the conservation law

where the divergence theorem has been used to convert the divergence term to a surface integral over the surface and is the unit outward normal to . The above equation is the starting point for the finite volume method.

## Finite volume formulation

The computational domain is divided into non-overlapping cells or finite volumes

Usually these cells are polygons (triangles, quadrilaterals) in 2-D and polyhedra (tetrahedron, hexahedron, prims, etc) in 3-D.

Introduce the cell average value over the cell

which is the basic unknown quantity in the finite volume method. Let

Then the integral conservation law for cell is

It now remains to approximate the flux integral. This can be achieved using gaussian quadrature. Taking gaussian points we get

where are the gaussian weights and is an approximation to .

## External links

- The Finite Volume Method (FVM) - An introduction by Oliver Rübenkönig of Albert Ludwigs University of Freiburg, available under the GNU Free Document License|GFDL.

* Return to Numerical Methods *