Generic scalar transport equation

From Wikipedia, the free encyclopedia

The generic scalar transport equation is a general partial differential equation that describes transport phenomena such as heat transfer, mass transfer, fluid dynamics (momentum transfer), etc.

All the transfer processes express a certain conservation principle. In this respect, any differential equation addresses a certain quantity as its dependent variable and thus expresses the balance between the phenomena affecting the evolution of this quantity. 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 express a balance between convection and diffusion.

A brief inspection of the equations governing various transport phenomena 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 concentration of a substance ci is

\frac{\partial{\rho c_i}}{\partial{t}} + \nabla \cdot (\rho \vec u c_i + \vec J) = R_i

where \vec u denotes the velocity field, \vec J denotes the diffusion flux the of the chemical species, and Ri denotes the rate of generation of Ri caused by the chemical reaction.

The x-momentum equation for a Newtonian fluid can be written as

\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

where Bx is the body force in the x-direction and Vx includes the viscous terms that are not expressed by \nabla \cdot (\mu \nabla u ).

Upon inspection of the above equations, it can be inferred 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

\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}

where Γ is the diffusion coefficient, or diffusivity.

  • The transient term, \frac{\partial{\rho \phi}}{\partial t}, accounts for the accumulation of φ in the concerned control volume.
  • The convection term, \nabla \cdot (\rho \vec u \phi ), accounts for the transport of φ due to the existence of the velocity field (note the velocity \vec u multiplying φ).
  • The diffusion term, \nabla \cdot (\Gamma \nabla \phi ), accounts for the transport of φ due to its gradients.
  • The source term, Sφ, 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, etc.) 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.

[edit] External links