Convection-Diffusion Equations (CDEs)

The equations

The convection--diffusion equations (CDEs) is a set of conservation laws, constituted by the continuity equation

$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$

and the general convection-diffusion equation (GCDE)

$\rho \frac{\partial \phi}{\partial t} + \rho \mathbf{v} \cdot \nabla \phi = \nabla \cdot \left( \Gamma_\phi \nabla \phi \right) + \dot{s}_\phi$

Here, $\rho$ is the density of the fluid, $\mathbf{v}$ is the velocity field, $\phi$ is a scalar magnitude of the fluid (such as the temperature or the concentration of a pollutant), $\Gamma_\phi$ is the diffusion coefficient and $\dot{s}_\phi$ is the source term. When $\phi$ is a vector magnitude (for instance the velocity field $\mathbf{v}$ ), the equation is

$\frac{\partial(\rho \phi)}{\partial t} + \nabla \cdot (\rho \mathbf{v} \otimes \mathbf{\phi}) = \nabla \cdot \left( \Gamma_\phi \nabla \phi \right) + \dot{s} _\phi$

where $\otimes$ denotes the exterior product of two vectors.

Notice that the GCDE for a scalar magnitude $\phi$ is ''combination'' of the linear transport equation

$\frac{\partial \phi}{\partial t} + \mathbf{v} \cdot \nabla \phi = \dot{s}_\phi$

and the diffusion/heat equation

$\frac{\partial \phi}{\partial t} = \Delta \phi + \dot{s}_\phi$

Therefore it is expected that the solution to a boundary-value problem involving the GCDE, when it exists, shares properties with the solution to a linear transport problem and a diffusion problem.

Project outline

The project focuses on solving the CDEs numerically following a finite volume approach. The problems considered always take place in rectangular domains, hence a cartesian mesh is suitable in order to solve them. The report is structured in the following way:

Introduction: brief summary of the project.
Convection-diffusion equations: rigorous derivation of the CDEs.
Numerical study of the convection-diffusion equations: discretization of the CDEs in a rectangular domain discretized by means of a cartesian mesh. The algorithm to solve a general transient problem is also given.
Diagonal flow case: numerical solution to a CDEs steady state problem, taking place in a square domain with a flow in the diagonal diretion.
Smith-Hutton case: numerical solution to a CDEs steady state problem, taking place in a rectangular domain with a ''circular'' flow.
Appendices: quick reference for some facts in Measure Theory, Ordinary Differential Equations and Numerical resolution of linear systems.

A C++ code was developed in order to solve these problems numerically.

Hereinafter, $\rho$ , $\Gamma$ are known constants, $\dot{s}_\phi = 0$ and $\mathbf{v}$ does not depend upon time. Under these hypothesis, the CDEs are

$\nabla \cdot \mathbf{v} = 0$

$\frac{\rho}{\Gamma} \mathbf{v} \cdot \nabla \phi = \Delta \phi$

Diagonal flow case

Let $\phi_\text{low} < \phi_\text{high}$ be two given constants, and let $L > 0$ be a fixed length giving the square domain $\Omega = (0,L) \times (0,L)$ . The velocity field is given by $\mathbf{v} = \frac{v_0}{\sqrt{2}} \mathbf{i} + \frac{v_0}{\sqrt{2}} \mathbf{j}$ , with $v_0 > 0$ a known constant. The the CDE is

$\frac{\rho}{\Gamma} \mathbf{v} \cdot \nabla \phi = \frac{1}{\sqrt{2} L} \frac{\rho v_0 L}{\Gamma} \left( \frac{\partial \phi}{\partial x} + \frac{\partial \phi}{\partial y} \right) = \beta \, \mathrm{Pe} \left( \frac{\partial \phi}{\partial x} + \frac{\partial \phi}{\partial y} \right) = \Delta \phi$

where $\mathrm{Pe}$ is Péclet's number. Let $C_1 = ([0,L) \times \{0\}) \cup ({L} \times [0,L))$ and $C_2 = (\{0\} \times (0,L]) \cup ((0,L] \times {L})$ be two curves in $\partial \Omega$ and consider the function

$g(x,y) = \left\{ \begin{aligned} &\phi_\text{low} & &\text{if } (x,y) \in C_1 \\ &\phi_\text{high} & &\text{if } (x,y) \in C_2 \\ &0 & &\text{otherwise} \end{aligned} \right.$

The diagonal flow case problem is the following boundary-value problem:

$\left\{ \begin{aligned} \Delta \phi - \left( \frac{\partial \phi}{\partial x} + \frac{\partial \phi}{\partial y} \right) \beta \, \mathrm{Pe} &= 0 & &\text{in } \Omega \\ \phi &= g & &\text{on } \partial \Omega \end{aligned} \right.$

Analytical study

The boundary-value problem is studied for two values of $\mathrm{Pe}$ :

For $\mathrm{Pe} = \infty$ , the previous boundary-value problem becomes a transport problem. The analytical solution is found via the method of characteristics. However, it is not a solution in the classical sense. The weak solution is not studied.
For $\mathrm{Pe} \in [0,\infty)$ , a second-order elliptic PDE is obtained. By means of energy methods in Sobolev spaces (explained in Chapter 6 of Lawrence C. Evan's excellent book ''Partial Differential Equations''), the existence of weak solution is proved, although it cannot be said whether it is unique or not. This solution turns out to be a $\mathcal{C}^\infty$ function in $\Omega$ .

Numerical study

For several values of Péclet's number in the range $[10^{-9}, 10^9]$ , the numerical solution is computed using the aforementioned C++ code.

Smith-Hutton case

Let $L > 0$ be a fixed length and the domain $\Omega = (-L,L) \times (0,L)$ . Take the velocity field given by

$\mathbf{v} = 2 y (1 - x^2) \mathbf{i} - 2 x (1 - y^2) \mathbf{j}$

defined in $\overline{\Omega}$ . Consider the curves $C_1 = [-L,0] \times \{ 0 \}$ , $C_2 = (\{-L\} \times (0,L)) \cup ([-L,L] \times \{ L \}) \cup (\{ L \} \times [0,L))$ and $C_3 = (0,L) \times \{0\}$ , which give a partition of $\partial \Omega$ . Now let $g \colon C_1 \cup C_2 \to \mathbb{R}$ be the function given by

$g(x,y) = \left\{ \begin{aligned} &1 + \tanh{(10(2x+1))} & &\text{if } (x,y) \in C_1 \\ &1 - \tanh{(10)} & &\text{if } (x,y) \in C_2 \end{aligned} \right.$

The boundary-value problem for the steady state Smith-Hutton problem is

$\left\{ \begin{aligned} \Delta \phi - \frac{\rho}{\Gamma} \mathbf{v} \cdot \nabla \phi &= 0 & &\text{in } \Omega \\ \phi &= g & &\text{on } C_1 \cup C_2 \\ \frac{\partial \phi}{\partial y} &= 0 & &\text{on } C_3 \end{aligned} \right.$

Hereinafter, we shall assume $L = 1 \ \mathrm{m}$ .

Velocity field

Smith-Hutton's problem velocity field studied, in particular, the existence and uniqueness of streamlines is proved. In addition, these lines are computed numerically for some initial conditions, yielding the following plot:

Analytical solution

Since the value $\rho / \Gamma$ appears in the PDE, it seems reasonable to assume that the solution will depend on this constant:

For $\rho / \Gamma = \infty$ , the boundary-value problem turns into a transport-like problem. The characteristics of this problem are the streamlines found in the previous section. Due to the nature of these curves, the analytical solution is harder to find. The expression for this solution is given in terms of another function, although it is not proved whether it is unique. Nonetheless, it is demonstrated that if two classical solutions to the transport-like problem exist, then they must be equal.
For $\rho / \Gamma \in [0,\infty)$ , a second-order elliptic PDE is obtained. The theory studied for the diagonal-flow case cannot be applied here (at least in a straightforward fashion), since on $C_3$ a Neumann boundary condition is prescribed (instead of a Dirichlet boundary condition).

Numerical solution

For several values in $\rho / \Gamma$ in the range $[10^{-9}, 10^9]$ , the numerical solution to the Smith-Hutton case is computed using the C++ code developed. Below some examples are shown:

Name		Name	Last commit message	Last commit date
Latest commit History 130 Commits
code		code
gnuplot		gnuplot
notes		notes
plots		plots
readme_images		readme_images
report		report
svg		svg
.gitattributes		.gitattributes
.gitignore		.gitignore
README.md		README.md
README.md.tex		README.md.tex
convection_diffusion_equations.pdf		convection_diffusion_equations.pdf

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Repository files navigation

Convection-Diffusion Equations (CDEs)

The equations

Project outline