Please wait a minute...

过刊目录

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

    论文   
    全选选: 隐藏/显示图片
    论文
    物态方程对恒星星核区坍缩的影响
    王贻仁, 李鸿, 姚进, 汪惟中
    1989, 6(3): 257-267. 
    摘要 ( )   PDF(707KB) ( )  
    本文对物态方程BBAL和EOS(1)的特性作了对比分析,在低密度区(ρ<2×1010gcm-8),二物态方程的中子,质子化学势有显著差别,这导致了它们的自由核子丰度的巨大差别。考察了有限温度修正与"零温近似"平衡关系式的差异,有限温度修正(热效应)对自由中子、质子化学势影响较小,但对平均重核质量数、电荷数有较大影响。以质量为15M恒星演化模型为初始模型,用上述二物态方程从星核坍缩起动至反弹前过程作了计算。虽然低密度区二物态方程之自由核子丰度相差很大,但星核坍缩计算结果相差不大,而将每一物态方程的有限温度修正和"零温近似"用于坍缩计算,其结果相差较显著。
    有限球中14MeV中子点源问题的数值模拟
    符尚武
    1989, 6(3): 268-276. 
    摘要 ( )   PDF(469KB) ( )  
    点源问题的物理实验与数值计算是研究反应堆防护和辐射效应,考查中子参数和计算方法的重要手段。本工作解决了有限元节点上的加源问题,并对几个典型的14MeV中子点源问题的物理实验,用双向间断有限元法进行了数值模拟,并和其它方法的计算结果及物理实验的测量结果进行了比较,取得了较为满意的结果。
    二维翼型跨音速绕流的数值模拟
    陈让福
    1989, 6(3): 277-288. 
    摘要 ( )   PDF(697KB) ( )  
    本文结合高精度TVD格式的数值通量和时间进展多步法给出了一种求解定常流问题的数值方法。同时给出一些特殊处理来加快数值解的收敛速度。本文用以上方法计算了翼型跨音速绕流问题,结果表明此方法具有分辨率高,收敛速度较快之性质。
    用大气环流模式研究两个月的气候平均环流演变
    姜达雍, 杨梅玉, 闵伟
    1989, 6(3): 289-301. 
    摘要 ( )   PDF(856KB) ( )  
    地球上的气候演变,是在一个包括大气、海洋、冰、土壤和植物的多相气候系统中,由许多物理、化学、生物过程相互作用的综合结果。对这些过程相互作用的了解,是基于有能力定量地、以相容的表达式将它们结合起来考虑、这种模试依靠一套方程组的数值积分和构成气候数值模式的经验关系。本文给出了一个全球经纬度网格的大气环流模式(AMS GCM),考虑了各种物理过程,并利用实际资料进行了两个月的积分,初步结果分析,效果是令人鼓舞的。
    H在Ni表面的扩散性质及表面缺陷的影响-嵌入原子法(EAM)计算研究
    曹培林, 施丹华
    1989, 6(3): 302-308. 
    摘要 ( )   PDF(492KB) ( )  
    本文用原子集团模型和嵌入原子方法,计算研究了H在Ni(100),(110)和(111)表面上的扩散特性,其扩散激活能分别是0.152eV;0.343eV和0.142eV。Ni(100)面上的吸附Ni原子和Ni原子空位,分别是在其表面上扩散的H原子的一种陷阱和位垒;Ni(100)面上的台阶,将使通过此台阶的H原子的扩散势垒增高,激活能增大,且引起扩散的各向异性。
    计算二维理论地震图的分步方法
    何柏荣, 冯德益, 王振彪, 葛晓斌
    1989, 6(3): 309-320. 
    摘要 ( )   PDF(672KB) ( )  
    本文以二维半无限空间内压缩线源激发的P-SV波为对象,讨论了求解具有间断系数的二维弹性波方程组初边值问题的几种分步方法。震源的贡献是以一个边界条件给出的。为了检验本方法的功能,我们给出了两个计算实例。
    定常不可压缩粘性流体流动Navier-Stokes方程的推进迭代法
    苏铭德, 陈霜立
    1989, 6(3): 321-334. 
    摘要 ( )   PDF(664KB) ( )  
    本文将推进迭代法推广到不可压缩流的N-S方程求解。给出了曲线坐标下的一般表达式,将推进迭代法应用于方框流及圆柱绕流,得到的结果表明该方法具有较快的收敛速度和能够得到满意的结果。
    运动边界产生的孤立波解
    吕秋强
    1989, 6(3): 335-339. 
    摘要 ( )   PDF(266KB) ( )  
    本文给出了求解劝边界Boussinesq方程的差分格式,计算了动边界以不同形式运动时所产生的波。数值结果表明了孤立波解的存在。
    Petrov-Galerkin有限元法求正则化长波方程的数值解
    唐世敏, 王唯
    1989, 6(3): 340-346. 
    摘要 ( )   PDF(423KB) ( )  
    本文给出数值求解正则化长波方程(简称为KLW方程)初值问题的一种Petrov-Galerkin有限无法。由于它对时间和空间变量分别有二阶和四阶精度,所以对一个孤立波的进波解,数值结果与精确解极为吻合。对于两个孤立波的相互碰撞,则发现碰撞后有微小振荡尾波存在,从而严格说来,它们不具备孤立子的性质。
    毫米波自由电子激光的空间电荷效应
    杨震华, 田世洪
    1989, 6(3): 347-354. 
    摘要 ( )   PDF(501KB) ( )  
    本文编制了三维数值模拟程序,考虑了电子在摇摆器(Wiggler)内纵向速度βz=υz/c<1的影响,考虑了电子在有质动力势阱中因"群聚"而引起的纵向空间电荷力的影响,考虑了电子束内静电波和激光场、Wiggler静磁波的三波相互作用对电子位相的影响等。对因"群聚"而引起的纵向电场满足的方程求出了近似的二维解析解,用4096个电子模拟求解了电子运动方程和光场方程。数值模拟的结果和Livermore实验室的实验结果进行了比较。结果表明,由于上述几项因素的考虑,我们的数值模拟结果和实验位基本一致,消除了文献[1,2,3]中出现的失谐曲线(detuning curve)峰值磁场漂移300-500高斯的差距。
    一般表载作用下的层状弹性系统的变形
    潘尔年
    1989, 6(3): 355-365. 
    摘要 ( )   PDF(685KB) ( )  
    本文提出了一种适应于求解多层各向同性弹性静力学问題的一般的且简洁的处理方法——即基于笛卡尔坐标与柱坐标两组向量函数系的传播矩阵法,从而避免了此类问题研究中的复杂性。在两向量函数系下分别得到了在一般表载作用下介质内任意点处的位移及应力分量的解析表达式。文中给出的算例既证实了所列公式的正确性,又表明:在数字积分中结合改进的Padé近似求和法后,本文所提出的方法是非常简单且易于实际应用的。
    人造软X射线单色器结构与性能的模拟计算
    刘文, 刘文汉, 吴自勤
    1989, 6(3): 366-370. 
    摘要 ( )   PDF(317KB) ( )  
    本文用计算机模拟计算的方法较详细地研究了W/C和Mo/Si人工多层膜x射线单色器结构与性能的关系。得到了多层膜层数和多层膜膜厚偏差对x射线衍射峰峰度、峰宽影响的一系列数据。本文还比较了膜厚无规涨落与系统线性偏差这两种偏差对单色器能量分辨率的影响,计算的结果表明后者的危害较大。
    随机数检验程序系统SUTEST
    张建中
    1989, 6(3): 371-377. 
    摘要 ( )   PDF(404KB) ( )  
    对计算机上给出的(0,1)上的数字序列r1,r2,…,ri,…,rN能否当做标准均匀分布的随机数用于Monte Carlo模拟计算,是Monte Carlo方法研究及其实际应用中的一个极其重要而又十分难于解决的问題。统计软件包SASD提供的随机数检验程序系统SUTEST,对数字序列{ri|0< ri< 1,i=1,2,…,N}给出十二类、二十七种不同的统计检验方法,由六十一个统计量从均匀性、随机性、独立性和统计模拟计算等许多不同的方面,为随机数检验提供了完整、配套的统计检验方法和有效可行的算法。给出了便于使、易于移植的FOR-TRAN源程序。在上述意义下,从统计上和实际应用上解决了这一问题。
    GNASH程序及其在中子核反应研究中的应用
    宓詠, 霍裕昆, 李向阳, 袁竹书
    1989, 6(3): 378-384. 
    摘要 ( )   PDF(476KB) ( )  
    GNASH是一个功能很强,基于Hauser-Feshbach统计理论并考虑预平衡等非统计效应修正的计算核反应截面及发射粒子能谱的大型计算机程序。本文简要地介绍了这一程序的功能、理论模型、特点,以及我们对这一程序的应用和开发的情况。
分享到: