Please wait a minute...

过刊目录

    1985年, 第2卷, 第1期
    刊出日期:1985-03-25
    上一期    下一期

    论文   
    全选选: 隐藏/显示图片
    论文
    单调二次样条插值及其在状态方程计算中的应用
    张宝琳
    1985, 2(1): 1-7. 
    摘要 ( )   PDF(416KB) ( )  
    本文证明了如下结果:设a,b,ya,yb,ya',yb',均为实数,满足aab;ya'≥0,yb'≥0;则当ya'+yb'≤4(yb-ya)/(b-a)时,满足s(a)=ya,s(b)=yb,S'(a)=ya',S'(b)=yb'的以a,(a+c)/(2),b为插值节点的单调增加二次样条函数S(x)εC'[a,b]存在且唯一,而且当ya'+yb'>4((yb-ya)/(b-a))时,满足上述要求的二次插值样条函数不存在。根据这个结果讨论了多节点单调二次插值样条函数的存在性,并且给出了在不同区域间状态方程的能量与压强曲面光滑衔接的一个计算方法。
    等离子体平衡逆问题的一种数值解法
    冯桂云
    1985, 2(1): 8-16. 
    摘要 ( )   PDF(650KB) ( )  
    本文借助离散化算子,把滤波方法推广到非圆截面,对具有非圆截面的等离子体平衡逆问题建立了一个稳定的显式逆推滤波算法;给出了部分非圆截面的计算结果;对圆截面情形的数值解与解析解作了比较。
    电子与固体相互作用的Monte Carlo模拟
    何延才
    1985, 2(1): 17-23. 
    摘要 ( )   PDF(481KB) ( )  
    本文提出了一个描述电子与固体相互作用的物理模型及Monte Carlo模拟方法。其计算结果较为简单,对背散射系数、透射系数及x射线深度分布函数的计算结果与文献报道的实验结果相符。本模型己成功地应用于薄膜、微颗粒体x射线定量分析的工作中。
    低能强流离子束在加速-减速系统中传输特性的数值计算方法
    华大平, 方锦清
    1985, 2(1): 24-37. 
    摘要 ( )   PDF(926KB) ( )  
    本文介绍了对于强流离子束在加速-减速系统中传输特性的有效的二维数值解方法。我们要求如下的数值解:(1)非线性泊松方程。(2)非相对论的离子运动方程。(3)离子流连续性方程。因为几何边界的复杂性。采用有限元方法求解方程(1),计算机程序己编制好。给出空间电位分布,离子轨迹和空间电荷密度分配的自洽解。方法对于研究强流离子束在加速-减速系统中传输问题是适用的,特别是对于高强流离子注入设备、强流离子束外延机、强流离子发生源等等都是适用的。
    源凸域正贡献蒙特卡罗方法
    裴鹿成
    1985, 2(1): 38-46. 
    摘要 ( )   PDF(611KB) ( )  
    在Duoi和Gerstl的贡献蒙特卡罗方法中,不少次级粒子的贡献是负的[3,10]。负贡献粒子常使蒙特卡罗估计的统计涨落变大,浪费较多的计算时间,因此,解决次级粒子的负贡献问题是很必要的。在本文中,一个新的称为源凸域正贡献方法的贡献蒙特卡罗方法被提出来了。在任何情况下,源凸域正贡献方法一定比Dubi和Gerstl的贡献方法好,因此,前者完全可以替代后者。对于源凸域R解决了次级粒子的负贡献问题,通过例子的实际计算结果表明,本方法的效率可比Dubi和Gerstl的方法提高1.5至3.4倍。
    吸引子维数的数值计算
    王光瑞, 陈式刚, 郝柏林
    1985, 2(1): 47-54. 
    摘要 ( )   PDF(554KB) ( )  
    对自冶的强迫布鲁塞尔振子采用Poincare截面方法,把吸引子维数计算中涉及的相空间降低一维。这样就可能在内存较小的计算机上计算柯尔莫哥洛夫容量dc,并且也方便了李雅普诺夫维数dL的计算。在计算dc时,使用不断含并小格子的办法,在只解一次微分方程的情况下,可同时得到一组不同ε的含有点集元素的小格子数N(ε)。这样做,既节省机器时间,也节省了内存。本文着重讨论数值计算中遇到的方法和技巧问题。
    散心爆轰的数值模拟
    李智伟, 刘邦弟
    1985, 2(1): 55-66. 
    摘要 ( )   PDF(783KB) ( )  
    本文在瞬时稳态爆轰的假定下,简要叙述了散心爆轰的不定常性质。在散心爆轰的数值模拟方法中,我们改进了人造粘性并用Cochran反应率和JWL状态方程来描述炸药的引爆过程。计了几个典型的散心爆轰问题,得到了比较满意的结果。
    广延X射线吸收精细结构的一种数据分析方法
    汪西川
    1985, 2(1): 67-75. 
    摘要 ( )   PDF(690KB) ( )  
    本文较详细地讨论了广延X射线吸收精细结构(Extended X-ray Absorption Fine Structu-re,简称EXAFS)的数据分析的整个过程。包括背景扣除、富里叶滤波、能量阈值的选取、以及结构参数的确定等问题。引入曲线拟合方法,可使分析结果的精度得到提高。利用所述方法我们分别对Cu、CuNi和CuNiMn等作了计算和分析,获得了较为满意的结果。
    中子慢化积分方程的数值解
    胡理清, 陈文韬
    1985, 2(1): 76-82. 
    摘要 ( )   PDF(430KB) ( )  
    本文对同心环几何条件,用Lobatto求积公式近似求解了中子慢化积分方程,得到了比较精确的超热能谱,并由此算出了共振积分。用逐步递推方法计算慢化源,用予制的Pij表内插求首次碰撞几率,大量节省了机器时间。同时,本文用多能级公式计算235U、239Pu核的共振截面,从而使共振峰间的重叠和干涉也得到了考虑。用文献提供的模型做了校核计算,取得比较满意的结果。
    矩阵偶的特征值与线性自动调节系统稳定性的计算程序
    董浩喆
    1985, 2(1): 83-90. 
    摘要 ( )   PDF(541KB) ( )  
    线性自动调节系统的稳定性分析,可以借助把问题转化为λB-A的特征值问题来解决。本文提出采用QZ算法求解特征值,并且给出实现这一算法的计算程序。它不仅能够解出特征值,而且给出将矩阵偶(A,B)变换成上三角形矩阵偶的变换矩阵。
    推广的Lindemann熔化定律
    张万箱, 张世泽
    1985, 2(1): 91-98. 
    摘要 ( )   PDF(496KB) ( )  
    本文将Lindemann理论应用于非谐振固体的情况,推广了Lindemann熔化定律,并且计算了Fe,Cu,Al的高压熔化曲线。计算结果与实验数据的符合是满意的。
    圆柱绕流问题的初期流动数值模拟
    卢建群
    1985, 2(1): 99-108. 
    摘要 ( )   PDF(607KB) ( )  
    本文继Zheng-huan Teng[1]中介绍的解Navier-stokes方程的椭圆涡团法,研究了一种新的变形涡团法,用以模拟不可压粘性流体绕圆柱的不定常流动。圆柱在静止流体中突然起动并做匀速直线运动。对整个流动区域构造完全Navier-Stokes方程的解并不容易,近十年出现很多数值研究,本文对算法有所推进。把本文的方法称作变形涡团法是因为圆柱边界附近的流体中用椭圆涡团,远离边界时用圆形涡团。计算圆柱绕流比平板绕流在满足附着条件上更为困难,本文分析了怎样在圆柱边界上给出适当的附着条件的数值方法。在算例中雷诺数分别取200、550、3000,得到了不定常边界层分离,二次涡等复杂的物理现象,这些数值结果与近年实验结果[2]是一致的。
    交替方向法在地下水渗流计算中的应用
    雷光耀, 张锁春, 高骥
    1985, 2(1): 109-115. 
    摘要 ( )   PDF(532KB) ( )  
    本文应用交替方向显式法、交替方向隐式法和隐式交替方向迭代法进行了地下水饱和-不饱和渗流的数值计算和比较,研究了计算中出现的迭代不收敛和结果失真等问题。在采用特定形式的附加项的基础上,给出了全隐式交替方向迭代格式,并对不同尺度的计算模型得到了良好的计算结果。文中对附加项的形式和有关参数的选取作了分析比较,给出了较为理想的结果。
    在核物质理论中一个15重积分的化简和计算
    申庆彪, 张竞上
    1985, 2(1): 116-123. 
    摘要 ( )   PDF(442KB) ( )  
    本文讨论了在核物质理论中所出现的受Pauli原理限制的对五个动量矢量的15重积分。首先利用能量守恒和动量守恒条件以及通过适当的变量变换和坐标系的选择,将15重积分化成5重积分,然后再用计算高重积分的数论网格法进行数值计算,并以表格形式在很宽的范围内给出了计算结果。
    孤子的碰撞——带外场的非线性Shrödinger方程的孤子之间的碰撞
    王红鹰
    1985, 2(1): 124-127. 
    摘要 ( )   PDF(209KB) ( )  
    本文用计算机考察了带外场(即变系数)非线性Shrödinger方程的孤子间碰撞机制。发现这种作变速运动的孤子碰后仍能恢复各自的形状。近年来,对孤立子的研究已深入到带外场的非线性Shrödinger方程。
分享到: