Matrix-Free Parallel Preconditioned Iterative Solvers for the 2D Helmholtz Equation Discretized with Finite Differences

We present a matrix-free parallel iterative solver for the Helmholtz equation related to applications in seismic problems and study its parallel performance. We apply Krylov subspace methods, GMRES, Bi-CGSTAB and IDR(s), to solve the linear system obtained from a second-order finite difference discretization. The Complex Shifted Laplace...

Hybrid-dimensional modeling for fluid flow in heterogeneous porous media using dual fracture-pore model with flux interaction of fracture–cavity network

Developing numerical method of fractured porous media is of paramount importance in geoscience applications. Previous studies have revealed that the discrete fractures and cavities as well as the heterogeneity have considerable influences on hydraulic property of porous media. This work presents a numerical investigation on fluid flow in...

Investigating Effects of Heterogeneity and Fracture Distribution on Two-Phase Flow in Fractured Reservoir with adaptive time strategy

Modeling of fluid flow in porous media is a pillar in geoscience applications. Previous studies have revealed that heterogeneity and fracture distribution have considerable influence on fluid flow. In this work, a numerical investigation of two-phase flow in heterogeneous fractured reservoir is presented. First, the discrete fracture model is...

Composite Anderson acceleration method with two window sizes and optimized damping

In this article, we propose and analyze a set of fully nonstationary Anderson acceleration (AA) algorithms with two window sizes and optimized damping. Although AA has been used for decades to speed up nonlinear solvers in many applications, most authors are simply using and analyzing the stationary version of AA (sAA) with fixed window size...

On an improved PDE-based elliptic parameterization method for isogeometric analysis using preconditioned Anderson acceleration

Constructing an analysis-suitable parameterization for the computational domain from its boundary representation plays a crucial role in the isogeometric design-through-analysis pipeline. PDE-based elliptic grid generation is an effective method for generating high-quality parameterizations with rapid convergence properties for the planar...

Asymptotic convergence analysis and influence of initial guesses on composite Anderson acceleration

Although Anderson acceleration AA(m) has been widely used to speed up nonlinear solvers, most authors are simply using and studying the stationary version of Anderson acceleration. The behavior and full potential of the non-stationary version of Anderson acceleration methods remain an open question. Motivated by the hybrid linear solver...

A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation

The Helmholtz equation is related to seismic exploration, sonar, antennas, and medical imaging applications. It is one of the most challenging problems to solve in terms of accuracy and convergence due to the scalability issues of the numerical solvers. For 3D large-scale applications, high-performance parallel solvers are also needed. In this...

An efficient adaptive implicit scheme with equivalent continuum approach for two-phase flow in fractured vuggy porous media

This work investigates numerical method and equivalent continuum approach (ECA) of fluid flow in fractured porous media. The commonly used discrete fracture model (DFM) without upscaling needs full discretization of all fractures. It enjoys the merit of capturing each fracture accurately but will get in trouble with mesh partition and low...

Globalization technique for projected Newton–Krylov methods

Large-scale systems of nonlinear equations appear in many applications. In various applications, the solution of the nonlinear equations should also be in a certain interval. A typical application is a discretized system of reaction diffusion equations. It is well known that chemical species should be positive otherwise the solution is not...