CFD Online URL
[Sponsors]
Home > Wiki > Source term linearization

Source term linearization

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
(Picard's Method)
 
(4 intermediate revisions not shown)
Line 9: Line 9:
where <math>S_C</math> denotes the '''constant''' part of S and <math>S_P</math> denotes the coefficient of <math>\phi_P</math> (not the value of S at P). This allows us to place <math>S_P</math> in the coefficients for <math>\phi_P</math>. <br>
where <math>S_C</math> denotes the '''constant''' part of S and <math>S_P</math> denotes the coefficient of <math>\phi_P</math> (not the value of S at P). This allows us to place <math>S_P</math> in the coefficients for <math>\phi_P</math>. <br>
-
Let <math>\phi_P^*</math> denote the value of <math>\phi_P</math>at the current itertaion. We now write a Taylor series expansion of S about <math>\phi_P^*</math> as <br>
+
Let <math>\phi_P^*</math> denote the value of <math>\phi_P</math>at the current iteration. We now write a Taylor series expansion of S about <math>\phi_P^*</math> as <br>
<math> S = S^* + \left (  \frac {\partial S}{\partial \phi} \right ) ^* \left ( \phi_P - \phi_P^* \right ) </math> <br>
<math> S = S^* + \left (  \frac {\partial S}{\partial \phi} \right ) ^* \left ( \phi_P - \phi_P^* \right ) </math> <br>
therefore<br>
therefore<br>
Line 15: Line 15:
<math> S_P = \left(  \frac {\partial S}{\partial \phi} \right ) ^* </math> <br>
<math> S_P = \left(  \frac {\partial S}{\partial \phi} \right ) ^* </math> <br>
where <math>\left  (  \frac {\partial S}{\partial \phi} \right ) ^* </math> is the gradient of S evaluated at <math>\phi_P^*</math>. <br>
where <math>\left  (  \frac {\partial S}{\partial \phi} \right ) ^* </math> is the gradient of S evaluated at <math>\phi_P^*</math>. <br>
-
<br>
+
 
 +
==Example==
As an illustrative example, consider <math> S = -T^3 + 10 \,</math>. Following Picard's method, we have <br>
As an illustrative example, consider <math> S = -T^3 + 10 \,</math>. Following Picard's method, we have <br>
<math>\left( \frac {\partial S}{\partial \phi} \right ) = -3T^2 </math> <br>
<math>\left( \frac {\partial S}{\partial \phi} \right ) = -3T^2 </math> <br>
<math> S_C = -T_P^{*3} +10  + 3T_P^{*2}T_P^* = 2T_P^{*3} +10 </math> <br>
<math> S_C = -T_P^{*3} +10  + 3T_P^{*2}T_P^* = 2T_P^{*3} +10 </math> <br>
<math> S_P = -3T_P^{*2} </math> <br>
<math> S_P = -3T_P^{*2} </math> <br>
 +
 +
==References==
 +
#{{reference-book|author=Patankar, S.V.|year=1980|title=Numerical Heat Transfer and Fluid Flow|rest=ISBN 0070487405, Hemisphere Publishing Corporation, USA.}}
 +
#{{reference-paper|author=Murthy, Jayathi Y.|year=1998|title=Numerical Methods in Heat, Mass, and Momentum Transfer|rest=Draft Notes, Purdue University ([http://widget.ecn.purdue.edu/%7Ejmurthy/me608/main.pdf/ download])}}
 +
#{{reference-paper|author=[http://webfea-lb.fea.aub.edu.lb/fea/me/CFD/ Darwish, Marwan]|year=2003|title=CFD Course Notes|rest=Notes, American University of Beirut}}

Latest revision as of 10:25, 22 February 2007

Contents

Introduction

In seeking the solution of the general transport equation for a scalar \phi, the main objective is to correctly handle the non-linearities by transforming them into linear form and then iteratively account for the non-linearity. The source term plays a central role in this respect when it is non-linear. For example, in radiation heat transfer, the source term in energy equation is expressed as fourth powers in the temperature.
When the source is constant and independent of the conserved scalar, the finite volume method assumes that the value of S prevails of the control volume and thus can be easily integrated. For a given control volume P, we obtain
\int_{\Omega} S d\Omega = S\Omega \,

Picard's Method

Picard's method is the most popular method used in conjunction with the finite volume method. For a given control volume P, we start by writing the source term as
 S = S_C + S_PT_P \,
where S_C denotes the constant part of S and S_P denotes the coefficient of \phi_P (not the value of S at P). This allows us to place S_P in the coefficients for \phi_P.

Let \phi_P^* denote the value of \phi_Pat the current iteration. We now write a Taylor series expansion of S about \phi_P^* as
 S = S^* + \left (  \frac {\partial S}{\partial \phi} \right ) ^* \left ( \phi_P - \phi_P^* \right )
therefore
 S_C = S^* - \left (  \frac {\partial S}{\partial \phi} \right ) ^* \phi_P^*
 S_P = \left(  \frac {\partial S}{\partial \phi} \right ) ^*
where \left  (  \frac {\partial S}{\partial \phi} \right ) ^* is the gradient of S evaluated at \phi_P^*.

Example

As an illustrative example, consider  S = -T^3 + 10 \,. Following Picard's method, we have
\left( \frac {\partial S}{\partial \phi} \right ) = -3T^2
 S_C = -T_P^{*3} +10  + 3T_P^{*2}T_P^* = 2T_P^{*3} +10
 S_P = -3T_P^{*2}

References

  1. Patankar, S.V. (1980), Numerical Heat Transfer and Fluid Flow, ISBN 0070487405, Hemisphere Publishing Corporation, USA..
  2. Murthy, Jayathi Y. (1998), "Numerical Methods in Heat, Mass, and Momentum Transfer", Draft Notes, Purdue University (download).
  3. Darwish, Marwan (2003), "CFD Course Notes", Notes, American University of Beirut.
My wiki