Convergence of nonconforming multigrid methods without full elliptic regularity

We consider nonconforming multigrid methods for symmetric positive definite second and fourth order elliptic boundary value problems which do not have full elliptic regularity. We prove that there is a bound () for the contraction number of the -cycle algorithm which is independent of mesh level, provided that the number of smoothing steps is sufficiently large. We also show that the symmetric variable -cycle algorithm is an optimal preconditioner.

     

抛物问题的拉格朗日乘子区域分解法

１．引言 近年来,一类新的非重叠区域分解方法一非匹配网格区域分解法,日益引起人们的广泛兴趣,并已成为当今区域分解方法研究的热门课题。这类区域分解方法的特点是：相邻子区域在公共边（或面）上的结点可以不重合,从而能解决许多传统区域分解方法不便解决的问题（如变动网格问题）．目前主要有两类方法来处理这种区域分解的强非协调性：Ｍｏｒｔａｒ无法（见［１－２］和［９－１０］）和拉格朗日乘子法（见［５］,［８］,［１１］和［１２］）．拉格朗日乘子法比Ｍｏｒｔａｒ无法有明显的优点：（１）界面变量（即拉格朗日乘子）…     

Solution of large‐scale weighted least‐squares problems

A sequence of least‐squares problems of the form min_y∥G^1/2(A^T y?h)∥₂, where G is an n×n positive‐definite diagonal weight matrix, and A an m×n (m?n) sparse matrix with some dense columns; has many applications in linear programming, electrical networks, elliptic boundary value problems, and structural analysis. We suggest low‐rank correction preconditioners for such problems, and a mixed solver (a combination of a direct solver and an iterative solver). The numerical results show that our technique for selecting the low‐rank correction matrix is very effective. Copyright © 2002 John Wiley & Sons, Ltd.     

An overlapping additive Schwarz preconditioner for interpolation on the unit sphere with spherical radial basis functions

The problem of interpolation of scattered data on the unit sphere has many applications in geodesy and Earth science in which the sphere is taken as a model for the Earth. Spherical radial basis functions provide a convenient tool for constructing the interpolant. However, the underlying linear systems tend to be ill-conditioned. In this paper, we present an additive Schwarz preconditioner for accelerating the solution process. An estimate for the condition number of the preconditioned system will be discussed. Numerical experiments using MAGSAT satellite data will be presented.     

On the optimal convergence factor of the accelerated parameterized inexact Uzawa method with three parameters for augmented systems

Under the assumption that all eigenvalues of the preconditioned Schur complement are real, we present an analytical proof for obtaining the optimal convergence factor of the real accelerated parameterized inexact Uzawa (APIU) method when P=A. It is proved that the optimal convergence factor is the same as that of the generalized successive overrelaxation method, which was published at the same time, and that it can be attained only at the unique optimum point of parameters, regardless of whether m>n or m=n. In addition, we generalize the APIU method and analyze the relationship between the APIU method and 10 additional Uzawa‐like methods.     

A note on spectrum distribution of constraint preconditioned generalized saddle point matrices

In this note, results concerning the eigenvalue distribution and form of the eigenvectors of the constraint preconditioned generalized saddle point matrix and its minimal polynomial are given. These results extend previous ones that appeared in the literature. Copyright © 2009 John Wiley & Sons, Ltd.     

An Uzawa-type algorithm for the coupled Stokes equations

An Uzawa-type algorithm is designed for the coupled Stokes equations discretized by the mixed finite element method. The velocity solved by the presented algorithm is weakly divergence-free, which is different from the one solved by the common Uzawa method. Besides, an optimal relaxation parameter of the presented algorithm is provided.     

A characteristic nonoverlapping domain decomposition method for multidimensional convection‐diffusion equations

We develop a quasi‐two‐level, coarse‐mesh‐free characteristic nonoverlapping domain decomposition method for unsteady‐state convection‐diffusion partial differential equations in multidimensional spaces. The development of the domain decomposition method is carried out by utilizing an additive Schwarz domain decomposition preconditioner, by using an Eulerian‐Lagrangian method for convection‐diffusion equations and by delicately choosing appropriate interface conditions that fully respect and utilize the hyperbolic nature of the governing equations. Numerical experiments are presented to illustrate the method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

An improvement to the Kunz preconditioner and numerical investigation of hydrofoil interactions in tandem

ABSTRACT

Numerical simulations of the tandem-system flow are widely conducted because of the complex interactions of the configuration. The cavitation phenomenon is an important factor that affects the hydrofoil in tandem. In this paper, we developed a new parameter for the Kunz preconditioner based on the local cavitation volume fraction. To assess this parameter, a type of Fourier footprint analysis and numerical test of the hydrofoil are adopted. The preconditioning method is also conducted in hydrofoil turbulent cavitation flows with RANS turbulence models, to prove that this method has good stability and convergence. Based on this, a tandem Clark-Y hydrofoil configuration is investigated. The results show that the distance of components has a strong influence on the cavitation and lift coefficient of the tandem hydrofoils. Therefore, the purpose of this study is to provide guidance on the design of marine vehicles that have tandem configurations.