Assionvi H. Kouevi,Gabriel J. Lord
Abstract:This paper is focussed on the numerical resolution of diffusion advection and reaction equations (DAREs) with special features (such as fractures, walls, corners, obstacles or point loads) which globally, as well as locally, have important effects on the solution. We introduce a multilevel and local time solver of DAREs based on the discontinuous Galerkin (DG) method for the spatial discreization and time stepping methods such as exponential time differencing (ETD), exponential Rosenbrock (EXPR) and implicit Euler (Impl) methods. The efficiency of our solvers is shown with several experiments on cyclic voltammetry models and fluid flows through domains with fractures.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to solve the problems encountered in the numerical solution of diffusion - advection - reaction equations (DAREs) with special features such as fractures, walls, corners, obstacles or point loads. These special features have important impacts on both the global and local behaviors of the solution. Specifically, the paper focuses on the following aspects:
1. **The influence of complex geometric structures**: Special features (such as fractures, walls, etc.) will cause changes in the flow and transport of solutes in local areas, which poses a challenge to numerical simulation.
2. **The choice of time step**: In order to accurately capture these local behaviors, spatial local refinement is required. However, this requires reducing the time step (\(\Delta t\)). To ensure the stability of explicit time integrators and the accuracy of implicit time integrators. If a smaller time step is uniformly applied throughout the simulation domain, it will lead to excessive computation time. Therefore, the Local Time Stepping (LTS) method needs to be used to improve efficiency.
3. **Efficient solution methods**: The paper proposes a multi - level and local - time - step solver, which is based on the Discontinuous Galerkin (DG) method for spatial discretization and combines time - integration methods such as Exponential Time Differencing (ETD), Exponential Rosenbrock (EXPR) and Implicit Euler (Impl). Through this method, the CPU time can be significantly reduced while ensuring accuracy.
4. **Verification in practical applications**: The paper verifies the effectiveness of the proposed method through multiple experiments, including cyclic voltammetry models and fluid - flow problems with fractures.
### Summary
The core problem of the paper is to develop an efficient numerical method for solving diffusion - advection - reaction equations with complex geometric features. By introducing local time steps and the Discontinuous Galerkin method, the computational efficiency can be significantly improved while ensuring accuracy.