Rui M.P. Almeida,José C.M. Duque,Belchior C.X. Mário
Abstract:In this paper we make a study of a partial integral differential equation with $p$-Laplacian using a mixed finite element method. Two stable and convergent fixed point schemes are proposed to solve the nonlinear algebraic system. Using the implementation of the method in Matlab environment, we numerically analyse the convergence with an example. Some other examples are presented in order to illustrate several asymptotic behaviours and some localization effects of the solutions.
What problem does this paper attempt to address?
This paper attempts to solve the problem of numerical solutions of partial integro - differential equations (PIDE) with p - Laplacian and memory terms. Specifically, the author focuses on the following evolutionary integro - differential equations:
\[
\begin{cases}
u_t(x,t)-\Delta_p u(x,t)=\int_0^t g(t - s)\Delta_p u(x,s)\,ds + f(x,t),&(x,t)\in\Omega\times(0,T],\\
u(x,t) = 0,&(x,t)\in\partial\Omega\times(0,T],\\
u(x,0)=u_0(x),&x\in\Omega,
\end{cases}
\]
where:
- $\Delta_p u=\text{div}(|\nabla u|^{p - 2}\nabla u)$ is the p - Laplacian operator;
- $g(t - s)$ represents the memory term;
- $f(x,t)$ and $u_0(x)$ are given functions;
- $\Omega\subset\mathbb{R}^n$ is a bounded region with Lipschitz - continuous boundary.
### Research Background and Motivation
1. **Applications of p - Laplacian**: The p - Laplacian appears widely in problems such as porous media flow, power - law materials, nonlinear diffusion and filtration, non - Newtonian fluids, and elastic - plastic torsional creep.
2. **Influence of Memory Terms**: The introduction of memory terms is to reflect the influence of the system's past history on its current state, such as the material memory effect in heat conduction, path - dependence in financial mathematics, etc.
3. **Challenges of Nonlinear Problems**: Due to the existence of nonlinear terms (such as p - Laplacian and memory terms), it is very difficult to solve these equations analytically. Therefore, it is necessary to develop effective numerical methods to study the dynamic behavior of such equations.
### Main Contributions
1. **Mixed Finite Element Method**: The author proposes a numerical solution method based on the mixed finite element method, combined with the Crank - Nicolson method and the trapezoidal integration formula for time discretization.
2. **Stability and Convergence Analysis**: Two fixed - point iteration schemes are proposed to solve the nonlinear algebraic system, and their stability and convergence are proved.
3. **Numerical Simulation and Verification**: Through numerical experiments in the Matlab environment, the convergence of the method is verified, and the influence of different parameters and initial conditions on the solution behavior is shown, including asymptotic behavior, localization effects, etc.
### Conclusion
The research in this paper is a supplement to the theoretical work of Antontsev et al. and Almeida et al., providing a simple and robust numerical method that can effectively simulate partial integro - differential equations with p - Laplacian and memory terms and observe their asymptotic behavior and localization characteristics.