A convection-diffusion problem with a large shift on Duran meshes

Mirjana Brdar,Sebastian Franz,Hans-Goerg Roos

2023-04-25

Abstract:A convection-diffusion problem with a large shift in space is considered. Numerical analysis of high order finite element methods on layer-adapted Duran type meshes, as well as on coarser Duran type meshes in places where weak layers appear, is provided. The theoretical results are confirmed by numerical experiments.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the convection - diffusion problem with large - space displacement. Specifically, when the small perturbation parameter \(\varepsilon\) is very small, how to solve this type of singularly perturbed problem efficiently and accurately. Such problems are widely used in fields such as control theory and biological sciences, for example, chemical medium models, circadian rhythms, epidemiology, respiratory systems, tumor growth, and neural networks.

### Problem Description

The specific problem studied in the paper is the convection - diffusion equation in the following form:

\[
- \varepsilon u''(x) - b(x) u'(x) + c(x) u(x) + d(x) u(x - 1)=f(x), \quad x\in\Omega=(0,2)
\]

The boundary condition is:

\[
u(2) = 0
\]

The initial condition is:

\[
u(x)=\Phi(x), \quad x\in(- 1,0]
\]

where \(b\geq\beta>0\), \(d\geq0\) and \(0 < \varepsilon\ll1\) is a small perturbation parameter. In addition, it is assumed that \(c-\frac{b'}{2}-\frac{\|d\|_{L^\infty(1,2)}}{2}\geq\gamma>0\) and \(\Phi(0) = 0\).

### Main Challenges

1. **Large Displacement**: The problem contains a large - space displacement \(u(x - 1)\), which makes the problem more complex.
2. **Singular Perturbation**: Since \(\varepsilon\) is very small, the solution has boundary layers or internal layers in some regions (such as near \(x = 0\) and \(x = 1\)), which imposes higher requirements on numerical methods.
3. **High - Order Finite Element Method**: It is necessary to perform numerical analysis of the high - order finite element method on the Durán - type mesh to ensure the accuracy and stability of the solution.

### Solutions

The paper solves the above problems through the following methods:

1. **Durán - type Mesh**: The Durán - type hierarchical mesh is introduced. This mesh can generate mesh points without using transition points and has a certain robustness, that is, the numerical results are not strongly affected by parameter changes.
2. **High - Order Finite Element Method**: The high - order finite element method is applied on the Durán - type mesh to deal with the boundary layers and internal layers in the singularly perturbed problem.
3. **Error Analysis**: A detailed error estimate is provided to prove the effectiveness and convergence of the numerical method.

### Result Verification

The theoretical results are verified through numerical experiments, indicating that using the high - order finite element method on the Durán - type mesh can effectively solve the convection - diffusion problem with large - space displacement and has good convergence.

In summary, the main objective of this paper is to propose an effective numerical method to solve the singularly perturbed convection - diffusion problem with large - space displacement, and the effectiveness of this method is verified through theoretical analysis and numerical experiments.