Calculation on non-orthogonal curvelinear structured grids, finite-volume method

From CFD-Wiki

(Redirected from Calculation on non-orthogonal curvelinear grids)

2D case

For calculations in complex geometries boundary-fitted non-orthogonal curvlinear grids is usually used.

General transport equation is transformed from the physical domain $(x,y)$ into the computational domain $\left( \xi , \eta \right)$ as the following equation

$\frac{\partial}{\partial \xi} \left( \rho U \phi \right) + \frac{\partial }{ \partial \eta } \left( \rho V \phi \right) = \frac{\partial}{\partial \xi} \left[ \frac{\Gamma }{J} \left( \alpha \frac{\partial \phi}{\partial \xi} - \beta \frac{\partial \phi}{ \partial \eta} \right) \right] + \frac{\partial}{\partial \eta} \left[ \frac{\Gamma}{J} \left( \gamma \frac{\partial \phi}{\partial \eta} - \beta \frac{\partial \phi}{ \partial \xi} \right) \right] + J S^{\phi}$

(2)

where

$U = \overline{u} \frac{\partial y}{\partial \eta} - \overline{v} \frac{\partial x}{\partial \eta}$

(3)

$V = \overline{v} \frac{\partial x}{ \partial \xi} - \overline{u} \frac{\partial y}{ \partial \xi}$

(4)

$\alpha = \left( \frac{\partial x}{\partial \eta } \right)^2 + \left( \frac{\partial y}{\partial \eta } \right)^2$

(5)

$\gamma = \left( \frac{\partial x}{ \partial \xi } \right)^2 + \left( \frac{\partial y}{ \partial \xi } \right)^2$

(6)

$\beta = \frac{\partial x}{ \partial \xi} \frac{\partial x}{ \partial \eta} + \frac{\partial y}{ \partial \xi} \frac{\partial y}{ \partial \eta}$

(7)

$J = \frac{\partial x}{ \partial \xi} \frac{\partial y}{ \partial \eta} - \frac{\partial y}{ \partial \xi} \frac{\partial x}{ \partial \eta}$

(8)

Using the finite volume method the trnsformed equations can be integrated as follows:

$\left[ \left( \rho U \Delta \eta \right) \right]^{e}_{w} + \left[ \left( \rho V \Delta \xi \right) \right]^{n}_{s} = \left[ \frac{\Gamma \Delta \eta}{J} \left( \alpha \frac{\partial \phi}{\partial \xi} - \beta \frac{\partial \phi}{ \partial \eta} \right) \right]^{e}_{w} + \left[ \frac{\Gamma \Delta \xi}{J} \left( \gamma \frac{\partial \phi}{\partial \eta } - \beta \frac{\partial \phi}{\partial \xi} \right) \right]^{n}_{s} + \left( J \Delta \xi \Delta \eta \right) \overline{S}^{\phi}_{P}$

(9)

The convection terms are approximated as described in section http://www.cfd-online.com/Wiki/Discretization_of_the_convection_term .

Diffusion terms are approximated by the second-oder central differencing scheme.

The standard form of the finite volume equation can be obtained as

$A^{\phi}_{P} \phi_{P} = A^{\phi}_{E} \phi_{E} + A^{\phi}_{W} \phi_{W} + A^{\phi}_{N} \phi_{N} + A^{\phi}_{S} \phi_{S} + b^{\phi}$

(10)

where

$A^{\phi}_{E} = \left(\frac{\Gamma}{J} \alpha \frac{\Delta \eta}{\Delta \xi} \right)_{e} + max \left[ 0, - \left( \rho U \Delta \eta \right)_{e} \right]$

(11)

$A^{\phi}_{W} = \left(\frac{\Gamma}{J} \alpha \frac{\Delta \eta}{\Delta \xi} \right)_{w} + max \left[ 0, \left( \rho U \Delta \eta \right)_{w} \right]$

(12)

$A^{\phi}_{N} = \left( \frac{ \Gamma }{J} \gamma \frac{\Delta \xi}{\Delta \eta} \right)_{n} + max \left[ 0, - \left( \rho V \Delta \xi \right)_{n} \right]$

(13)

$A^{\phi}_{S} = \left( \frac{ \Gamma }{J} \gamma \frac{\Delta \xi}{\Delta \eta} \right)_{s} + max \left[ 0, \left( \rho V \Delta \xi \right)_{s} \right]$

(14)

$b^{\phi} = \left( J \Delta \xi \Delta \eta \right) \overline{S}^{\phi}_{P} - \left[ \frac{\Gamma \Delta \eta}{J} \left( \beta \frac{\partial \phi}{\partial \eta} \right) \right]^{e}_{w} - \left[ \frac{\Gamma \Delta \xi}{J} \left( \beta \frac{\partial \phi}{\partial \xi} \right) \right]^{n}_{s}$

(15)

Calculation on non-orthogonal curvelinear structured grids, finite-volume method

From CFD-Wiki

2D case

Views

My wiki

wiki navigation

wiki search

wiki toolbox