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.     

弱Orlicz鞅空间的弱原子分解及其应用

对3类由凹函数生成的弱Orlicz鞅空间建立了相应的弱原子分解.作为应用,首先给出了这些弱Orlicz鞅空间上次线性算子有界的一个充分条件,并在此基础上证明了一些弱型鞅不等式,然后证明了关于这些弱Orlicz鞅空间的Marcinkiewicz型插值定理.     

Interpolating solutions of the Poisson equation in the disk based on Radon projections

We consider an algebraic method for reconstruction of a function satisfying the Poisson equation with a polynomial right-hand side in the unit disk. The given data, besides the right-hand side, is assumed to be in the form of a finite number of values of Radon projections of the unknown function. We first homogenize the problem by finding a polynomial which satisfies the given Poisson equation. This leads to an interpolation problem for a harmonic function, which we solve in the space of harmonic polynomials using a previously established method. For the special case where the Radon projections are taken along chords that form a regular convex polygon, we extend the error estimates from the harmonic case to this Poisson problem. Finally we give some numerical examples.     

A new approach to constant term identities and Selberg-type integrals

Selberg-type integrals that can be turned into constant term identities for Laurent polynomials arise naturally in conjunction with random matrix models in statistical mechanics. Built on a recent idea of Karasev and Petrov we develop a general interpolation based method that is powerful enough to establish many such identities in a simple manner. The main consequence is the proof of a conjecture of Forrester related to the Calogero–Sutherland model. In fact we prove a more general theorem, which includes Aomoto's constant term identity at the same time. We also demonstrate the relevance of the method in additive combinatorics.     

Convergence of FFT‐based homogenization for strongly heterogeneous media

The FFT‐based homogenization method of Moulinec–Suquet has recently attracted attention because of its wide range of applicability and short computational time. In this article, we deduce an optimal a priori error estimate for the homogenization method of Moulinec–Suquet, which can be interpreted as a spectral collocation method. Such methods are well‐known to converge for sufficiently smooth coefficients. We extend this result to rough coefficients. More precisely, we prove convergence of the fields involved for Riemann‐integrable coercive coefficients without the need for an a priori regularization. We show that our L² estimates are optimal and extend to mildly nonlinear situations and L^p estimates for p in the vicinity of 2. The results carry over to the case of scalar elliptic and curl ? curl‐type equations, encountered, for instance, in stationary electromagnetism. Copyright © 2014 John Wiley & Sons, Ltd.     

Acceleration of self‐consistent field convergence in ab initio molecular dynamics simulation with multiconfigurational wave function

The Lagrange interpolation of molecular orbital (LIMO) method, which reduces the number of self‐consistent field iterations in ab initio molecular dynamics simulations with the Hartree–Fock method and the Kohn–Sham density functional theories, is extended to the theory of multiconfigurational wave functions. We examine two types of treatments for the active orbitals that are partially occupied. The first treatment, as denoted by LIMO(C), is a simple application of the conventional LIMO method to the union of the inactive core and the active orbitals. The second, as denoted by LIMO(S), separately treats the inactive core and the active orbitals. Numerical tests to compare the two treatments clarify that LIMO(S) is superior to LIMO(C). Further applications of LIMO(S) to various systems demonstrate its effectiveness and robustness. © 2014 Wiley Periodicals, Inc.     

Wavelet multiresolution interpolation Galerkin method for nonlinear boundary value problems with localized steep gradients

The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix. The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities, including transcendental ones, in which the discretization process is as simple as that in solving linear problems, and only common two-term connection coefficients are needed. All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method, which does not require numerical integration in the resulting nonlinear discrete system. The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers. The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids, and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids. In addition, Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method, including the initial guess far from real solutions.

     

构造低次有理插值函数的一种方法

关于有理插值的算法已有很多[1,4,5],受二元多项式插值迭加算法[6]的启发,我们给出一种简便的求低次有理插值函数的方法,同时给出有理插值函数存在的充分条件,便于检验.所给方法具有可操作性和实际应用价值,且具有较好的灵活性.     

优化端点条件的平面二次均匀B样条插值曲线

在利用反求法构造B样条插值曲线时,往往需要选取端点条件。 因此,可对端点条件进行优化选取,使得构造的B样条插值曲线满足特定要求。提出了一种利用曲线内能极小选取平面二次均匀B样条插值曲线端点条件的算法。首先给出了二次均匀B样条插值曲线分控制顶点与首个控制顶点（即端点条件）的递推关系式;然后给出了利用曲线内能极小优化选取首个控制顶点的算法,证明了利用该算法构造的C¹连续二次均匀B样条插值曲线为保形插值,并通过数值算例证明了算法的有效性;最后,为便于实际应用,基于MATLAB平台设计了算法所对应的图形用户界面,用户通过简单的操作即可获得光顺的C¹连续二次均匀B样条保形插值曲线。     

有耗介质层上多导体传输线的电磁耦合时域分析方法

目前,针对空间电磁场作用有耗介质层上传输线的电磁耦合,仍缺乏有效的数值分析方法.因此,本文提出一种高效的时域混合算法,很好地解决了有耗介质层上传输线电磁耦合建模难的问题.首先,对经典传输线方程进行改进,推导了适用于有耗介质层上多导体传输线电磁耦合分析的修正传输线方程.然后,结合时域有限差分方法和相应插值技术,求解修正传输线方程,获得多导线及其端接负载上的电压和电流响应,并实现空间电磁场辐射与多导线瞬态响应的同步计算.最后,通过相应计算实例的数值模拟,与CST软件的仿真结果进行对比,验证了时域混合算法的正确性和高效性.