Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

线性互补问题的并行多分裂松弛迭代算法

运用矩阵多重分裂理论,同时考虑并行计算与松弛迭代法,得到一类求解线性互补问题的高效数值算法．当问题的系数矩阵为对角元为正的H-矩阵或对称半正定矩阵时,证明了算法的全局收敛性；该算法与已有算法相比,具有计算量小、计算速度快等特点,因而特别适于求解大规模问题．数值试验的结果说明了算法的有效性．     

�ɲ����˶�����ά��������������һ�������ڵ����Զ������������������

??We study the linear quadratic optimal stochastic control problem which is jointly driven by Brownian motion and L\'{e}vy processes. We prove that the new affine stochastic differential adjoint equation exists an inverse process by applying the profound section theorem. Applying for the Bellman's principle of quasilinearization and a monotone iterative convergence method, we prove the existence and uniqueness of the solution of the backward Riccati differential equation. Finally, we prove that the optimal feedback control exists, and the value function is composed of the initial value of the solution of the related backward Riccati differential equation and the related adjoint equation.     

Iterative solutions for low lying excited states of a class of SchrSdinger equation

The convergent iterative procedure for solving the groundstate Schr?dinger equation is extended to derive the excitation energy and the wavefunction of the low-lying excited states. The method is applied to the one-dimensional quartic potential problem. The results show that the iterative solution converges rapidly when the coupling g is not too small.     

响应矩阵分析方法在BEPC储存环上的应用

由于储存环中各种元件误差的存在, 机器的实际运行模式与设计模式有一定的偏差. 目前广泛开展的响应矩阵方法研究, 可以分析出磁铁元件以及束流位置测量元件的误差, 使束流基本参数得到校正. 介绍了用响应矩阵分析方法, 在BEPC储存环上进行的局部轨道校正的实验研究, 以及BEPC储存环束流参数校正的模拟研究.     

一类二阶两点边值问题的单调迭代方法

通过改进经典的单调迭代方法对于一类二阶两点边值的问题的正解建立了单调迭代程序。这些迭代程序都是从常值函数开始的，因而是可行并且有效的。     

X光底片在位相对标定

本文描述了X光底片在位相对标定技术。其原理是使X射线谱经阶梯形吸收滤片透射后,对X光底片曝光,测量底片的曝光量。文中给出了标定方法和数据处理方法,而且也给出了在X光激光实验中得到的Kodak AA5底片的特性曲线。     

用最小电离的质子刻度铅–闪烁光纤夹层电磁量能器

研究利用最小电离的质子和已知能量的电子来刻度铅–闪烁光纤夹层电磁量能器.对最小电离粒子输出响应进行光纤衰减修正后,利用最小电离质子的输出响应对量能器测量单元进行标定,然后对不同入射能量的电子事例进行能量重建,重建能量与入射能量有很好的线性关系.本文还对最小电离的质子与μ子对量能器测量单元的输出响应进行了比较,发现二者在误差范围内是一致的.