From charlesreid1

No edit summary
No edit summary
Line 1: Line 1:
Reynolds Transport Theorem is always a safe starting point
Reynolds Transport Theorem is a good starting point for deriving governing equations of quantities in a fluid.
 
In a way, it is the most fundamental starting point.


Fundamental balance equation:
Fundamental balance equation:
Line 11: Line 9:
Reynolds Transport Theorem is a formal way of performing this balance over a fluid control volume, or fluid material volume, of arbitrary shape, moving at an arbitrary velocity.
Reynolds Transport Theorem is a formal way of performing this balance over a fluid control volume, or fluid material volume, of arbitrary shape, moving at an arbitrary velocity.


A derivation of Reynolds Transport Theorem is given here: [[Reynolds Transport Theorem Derivation]]
Definitions and derivations of some important terms and identities used to construct Reynolds Transport Theorem are given here: [[Reynolds Transport Theorem Derivation]]
 
= Definition =
 
For any extensive fluid property <math>B</math>, there is a corresponding intensive property <math>b = \frac{\partial B}{\partial m}</math>.


Given an extensive property <math>F</math> and a corresponding intensive property <math>\mathcal{F} = \frac{\partial F}{\partial t}</math>, <math>F</math> can be defined in terms of <math>\mathcal{F}</math>:
A material volume corresponding to a mass of fluid with the property <math>B</math> can be written as <math>V_{B}^{M}(t)</math>


<math>
<math>
F = \iiint_{V_{\mathcal{F}}(t)} \rho \mathcal{F} dV
B(t) = \iiint_{V_{B}^{M}} \rho b(\boldsymbol{x}, t) dV
</math>
</math>


where <math>\rho</math> is the fluid density and <math>V_{\mathcal{F}}</math> is the material volume that corresponds to the property <math>F</math> (which moves at a velocity <math>\boldsymbol{v}_{\mathcal{F}}</math>.
where <math>\rho</math> is the fluid density.


The partial derivative with respect to volume for the intensive property <math>\mathcal{F}</math> is,
{{ambox
| type=notice
| text =  '''Units analysis:'''
 
The presence of <math>\rho</math> in the above integral definition of <math>B(t)</math> can be verified by analyzing the units of some quantities.
 
The quantity <math>b = \frac{\partial B}{\partial m}</math> has units <math>[b] = \frac{ [B] }{ [m] }</math>
 
The quantity <math>\iiint_{V} b dV</math> has units <math>\frac{ [B] }{ [m] } [vol]</math>
 
The quantity <math>\iiint_{V} \rho b dV</math> has units <math>\frac{ [B] }{ [m] } \frac{ [m] }{ [vol] } [vol] = [B] </math>
}}
 
The quantity of interest for the material volume is material derivative of the extensive property, <math>\frac{dB}{dt}</math> - that is, the rate of change of <math>B</math> in the material volume <math>V_{B}^{M}</math> (which is a volume whose boundaries move such that the flux of <math>B</math> across the boundaries is 0, see [[Reynolds Transport Theorem Derivation]] page for details).
 
This quantity can be written as:


<math>
<math>
\frac{ \partial \mathcal{F} }{ \partial V } = 0
\frac{dB}{dt} = \frac{d}{dt} \iiint_{V_{B}^{M}(t)} \rho b(\boldsymbol{x}, t) dV
</math>
</math>


and for the extensive property <math>F</math>,
Because the integral is with time, the integral <math>\iiint_{V_{B}^{M}(t)} and the derivative <math>\frac{d}{dt}</math> do not commute.
 
However, if the integration were over a fixed volume that did not change in time (for example, <math>V_{B,0}</math>, the integral and the derivative ''would'' commute.
 
The Jacobian <math>J = \frac{ d V_{B}^{M} }{ d V_{0,B} }</math> allows us to do this.
 
The volume integral can be converted as follows:


<math>
<math>
\frac{ \partial F }{ \partial V }  
\begin{array}{rcl}
= \frac{ \partial F }{ \partial m } \frac{ \partial m }{ \partial V }  
\frac{ dB }{ dt } &=& \frac{d}{dt} \iiint_{V_{B}^{M}(t)} \rho b(\boldsymbol{x},t) dV \\
= \frac{1}{\rho} \frac{ \partial F}{\partial m}
&=& \frac{d}{dt} \iiint_{V_{B,0}} \rho b(\boldsymbol{x},t) J dV
\end{array}
</math>
</math>


Then Reynolds Transport Theorem is given by:
Next, the temporal derivative can be taken inside the integral (since they now commute):


<math>
<math>
\frac{d}{dt} \iiint_{V_{\mathcal{F}}(t)}  
\begin{array}{rcl}
= \iiint_{V_{\mathcal{F}}(t)} \frac{ \partial \rho \mathcal{F} }{\partial t} dV  
\frac{d}{dt} \iiint_{V_{B,0}} \rho b(\boldsymbol{x},t) J dV
+ \iint_{S_{\mathcal{F}}(t)} \rho \mathcal{F} \boldsymbol{v}_{\mathcal{F}} \cdot \boldsymbol{n} dS
&=& \iiint_{V_{B,0}} \left[ \frac{ d \rho b }{ dt } J + \rho b \frac{ dJ }{ dt } \right] dV \\
&=& \iiint_{V_{B,0}} \left[ \frac{ d \rho b }{ dt } J + \rho b \frac{ 1 }{ J } \frac{dJ}{dt} \right] J
\end{array}
</math>
 
Next, the Euler expansion formula can be used to simplify the second term in the integral to yield:
 
<math>
\frac{dB}{dt} = \iiint_{V_{B,0}} \left[ \frac{ d \rho b }{ dt } + \rho b \left( \nabla \cdot \boldsymbol{v}_{B} \right] J dV
</math>
 
and this can now be transformed back to the integral over the material volume <math>V_{B}^{M}</math> to yield Reynolds' Transport Theorem:
 
{|style="border:1px solid black;"
|<math>
\frac{dB}{dt} = \frac{d}{dt} \iiint_{V_{B}^{M}(t)} \rho b dV = \iiint_{V_{B}^{M}(t)} \left[ \frac{d \rho b}{ dt } + \rho b \left( \nabla \cdot \boldsymbol{v}_{B} \right) \right] dV
</math>
|}
 
Now, the definition of the substantial derivative can be used to put this into a different (more familiar) form:
 
{|style="border:1px black solid;"
|<math>
\frac{dB}{dt} = \iiint_{V_{B}^{M}(t)} \left[ \frac{ \partial \rho b }{\partial t} + \nabla \cdot \left( \rho b \boldsymbol{v}_{B} \right) \right] dV
</math>
|}
 
Green's Theorem can now be applied to this last divergence term, so that
 
{|style="border:1px black solid;"
|<math>
\frac{dB}{dt} = \iiint_{V_{B}^{M}(t)} \frac{ \partial \rho b }{\partial t} dV + \iint_{S_{B}^{M}(t)} \rho b \boldsymbol{v}_{B} \cdot \boldsymbol{n} dS
</math>
</math>
|}


where <math>\boldsymbol{v}_{\mathcal{F}}</math> is the velocity vector of the material volume containing <math>F</math>.
where <math>S_{B}^{M}</math> is the surface of the material volume <math>V_{B}^{M}</math>, <math>\boldsymbol{n}</math> is the vector normal to this surface, and <math>dS</math> is a differential area on the surface.

Revision as of 00:45, 4 November 2010

Reynolds Transport Theorem is a good starting point for deriving governing equations of quantities in a fluid.

Fundamental balance equation:

$ \mbox{in} - \mbox{out} + \mbox{generation} - \mbox{consumption} = \mbox{accumulation} $

Reynolds Transport Theorem is a formal way of performing this balance over a fluid control volume, or fluid material volume, of arbitrary shape, moving at an arbitrary velocity.

Definitions and derivations of some important terms and identities used to construct Reynolds Transport Theorem are given here: Reynolds Transport Theorem Derivation

Definition

For any extensive fluid property $ B $, there is a corresponding intensive property $ b = \frac{\partial B}{\partial m} $.

A material volume corresponding to a mass of fluid with the property $ B $ can be written as $ V_{B}^{M}(t) $

$ B(t) = \iiint_{V_{B}^{M}} \rho b(\boldsymbol{x}, t) dV $

where $ \rho $ is the fluid density.

The quantity of interest for the material volume is material derivative of the extensive property, $ \frac{dB}{dt} $ - that is, the rate of change of $ B $ in the material volume $ V_{B}^{M} $ (which is a volume whose boundaries move such that the flux of $ B $ across the boundaries is 0, see Reynolds Transport Theorem Derivation page for details).

This quantity can be written as:

$ \frac{dB}{dt} = \frac{d}{dt} \iiint_{V_{B}^{M}(t)} \rho b(\boldsymbol{x}, t) dV $

Because the integral is with time, the integral $ \iiint_{V_{B}^{M}(t)} and the derivative <math>\frac{d}{dt} $ do not commute.

However, if the integration were over a fixed volume that did not change in time (for example, $ V_{B,0} $, the integral and the derivative would commute.

The Jacobian $ J = \frac{ d V_{B}^{M} }{ d V_{0,B} } $ allows us to do this.

The volume integral can be converted as follows:

$ \begin{array}{rcl} \frac{ dB }{ dt } &=& \frac{d}{dt} \iiint_{V_{B}^{M}(t)} \rho b(\boldsymbol{x},t) dV \\ &=& \frac{d}{dt} \iiint_{V_{B,0}} \rho b(\boldsymbol{x},t) J dV \end{array} $

Next, the temporal derivative can be taken inside the integral (since they now commute):

$ \begin{array}{rcl} \frac{d}{dt} \iiint_{V_{B,0}} \rho b(\boldsymbol{x},t) J dV &=& \iiint_{V_{B,0}} \left[ \frac{ d \rho b }{ dt } J + \rho b \frac{ dJ }{ dt } \right] dV \\ &=& \iiint_{V_{B,0}} \left[ \frac{ d \rho b }{ dt } J + \rho b \frac{ 1 }{ J } \frac{dJ}{dt} \right] J \end{array} $

Next, the Euler expansion formula can be used to simplify the second term in the integral to yield:

$ \frac{dB}{dt} = \iiint_{V_{B,0}} \left[ \frac{ d \rho b }{ dt } + \rho b \left( \nabla \cdot \boldsymbol{v}_{B} \right] J dV $

and this can now be transformed back to the integral over the material volume $ V_{B}^{M} $ to yield Reynolds' Transport Theorem:

$ \frac{dB}{dt} = \frac{d}{dt} \iiint_{V_{B}^{M}(t)} \rho b dV = \iiint_{V_{B}^{M}(t)} \left[ \frac{d \rho b}{ dt } + \rho b \left( \nabla \cdot \boldsymbol{v}_{B} \right) \right] dV $

Now, the definition of the substantial derivative can be used to put this into a different (more familiar) form:

$ \frac{dB}{dt} = \iiint_{V_{B}^{M}(t)} \left[ \frac{ \partial \rho b }{\partial t} + \nabla \cdot \left( \rho b \boldsymbol{v}_{B} \right) \right] dV $

Green's Theorem can now be applied to this last divergence term, so that

$ \frac{dB}{dt} = \iiint_{V_{B}^{M}(t)} \frac{ \partial \rho b }{\partial t} dV + \iint_{S_{B}^{M}(t)} \rho b \boldsymbol{v}_{B} \cdot \boldsymbol{n} dS $

where $ S_{B}^{M} $ is the surface of the material volume $ V_{B}^{M} $, $ \boldsymbol{n} $ is the vector normal to this surface, and $ dS $ is a differential area on the surface.