A numerical experimental study of inverse preconditioning for the parallel iterative solution to 3D finite element flow equations

Integration of the subsurface flow equation by finite elements (FE) in space and finite differences (FD) in time requires the repeated solution to sparse symmetric positive definite systems of linear equations. Iterative techniques based on preconditioned conjugate gradients (PCG) are one of the most attractive tool to solve the problem on sequential computers. A present challenge is to make PCG attractive in a parallel computing environment as well. To this aim a key factor is the development of an efficient parallel preconditioner. FSAI (factorized sparse approximate inverse) and enlarged FSAI relying on the approximate inverse of the coefficient matrix appears to be a most promising parallel preconditioner. In the present paper PCG using FSAI, diagonal and pARMS (parallel algebraic recursive multilevel solvers) preconditioners is implemented on the IBM SP4/512 and CLX/768 supercomputers with up to 32 processors to solve underground flow problems of a large size. The results show that FSAI may allow for a parallel relative efficiency larger than 50% on the largest problems with p=32 processors. Moreover, FSAI turns out to be significantly less expensive and more robust than pARMS. Finally, it is shown that for p in the upper range may be much improved if PCG–FSAI is implemented on CLX.     

等腰三角形Mindlin板的自由振动分析

提出了一种新方法来对基于 Mindlin剪切变形理论的等腰三角形板进行自由振动分析 .此方法采用了一种新的基函数并利用 pb-2 Rayleigh-Ritz边界函数得到了一种新型的 Ritz方法 .这种方法的有效性通过收敛性和对比性分析得到了证实 .数值结果表明此方法相当精确有效 .     

Finite element modal solutions of planar photonic crystal fibers with rectangular air-holes

Modal solutions of planar photonic crystal waveguides with rectangular air-holes are presented by using a rigorous full-vectorial finite element-based approach. The effective indices, mode field profiles, spot-sizes, power confinements, modal hybridness, beat lengths and group velocity dispersions are shown for the fundamental and higher order modes of the quasi-TE and TM polarizations.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

On the Modified Crystalline Stefan Problem with Singular Data

We study the behavior of solutions of the modified Stefan problem in the plane for polygonal interfaces. We are particularly interested in a solution near a singularity of either the loss of a facet or the breaking of a facet. We establish precise regularity results if a facet disappears. We use them to establish the existence of a weak solution with singular data, i.e., when some of the zero-crystalline-curvature facets have zero length.     

On Pairs of Diagonal Quintic Forms

We demonstrate that a pair of additive quintic equations in at least 34 variables has a nontrivial integral solution, subject only to an 11-adic solubility hypothesis. This is achieved by an application of the Hardy–Littlewood method, for which we require a sharp estimate for a 33.998th moment of quintic exponential sums. We are able to employ p-adic iteration in a form that allows the estimation of such a mean value over a complete unit square, thereby providing an approach that is technically simpler than those of previous workers and flexible enough to be applied to related problems.     

Optimal Investment in a Levy Market

In this paper we consider the optimal investment problem in a market where the stock price process is modeled by a geometric Levy process (taking into account jumps). Except for the geometric Brownian model and the geometric Poissonian model, the resulting models are incomplete and there are many equivalent martingale measures. However, the model can be completed by the so-called power-jump assets. By doing this we allow investment in these new assets and we can try to maximize the expected utility of these portfolios. As particular cases we obtain the optimal portfolios based in stocks and bonds, showing that the new assets are superfluous for certain martingale measures that depend on the utility function we use.     

A Verified Optimization Technique to Locate Chaotic Regions of Hénon Systems

We present a new verified optimization method to find regions for Hénon systems where the conditions of chaotic behaviour hold. The present paper provides a methodology to verify chaos for certain mappings and regions. We discuss first how to check the set theoretical conditions of a respective theorem in a reliable way by computer programs. Then we introduce optimization problems that provide a model to locate chaotic regions. We prove the correctness of the underlying checking algorithms and the optimization model. We have verified an earlier published chaotic region, and we also give new chaotic places located by the new technique.     

中物院远红外自由电子激光实验研究

中国工程物理研究院基于射频直线加速器技术的远红外自由电子激光(FIR-FEL)实验取得阶段性进展,于2005年3月24晚8时30分首次出光, 并多次重复. 中心波长115μm, 谱宽1%.介绍了实验系统的主要组成部分和主要实验结果.     

牛顿弦截法预估校正迭代格式的收敛阶

研究如下形式的牛顿弦截法的预估校正(P.C.)格式:P(预估):~xk+1=xk-(xk-xk-1)f(xk)f(xk)-f(xk-1)C(校正):xk+1=xk-(~xk+1-xk)f(xk)f~(xk+1)-f(xk)证明了它的收敛阶为2.618.