High‐order implicit time integration for unsteady incompressible flows

The spatial discretization of unsteady incompressible Navier–Stokes equations is stated as a system of differential algebraic equations, corresponding to the conservation of momentum equation plus the constraint due to the incompressibility condition. Asymptotic stability of Runge–Kutta and Rosenbrock methods applied to the solution of the resulting index‐2 differential algebraic equations system is analyzed. A critical comparison of Rosenbrock, semi‐implicit, and fully implicit Runge–Kutta methods is performed in terms of order of convergence and stability. Numerical examples, considering a discontinuous Galerkin formulation with piecewise solenoidal approximation, demonstrate the applicability of the approaches and compare their performance with classical methods for incompressible flows. Copyright © 2011 John Wiley & Sons, Ltd.     

Mixed‐order interpolation for the Galerkin coarse‐grid approximations in algebraic multigrid solvers

An empirical investigation is made of AMG solver performance for the fully coupled set of Navier–Stokes equations. The investigation focuses on two different FV discretizations for the standard driven cavity test problem. One is a collocated vertex‐based discretization; the other is a cell‐centred staggered‐grid discretization. Both employ otherwise identical orthogonal Cartesian meshes. It is found that if mixed‐order interpolation is used in the construction of the Galerkin coarse‐grid approximation (CGA), a close‐to‐optimum mesh‐independent scaling of the AMG convergence is observed with similar convergence rates for both discretizations. If, on the other hand, an equal‐order interpolation is used, convergence rates are mesh‐dependent but the scaling differs in each case. For the collocated‐grid case, it depends both on the mesh size, h (or bandwidth Q～h^?1) and on the total number of grids, G, whereas for the staggered‐grid case it depends only on Q. Comparing the two characteristics reveals that the Q‐dependent parts are very similar; it is only in the G‐dependent convergence for the collocated‐grid case that they differ. This takes the form of stepped reductions in the AMG convergence rate (implying step reductions in the quality of the Galerkin CGA that correlate exactly with step increases in G). These findings reinforce previous evidence that, for optimum mesh‐independent performance, mixed‐order interpolations should be used in forming Galerkin CGAs for coupled Navier–Stokes problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of an aggregation‐based algebraic two‐grid method for a rotated anisotropic diffusion problem

A two‐grid convergence analysis based on the paper [Algebraic analysis of aggregation‐based multigrid, by A. Napov and Y. Notay, Numer. Lin. Alg. Appl. 18 (2011), pp. 539–564] is derived for various aggregation schemes applied to a finite element discretization of a rotated anisotropic diffusion equation. As expected, it is shown that the best aggregation scheme is one in which aggregates are aligned with the anisotropy. In practice, however, this is not what automatic aggregation procedures do. We suggest approaches for determining appropriate aggregates based on eigenvectors associated with small eigenvalues of a block splitting matrix or based on minimizing a quantity related to the spectral radius of the iteration matrix. Copyright © 2015 John Wiley & Sons, Ltd.     

一类幂等半环的性质和结构

研究了加法半群为半格的半环类S+l中的乘法带半环和矩形带半环类BR中的乘法带半环;给出了ID半环中乘法带半环的结构定理,即ID∩.■°D=.■z∨.■z∨D.     

NO分子宏观气体热力学性质的理论研究

采用量子统计系综理论,研究了基态NO分子宏观气体摩尔熵、摩尔内能、摩尔热容等热力学性质.首先应用课题组前期建立的变分代数法(variational algebraic method, VAM)计算获得了基态NO分子的完全振动能级,得到的VAM振动能级作为振动部分,结合欧拉-麦克劳林渐进展开公式的转动贡献,应用于经典的热力学与统计物理公式中,从而计算得到了1000-5000 K温度范围内NO宏观气体的摩尔内能、摩尔熵和摩尔热容.将不同方法计算得到的摩尔热容结果分别与实验值进行比较,结果表明基于VAM完全振动能级获得的结果优于其他方法获得的理论结果.振动部分采用谐振子模型对无限能级求和计算热力学性质的方法有一定的局限性,应当使用有限的完全振动能级进行统计求和.     

ALGEBRAIC EXTENSION OF ＊-A OPERATOR

In this paper, we study various properties of algebraic extension of *-A operator.Specifically, we show that every algebraic extension of *-A operator has SVEP and is isoloid.And if T is an algebraic extension of *-A operator, then Weyl's theorem holds for f(T), where f is an analytic functions on some neighborhood of σ(T) and not constant on each of the components of its domain.     

一类随机延迟微分代数系统的Euler-Maruyama方法

我们主要构造了数值求解一类1指标随机延迟微分代数系统的Euler-Maruyama方法,并且证明用该方法求解此类问题可达到1/2阶均方收敛.最后的效值试验验证了方法的有效性及所获结论的正确性.     

解强刚性块线代数方程组的两类L-收敛迭代法

对解强刚性块线代数方程组X=(A(?)J)X φ,本文提出了L-收敛的最佳单参数迭代法(L-OOPI)和L-收敛的多参数迭代直接法(L-MPID),并给出了数值例子．数例表明,对于强刚性块线代数方程组,该二迭代法是有效的．     

人造地球卫星运动方程的代数动力学算法数值解

代数动力学算法首次被用于求人造地球卫星运动方程的数值解，在四阶算法下，与Runge-Kutta算法和辛算法的计算结果作了比较.结果表明，代数动力学算法对于人造地球卫星长期轨道的计算有较高的精度.并讨论了地球四极和八极带谐项对卫星轨道的影响. 关键词： 人造地球卫星运动方程 代数动力学算法数值解 地球四极和八极带谐项对卫星轨道的影响