Generic scalar transport equation
From CFD-Wiki
(+wikilink) |
|||
(3 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
- | A differential equation | + | A differential equation expresses a certain conservation principle. Whether be it in electromagnetics, fluid dynamics, heat transfer, radiation, electronics... conservation principles are the basis for the derivation of differential or integro-differential equations. In this respect, any differential equation addresses a certain quantity as it dependent variable and thus expresses the balance between the phenomena affecting the evolution of this quanitity. For example, the temperature of a fluid in a heated pipe is affected by convection due to the solid-fluid interface, and due to the fluid-fluid interaction. Furthermore, temperature is also diffused inside the fluid. For a state state problem, with the absence of sources, a differential equation governing the temperature will definetely express a balance between convection and diffusion. |
- | A brief inspection of the equations governing various physical phenomena will reveal that all of these equations can be put into a generic form thus allowing a systematic approach for a computer simulation. | + | |
- | For example, the conservation equation of a chemical species <math> c_i </math> is | + | A brief inspection of the equations governing various physical phenomena will reveal that all of these equations can be put into a generic form thus allowing a systematic approach for a computer simulation. For example, the conservation equation of a chemical species <math> c_i </math> is |
- | <math> \frac{\partial{\rho c_i}}{\partial t} + \nabla \cdot (\rho \vec u c_i + \vec J) = R_i </math> | + | |
- | where <math> \vec u </math> denotes the velocity field, <math> \vec J </math> denotes the diffusion flux the of the chemical species, and <math> R_i </math> denotes the rate of generation of <math> R_i </math> caused by the chemical reaction. | + | :<math> \frac{\partial{\rho c_i}}{\partial t} + \nabla \cdot (\rho \vec u c_i + \vec J) = R_i </math> |
- | The '''x-momentum''' equation for a Newtonian fluid can be written as | + | |
- | <math> \frac{\partial{\rho u}}{\partial t} + \nabla \cdot (\rho \vec u u ) =\nabla \cdot (\mu \nabla u ) - \frac {\partial p}{\partial x} + B_x + V_x </math> | + | where <math> \vec u </math> denotes the velocity field, <math> \vec J </math> denotes the diffusion flux the of the chemical species, and <math> R_i </math> denotes the rate of generation of <math> R_i </math> caused by the chemical reaction. |
- | where <math> B_x </math> is the body force in the x-direction and <math>V_x</math> includes the viscous terms that are not expressed by <math>\nabla \cdot (\mu \nabla u )</math | + | |
- | + | The '''x-momentum''' equation for a Newtonian fluid can be written as | |
- | Upon inspection of the above equations, it can be infered that all the dependent variables seem to obey a generalized conservation principle. If the dependent variable (scalar or vector) is denoted by <math>\phi</math>, the '''generic''' differential equation is | + | |
- | <math> \underbrace{ \frac{\partial{\rho \phi}}{\partial t}}_{Transient \ term} + \underbrace{ \nabla \cdot (\rho \vec u \phi )}_{Convection \ term} =\underbrace {\nabla \cdot (\Gamma \nabla \phi )}_{Diffusion \ term} + \underbrace {S_{\phi}}_{Source \ term}</math | + | :<math> \frac{\partial{\rho u}}{\partial t} + \nabla \cdot (\rho \vec u u ) =\nabla \cdot (\mu \nabla u ) - \frac {\partial p}{\partial x} + B_x + V_x </math> |
- | + | ||
- | where <math> \Gamma </math> is the '''diffusion''' coefficient or '''diffusivity''' <br>. | + | where <math> B_x </math> is the body force in the x-direction and <math>V_x</math> includes the viscous terms that are not expressed by <math>\nabla \cdot (\mu \nabla u )</math> |
+ | |||
+ | Upon inspection of the above equations, it can be infered that all the dependent variables seem to obey a generalized conservation principle. If the dependent variable (scalar or vector) is denoted by <math>\phi</math>, the '''generic''' differential equation is | ||
+ | |||
+ | :<math> \underbrace{ \frac{\partial{\rho \phi}}{\partial t}}_{Transient \ term} + \underbrace{ \nabla \cdot (\rho \vec u \phi )}_{Convection \ term} =\underbrace {\nabla \cdot (\Gamma \nabla \phi )}_{Diffusion \ term} + \underbrace {S_{\phi}}_{Source \ term}</math> | ||
+ | |||
+ | where <math> \Gamma </math> is the '''diffusion''' coefficient or '''diffusivity'''. <br> | ||
+ | * The transient term, <math> \frac{\partial{\rho \phi}}{\partial t} </math>, accounts for the accumulation of <math> \phi </math> in the concerned control volume <br> | ||
+ | * The convection term, <math>\nabla \cdot (\rho \vec u \phi )</math>, accounts for the transport of <math> \phi </math> due to the existence of the velocity field (note the velocity <math> \vec u </math> multiplying <math> \phi </math>) <br> | ||
+ | * The diffusion term, <math> \nabla \cdot (\Gamma \nabla \phi )</math>, accounts for the transport of <math> \phi </math> due to its gradients. <br> | ||
+ | * The source term, <math> S_{\phi} </math>, accounts for any sources or sinks that either create or destroy <math> \phi </math>. Any extra terms that cannot be cast into the convection or diffusion terms are considered as source terms. | ||
+ | |||
+ | The objective of all discretization techniques (Finite Difference, Finite Element, Finite Volume, Boundary Element...) is to devise a mathematical formulation to transform each of these terms into an algebraic equation. Once applied to all control volumes in a given mesh, we obtain a full linear system of equations that needs to be solved. | ||
+ | ==External link== | ||
+ | *[http://en.wikipedia.org/wiki/Generic_scalar_transport_equation Generic scalar transport equation article on Wikipedia] |
Latest revision as of 03:37, 14 February 2006
A differential equation expresses a certain conservation principle. Whether be it in electromagnetics, fluid dynamics, heat transfer, radiation, electronics... conservation principles are the basis for the derivation of differential or integro-differential equations. In this respect, any differential equation addresses a certain quantity as it dependent variable and thus expresses the balance between the phenomena affecting the evolution of this quanitity. For example, the temperature of a fluid in a heated pipe is affected by convection due to the solid-fluid interface, and due to the fluid-fluid interaction. Furthermore, temperature is also diffused inside the fluid. For a state state problem, with the absence of sources, a differential equation governing the temperature will definetely express a balance between convection and diffusion.
A brief inspection of the equations governing various physical phenomena will reveal that all of these equations can be put into a generic form thus allowing a systematic approach for a computer simulation. For example, the conservation equation of a chemical species is
where denotes the velocity field, denotes the diffusion flux the of the chemical species, and denotes the rate of generation of caused by the chemical reaction.
The x-momentum equation for a Newtonian fluid can be written as
where is the body force in the x-direction and includes the viscous terms that are not expressed by
Upon inspection of the above equations, it can be infered that all the dependent variables seem to obey a generalized conservation principle. If the dependent variable (scalar or vector) is denoted by , the generic differential equation is
where is the diffusion coefficient or diffusivity.
- The transient term, , accounts for the accumulation of in the concerned control volume
- The convection term, , accounts for the transport of due to the existence of the velocity field (note the velocity multiplying )
- The diffusion term, , accounts for the transport of due to its gradients.
- The source term, , accounts for any sources or sinks that either create or destroy . Any extra terms that cannot be cast into the convection or diffusion terms are considered as source terms.
The objective of all discretization techniques (Finite Difference, Finite Element, Finite Volume, Boundary Element...) is to devise a mathematical formulation to transform each of these terms into an algebraic equation. Once applied to all control volumes in a given mesh, we obtain a full linear system of equations that needs to be solved.