计算物理杂志 北大期刊 CSCD期刊 统计源期刊

自适应间断有限元方法求解双曲守恒律方程

作者：徐云 蔚喜军 刊期：2009年第02期

研究自适应Runge-Kutta间断Galerkin（RKDG）方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量.

一种抗单元逆转的软组织有限元模型

作者：张秋葵 罗立民 刊期：2009年第02期

由于Green应变具有单元翻转不变性,在大变形的仿真中有可能会出现单元的逆转,破坏模型网格的拓扑结构,导致仿真失败,因此提出一个新的有限单元模型,在采用Green应变度量变形同时把单元体积变化率和单元内压强关系引入到本构方程中,解决单元逆转问题;同时,模型导出的单元弹性势能是单元结点位移的多项式函数,有效地简化了共轭梯度法的求解,效率显...

非结构四边形网格下的一类保对称有限体元格式

作者：聂存云 舒适 盛志强 刊期：2009年第02期

针对定常扩散问题,在非结构四边形网格下,通过选取特殊的控制体和有限体元空间,建立两种保对称有限体元格式,在拟一致网格剖分下,当扩散系数光滑时,证明有限体元解函数在L2和H1范数下均具有饱和误差阶.数值实验验证理论结果的正确性,同时表明新格式对扭曲大变形四边形网格、间断系数问题具有较强的适应性.在正交网格下,第二种格式对流（flux）函...

基于内存约束的一维负载平衡方法及其应用

作者：刘旭 莫则尧 曹小林 刊期：2009年第02期

在内存约束的条件下,提出一种新的一维负载平衡方法.该方法是一个迭代法,严格的理论证明表明,它可以线性收敛到最优剖分.针对典型的负载不平衡模型和三维分子动力学问题,1 024个处理器上的实际应用表明,该方法可以取得较好的效果.

块三对角线性方程组的一类二维区域分解并行不完全分解预条件

作者：吴建平 宋君强 张卫民 李晓梅 刊期：2009年第02期

基于二维重叠区域分解,对每个子区域上局部不完全LU分解所得到的上、下三角因子分别进行组合,给出一类全局并行不完全分解型预条件.所给出的并行化方法适用于任何不完全LU分解型预条件.对采用二维区域分解与一维区域分解时所得并行预条件的并行计算性能进行分析比较.实验结果表明,提出的并行化方法普遍优于加性Schwarz并行化方法,且当处理器个数...

含运动变形物体流场的数值模拟

作者：邓树升 耿继辉 叶经方 谭俊杰 刊期：2009年第02期

基于含拓扑结构改变、可变形运动网格生成技术及动边界条件下的数值处理方法,对爆炸波通过有限入口进入建筑物内部的透射（infiltration）现象进行直接模拟.系统研究入口处门板不同运动、变形、脱落方式下,爆炸波在建筑物内部的传播过程和内壁面上的爆炸载荷,直接模拟门与挡板的碰撞、门的断裂及破片的抛洒过程.结果表明,所发展的用于含拓扑结构...

高阶精度CE/SE算法及其应用

作者：张德良 王景焘 王刚 刊期：2009年第02期

对时-空守恒元解元算法（CE/SE）的网格设置做较大改进,提出一种新的六面体解元和元定义;同时在解元中对物理量进行高阶Taylor展开,给出一种在时间和空间上均具有高阶精度CE/SE算法.在此基础上,把新型的高阶精度CE/SE算法推广应用于高速流动捕捉激波间断、气相化学反应流动、计及固体动态效应的流体-弹塑性流动和非稳态多相不可压缩粘性流动中....

用RANS/LES混合方法研究超声速底部流动

作者：肖志祥 符松 刊期：2009年第02期

在剪切应力输运和弱非线性k-ω模式中引入压缩性修正,并在此基础上构造RANS/LES混合方法,对包含丰富流动结构和复杂流动机理的导弹类超声速底部流场进行数值分析,发现压缩性修正可更好地描述流场特征,适于具有较强压缩性效应流动的数值分析;混合方法捕捉到丰富的流动现象及非定常特性.

二维地效翼及地效流动特性数值研究

作者：杨韡 杨志刚 刊期：2009年第02期

用数值模拟的方法,对二维NACA0012翼型在地面效应下的空气动力特性和地效流动特性进行研究,得到地效翼的升力、阻力和翼型表面压力分布随攻角及相对飞行高度的变化情况.通过对计算结果的分析,可以看出,在一定的攻角,靠近地面飞行,机翼的升力得到提高;随着飞行高度的降低,地面效应增强,机翼的失速攻角减小;地面附近的粘性流动对机翼的空气动力特...

第十届全国蒙特卡罗方法及其应用学术交流会第一轮通知

刊期：2009年第02期

氯气/BHP液滴化学反应流动一维数值模拟

作者：薛社生 刘全 李守先 束小建 刊期：2009年第02期

为研究横向均匀液滴发生器的机理,对Cl2/He气体横向流过BHP液滴场的物理化学行为建立一维化学反应流动模型.其中气流场按不可压缩假设近似为一维均匀流动;发生于液滴内的化学反应由四种物质的反应扩散方程描述;气相中各种物质浓度由对流扩散方程描述.液滴是这些物质的源或汇,两个模型方程通过此源项相联系.对模型方程作数值求解.结果显示,氯气利...

辐射与物质强耦合问题的二维输运数值模拟

作者：李双贵 杭旭登 李敬宏 刊期：2009年第02期

分析确定性输运数值方法模拟辐射与物质强耦合问题的难点,给出在任意四边形网格上求解该类问题的简单隅角平衡算法,并应用灰体输运综合加速算法提高输运方程的源迭代收敛速度,最后给出数值例子.

原子轨道强耦合方法研究He^2＋-H^-碰撞的电荷转移过程

作者：闫世杰 刘玲 王建国 刊期：2009年第02期

应用双中心原子轨道强耦合方法研究He^2＋-H-碰撞的单次电荷转移过程.计算中,对入射粒子He^2＋,包含n=1-7的所有束缚态,计算的能量本征值与NIST标准数据在百分之几的精度内符合很好;对靶H-,包括一个束缚态1s和五个连续态ns（n=2-6）,束缚态能量与他人理论结果一致.在4-400 keV的入射粒子能量范围,计算单电子俘获过程的总截面及到各个壳层上的态...

双原子分子电子能量损失谱的解谱方法

作者：刘雯静 朱林繁 葛敏 徐克尊 刊期：2009年第02期

介绍一种双原子分子电子能量损失谱的解谱方法,在利用已知的双原子分子光谱信息及理论计算的Franck-Condon（FC）因子的基础上,通过二次函数调整整个电子态振动强度的轮廓,减少拟合过程中的待定参数个数,降低对FC因子精确性的要求,提高拟合的精度和可信度.用Mathematica语言编写基于以上方法的数据处理程序,改进算法,提高了软件的运行效率.利用...

半透明梯度折射率介质内辐射熵传递方程及其数值模拟

作者：刘林华 刊期：2009年第02期

在非相干辐射条件下,基于Planck光谱辐射熵强度定义,导出半透明梯度折射率介质内光谱辐射熵传递方程,以及局部辐射熵产率理论表达式.基于离散坐标法对辐射熵传递方程进行数值求解.以一维半透明梯度介质平板为例,对辐射熵方程及其算法进行验证.平板整体无因次辐射熵产的计算结果与宏观热力学定律的结果一致.