间断问题的并行多重网格预条件共轭梯度法求解

首先比较了四种串行多重网格预条件共轭梯度法求解四类典型间断问题时,数值收敛因子与间断类型、间断程度、网格步长的关系;然后,在算法中引入Schwarz并行,探讨了并行算法收敛因子与子区域个数的关系,以及Schwarz并行的获益。给出了详细的性能比较结果。     

三维问题的局部块分解预条件

利用块三对角矩阵的嵌套局部块分解构造了一个不完全分解预条件子,并考虑了其修正型变种,分析了两者的存在性及若干性质.针对标准七点差分矩阵,给出了预条件后的实际条件数.结果表明,采用局部块分解预条件时条件数与矩阵阶数的2/3次幂成正比,而采用修正型预条件时条件数与矩阵阶数的立方根成正比.最后考虑了预条件的高效实现并在主频为550MHz、内存为256M的微机上作了若干数值实验,并与其它较有效的预条件方法进行了比较.     

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.     

Preconditioning the Poincaré-Steklov operator by using Green's function

This paper is concerned with the Poincaré-Steklov operator that is widely used in domain decomposition methods. It is proved that the inverse of the Poincaré-Steklov operator can be expressed explicitly by an integral operator with a kernel being the Green's function restricted to the interface. As an application, for the discrete Poincaré-Steklov operator with respect to either a line (edge) or a star-shaped web associated with a single vertex point, a preconditioner can be constructed by first imbedding the line as the diameter of a disk, or the web as a union of radii of a disk, and then using the Green's function on the disk. The proposed technique can be effectively used in conjunction with various existing domain decomposition techniques, especially with the methods based on vertex spaces (from multi-subdomain decomposition). Some numerical results are reported.

     

Iterative proximal regularization of the modified Lagrangian functional for solving the quasi-variational Signorini inequality

The iterative Uzawa method with a modified Lagrangian functional is used to numerically solve the semicoercive Signorini problem with friction (quasi-variational inequality).     

On the rotating Navier-Stokes equations with mixed boundary conditions

The stationary and nonstationary rotating Navier-Stokes equations with mixed boundary conditions are investigated in this paper. The existence and uniqueness of the solutions are obtained by the Galerkin approximation method. Next, θ-scheme of operator splitting algorithm is applied to rotating Navier-Stokes equations and two subproblems are derived. Finally, the computational algorithms for these subproblems are provided.     

Parameterized (in)approximability of subset problems

We discuss approximability and inapproximability in FPT-time for a large class of subset problems where a feasible solution $S$ is a subset of the input data. We introduce the notion of intersective approximability that generalizes the one of safe approximability introduced in Guo et al. (2011) and show strong parameterized inapproximability results for many of the subset problems handled.     

椭圆PDE-约束优化问题的一个预条件子

针对由Galerkin有限元离散椭圆PDE-约束优化问题产生的具有特殊结构的3×3块线性鞍点系统,提出了一个预条件子并给出了预处理矩阵特征值及特征向量的具体表达形式.数值结果表明了该预条件子能够有效地加速Krylov子空间方法的收敛速率,同时也验证了理论结果.     

Block ILU Preconditioners for a Nonsymmetric Block-Tridiagonal M-Matrix

We propose block ILU (incomplete LU) factorization preconditioners for a nonsymmetric block-tridiagonal M-matrix whose computation can be done in parallel based on matrix blocks. Some theoretical properties for these block ILU factorization preconditioners are studied and then we describe how to construct them effectively for a special type of matrix. We also discuss a parallelization of the preconditioner solver step used in nonstationary iterative methods with the block ILU preconditioners. Numerical results of the right preconditioned BiCGSTAB method using the block ILU preconditioners are compared with those of the right preconditioned BiCGSTAB using a standard ILU factorization preconditioner to see how effective the block ILU preconditioners are.     

Regularized DPSS Preconditioners for Singular Saddle Point Problems

Recently, Cao proposed a regularized deteriorated positive and skew-Hermitian splitting (RDPSS) preconditioner for the non-Hermitian nonsingular saddle point problem. In this paper, we consider applying RDPSS preconditioner to solve the singular saddle point problem. Moreover, we propose a two-parameter accelerated variant of the RDPSS (ARDPSS) preconditioner to further improve its efficiency. Theoretical analysis proves that the RDPSS and ARDPSS methods are semi-convergent unconditionally. Some spectral properties of the corresponding preconditioned matrices are analyzed. Numerical experiments indicate that better performance can be achieved when applying the ARDPSS preconditioner to accelerate the GMRES method for solving the singular saddle point problem.