NS方程在隧道通风工程中的应用

在隧道通风工程中,为了研究通风过程中风压风速度的变化情况,将隧道中空气的流动看成流体力学中所研究的粘性不可压缩流体的非定常流动。用NS方程作为隧道通风的数学模型,采用三次插值伪质点法求解NS方程。以浙江台州104国道黄土岭隧道风改建工程为研究对象,对不同长度和高度的引洞通风状况给出计算机数值模拟,并进一步分析通风过程中风压风速的变化情况。从达到通风要求、提高通风效率的角度出发,可根据本文的计算结果与数据分析,选择最优化的通风方案。     

弹性波方程有限元模拟的曲边地表处理

本文对比研究了关于弹性波模拟中的曲边地表形状处理的两种方法,一种是用给定的实际介质数值划定的地表形状,另一种是用样条插值逼近地表形状.本文采用有限元方法进行弹性波数值模拟,给出了基于这两种方法计算的数值例子,并对结果进行了分析比较.结果表明使用后一种方法对地表进行处理时,地表人工离散产生的干扰明显减少,优于前一种方法.     

线性Poisson-Boltzmann方程的Mortar有限元方法的数值计算

本文对分子生物物理学中产生的线性Poisson-Boltzmann方程（PBE），给出了Mortar有限元方法的计算过程，数值计算例子表明，与一般的协调有限元方法相比，用Mortar元方法求解此类有间断系数的问题有非常有效的。     

二阶椭圆问题带单位分解技巧的两重网格方法

标准的两重网格方法是一种求解二阶椭圆问题的局部并行方法,其计算所得数值解在整个求解区域上并不连续使用单位分解技术,将各个子区域上的局部解粘合在一起,从而得到全局连续解,并证明此解在H1范数意义下最优.更进一步,可以证明通过在粗网格上修正,能够改善其L2误差.数值例子验证了理论的正确性.     

基于透射边界条件的高阶离散型角点条件

对波动方程的数值模拟中，在有限区域建立吸收边界条件，其中对区域角点的处理是一个很重要的问题．随着吸收边界条件阶数的提高，与之匹配的角点条件也越难建立．MTF是一种离散型吸收边界条件．在此，对于二维问题，基于MTF建立离散型高阶角点条件，对计算区域角点处理时，在区域对角线方向上建立N阶MTF公式．问题也可推广到三维．数值结果证实了我们的猜测．     

Numerical computation of the flow around two square cylinders arranged side-by-side

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｗａｋｅｉｎｔｅｒｆｅｒｅｎｃｅｗｉｔｈｃｏｍｐｌｅｘｃｏｎｆｉｇｕｒａｔｉｏｎｓｃｏｎｓｉｓｔｉｎｇｏｆｍｕｌｔｉｐｌｅｂｌｕｆｆｂｏｄｉｅｓｉｓａｔｔｒａｃｔｉｎｇａｔｔｅｎｔｉｏｎｏｆａｌｏｔｏｆｒｅｓｅａｒｃｈｅｒｓ,ｂｅｃａｕｓｅｏｆｉｔｓｐｒａｃｔｉｃａｌａｐｐｌｉｃａｔｉｏｎｓｔｏｉｎｄｕｓｔｒｙ.Ｆｏｒｅｘａｍｐｌｅ,ｔｗｉｎｓｔｒｕｔｓｔｏｓｕｐｐｏｒｔｗｉｎｇｓｉｎｔｈｅｆｉｅｌｄｓｏｆａｅｒｏｎａｕｔｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ;ｔｗｉｎｃｈｉｍ…     

两并列方柱绕流相互干扰的数值研究

用数值方法研究了绕流两并列方柱的气动力及其相互干扰.使用了陈素琴等人提出的对流项用三阶逆风格式的带多重网格的改进的MAC方法,对两并列方柱的流场进行了模拟.结果表明:两并列方柱在不同的间距比时其干扰特征有很大的不同,小间距比时出现双稳态偏流,流动并不对称于双方柱间隙的中心轴线,而是偏向其中的一个方柱.     

A MULTI-GRID ALGORITHM FOR MIXED PROBLEMS WITH PENALTY BY C°-PIECEWISE LINEAR ELEMENT APPROXIMATION

In this paper we describe a multi-grid algorithm for mixed problems with penalty by the linear finite element approximation. It is proved that the convergence rate of the algorithm is bound ed away from 1 independently of the meshsize. For convenience, we only discuss Jacobi relaxation as smoothing operator in detail.     

Two-grid partition of unity method for second order elliptic problems

A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.