Abstract:We study random compressible viscous magnetohydrodynamic flows. Combining the Monte Carlo method with a deterministic finite volume method we solve the random system numerically. Quantitative error estimates including statistical and deterministic errors are analyzed up to a stopping time of the exact solution. On the life-span of an exact strong solution we prove the convergence of the numerical solutions. Numerical experiments illustrate rich dynamics of random viscous compressible magnetohydrodynamics.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is the error analysis of random solutions in the compressible magnetohydrodynamics (MHD) system. Specifically, the author combines the Monte Carlo method and the deterministic finite - volume method to numerically solve the stochastic system and analyzes the statistical error and the deterministic error until the stopping time of the exact solution. Within the existence time of the exact strong solution, the author proves the convergence of the numerical solution and shows the rich dynamics of stochastic viscous compressible magnetohydrodynamics through numerical experiments.
### Specific description of the problem
1. **Research background**:
- The compressible magnetohydrodynamics (MHD) system describes the time evolution of a conducting fluid with viscosity and compressibility in the presence of a magnetic field.
- The equations of this system include mass conservation, momentum conservation, and magnetic field equations, and require the divergence of the magnetic field to be zero (i.e., no magnetic monopoles).
2. **Research objectives**:
- Analyze the error of the Monte Carlo finite - volume method in solving the stochastic compressible MHD system.
- Provide quantitative error estimates, including statistical error and deterministic error.
- Prove the convergence of the numerical solution within the existence time of the exact strong solution.
3. **Main challenges**:
- Error estimates can only be obtained in the class of regular solutions, and these solutions may not exist globally.
- The norm of the regular solution will diverge when approaching the maximum existence time \(T_{\text{max}}\).
- Data, solutions, and \(T_{\text{max}}\) are random variables, and their stochastic properties need to be considered.
- Under non - homogeneous boundary conditions, the numerical stability results of the compressible MHD system have not been fully established, and appropriate numerical methods need to be developed.
4. **Main results**:
- The author proves the convergence of the numerical method within the existence time \([0, T_{\text{max}})\) of the exact strong solution.
- Provides quantitative error estimates on any interval \([0, T]\), where \(T < T_{\text{max}}\).
- By introducing the stopping time \(T_M\), the author can control the statistical error and the deterministic error in the error estimate.
### Mathematical model and method
- **Compressible MHD system**:
\[
\begin{aligned}
&\partial_t \varrho+\text{div}(\varrho \mathbf{u}) = 0, \\
&\partial_t (\varrho \mathbf{u})+\text{div}(\varrho \mathbf{u} \otimes \mathbf{u})+\nabla p=\text{div} \mathbf{S}+\text{curl} \mathbf{B} \times \mathbf{B}+\varrho \nabla G, \\
&\partial_t \mathbf{B}=\text{curl}(\mathbf{u} \times \mathbf{B})-\text{curl}(\zeta \text{curl} \mathbf{B}), \quad \text{div} \mathbf{B} = 0,
\end{aligned}
\]
where \(\varrho\) is density, \(\mathbf{u}\) is the velocity field, \(\mathbf{B}\) is the magnetic field, \(p\) is pressure, \(\mathbf{S}\) is the viscous stress tensor, \(\zeta\) is the magnetic diffusion coefficient, and \(G\) is the gravitational potential.
- **Numerical method**:
- Use the Monte Carlo method to handle stochastic data.
- Combine the finite - volume method for numerical discretization to ensure mass conservation, density positivity, and discrete energy balance.
- Control the divergence - free property of the magnetic field.
### Summary
This paper, through the stochastic compressible MHD system...