Hanyu Li,Wing Tat Leung,Mary F. Wheeler
Abstract:Convergence failure and slow convergence rates are among the biggest challenges with solving the system of non-linear equations numerically. Although mitigated, such issues still linger when using strictly small time steps and unconditionally stable fully implicit schemes. The price that comes with restricting time steps to small scales is the enormous computational load, especially in large-scale models. To address this problem, we introduce a sequential local mesh refinement framework to optimize convergence rate and prevent convergence failure, while not restricting the whole system to small time steps, thus improving computational efficiency. We test the algorithm with the non-linear two-phase flow model. Starting by solving the problem on the coarsest space-time mesh, the algorithm refines the domain in time at water saturation front to prevent convergence failure. It then deploys fine spatial grid in regions with large saturation variation to assure solution accuracy. After each refinement, the solution from the previous mesh is used to estimate the initial guess of unknowns on the current mesh for faster convergence. Numerical results are presented to confirm accuracy of our algorithm as compared to the uniformly fine time step and fine spatial discretization solution. We observe approximately 25 times speedup in the solution time by using our algorithm.
What problem does this paper attempt to address?
This paper attempts to solve the problems of convergence failure and slow convergence speed in nonlinear two - phase flow problems. Specifically, the article proposes a sequential local mesh refinement solver to optimize the convergence speed and prevent convergence failure, without restricting the entire system to small time steps, thereby improving computational efficiency.
### Problem Background
When numerically solving nonlinear equation systems, convergence failure and slow convergence speed are among the greatest challenges. Even when using strictly small time steps and unconditionally stable fully - implicit schemes, these problems still exist. Restricting the time step to a very small scale can alleviate these problems, but it will bring a huge computational burden, especially in large - scale models.
### Solution
To solve the above - mentioned problems, the author introduced a Sequential Local Mesh Refinement framework, which aims to:
1. **Optimize the Convergence Rate**: Optimize the convergence speed of the Newton method by locally refining the mesh.
2. **Prevent Convergence Failure**: Prevent convergence failure by performing temporal refinement at the water saturation front.
3. **Improve Computational Efficiency**: Avoid using overly small time steps in the entire system, thereby reducing the amount of computation.
### Specific Methods
- **Initial Coarse - grid Solution**: Start solving the problem from the coarsest space - time grid.
- **Temporal Refinement**: Perform temporal refinement at the water saturation front to prevent convergence failure.
- **Spatial Refinement**: Perform spatial refinement in areas with large saturation changes to ensure the accuracy of the solution.
- **Initial Guess Estimation**: After each refinement, use the solution of the previous grid as the initial guess for the unknowns of the current grid to accelerate the convergence speed.
### Results
Numerical results show that compared with uniform fine time steps and fine - space discretization, this algorithm improves the solution speed by about 25 times and maintains high accuracy.
### Keywords
Space - time domain decomposition, mixed finite element method, sequential local mesh refinement, iterative solver, nonlinear system
Through this method, the author has successfully solved the convergence and computational efficiency problems in nonlinear two - phase flow problems, especially showing significant advantages in large - scale models.