# Convection Diffusion Application

**The numerical solution of convection-diffusion transport problems arises in many important applications in science and engineering. These problems occur in many applications such as in the transport of air and ground water pollutants, oil reservoir flow, in the modeling of semiconductors, and so forth. This paper describes the Convection Diffusion Applications for solving this equation.**

## Contents |

## Theory

Under the assumption of incompressibility, the governing equations are given by

(1)

(2)

In the context of mass difussion within a fluid, (1) is is the mass
conservation equation and (2) is a contitutive law proposed by Fourier. The
notation is standard: ρ is the density, *C* the heat capacity, κ
the thermal conductivity, *T* is the temperature, *v* is the velocity field
and *q* is the diffusive flux per unit fluid density.

Remark: systen can be decoupled since we can plug (2) into (1) and solve the
scalar equation

(3)

**Problem statement**

Let us consider the transport by convection and diffusion in an open set
Ω (d=2 or 3) \ with piecewise smooth boundary Γ,
such that . The unit outward normal vector to
Γ is denoted *n*. The convection-diffusion initial-boundary value
problem can be stated as follows: given a divergence-free velocity field a,
the diffusion tensor κ and adequate initial and boundary conditions, find *T* :
such
that

(4)
*T*(*x*,0) = *T*_{0}(*x*)*o**n*Ω

#### Space discretization method

Multiplying Eq.(4) by a test function *W* and intehrating on the whole domain
Ω the equation reads

(5)

Integratin by parts the right term of Eq.(5) leads to

(6)

**Finite element discretization**

The temperature is discretized in the standard finite element method manner as
(7)
where *N*_{i} are the nodal shape functions. Substituting the finite element
approximation (7) into the variational equation () and choosing a Galerling formulation (*W*_{i} = *N*_{i}) leads to the following equation:

(7)

#### Time discretization method

Consider a first-order BDF (that is, the Euler implicit scheme)

or a second-order BDF

## How to analyse using the current application

We will use GiD as pre and post-processor link gid .

## Contact

For any additional information please contact with marti.jm@gmail.com