Abstract:This paper aims to present a local discontinuous Galerkin (LDG) method for solving backward stochastic partial differential equations (BSPDEs) with Neumann boundary conditions. We establish the $L^2$-stability and optimal error estimates of the proposed numerical scheme. Two numerical examples are provided to demonstrate the performance of the LDG method, where we incorporate a deep learning algorithm to address the challenge of the curse of dimensionality in backward stochastic differential equations (BSDEs). The results show the effectiveness and accuracy of the LDG method in tackling BSPDEs with Neumann boundary conditions.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to provide an effective numerical solution method for backward stochastic partial differential equations (BSPDEs) with Neumann boundary conditions. Specifically, the paper proposes a local discontinuous Galerkin (LDG) method to solve such equations and combines deep - learning algorithms to address the curse - of - dimensionality challenge in high - dimensional problems.
### Problem Background
1. **Research Object**: Backward stochastic partial differential equations (BSPDEs) with Neumann boundary conditions, which are in the form of:
\[
\begin{cases}
-du(x,t) = \left[\frac{1}{2}(|\sigma|^2 + |\bar{\sigma}|^2) u_{xx} + \sigma \psi_x + \Gamma(x,t,u,u_x,\psi)\right] dt - \psi(x,t) dW_t, & (x,t) \in [0,b] \times [0,T), \\
u_x(0,t) = g(0,t), \quad u_x(b,t) = g(b,t), & t \in [0,T], \\
u(x,T) = G(x), & x \in [0,b].
\end{cases}
\]
2. **Limitations of Existing Methods**:
- For BSPDEs, existing work mainly focuses on problems in the full space and under Dirichlet boundary conditions, while there is less research on Neumann boundary conditions.
- In high - dimensional cases, traditional numerical methods (such as the finite element method) face the curse - of - dimensionality problem, that is, the computational complexity rises sharply as the dimension increases.
### Proposed Method
1. **Local Discontinuous Galerkin (LDG) Method**:
- The LDG method is a numerical method suitable for solving nonlinear hyperbolic systems. It has the characteristics of being easy to design high - precision algorithms, adapting to different geometric shapes, and being suitable for efficient parallel computing.
- By discretizing BSPDEs into high - dimensional BSDEs and combining deep - learning algorithms (such as the deep BSDE method), the curse - of - dimensionality in high - dimensional problems can be effectively addressed.
2. **Application of Deep - Learning Algorithms**:
- Deep - learning algorithms are used to solve high - dimensional BSDEs. In particular, neural - network - based methods can transform PDEs into equivalent forward - backward stochastic differential equation systems (FBSDEs) and perform numerical solutions through neural networks.
- The deep backward dynamic programming method is introduced. Through the backward discretization scheme, the single - step dynamic programming principle is derived, the training tasks are decomposed, and the original problem is iteratively solved, thereby improving the algorithm stability.
### Research Results
1. **Stability Analysis**: The $L^2$ stability of the LDG method is established.
2. **Error Estimation**: The optimal error estimate is provided, and the existence and uniqueness of the numerical solution are proved.
3. **Numerical Experiments**: The effectiveness and accuracy of the LDG method in solving BSPDEs with Neumann boundary conditions are verified through two numerical examples.
In summary, this paper aims to provide an efficient and accurate numerical solution method for BSPDEs with Neumann boundary conditions by introducing the LDG method and deep - learning algorithms, especially suitable for high - dimensional problems.