抛物问题的mortar有限元逼近的瀑布型多重网格法

用瀑布型多重网格法解决椭圆、抛物问题，已有不少研究工作[1-2]，本文对抛物问题的mortar有限元的全离散格式提出瀑布型多重网格法，证明了该方法是最优的，即具有最优精确度和复杂度．     

不可压缩流动的移动粒子半隐式法数值模拟

介绍了一种新的移动粒子半隐式方法。该方法采用Lagrange描述,用一系列离散粒子代替流体,通过半隐式求解和时间推进法来进行计算,弥补了传统网格法对复杂形状、大变形、高速撞击等情况下网格划分和重构过于繁琐的缺点,在模拟不可压缩流动问题中得到了较好的效果。     

三维目标与粗糙面复合散射的广义稀疏矩阵平面迭代及规范网格算法

提出了快速计算二维导体粗糙面与面上金属目标复合散射的广义稀疏矩阵平面迭代及规范网格法(G-SMFSIA/CAG).推导了二维导体粗糙面与面上目标相互作用的耦合积分方程,用稀疏矩阵平面迭代及规范网格法(SMFSIA/CAG)求解粗糙面部分的表面积分方程,而用基于RWG基函数的矩量法(MOM)计算目标部分的表面积分方程,并通过更新方程的激励项迭代求解目标与粗糙面的相互耦合作用.结合Monte-Carlo方法产生具有PM(Pierson-Moskowitz)海浪谱的随机海洋粗糙面,数值分析了海面上不同形状导体目 关键词： 复合散射 广义稀疏矩阵平面迭代及规范网格法 随机海洋粗糙面 双站散射系数     

二次Lagrangian有限元方程的几何多重网格法

为了构造快速求解二次Lagrangian有限元方程的几何多重网格法,在选择二次Lagrangian有限元空间和一系列线性Lagrangian有限元空间分别作为最细网格层和其余粗网格层以及构造一种新限制算子的基础上,提出了一种新的几何多重网格法,并对它的计算量进行了估计.数值实验结果,与通常的几何多重网格法和AMG01法相比,表明了新算法计算量少且稳健性强.     

Using Digital Imaging to Characterize Threshold Dynamic Parameters in Porous Media Based on Lattice Boltzmann Method

Digital images (DI) and lattice Boltzmann method (LBM) are used to characterize the threshold dynamic parameters of porous media. Two-dimensional representations of the porous structure are reconstructed from segmentation of digital images obtained from a series of tiny samples. The threshold pressure gradients and threshold Péclet numbers are researched on seven test samples by using LBM. Numerical results are in agreement with that obtained by integrating Darcy’s law. The results also indicate that fluids can flow through porous media only if the fluid force is large enough to overcome threshold pressure gradient in porous media. One synthetic case is used to further illustrate the applicability of the proposed technique. In addition, the dynamical rules in our model are local, therefore it can be run on parallel computers with well computational efficiency.     

二维抛物型方程的高精度多重网格解法

提出了数值求解二维抛物型方程的一种新的高精度加权平均紧隐格式，利用Fourier分析方法证明了该格式是无条件稳定的，为了克服传统迭代法在求解隐格式是收敛速度慢的缺陷，利用了多重网格加速技术，大大加快了迭代收敛速度，提高了求解效率，数值实验结果验证了方法的精确性和可靠性。     

自适应网格加密和多重网格求解严重病态方程

本文探讨用自适应网格加密技术结合多重网格法求解由有限元离散导致的病态方程。计算结果表明；自适应多重网格有限元法是求解严重病态方程珠有效方法。     

ADAPTIVE MESHLESS METHOD BASED ON LOCAL FIT TECHNOLOGY

An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement.New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.     

RKPM形状函数的矩式显式表述及快速计算

给出了计算再生核质点法(RKPM)形状函数及其导数的矩式显式处理方法。其特点是在计算形状函数及其导数时不涉及矩阵的求逆或者线性方程组的求解，从而减少计算误差的产生并提高了计算速度。二维及三维形状函数计算算例表明该方法是提高RKPM计算效率的一种有效途径。