Documentation

The transport equation describes how a scalar quantity is transported in a space. Usually, it is applied to the transport of a scalar field (e.g. chemical concentration, material properties or temperature) inside an incompressible flow. From the mathematical point of view, the transport equation is also called the convection-diffusion equation, which is a first-order PDE (partial differential equation). The convection-diffusion equation is the basis for the most common transportation models.

The transport equation can be seen as the generalization of the continuity equation\(^1\). While the continuity equation (extensively described in the article on incompressible flow) usually describes the conservation of mass, the convection-diffusion equation describes the continuity/conservation of any scalar field in any space. Let’s consider an infinitesimal portion of space and its boundaries, as described in Figure 1.

The continuity principle states that the rate of change for a scalar quantity in any differential control volume is given by the flow and diffusion of the scalar into and out of the system, along with any generation or consumption inside the control volume. In practice, it means that the variation of concentration of a certain quantity in the volume is given by the balance of this quantity flow across the boundary and the amount of quantity produced or removed in the volume. From the mathematical point of view, this balance is expressed by the following equation:

$$ \frac{\partial c}{\partial t}+\nabla \cdot j=S \tag{1}$$

where \(c\) is the scalar field to be analyzed, \(j\) is the flux of \(c\) through the boundary, and \(S\) is the source/sink term inside \(\Omega\). Equation 1 is nothing more than the balance of a scalar quantity inside the volume. The first term (\(\frac{\partial c}{\partial t}\)) represents the time-dependent variation of the scalar quantity field inside the control volume, the second term (\(\nabla\cdot j\)) represents the net balance of the quantity which enters and exits the control volume, and the third term (\(S\)) represents the amount of the scalar quantity “created” or “destroyed” inside the volume.

Equation 1 can be further detailed by developing its terms:

- Flux \(j\) can be divided into two terms: the convective and the diffusive terms. The convection term is the quantity of the transported field which moves across the boundaries because of the flow; thus it is proportional to the velocity and can be written as \(j_{convection}=cu\), where \(c\) is the transported scalar quantity and \(u\) is the velocity of the means which transports this quantity. The diffusion term is the transportation of the scalar quantity according to its gradient, so \(j_{diffusion}=D\nabla c\), where \(D\) is the diffusivity.
- The source term can be divided into a pure source term and a reaction term. The pure source term (\(S_S\)) represents the creation/destruction rate of the field inside the volume. The reaction term (\(S_R\)) describes the creation/destruction of the transported quantity as a reaction to this quantity itself; it is, therefore, proportional to the transported field and can be written as \(S_R=f(c)\), where \(f(c)\) is a function of the transported scalar field. \(S_R\) is quite uncommon for engineering applications and is often neglected.

Equation 1 can now be re-written in its fully developed form as:

$$ \frac{\partial c}{\partial t}+\nabla\cdot(D\nabla c)+\nabla\cdot (uc)=S_S+S_R \tag{2}$$

During thermal simulations, the temperature field (which is scalar) is transported according to the convection-diffusion equation. In this specific case, the following notation is commonly used by the science community:

$$ \frac{\partial (\rho c_{p} T)}{\partial t}+\nabla\cdot (-k\nabla T)+\nabla\cdot (uT)=Q+\epsilon\sigma A\left(T_{hot}^4-T_{cold}^4\right) \tag{3}$$

where:

- \(\rho\) is the material density
- \(c_p\) is the heat capacity
- \(T\) is the temperature
- \(k\) is the thermal conductivity
- \(Q\) is the volumetric heat flux
- \(\epsilon\) is the emissivity
- \(\sigma\) is the Stefan-Boltzmann constant
- \(A\) is the boundary surface in which heat is exchanged by radiation.

Equations 2 and 3 differ only for the notation and for the complexity of the reaction term, coming from the physical modeling of heat transfer phenomena\(^3\).

Equation 3 is often simplified in common engineering applications thanks to the following hypotheses:

- Homogeneous material (i.e. constant and uniform material parameters):

$$ \rho c_p\frac{\partial(T)}{\partial t}-k\Delta T+\nabla\cdot(uT) = Q+\epsilon\sigma A\left(T_{hot}^4-T_{cold}^4\right) \tag{4}$$

- Incompressible flow (i.e. \(\nabla\cdot u=0\)):

$$ \rho c_p\frac{\partial(T)}{\partial t}-k\Delta T+u\cdot\nabla T= Q+\epsilon\sigma A\left(T_{hot}^4-T_{cold}^4\right) \tag{5}$$

- Neglection of heat exchange by radiation:

$$ \rho c_p\frac{\partial(T)}{\partial t}-k\Delta T+u\cdot\nabla T= Q \tag{6}$$

- Steady state (i.e. \(\frac{\partial T}{\partial t}=0\)):

$$ -k\Delta T+u\cdot\nabla T= Q \tag{7}$$

Transportation models are commonly used to analyze the dispersion of a certain chemical component in a fluid — some pollution particles, for instance. In this case, the transport is defined as “passive transport”, because the presence of the chemical concentration does not affect the fluid flow. Referring to equation 2, the convection term represents the transportation of the chemical component with the fluid, while the diffusive term represents a chemical reaction or molecular diffusion phenomena that could occur in the flow. When a simple recirculation of a certain amount of chemical components is modeled, no explicit boundary conditions are needed for the simulation and only the initial chemical concentration distribution is required (Figure 2).

On the other side, explicit boundary conditions are needed when the injection of a chemical component in a region is simulated, e.g. a chimney stack (Figure 3). In this case, a boundary condition on the inlet boundary is needed to be imposed (for instance \(c=1\), if pure smoke is emitted into the atmosphere from the chimney or \(c=0.5\) if the chimney output is a mixture of 50% smoke and 50% air).

The level set is a particular family of transportation models in which a distance function (called level-set function) is transported \(^{2,3}\). This distance function is computed with respect to an interface (surface for 3D problems, line for 2D problems, and point for 1D problems) and must have two characteristics:

- It must be an eulerian function, i.e. \(||\nabla\varphi||=1\) where \(\varphi\) is the level-set function
- It has a sign convention i.e. it is positive on one side of the interface and negative on the other side.

For the sake of clarity, examples of the 1D and 2D distance function to be transported are shown in Figure 4 and Figure 5 respectively.

Level-set methods follow the principles of the classic transport equation, but they require a further computation step (called re-initialization) in order to maintain the Eulerian distance function condition. The main advantage of the use of a level-set function is the possibility to solve the transport equation just once, but to obtain many results in the post-processing. For instance, for multiphase flows, it is possible to compute the transported viscosity & the density fields, the surface tension, and to impose a mixing law across the interface with very little computational effort. In addition, the use of the level-set function allows dealing with fields that show high gradients (e.g. the density field for the simulation of a water droplet falling through the air). The transport of the density field would lead to numerical instabilities across the interface where the density gradient is very high, while the transport of the level-set is more stable and enables computing the density as post-processing. The transportation of a distance function has several engineering applications; in the following sections, we will report the most common ones.

The purpose of level-set methods in image processing is to identify and export the contour of a certain profile in the image. This process is called image segmentation and is widely used in computer graphics and medical applications\(^4\). It is commonly used to analyze ultrasounds and tomographies’ outputs in the localization of tumors; Animation 1 shows the segmentation process for an ultrasound of lungs. With the segmentation process, it is possible to identify objects inside a digital image (i.e. lungs) through the identification of its boundaries and, eventually, use this information for analysis of the pathology.

In the case of image segmentation, convection and source terms in equation 3 are not considered; the only means of transport is diffusion according to the gradient of the pixel density of the image.

The level-set method is also commonly used to simulate multiphase flows. In this case, the transport model is defined as “active transport” because the convection/diffusion of the interface (through the distance function) affects the flow. The distance from the interface influences the flow because it is used to compute the material parameters in the computational domain. Given a flow composed of two fluids characterized by different material properties \(\eta_1\) and \(\eta_2\), the value in any given position (or element in the case of a numerical simulation) is given by:

$$ \eta=\eta_1 H(\varphi)+\eta_2(1-H(\varphi)) \tag{8}$$

where \(H(\phi)\) is the Heaviside function defined as

$$ H(\varphi)=1 \qquad \forall\varphi\geq0 \tag{9}$$

$$ H(\varphi)=0 \qquad \forall\varphi<0 \tag{10}$$

In this way, any complex multiphase flow with non-uniform material properties can be modeled through the simple transport of a scalar field (the level-set function) and the interface can be tracked implicitly as the zero Isocontour of this function. In Figure 6, the results of the simulation of a rising bubble are shown; in this case, the interface represents the \(\varphi=0\) isoline, while the region where \(\varphi>0\) is colored in blue and the region where \(\varphi<0\) is colored in red.

Similarly to the simulation of multiphase flows, the level-set method is used for the modeling of material crystallization. The transport of the level-set function is the starting point to analyze the phase change and crystal growth within a material. In Figure 7, a 2D representation of the crystal structure of a material through level-set method is presented. Each microstructural grain is identified through its boundaries (the zero-level of the level-set function), while the color is based on the level-set function in order to distinguish the grain.

References

- https://www.io-warnemuende.de/tl_files/project/amber/workshops/Lagoon%20Ecosystem%20Modelling/SSK-Lecture4.pdf
- https://math.berkeley.edu/~sethian/Papers/sethian.vallambrosa.87.pdf
- https://math.berkeley.edu/~sethian/Papers/sethian.osher.88.pdf
- http://physbam.stanford.edu/~fedkiw/
- https://www.youtube.com/watch?v=1R9oJ4OrTD0
- http://chaire-digimu.cemef.mines-paristech.fr/wp-content/uploads/2017/01/ManuscrittheseScholtes.pdf

Last updated: September 3rd, 2021

We appreciate and value your feedback.

What's Next

What is Joule Heating?