Fluid Mechanics for Mechanical Engineers/Transport Equations

< Fluid Mechanics for Mechanical Engineers

Diffusion Phenomena

The direction of diffusion of mass, energy and momentum

The fluid motion (convection) takes place with two mechanisms:

The molecular motion is responsible for the diffusion phenomena. In fluid mechannics, we consider mostly the diffusion of mass, momentum and energy.

The diffusion of a quantity takes place always from rich to poor regions. For example, when one puts a drop of ink in water, ink molecules moves (diffuses) from ink rich region to ink poor region. Heat diffuses from high temperature regions to colder regions. Similarly, momentum diffuses from momentum rich regions to momentum poor regions. The diffusive flux in one diemension is formulated by Fick's first law as follows for mass, energy and momentum respectively:

J_{m}=-k_m\frac{d C}{d x} [\frac{kg}{m^2s}] \quad \text{where} \quad C \quad \text{is the concentration with the unit} \quad [mol/m^3] .

J_{E}=-k\frac{d T}{d x} [\frac{J}{m^2s}]

J_{mom}=-\mu\frac{d U}{d x} [\frac{kg\,m/s}{m^2s}]

The negative sign is used since the gradient is negative opposite to the direction of diffusion, i.e. with the negative sign the direction of diffusion is correctly assesed. The constants of diffusion are the diffusion constant of one species into another, thermal conductivity and the dynamic viscosity of the fluid, respectively. The diffusion of energy in the form of heat is also called as conduction.

General Form of Transport Equations

The RTT theorem is given only with advection and without any sources. In the more general form of transport equation, the diffusion and the sources of the transported quantities have to be considered. For an extensive quantity B , the transport equation is:

{\displaystyle \displaystyle \frac{dB}{dt}|_{sys} = \frac{\partial}{\partial t}\int_{cv}{\rho b dV} + \int_{cs}{\rho b \vec{U}\cdot\vec{n}dA}+\int_{cs}{\vec{J}_B\cdot\vec{n}dA}-\int_{cv}{\dot{q}_B dV} }

where \dot{q}_b is the source of B per unit volume within CV and \vec{J}_b is the diffusive flux of B through the CS of CV.

Conservation of Mass Equation (Continuity Equation)

The integral equation for the conservation of mass law is {\displaystyle \displaystyle \frac{dM}{dt}|_{sys} = \frac{\partial}{\partial t}\int_{cv}{\rho dV} + \int_{cs}{\rho \vec{U}\cdot\vec{n}dA}+\int_{cs}{\vec{J}_{m}\cdot\vec{n}dA}-\int_{cv}{\dot{q}_{m} dV}=0 }

According to this equation, when adjective and diffusive flux are ignored, it can be seen that the rate of change of mass in the CV is due to the mass sources:

{\displaystyle \frac{\partial}{\partial t}\int_{cv}{\rho dV} =\int_{cv}{\dot{q}_{m} dV} }In case of no adjective flux and mass sources, {\displaystyle \frac{\partial}{\partial t}\int_{cv}{\rho dV}=-\int_{cs}{\vec{J}_{m}\cdot\vec{n}dA} }the mass in CV, can change only through diffusive flux. Note that, incoming diffusive flux will have a negative sign and together with the extra negative sign its contribution will be positive.

For a differential volume dV=dx_1dx_2dx_3, each term in the integral equation can be formulated as follows:

The time rate of change of mass in the CV:{\displaystyle \frac{\partial}{\partial t}\int_{cv}{\rho dV}\approx\frac{\partial \rho}{\partial t}dV }

The advective flow rate of mass through CS:

The advective flow rate can be decomposed into the surfaces having their surface normal vectors along x_1, x_2 and x_3 axis:{\displaystyle \int_{cs}{\rho \vec{U}\cdot\vec{n}dA}=\int_{cs_{x1}}{\rho \vec{U}\cdot\vec{n}dA}+\int_{cs_{x2}}{\rho \vec{U}\cdot\vec{n}dA}+\int_{cs_{x3}}{\rho \vec{U}\cdot\vec{n}dA} }

{\displaystyle \int_{cs_{x1}}{\rho \vec{U}\cdot\vec{n}dA}\approx -\rho U_1dx_2dx_3+\rho U_1dx_2dx_3+{\partial \rho U_1 \over \partial x_1}dx_1dx_2dx_3={\partial \rho U_1 \over \partial x_1}dV }Hence the total advective flux is:

{\displaystyle \int_{cs}{\rho \vec{U}\cdot\vec{n}dA}\approx{\partial \rho U_1 \over \partial x_1}dV+{\partial \rho U_2 \over \partial x_2}dV+{\partial \rho U_3 \over \partial x_3}dV= {\partial \rho U_i \over \partial x_i}dV } In fact, one can see clearly the Divergence theorem: {\displaystyle \int_{cs}{\rho \vec{U}\cdot\vec{n}dA}=\int_{cv}\vec{\nabla}\cdot \rho \vec{U} dV=\int_{cv}{\partial \rho U_i \over \partial x_i}dV\approx{\partial \rho U_i \over \partial x_i}dV } The diffusive flow rate of mass through CS:

By using the divergence theorem, the diffusive flow rate can be written for a differential volume as follows:

{\displaystyle \int_{cs}{\vec{J}_{m}\cdot\vec{n}dA}=\int_{cv}\vec{\nabla}\cdot \vec{J}_{m} dV=\int_{cv}{\partial J_{m_i} \over \partial x_i}dV\approx{\partial J_{m_i} \over \partial x_i}dV }For isotropic diffusion, the i^{th} component of the diffusion flux is

{\displaystyle J_{m_i}=-k_m\frac{\partial C}{\partial x_i}}Hence the diffusive flow rate can be written as

{\displaystyle \int_{cs}{\vec{J}_{m}\cdot\vec{n}dA}\approx-k_m{\partial^2 C \over \partial x_i^2}dV }

The net source of mass in CV:

{\displaystyle \int_{cv}{\dot{q}_{m} dV}\approx\dot{q}_{m} dV }The final differential form of the conservation of mass equation reads

{\displaystyle \frac{\partial \rho}{\partial t}dV+{\partial \rho U_i \over \partial x_i}dV -k_m{\partial^2 C \over \partial x_i^2}dV-\dot{q}_{m} dV=0}and per nunit volume

{\displaystyle \frac{\partial \rho}{\partial t}+{\partial \rho U_i \over \partial x_i} -k_m{\partial^2 C \over \partial x_i^2}-\dot{q}_{m}=0}

Conservation of Energy Equation

The conservation of energy equaiton is very interesting since it shows the dissipation of energy and reversible conversion of energy as the fluid flows. The transport equation of energy reads:

{\displaystyle \displaystyle \frac{dE}{dt}|_{sys} =[\dot{Q}+\dot{W}]_{on\ sys}= \frac{\partial}{\partial t}\int_{cv}{\rho e dV} + \int_{cs}{\rho e \vec{U}\cdot\vec{n}dA}+\int_{cs}{\vec{J}_E\cdot\vec{n}dA}-\int_{cv}{\dot{q}_E dV} }where e=u+\frac{U_iU_i}{2} is the sum of internal energy and kinetic energy per unit mass.

Rearranging the above equaiton for an instant at whşch CV and the system collapse on each other:

{\displaystyle \frac{\partial}{\partial t}\int_{cv}{\rho e dV} + \int_{cs}{\rho e \vec{U}\cdot\vec{n}dA}=[\dot{Q}+\dot{W}]_{on\ cv}-\int_{cs}{\vec{J}_E\cdot\vec{n}dA}+\int_{cv}{\dot{q}_E dV} }Heat \dot{Q} can be added to the CV in the form of radiation or a surface heater. The source of energy \dot{q}_E can be the energy released or absorbed during chemical reaction or the energy of the added mass from the source of mass. Therefore, \dot{Q} can be treated as a part of energy source \dot{q}_E if it is a volumetric addition and/or as a part of conduction \vec{J}_E term on the surface of the CV if it is a heating surface.

The rate of change of energy in a differential volume is

{\displaystyle \frac{\partial}{\partial t}\int_{cv}{\rho e dV} \approx \frac{\partial \rho e}{\partial t}dV }

The advection term for a differential volume can be formulated by using the divergence theorem

{\displaystyle \int_{cs}{\rho e \vec{U}\cdot\vec{n}dA}=\int_{cv}\vec{\nabla}\cdot({\rho e \vec{U})dV} \approx \frac{\partial \rho e U_i}{\partial x_i} dV }Body and surface forces work on CV (\dot{W})

{\displaystyle \dot{W}_{body}=\int_{cv} \vec{U}\cdot \rho \vec{g} dV }{\displaystyle \dot{W}_s=\int_{cs} \vec{U} \cdot \vec{\sigma} dA }where \vec{\sigma} is the sum of the stress caused by pressure and viscous stresses

{\displaystyle \vec{\sigma}=-P\vec{n}+\vec{\tau} }Hence, {\displaystyle \dot{W}=\int_{cv} \vec{U}\cdot \rho \vec{g} dV -\int_{cs} P \vec{U} \cdot \vec{n} dA + \int_{cs} \vec{U} \cdot \vec{\tau} dA }Since \vec{\tau}=\bar{\bar{\tau}}\cdot \vec{n} is the local stress vector and \bar{\bar{\tau}} is the stress tensor.

{\displaystyle \dot{W}=\int_{cv} \vec{U}\cdot \rho \vec{g} dV -\int_{cs} P \vec{U} \cdot \vec{n} dA + \int_{cs} (\vec{U} \cdot \bar{\bar{\tau}}) \cdot \vec{n} dA }

For a diffrential volume, \vec{\tau}=\bar{\bar{\tau}}\cdot \vec{n}=\tau_{ij}n_i is the stress vector on a face having its nomal in i^{th}-direction.

{\displaystyle \dot{W}\approx \rho U_i g_i dV - \frac{\partial P U_i}{\partial x_i}dV + \frac{\partial U_j \tau_{ij}}{\partial x_i}dV }In order to have a better decomposition later, one can switch i and j index since \tau_{ij}=\tau_{ji} and use the equality \frac{\partial U_j \tau_{ij}}{\partial x_i}=\frac{\partial U_i \tau_{ji}}{\partial x_j}

{\displaystyle \dot{W}\approx \rho U_i g_i dV - \frac{\partial P U_i}{\partial x_i}dV + \frac{\partial U_i \tau_{ji}}{\partial x_j}dV }

The heat diffusion (conduction) term for a differential volume can be written by using the divergence term as:

{\displaystyle \int_{cs}{\vec{J}_E\cdot\vec{n}dA}=\int_{cv} \vec{\nabla} \cdot \vec{J}_E dV= -\int_{cv} k{\partial^2 T \over \partial x_i^2 } dV \approx -k{\partial^2 T \over \partial x_i^2 } dV }Rewriting the energy equation, using the above forms derived for a differential volume:

{\displaystyle \frac{\partial \rho e}{\partial t}dV+\frac{\partial \rho e U_i}{\partial x_i} dV=\rho U_i g_i dV - \frac{\partial P U_i}{\partial x_i}dV + \frac{\partial U_i \tau_{ji}}{\partial x_j}dV +k{\partial^2 T \over \partial x_i^2 } dV+\dot{q}_EdV }or per unit volume

{\displaystyle \frac{\partial \rho e}{\partial t}+\frac{\partial \rho e U_i}{\partial x_i}=\rho U_i g_i - \frac{\partial P U_i}{\partial x_i} + \frac{\partial U_i \tau_{ji}}{\partial x_j} +k{\partial^2 T \over \partial x_i^2 } + \dot{q}_E }

When left side of the equation is expanded,

{\displaystyle e \left(\frac{\partial \rho}{\partial t}+\frac{\partial \rho U_i}{\partial x_i}\right)+ \rho \left(\frac{\partial e}{\partial t}+U_i\frac{\partial e }{\partial x_i}\right) =\rho U_i g_i - \frac{\partial P U_i}{\partial x_i} + \frac{\partial U_i \tau_{ji}}{\partial x_j} +k{\partial^2 T \over \partial x_i^2 } + \dot{q}_E }it can be seen that the first term involves one part of the conservation of mass equation. Neglecting the diffusion and source of mass, the continuity equation can be set to zero and the energy equation can be reduced to:

{\displaystyle \rho \left(\frac{\partial e}{\partial t}+U_i\frac{\partial e }{\partial x_i}\right) =\rho U_i g_i - \frac{\partial P U_i}{\partial x_i} + \frac{\partial U_i \tau_{ji}}{\partial x_j} +k{\partial^2 T \over \partial x_i^2 } + \dot{q}_E }The left side becomes the substantial derivative of the energy per unit time:

{\displaystyle \frac{D e}{D t} =U_i g_i - \frac{1}{\rho}\frac{\partial P U_i}{\partial x_i} + \frac{1}{\rho}\frac{\partial U_i \tau_{ji}}{\partial x_j} +\frac{1}{\rho}k{\partial^2 T \over \partial x_i^2 } + \frac{1}{\rho}\dot{q}_E }

Decomposition of Energy Equation

Inserting the definition of e=u+U_iU_i/2 and expanding the work done by pressure and viscous stresses, {\displaystyle \frac{D u}{D t}+\frac{D U_iU_i/2}{D t} =U_i g_i - \frac{P}{\rho}\frac{\partial U_i}{\partial x_i} - \frac{U_i}{\rho}\frac{\partial P}{\partial x_i} + \frac{U_i}{\rho}\frac{\partial \tau_{ji}}{\partial x_j} + \frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} +\frac{1}{\rho}k{\partial^2 T \over \partial x_i^2 } + \frac{1}{\rho}\dot{q}_E }Since \frac{D U_iU_i/2}{D t}=U_i\frac{D U_i}{D t} ,

{\displaystyle \frac{D u}{D t}+\underline{U_i\frac{D U_i}{D t}} =\underline{U_i g_i} - \frac{P}{\rho}\frac{\partial U_i}{\partial x_i} \underline{- \frac{U_i}{\rho}\frac{\partial P}{\partial x_i}} + \underline{\frac{U_i}{\rho}\frac{\partial \tau_{ji}}{\partial x_j}} + \frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} +\frac{1}{\rho}k{\partial^2 T \over \partial x_i^2 } + \frac{1}{\rho}\dot{q}_E }Now we can decompose this equation into the mechanical and thermal energy equations. The underlined terms are nothing but the scalar product of the velocity vector and the momentum equation \left(\vec{U}\cdot \frac{D\vec{U}}{Dt}\right) and form the mechanical energy equation, the remaining terms form the thermal energy equation.

The mechanical energy equation:{\displaystyle U_i\frac{D U_i}{D t} =U_i g_i - \frac{U_i}{\rho}\frac{\partial P}{\partial x_i} +\frac{U_i}{\rho}\frac{\partial \tau_{ji}}{\partial x_j} }Since {\displaystyle \frac{U_i}{\rho}\frac{\partial \tau_{ji}}{\partial x_j} =\frac{1}{\rho}\frac{\partial U_i \tau_{ji}}{\partial x_j} -\frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} }and{\displaystyle \frac{U_i}{\rho}\frac{\partial P}{\partial x_i} =\frac{1}{\rho}\frac{\partial U_i P}{\partial x_i} - \frac{P}{\rho}\frac{\partial U_i}{\partial x_i} }The mechanical energy equation can be written in a more interpretable manner:

{\displaystyle U_i\frac{D U_i}{D t} =U_i g_i - \frac{1}{\rho}\frac{\partial U_i P}{\partial x_i} + \frac{P}{\rho}\frac{\partial U_i}{\partial x_i} + \frac{1}{\rho}\frac{\partial U_i \tau_{ji}}{\partial x_j} -\frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} }The thermal energy equation:{\displaystyle \frac{D u}{D t} =- \frac{P}{\rho}\frac{\partial U_i}{\partial x_i} + \frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} +\frac{1}{\rho}k{\partial^2 T \over \partial x_i^2 } + \frac{1}{\rho}\dot{q}_E }The meaning of each term is as follows:

Term Physical meaning
{\displaystyle - \frac{1}{\rho}\frac{\partial U_i P}{\partial x_i} } Rate of work done by pressure.
{\displaystyle \frac{P}{\rho}\frac{\partial U_i}{\partial x_i} } Rate of reversible conversion of kinetic energy into internal energy. (-) in the case of compression and (+) in the case of expansion, because \frac{\partial U_i}{\partial x_i}<0 for compression and larger than 0 for expansion. It appears with a negative sign in the internal energy equation.
{\displaystyle \frac{1}{\rho}\frac{\partial U_i \tau_{ji}}{\partial x_j} } Rate of work done by viscous stresses.
{\displaystyle -\frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} } Rate of dissipation of mechanical energy, in other words irreversible conversion of kinetic energy to internal energy. It is always negative. In other words, mechanical energy drops due to this term. In the internaş energy equation it appears with a possitive sign.
{\displaystyle -\frac{P}{\rho}\frac{\partial U_i}{\partial x_i} } Rate of reversible conversion of kinetic energy into internal energy. (+) in the case of compresion and (-) in the case of expansion, because \frac{\partial U_i}{\partial x_i}<0 for compression and larger than 0 for expansion.
{\displaystyle \frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} } Rate of increase of internal energy by irreversible viscous dissipation. It is always positive. In other words, internal energy increases due to this term.
{\displaystyle \frac{1}{\rho}k{\partial^2 T \over \partial x_i^2 } } Rate of heat added or extracted via conduction (thermal diffusion). It can be (+) or (-).
{\displaystyle \frac{1}{\rho}\dot{q}_E } Rate of head added via combustion or radiation.

Comment on the dissipation term

The dissipation term \frac{\tau_{ji}}{\rho}\frac{\partial U_i}{\partial x_j} is always positive and extracts energy from the mechanical energy equation and inreases the internal energy. As can be seen, without the viscous stresses or velocity gradients there is no dissipation. In fact, viscous stresses do occur only when there is a velocişty gradient. Hence, in order to avoid dissipation one should avoid or reduce unnecessary velocity gradients.

This article is issued from Wikiversity - version of the Tuesday, April 05, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.