Please wait a minute...

过刊目录

    1986年, 第3卷, 第2期
    刊出日期:1986-06-25
    上一期    下一期

    论文   
    全选选: 隐藏/显示图片
    论文
    解二维柱几何输运问题的离散——配置有限元方法
    杜明笙, 游腾高
    1986, 3(2): 129-141. 
    摘要 ( )   PDF(588KB) ( )  
    本文研究了求解多群多种介质的二维柱几何中子输运问题的离散一配置方法(TCC)。该方法对时间t用Crank-Nicholson中心差分,源项采用离散纵标近似和源迭代技巧,对空间变量采用配置有限元方法。本文讨论了算法。证明了守恒性。并且证明了本方法在K=1时等价于DSN[2]方法,因此本方法是DSN方法的高阶推广。本文还给出了数值计算结果。我们对精确解以及多个二维柱几何的临介物理问题进行了数值计算,由分析比较可见,取得了较好的结果。
    一维冲击波和爆轰波计算程序SSS
    孙承纬
    1986, 3(2): 142-154. 
    摘要 ( )   PDF(731KB) ( )  
    本文叙述一维程序SSS的理论和结构。介质采用流体弹塑性模型。爆轰反应有四种可选方式:Arrhenius律、C-J比容法、方波法和Forest Fire方程。爆轰产物可选用HOM或JWL状态方程。SSS程序的特殊功能有,激光效应、爆轰增长、强爆炸波、散心爆轰和瞬时爆轰等,初始起爆或能量释放位置可任意设定。本程序部分参考了美国洛斯阿拉莫斯实验室的一维程序SIN。本文列举的算例把这两个程序作了比较。
    低Re数非线性绕流问题的一种算法
    张峻岫, 雷旭明
    1986, 3(2): 155-160. 
    摘要 ( )   PDF(370KB) ( )  
    本文对低Re数非线性绕流问題提出了一种算法。做法是以Oseen线化方程为基础通过迭代修正来考虑非线性的惯性项的影响,即先将Navier-Stokes方程写作Oseen方程外加一"强迫函数"项的形式,然后用基本解‘Oseenlet’的积分形式给出其解式,它是一个含物面积分和流动空间积分的非线性积分微分方程。设沿边界上的‘Oseenlet’分布强度为待求量,由边条件加以确定。空间体积分由于有赖于流场,而它在求出解前是未知的,故采用迭代修正的做法进行处理。迭代过程从线性流场开始进行,直到算出达到规定精度要求的收敛解为止。作为算例,对圆柱绕流问題进行了计算,给出了圆柱阻力随Re数的变化规律,并同实验及有关计算资料做了比较。结果表明,本文方法是令人满意的。
    用有限元法解龟磁波穿透非均匀有耗等离子体鞘套时的损耗
    刘铁军, 石中岳
    1986, 3(2): 161-170. 
    摘要 ( )   PDF(596KB) ( )  
    本文用有限无法求解了电磁波穿透具有任意密度剖面的非均匀等离子体鞘套时的反射和透射问題。鞘套的剖面参数选自典型再入飞行器的粘性激波层平衡流场数据[1]。利用虚功原理和线性插值函数从Helmholtz方程导出有限元方程。鞘套按对数剖分为31个单元。然后,分别对频率f=400MHZ,4000MHZ和10GHZ时,计算了鞘套内的场值及反射,透射和吸收系数。并估算了信号的总功耗。误差分析表明,计算结果是令人满意的。与其它方法相比较,本文所给出的有限无法是简捷而精确的,所以该方程适于计算波穿透具有任意的而且变化剧烈的剖面的等离子体鞘套时的损耗
    高空核爆炸产生的电子数密度的近似估计
    王泰春, 王玉芝
    1986, 3(2): 171-182. 
    摘要 ( )   PDF(578KB) ( )  
    在高空核爆炸的条件下,由γ-康普顿散射产生的高能电子离化空气产生次级电子。本文推导了描述次级电子的速率方程组的解析解。并把解析解的结果与数值计算的结果作了比较,离地面的高度越高,两种方法的计算结果越接近。
    马赫反射的一种简便的数值模拟方法——分裂拟特征线法
    刘儒勋, 魏勇
    1986, 3(2): 183-193. 
    摘要 ( )   PDF(553KB) ( )  
    本文提出马赫反射的一种简易、有效的数值模拟方法——分裂拟特征线法,并对平面激波碰到斜劈时所形成的马赫反射进行了模拟计算。具体格式设计灵活,程序实现容易,运算量和信息量比已有的其它方法都有很大的节省。计算结果表明,马赫反射图象和流场基本特征与已有的实验或计算结果比较,是比较吻合的。
    解线性代数方程组的TCG迭代法
    陈创飞
    1986, 3(2): 194-206. 
    摘要 ( )   PDF(672KB) ( )  
    ICCG方法是解线性代数方程组较为理想的方法,但它仅适用于具有正定对称的系数阵。本文介绍的TCG方法便足改造过的ICCG方法,它适用于一般非奇异的非正定对称阵。TCG方法比常用的JLUCG方法,对于非定常问题,可提高效率18%,特别是取用SIP不完全L、U分解作预条件时,可提高效率40%,是计算非正定对称阵较为理想的迭代法之一。本文推导出在消去法不完全L,U分解下的TCG方法,并用数值结果论证出它比ILUCG方法加速收敛的所在。
    QZ算法误差分析
    董浩喆
    1986, 3(2): 207-216. 
    摘要 ( )   PDF(460KB) ( )  
    在这编文章中,我们给出求解λB-A型矩阵偶特征值的QZ算法的误差分析。全文分两部分。在第一部分中,我们估计QZ算法的变换次数。在第二部分中,我们给出误差估计式。
    解不可压缩流体力学问题的降阶法——Ⅱ.空间Vh的基函数和误差估计
    于欣
    1986, 3(2): 217-226. 
    摘要 ( )   PDF(589KB) ( )  
    本文对于一大类数值求解二维Navier-Stokes方程边值问题的有限元格式给出了零散度空间Vh的一组简单基函数,讨论了速度的数值误差对压力的数值解的影响,并提出一个改进算法。
    求解三维非定常流的一种拉氏方法
    刘邦弟
    1986, 3(2): 227-233. 
    摘要 ( )   PDF(396KB) ( )  
    本文提供了一种三维非定常流的有限差分方法。是一种拉格朗日方法。本文构造了一个二十四面体。将动量守恒方程在这个二十四面体上积分,建立的差分方程,形式简单,容易计算。
    定态对流扩散方程的修正Dennis格式
    陆金甫, 关治
    1986, 3(2): 234-240. 
    摘要 ( )   PDF(358KB) ( )  
    本文提出了求解定态对流扩散方程的修正Dennis格式。对模型问题的初步分析和数值计算表明在相当程度上克服了Dennis格式的异常不精确性质(Peclet数大时)。
    用质点源方法计算三维内电磁脉冲
    王玉芝, 王泰春
    1986, 3(2): 241-248. 
    摘要 ( )   PDF(412KB) ( )  
    本文给出用质点法求解三维长方腔体中的初级电流及其相应的内电磁脉冲的数值结果。井把这种方法得到的结果与三维欧拉方法的计算结果作了比较。在相同的计算条件下,两种方法得到的结果符合较好。
    EXAFS数据分析中振幅衰减因子标度常数近似影响的模拟研究
    郑宇, 华中一
    1986, 3(2): 249-254. 
    摘要 ( )   PDF(398KB) ( )  
    本文就EXAFS曲线拟合数据分析方法中振幅衰减因子的标度常数近似对晶态和非晶态材料结构参量分析精度的影响进行了模拟研究,并且讨论了解决问题的方法。
分享到: