Meshfree Generalized Multiscale Exponential Integration Method for Parabolic Problems

Djulustan Nikiforov,Leonardo A. Poveda,Dmitry Ammosov,Yesy Sarmiento,Juan Galvis

2024-10-15

Abstract:This paper considers flow problems in multiscale heterogeneous porous media. The multiscale nature of the modeled process significantly complicates numerical simulations due to the need to compute huge and ill-conditioned sparse matrices, which negatively affect both the computational cost and the stability of the numerical solution. We propose a novel combined approach of the meshfree Generalized Multiscale Finite Element Method (MFGMsFEM) and exponential time integration for solving such problems. MFGMsFEM provides a robust and efficient spatial approximation, allowing us to consider complex heterogeneities without constructing a coarse computational grid. At the same time, exponential integration, using the cost-effective MFGMsFEM matrix, provides a robust temporal approximation for stiff multiscale problems, allowing larger time steps. For the proposed multiscale approach, we provide a rigorous convergence analysis, including the new analysis of the MFGMsFEM spatial approximation. We conduct numerical experiments to computationally verify the proposed approach by solving linear and semi-linear flow problems in multiscale media. Numerical results demonstrate that the proposed multiscale method achieves significant reductions in computational cost and improved stability, even with larger time steps, confirming the theoretical analysis.

Numerical Analysis

What problem does this paper attempt to address?

The paper attempts to address the issues encountered in flow simulation within multi-scale heterogeneous porous media. Specifically, these problems become numerically complex due to the need to handle large-scale and ill-conditioned sparse matrices, which affects computational cost and the stability of numerical solutions.

The paper proposes a novel method that combines the Meshfree Generalized Multiscale Finite Element Method (MFGMsFEM) and the exponential time integration method. This approach addresses the issues through the following aspects:

1. **Spatial Approximation**: MFGMsFEM provides a robust and efficient spatial approximation method, allowing for the consideration of complex heterogeneities without the need to construct a coarse grid.
2. **Temporal Approximation**: Utilizing the efficient matrices of MFGMsFEM, the exponential integration method can offer robust temporal approximation for stiff multi-scale problems and allows for the use of larger time steps.

Through theoretical analysis and numerical experiments, the method significantly reduces computational cost and improves stability, maintaining good performance even with larger time steps. This indicates that the method has high potential in handling problems with high-contrast multi-scale coefficients.