Please wait a minute...

过刊目录

    1991年, 第8卷, 第3期
    刊出日期:1991-09-25
    上一期    下一期

    论文   
    全选选: 隐藏/显示图片
    论文
    声波方程数值反演的FR法
    张文飞
    1991, 8(3): 225-229. 
    摘要 ( )   PDF(312KB) ( )  
    本文提出了声波方程反演的FR方法,FR法是一个应用范围很广的波动方程反演方法,它不受空间维数与边界条件的限制,适应于各类波动反问题。文中推导了算法的基本公式,并给出了一维与二维声波反问题的数值算例。
    等离子体粒子模拟双时标算法
    聂小波, 张家泰, 许林宝
    1991, 8(3): 230-234. 
    摘要 ( )   PDF(300KB) ( )  
    给出能节省1/2-1/4费用的等离子体粒子模拟双时标算法,方法是用小步长积分电子运动方程,用大步长积分离子运动方程。已在一维半电磁粒子云模拟程序CIC-R上实现双时标算法。理论分析和计算机模拟结果肯定了这个算法的稳定性和正确性。
    模拟无粘超音速射流/绕流干扰的高分辨率自适应网格数值方法
    吴颂平, 李椿萱
    1991, 8(3): 235-242. 
    摘要 ( )   PDF(775KB) ( )  
    本文提出一种模拟飞行器超音速尾喷流与飞行器超音速绕流相互作用所形成的复杂流场的差分方法。本方法对流场中的光滑区域使用传统的二阶格式,而在流动参量变化剧烈的区域中采用细网格的、具有高分辨率的TVD格式。为此,文中设计了一种逐维自适应网格加密技术,并与Strang型分裂格式相结合,以求解多维的气体动力学方程。作为算例,本文给出了以Euler方程为模型的尾喷流/绕流干扰流动的计算结果。数值计算表明,本文所提出的自适应网格加密及局部使用TVD格式的算法不仅可以取得高清晰度的流场显示,同时还可使整个计算的效率得到大幅度提高。
    用蒙特卡罗方法计算托卡马克中高能α粒子的能量损失及分布
    齐昌玮
    1991, 8(3): 243-248. 
    摘要 ( )   PDF(380KB) ( )  
    本文采用蒙特卡罗方法计算了托卡马克工程试验堆中高能逃逸α粒子的能量损失及分布。发现慢化过程中由于碰撞扩散而引起的能量损失远比过去预料的大。对抛物分布的等离子体,慢化过程的能量损失占总能量损失的绝大部分。而逃逸α粒子仅打在壁的上半部分或下半部分(由磁场方向决定),在极向和纵向上分布不均匀。在纵向上损失能量峰值偏离纵向磁场最小处,偏移方向由电流方向决定。
    对流扩散方程的涡区分离解法
    张世雄
    1991, 8(3): 249-256. 
    摘要 ( )   PDF(561KB) ( )  
    对流扩散方程是流体计算中一个基本方程,常用的数值方法导至解一个高阶的代数方程组,要求较大的存贮量和较长的计算时间。本文提出一种涡区分离解法,它利用对流扩散方程的迎风性质,把涡区从对流支配区分离出来,仅在各个涡区建立代数方程组并求解。而在对流支配区,则充分利用其抛物性,只需采用显式格式进行计算。由于在各涡区建立的这些方程组阶数和带宽都较小,因此要求存贮量较小,计算速度较快。对于雷诺数较大,涡区范围较小的问题,该方法特别有效。
    正规三角形网格上二维非线性守恒律组的二阶精度MmB差分格式
    杨树礼
    1991, 8(3): 257-263. 
    摘要 ( )   PDF(367KB) ( )  
    本文在正规三角形网格上对二维非线性守恒双曲方程组的初值问题构造一类二阶精度MmB(locally Maximum-minimum Bounds preserving)差分格式,并用此格式对二维无粘Bergers方程的Riemann初值问题进行了计算,结果表明此格式具有高分辨率和非振荡等性质。
    部分潜入水中任意形状等截面柱体的振动分析
    周叮
    1991, 8(3): 264-270. 
    摘要 ( )   PDF(401KB) ( )  
    本文研究部分潜入水中的任意形状等截面柱体的弯曲自由振动问题,应用Fourier级数匹配法,提供了一个精确计算这类柱水耦联系统振动特性的解析方法,结果可由计算机解出。本文最后取水中椭圆柱体为例,给出了一些具体的数值计算结果。
    气藏两相渗流问题全隐式模拟的改进Newton方法及其数值比较
    林刚, 石济民, 林振宝, 吕涛, 林爱民
    1991, 8(3): 271-278. 
    摘要 ( )   PDF(445KB) ( )  
    本文分析含水气藏全隐式数值模似的计算问题,阐述了数值解偏微分方程的各种改进牛顿法,数值计算结果表明:采用萨曼斯基技巧改进牛顿方法可以提高气藏模拟的效率。
    Kaczmarz迭代与流场计算
    史金松, 庄卫华, 章卫军
    1991, 8(3): 279-286. 
    摘要 ( )   PDF(423KB) ( )  
    本文提出基于有限元法及Kaczmarz迭代的一种求解二维恒定流的数值方法,给出了用此方法求解的数值实例。
    气动计算中色散可控的迎风紧致格式
    马延文, 傅德薰
    1991, 8(3): 287-294. 
    摘要 ( )   PDF(480KB) ( )  
    文中通过对修正方程色散项的耗散类比方法,指出该项在改善数值解中非物理振荡的重要作用,给出了一类依赖于三个自由参量的色散可控迎风紧致格式。通过这三个参量可控制耗散量的大小,也可控制色散量的大小及方向,并给出了一个具体的色散协调因子。文中给出的格式有着精度高、方法简单、计算量小和有着强的对激波的捕捉能力等优点。对二维激波反射问题进行了数值实验。计算结果非常令人满意。
    一维无序系统本征矢的解法
    徐慧
    1991, 8(3): 295-304. 
    摘要 ( )   PDF(604KB) ( )  
    本文提出了一种直接求解五对角厄米特矩阵本征矢的方法,即通过代数运算求出本征矢,而不是用三对角块矩阵来解。该方法对计算一维无序系统的本征矢特别适用。计算结果表明,此方法的存储量小,计算精度高,速度快,并可计算上十万个点阵的无序系统本征矢。
    二维二相不规则网格的油藏数值模拟
    马峻, 雷光耀
    1991, 8(3): 305-311. 
    摘要 ( )   PDF(429KB) ( )  
    本文采用不规则网格的差分方法,对平面油水模型进行了数值模拟,并用IMPES方法计算了两个实例。例1的计算结果与精确解吻合良好,比Pedrosa的局部加密网格法使用的网格数少而且计算更简便。例2为矩形网格难以计算的问题,应用本文方法也得出了较好的结果。这表明本文的不规则网格差分方法是可靠而简便的,它可以有效地应用于复杂外边界和复杂地层结构的油藏模拟问题。
    四边形单元中的沙漏控制
    王肖钧, 胡秀章, 刘胜求
    1991, 8(3): 312-318. 
    摘要 ( )   PDF(435KB) ( )  
    本文对四边形单元的沙漏现象作了分析。通过引入沙漏应变和沙漏应力的方法来控制弹塑性大变形中的沙漏现象。采用二维有限元Lagrange数值计算程序计算了平头弹对平板的高速撞击问题。计算结果表明本文提出的方法是有效的。
    用TVD格式求解一维滑坡问题
    郭文海
    1991, 8(3): 319-324. 
    摘要 ( )   PDF(324KB) ( )  
    本文用二阶精度显式TVD格式计算了一维滑坡问题,得到了和实际值基本一致的结果。
    ORIGEN-2程序的移植、扩充及其应用
    李星月, 归寿造
    1991, 8(3): 325-329. 
    摘要 ( )   PDF(327KB) ( )  
    ORIGEN-2是一个功能很强,用途十分广泛的点燃耗及放射性衰变计算机程序,本文简要介绍它的功能、理论模型及主要特点,描述了程序的移植、扩充和计算裂变产物γ谱的应用。
    一个模拟光子照射人体的蒙特卡罗程序
    李星月, 蔡建华
    1991, 8(3): 330-336. 
    摘要 ( )   PDF(461KB) ( )  
    本文分绍了一个用于辐射防护和医疗目的的蒙特卡罗程序。它能计算出受到光子照射时人体20种器官(组织)的吸收剂量和人体有效剂量当量。程序中使用了一种特殊的方法以加快计算结果的收敛。本文所介绍的方法也可用于其他复杂几何系统的剂量计算。
分享到: