Stephen Metcalfe,Thomas P. Wihler
Abstract:This work is concerned with the development of an adaptive numerical method for semilinear heat flow models featuring a general (possibly) nonlinear reaction term that may cause the solution to blow up in finite time. The fully discrete scheme consists of a high order discontinuous Galerkin (dG) time stepping method and a conforming finite element discretisation (cG) in space. The proposed adaptive procedure is based on rigorously devised conditional a posteriori error bounds in the $L^{\infty}(L^{\infty})$ norm. Numerical experiments complement the theoretical results.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an adaptive numerical method for solving the semilinear heat flow model with a nonlinear reaction term that may blow - up in a finite time. Specifically, this method aims to provide an efficient and accurate numerical framework to handle cases where the solution tends to infinity in a finite time due to the growth of the reaction term and the sufficiently large energy of the initial data.
### Specific description of the problem
Consider the following initial - boundary value problem:
\[
u_t-\kappa\Delta u = f(u)\quad\text{in }\Omega,\;t > 0,
\]
\[
u = 0\quad\text{on }\partial\Omega,\;t > 0,
\]
\[
u(\cdot,0)=u_0\quad\text{in }\Omega,
\]
where:
- $\Omega\subset\mathbb{R}^d$ is a bounded polyhedral region, $d = 1,2,3$.
- $\kappa>0$ is a constant diffusion coefficient.
- $u_0\in W^{2,\infty}(\Omega)$ is the initial condition, and satisfies $u_0|_{\partial\Omega}=0$.
- The reaction term $f:\Omega\times\mathbb{R}_0^+\times\mathbb{R}\to\mathbb{R}$ is continuously differentiable and satisfies the local Lipschitz estimate:
\[
|f(x,t,v)-f(x,t,w)|\leq L(t,|v|,|w|)|v - w|\quad\forall x\in\Omega,\;\forall t\in\mathbb{R}_0^+,\;\forall v,w\in\mathbb{R},
\]
where $L:\mathbb{R}_0^+\times\mathbb{R}_0^+\times\mathbb{R}_0^+\to\mathbb{R}_0^+$ is a known function that satisfies certain regularity conditions.
If the reaction term $f$ has sufficient growth and the initial data $u_0$ has sufficient energy, then the solution may blow - up in a finite time, that is, there exists a maximum existence time $T_\infty<\infty$ such that
\[
\|u(t)\|_{L^\infty(\Omega)}<\infty\quad\text{for }0 < t < T_\infty,\quad\lim_{t\nearrow T_\infty}\|u(t)\|_{L^\infty(\Omega)}=\infty.
\]
### Main objectives of the paper
The main objective of the paper is to develop an adaptive numerical method by introducing the high - order discontinuous Galerkin (dG) time discretization method and the conforming finite - element space discretization method, combined with conditional a posteriori error estimates. This method can automatically adjust the time and space step sizes when approaching the blow - up time, thereby effectively capturing the behavior of the solution and accurately approximating the blow - up time.
### Key features of the method
1. **High - order time discretization**: Using the high - order dG method for time discretization, it is possible to use larger time steps and a higher order of approximation in the pre - blow - up phase, and switch to a low - order approximation when approaching the blow - up time.
2. **Conditional a posteriori error estimate**: Based on the conditional a posteriori error estimate, ensure that the error estimate is only valid in the pre - blow - up phase and is robust with respect to the distance from the blow - up time.
3. **Adaptive strategy**: Dynamically adjust the time and space step sizes through an adaptive algorithm to guide the numerical solution process to accurately approach the blow - up time.
In summary, this paper aims to provide an efficient numerical method that can handle the semilinear heat flow model with complex nonlinear reaction terms and blow - up in a finite time.