变度量各向异性网格生成算法和各向异性椭圆偏微分方程匹配网格

对于各向异性问题而言,在进行数值求解时,一个相匹配的各向异性网格至关重要.本文针对二阶椭圆偏微分方程,当其系数矩阵为各向异性时,给出了系数矩阵的逆作为相匹配的各向异性度量,在此度量下生成的各向异性网格即为匹配网格.本文还提出了新的变度量各向异性网格生成算法,该算法通过在各向异性背景网格上结合各向异性Delaunay准则和基于力平衡的结点优化函数生成高质量的各向异性三角化网格.通过数值算例可知,在匹配的各向异性三角化网格上,不仅有条件数小的代数离散系统,还有高精度的数值解,更甚者,在匹配网格结点上的l2范数误差有超收敛现象.     

非结构四边形二次Lagrangian有限元方程的代数多重网格法及收敛性分析

1 引言代数多重网格(AMG)法是求解有结构、非结构冈格下的PDEs离散化系统,以及其它大型稀疏矩阵方程的最为有效的方法之一,是当前多重网格法研究领域的热点．它主要侧重于粗化技术的研究,目前这方面已有大量的工作(见文[1～4]),其中基于部分几何信息和分析信息的AMG法,是当前AMG法发展的一个新趋势(见文[5～7])．     

电缆方程的多重网格并行代数法（英文）

致力于研究求解线性代数方程组的多重网格并行算法,该算法是基于构建矩阵序列的经典Runge-Stuben(RS)方法及其改进的并行修正独立集合(PMIS)方法的.展示了求解离散电缆方程式所得到的线性代数方程组的结果,而电缆方程是用作描述电信号传播的.在求解中用到了GPUPU技术.展示了模型问题在不同尺度的模拟区域上的数值结果.     

杆件轴向受迫振动的Galerkin有限元EEP法自适应求解

基于单元能量投影(element energy projection,EEP)法自适应分析在杆件静力问题以及离散系统运动方程组中所取得的成果,以直杆轴向受迫振动为例,研究并建立了一种在时间域和一维空间域同时实现自适应分析的方法.该方法在时间和空间两个维度都采用连续的Galerkin有限元法(finite element method,FEM)进行求解,根据半离散的思想,由空间有限元离散将模型问题的偏微分控制方程转化为离散系统运动方程组,对该方程组进行时域有限元自适应求解;然后再基于空间域超收敛计算的EEP解对空间域进行自适应,直至最终的时空网格下动位移解答的精度逐点均满足给定误差限要求.文中对其基本思想、关键技术和实施策略进行了阐述,并给出了包括地震波输入下的典型算例以展示该法有效可靠.     

高波数波动问题的多水平方法

本文主要讨论求解高波数Helmholtz方程的多水平方法,主要回顾了一些具有代表性的多重网格方法.如Erlangga等人的shifted Laplacian预处理的多重网格法;Elman等提出的修正的多重网格方法;以及我们的基于连续内罚有限元(CIP-FEM)离散代数系统的多水平算法.最后还介绍了求解高波数时谐Maxwell方程的CIP-FEM离散代数系统的多水平算法.     

三维薄结构热传导问题分层二次元方程的多水平方法

在利用有限元法对三维薄结构进行分析时,为了减少单元数目,常采用六面体薄单元,相应的高阶单元在计算精度、抗畸变程度等方面具有明显优势.但与低阶元相比,高阶单元需要更多的计算机存储空间,离散化线性系统具有更高的计算复杂性,并且系数矩阵是严重病态的,采用通常的求解方法其效率将大大降低.该文针对三维薄结构稳态热传导问题,利用局部块Gauss-Seidel光滑子和基于"距离矩阵"的DAMG法,为其分层二次元离散系统设计了一种具有更好计算效率和鲁棒性(robustness)的多水平方法.由于采用了分层基,程序实现中不再需要建立判定未知数变量指标与所属几何节点类型对应关系的代数判据,网格转换算子的构造也变得非常简单,从而大大提高了运算效率.数值实验结果验证了该方法的有效性和鲁棒性.     

一维强振荡变系数奇异摄动问题的多尺度有限元计算

针对含小参数的一维强振荡变系数奇异摄动扩散方程,提出了多尺度有限元的有效计算方法.通过求解微分算子的子问题得到多尺度基函数,能刻画出强振荡的微观信息.利用节点基函数映射,形成约化降阶矩阵,仅求解小规模代数方程组,在粗网格即能获得高精度结果,从而得到不依赖于摄动系数大小的、一致收敛的数值模拟.     

离散对偶代数Riccati方程异类约束解的双迭代算法

利用逆矩阵的Neumann级数形式,将在离散时间跳跃线性二次控制问题中遇到的含未知矩阵之逆的离散对偶代数Riccati方程(DCARE)转化为高次多项式矩阵方程组,然后采用牛顿算法求高次多项式矩阵方程组的异类约束解,并采用修正共轭梯度法求由牛顿算法每一步迭代计算导出的线性矩阵方程组的异类约束解或者异类约束最小二乘解,建立求DCARE的异类约束解的双迭代算法.双迭代算法仅要求DCARE有异类约束解,不要求它的异类约束解唯一,也不对它的系数矩阵做附加限定.数值算例表明,双迭代算法是有效的.     

Legendre多项式法求一类变阶数分数阶微分方程数值解

本文利用Legendre多项式求解一类变分数阶微分方程.结合Legendre多项式,给出三种不同类型的微分算子矩阵.通过微分算子矩阵,将原方程转化一系列矩阵的乘积.最后离散变量,将矩阵的乘积转化为代数方程组,通过求解方程组,从而得到原方程的数值解.数值算例验证了本方法的高度可行性和准确性.     

L-曲线估计确定正则参数的双网格迭代法

本文考虑对不适定问题离散化得到的大规模不适定线性方程组进行Tiknonov正则化,然后用双网格迭代法求解得到的Tikhonov正则化方程组,并用L-曲线估计法来确定正则参数.试验问题的数值结果表明双网格迭代法求解正则化后的对称正定线性方程组效果很好,且L-曲线估计法确定正则参数计算量很小.     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP ——In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

A unified approach to certain problems of best local approximation

In this paper we study best local quasi-rational approximation and best local approximation from finite dimensional subspaces of vectorial functions of several variables. Our approach extends and unifies several problems concerning best local multi-point approximation in different norms.     

Boundedness for commutators of multipliers on Hardy type spaces

In this paper, we study the commutators generalized by multipliers and a BMO function. Under some assumptions, we establish its boundedness properties from certain atomic Hardy space Hb^p（R^n） into the Lebesgue space L^p with p 〈 1.     

THE 8~(th) INTERNATIONAL CONGRESS ON INDUSTRIAL AND APPLIED MATHEMATICS

<正>August 10-14,2015Beijing,ChinaThe International Congress on Industrial and Applied Mathematics(ICIAM)is the premier international congress in the field of applied mathematics held every four years under the auspices of the International Council for Industrial and Applied Mathematics.From August 10 to 14,2015,mathematicians,scientists     

Towards Excellence in Science Chinese Academy of Sciences

<正>May 26,2014,Beijing Science is a human enterprise in the pursuit of knowledge.The scientific revolution that occurred in the 17th Century initiated the advances of modern science.The scientific knowledge system created by     

Bounds for the zeros of polynomials

Let P(z)=∑↓j=0↑n ajx^j be a polynomial of degree n. In this paper we prove a more general result which interalia improves upon the bounds of a class of polynomials. We also prove a result which includes some extensions and generalizations of Enestrǒm-Kakeya theorem.     

Boundedness of commutators related to Marcinkiewicz integrals on hardy type spaces

In this paper, the authors study the boundedness of the operator [μΩ, b], the commutator generated by a function b ∈ Lipβ(Rn)(0 ＜β≤ 1) and the Marcinkiewicz integrals μΩ, on the classical Hardy spaces and the Herz-type Hardy spaces in the case Ω∈ Lipα(Sn-1)(0 ＜α≤ 1).     

Jacobi polynomials used to invert the laplace transform

Given the Laplace transform F(s) of a function f(t), we develop a new algorithm to find an approximation to f(t) by the use of the classical Jacobi polynomials. The main contribution of our work is the development of a new and very effective method to determine the coefficients in the finite series expansion that approximation f(t) in terms of Jacobi polynomials. Some numerical examples are illustrated.     

Generalization of the interaction between Haar approximation and polynomial operators to higher order methods

In applications it is useful to compute the local average empirical statistics on u. A very simple relation exists when of a function f（u） of an input u from the local averages are given by a Haar approximation. The question is to know if it holds for higher order approximation methods. To do so, it is necessary to use approximate product operators defined over linear approximation spaces. These products are characterized by a Strang and Fix like condition. An explicit construction of these product operators is exhibited for piecewise polynomial functions, using Hermite interpolation. The averaging relation which holds for the Haar approximation is then recovered when the product is defined by a two point Hermite interpolation.