Abstract:The regularity of solutions to the stochastic nonlinear wave equation plays a critical role in the accuracy and efficiency of numerical algorithms. Rough or discontinuous initial conditions pose significant challenges, often leading to a loss of accuracy and reduced computational efficiency in existing methods. In this study, we address these challenges by developing a novel and efficient numerical algorithm specifically designed for computing rough solutions of the stochastic nonlinear wave equation, while significantly relaxing the regularity requirements on the initial data. By leveraging the intrinsic structure of the stochastic nonlinear wave equation and employing advanced tools from harmonic analysis, we construct a time discretization method that achieves robust convergence for initial values \((u^{0}, v^{0}) \in H^{\gamma} \times H^{\gamma-1}\) for all \(\gamma > 0\). Notably, our method attains an improved error rate of \(O(\tau^{2\gamma-})\) in one and two dimensions for \(\gamma \in (0, \frac{1}{2}]\), and \(O(\tau^{\max(\gamma, 2\gamma - \frac{1}{2}-)})\) in three dimensions for \(\gamma \in (0, \frac{3}{4}]\), where \(\tau\) denotes the time step size. These convergence rates surpass those of existing numerical methods under the same regularity conditions, underscoring the advantage of our approach. To validate the performance of our method, we present extensive numerical experiments that demonstrate its superior accuracy and computational efficiency compared to state-of-the-art methods. These results highlight the potential of our approach to enable accurate and efficient simulations of stochastic wave phenomena even in the presence of challenging initial conditions.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of calculating rough solutions to **stochastic nonlinear wave equations**. Specifically, the article focuses on how to effectively solve stochastic nonlinear wave equations numerically when the initial conditions are rough or discontinuous, and ensure the accuracy and computational efficiency of the solutions.
#### Main challenges:
1. **Low regularity of initial conditions**: Traditional numerical methods usually require solutions to have high regularity (such as \(C([0, T]; H^1\times L^2)\)), but in practical applications, the initial conditions may be rough or even contain discontinuities, which causes the solutions to fail to meet these high - regularity requirements.
2. **Decrease in computational efficiency and accuracy**: Rough initial conditions can lead to a significant reduction in the convergence order of existing numerical methods, thus affecting computational efficiency and the accuracy of results.
#### Solutions:
To solve the above problems, the author proposes a new and efficient numerical algorithm with the following characteristics:
- **Relaxing the regularity requirements of initial data**: The new method can still maintain good convergence when the initial data belongs to \(H^\gamma\times H^{\gamma - 1}\) (for all \(\gamma>0\)).
- **Improved error rate**: In different spatial dimensions, the new method can achieve a higher convergence order. For example, in one - dimensional and two - dimensional cases, for \(\gamma\in(0,\frac{1}{2}]\), the error rate is \(O(\tau^{2\gamma})\); in three - dimensional cases, for \(\gamma\in(0,\frac{3}{4}]\), the error rate is \(O(\tau^{\max(\gamma,2\gamma-\frac{1}{2})})\).
- **Using the solution operator of the linear wave equation**: By introducing the solution operator \(e^{(t - t_n)L}\) of the linear wave equation, the new method can estimate the error more accurately and does not need to rely on the Hölder continuity in the time direction of the solution.
#### Verification and advantages:
To verify the effectiveness of the new method, the author conducted extensive numerical experiments. The results show that the new method has higher accuracy and computational efficiency compared with existing methods, especially when dealing with low - regularity or discontinuous initial conditions.
In summary, the main goal of this paper is to develop a new numerical algorithm that can effectively solve stochastic nonlinear wave equations under low - regularity conditions, thereby improving the accuracy and efficiency of numerical simulations.