Abstract:We address in this paper a nonlinear parabolic system, which is built to retain the main mathematical difficulties of the P1 radiative diffusion physical model. We propose a finite volume fractional-step scheme for this problem enjoying the following properties. First, we show that each discrete solution satisfies a priori L -estimates, through a discrete maxi- mum principle; by a topological degree argument, this yields the existence of a solution, which is proven to be unique. Second, we establish uniform (with respect to the size of the meshes and the time step) L2 -bounds for the space and time translates; this proves, by the Kolmogorov theorem, the relative compactness of any sequence of solutions obtained through a sequence of discretizations the time and space steps of which tend to zero; the limits of converging subsequences are then shown to be a solution to the continuous problem. Estimates of time translates of the discrete solutions are obtained through the formalization of a generic argument, interesting for its own sake.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is how to construct and analyze a fractional - step finite - volume scheme to deal with the nonlinear parabolic system extracted from the simplified radiative transfer physical model (the so - called P1 model). Specifically, this system aims to retain the main mathematical difficulties of the initial physical model and prove the existence and uniqueness of the solutions of this discrete scheme, as well as its convergence to the solutions of the continuous problem.
### System Description
The nonlinear parabolic system studied in the paper is as follows:
\[
\begin{cases}
\frac{\partial u}{\partial t}-\Delta u + u^{4}-\phi = 0 & \text{in }(x,t)\in\Omega\times(0,T),\\
\phi-\Delta\phi - u^{4}= 0 & \text{in }(x,t)\in\Omega\times(0,T),\\
u(x,0)=u_{0}(x) & \text{in }x\in\Omega,\\
u(x,t) = 0 & \text{in }(x,t)\in\partial\Omega\times(0,T),\\
\nabla\phi\cdot n = 0 & \text{in }(x,t)\in\partial\Omega\times(0,T),
\end{cases}
\]
where $\Omega$ is a bounded connected subset in $\mathbb{R}^{d}$, $d = 2$ or $3$, and is assumed to be a polygonal or polyhedral region; $T<\infty$ is the final time; $u$ and $\phi$ are real - valued functions defined on $\Omega\times[0,T)$, $\partial\Omega$ represents the boundary of $\Omega$, and $n$ is the outer normal vector.
### Main Challenges
1. **Mathematical Difficulty**: This system contains nonlinear terms (such as $u^{4}$), which increases the complexity of the solution.
2. **Discretization Method**: It is necessary to design an effective discretization method to ensure the existence, uniqueness and convergence of the solution.
3. **Numerical Stability**: Ensure the stability and non - negativity of the numerical solution, which is very important for practical applications.
### Solution
The paper proposes a fractional - step finite - volume scheme, which has the following characteristics:
- **Discrete Maximum Principle**: Each discrete solution satisfies an a priori $L^{\infty}$ estimate, and the non - negativity of the solution is guaranteed by the discrete maximum principle.
- **Existence and Uniqueness**: Through topological degree arguments, the existence and uniqueness of the discrete solution are proved.
- **Convergence**: Through Kolmogorov's theorem, the relative compactness of the discrete solution sequence is proved, thereby proving that the discrete solution converges to the weak solution of the continuous problem.
### Conclusion
The paper successfully constructs a finite - volume scheme, solves the discretization problem of the above - mentioned nonlinear parabolic system, and proves the existence, uniqueness and convergence of the solutions of this scheme. This not only verifies that there indeed exists a weak solution to this problem, but also provides a theoretical basis and numerical tools for subsequent research.