适合于求解边界元方程组的GMRES算法的实用化和并行化研究

为了将GMRES算法应用于大型边界元方程组的求解,采用预条件技术和重正交技术相结合的方法实现了该算法的实用化,然后在实用化的基础上针对迭代算法具有良好并行性的特点,研究了该算法在网络机群环境下的并行化技术。数值试验和分析表明所用的这些技术是行之有效的,对于提高求解速度和增大求解问题的规模是有意义的。     

A Fast,Matrix-free Implicit Method for Computing Low Mach Number Flows on Unstructured Grids

A fast, matrix-free implicit method has been developed to solve low Mach number flow problems on unstructured grids. The preconditioned compressible Euler and Navier-Stokes equations are integrated in time using a linearized implicit scheme. A newly developed fast, matrix-free implicit method, GMRES + LU?SGS, is then applied to solve the resultant system of linear equations. A variety of computations has been made for a wide range of flow conditions, for both in viscid and viscous flows, in both 2D and 3D to validate the developed method and to evaluate the effectiveness of the GMRES + LU?SGS method. The numerical results obtained indicate that the use of the GMRES + LU?SGS method leads to a significant increase in performance over the LU?SGS method, while maintaining memory requirements similar to its explicit counterpart. An overall speedup factor from one to more than two order of magnitude for all test cases in comparison with the explicit method is demonstrated.     

Comparison of GMRES and ORTHOMIN for the black oil model on unstructured grids

This paper addresses an application of ORTHOMIN and GMRES to petroleum reservoir simulation using the black oil model on unstructured grids. Comparisons between these two algorithms are presented in terms of storage and total flops per restart step. Numerical results indicate that GMRES is faster than ORTHOMIN for all tested petroleum reservoir problems, particularly for large scale problems. The control volume function approximation method is utilized in the discretization of the governing equations of the black oil model. This method can accurately approximate both the pressure and velocity in the simulation of multiphase flow in porous media, effectively reduce grid orientation effects, and be easily applied to arbitrarily shaped control volumes. It is particularly suitable for hybrid grid reservoir simulation. Copyright © 2006 John Wiley & Sons, Ltd.     

GMRES and the minimal polynomial

We present a qualitative model for the convergence behaviour of the Generalised Minimal Residual (GMRES) method for solving nonsingular systems of linear equationsAx =b in finite and infinite dimensional spaces. One application of our methods is the solution of discretised infinite dimensional problems, such as integral equations, where the constants in the asymptotic bounds are independent of the mesh size.Our model provides simple, general bounds that explain the convergence of GMRES as follows: If the eigenvalues ofA consist of a single cluster plus outliers then the convergence factor is bounded by the cluster radius, while the asymptotic error constant reflects the non-normality ofA and the distance of the outliers from the cluster. If the eigenvalues ofA consist of several close clusters, then GMRES treats the clusters as a single big cluster, and the convergence factor is the radius of this big cluster. We exhibit matrices for which these bounds are tight.Our bounds also lead to a simpler proof of existing r-superlinear convergence results in Hilbert space.This research was partially supported by National Science Foundation grants DMS-9122745, DMS-9423705, CCR-9102853, CCR-9400921, DMS-9321938, DMS-9020915, and DMS-9403224.     

A simpler GMRES

The generalized minimal residual (GMRES) method is widely used for solving very large, nonsymmetric linear systems, particularly those that arise through discretization of continuous mathematical models in science and engineering. By shifting the Arnoldi process to begin with Ar₀ instead of r₀, we obtain simpler Gram–Schmidt and Householder implementations of the GMRES method that do not require upper Hessenberg factorization. The Gram–Schmidt implementation also maintains the residual vector at each iteration, which allows cheaper restarts of GMRES(m) and may otherwise be useful.     

A MODIFIED GMRES METHOD FOR SOLVING LARGE NONSYMMETRIC LINEAR SYSTEMS

A modified GMRES method is proposed in this paper, the method replaces the approximation xm obtained by the GMRES method with a new approximation xm which is a linear combination of xm and the wasted basis vector vm 1. The residual norm of the new approximation satisfies a small one-dimensional minimization problem. Relationships between the residual norms of xm and xm are given. We show that the resulting m-step modified GMRES method is better than the original m-step GMRES method in theory and is consi...     

Weighted FOM and GMRES for solving nonsymmetric linear systems

This paper presents two new methods called WFOM and WGMRES, which are variants of FOM and GMRES, for solving large and sparse nonsymmetric linear systems. To accelerate the convergence, these new methods use a different inner product instead of the Euclidean one. Furthermore, at each restart, a different inner product is chosen. The weighted Arnoldi process is introduced for implementing these methods. After describing the weighted methods, we give the relations that link them to FOM and GMRES. Experimental results are presented to show the good performances of the new methods compared to FOM(m) and GMRES(m). This revised version was published online in August 2006 with corrections to the Cover Date.     

Application of a new fast multipole BEM for simulation of 2D elastic solid with large number of inclusions

A fast multipole method (FMM) is applied for BEM to reduce both the operation and memory requirement in dealing with very large scale problems. In this paper, a new version of fast multipole BEM for 2D elastostatics is presented and used for simulation of 2D elastic solid with a large number of randomly distributed inclusions combined with a similar subregion approach. Generalized minimum residual method (GMRES) is used as an iterative solver to solve the equation system formed by BEM iteratively. The numerical results show that the scheme presented is applicable to certain large scale problems. The project supported by the National Nature Science Foundation of China (10172053) and the Ministry of Education     

On the block GMRES method with deflated restarting

In this paper, we consider the block-GMRES method with deflated restarting for solving nonsymmetric linear systems with multiple right-hand sides. We modify slightly the restarted block-GMRES method with deflation of eigenvalues proposed by Morgan to obtain a new one. It is shown that the modified method is mathematically equivalent to the Morgan’s original method. However, the error analysis shows the modified version minimizes the numerical errors during a restart and therefore is better suited if the linear systems have to be solved with high precision. Numerical experiments report the effectiveness of the modified method.