Please wait a minute...

过刊目录

    1992年, 第9卷, 第S1期
    刊出日期:1992-12-31
    上一期    下一期

    论文   
    全选选: 隐藏/显示图片
    论文
    差分格式的余项效应研究
    刘儒勋
    1992, 9(S1): 479-484. 
    摘要 ( )   PDF(437KB) ( )  
    采用差分格式的Modified PDE,对其余项进行了统一的、系统的分析和研究,定义了色散关系式和耗散关系式,根据这两个关系式对差分格式的稳定、耗散、色散,以及群速度效应作了较全面的讨论,并指出了余项分析的意义和实用价值。
    一种散乱位势数据网格化的方法
    蔡宗熹, 蔡越虹
    1992, 9(S1): 485-486. 
    摘要 ( )   PDF(160KB) ( )  
    通过偶层场强分量表达式作为第一类Fredholm积分方程,求解虚拟磁偶极子的层磁化强度J(α,β)。再将求得的磁化强度J(α,β)代入偶层场强分量表达式算出矩形网格节点上的位势值。它实现了保留位势数据物理特性的散乱位势数据网格化。
    济南岩溶区地下水水质模型及其数值解法研究
    杨天行, 王运国, 姚磊华
    1992, 9(S1): 487-491. 
    摘要 ( )   PDF(320KB) ( )  
    本文建立了济南岩溶区地下水SO42-及总硬度的双重介质地下水水质运移模型,并提出了一种Galerkin有限元法与广义迎风格式及多单元均衡格式相结合的数值方法。该法克服了数值弥散与数值解的振动现象,结果很好。
    标记零算法与程序软件
    赵玉鹏, 马邃, 徐玉民, 邓飞其
    1992, 9(S1): 492-494. 
    摘要 ( )   PDF(220KB) ( )  
    本文指出近年来关于指派问题若干解法的错误,并提出一种独立的新算法-标记零算法。我们对该算法进行了理论分析。给出了程序软件,并用实例验证了算法的有效性。
    T(q,r)阵BPSD迭代的收敛性
    胡家赣
    1992, 9(S1): 495-497. 
    摘要 ( )   PDF(224KB) ( )  
    本文由方程组Ax=f的系数矩阵AT(q,r)阵建立了BSSOR阵和块Jacobi阵特征值的关系式,从而对讨论了T(1,1)阵和T(1,2)阵BPSD迭代的收敛性和T(1,1)阵的最佳参数以及A为其它矩阵时PSD迭代的收敛性。推导简单有趣。
    利用物理原理研究非线性对流项的差分格式
    蒋锦良
    1992, 9(S1): 497-497. 
    摘要 ( )   PDF(88KB) ( )  
    在流体力学与对流传热等物理问题中的基本方程可以归结为非线性对流扩散方程,在理论分析与数值求解这类物理问题时所遇到的主要困难是对流项的非线性。目前非线性对流扩散方程的数学理论还不完善,为了求解实际问题的需要,不得不采用一些近似的方法,其中最简单而常用的方法就是采用局部线性化,这样就可以把常系数线性对流扩散方程的数学理论与数值计算方法应用到非线性方程之中。但是,从计算实践发现,这种局部线性化方法的计算精度会受到影响,尤其是在对流占优问题中的计算精度会发生明显的下降。
    轴对称体热传导问题边界元分析中奇异积分的一种统一处理方法
    姚寿广, 朱德书
    1992, 9(S1): 498-500. 
    摘要 ( )   PDF(240KB) ( )  
    针对轴对称热传导问题边界元分析中奇异积分的处理提出了一种新的统一处理方法,在此基础上通过数值试验进一步给出了一种更为简洁的处理方法。数值计算及应用结果表明,本文提出的方法计算工作量小,便于形成具有良好结构性能和高效率的通用分析程序。
    神经网络方法在高能物理数据分析中的应用
    张子平
    1992, 9(S1): 500-500. 
    摘要 ( )   PDF(85KB) ( )  
    神经网络方法经过了七十年代和八十年代的沉寂之后,近几年的重新兴起,已越来越引起人们的广泛兴趣和研究,特别是向后传播法(back-propagation)训练模式的成功,被应用于高能物理实验的物理分析中,取得了较之别的方法好得多的结果,使用该方法,以喷注的形状来鉴别b夸克喷注,以及区分胶子喷注和夸克喷注的工作已见诸于报导。我们在LEP/L3实验的τ极化分析中,设计了一个三层的网络,用通过L3探测器模拟的M.C.事例作为输入,以向后传播法训练之,用以筛选τ的衰变模式,得到了相当满意的选择效率和本底水平,以L3的实验数据检验之,实验和M.C,的结果符合得相当好。
    物理问题的数值模拟及高精度差分格式
    傅德薰, 马延文
    1992, 9(S1): 501-505. 
    摘要 ( )   PDF(341KB) ( )  
    本文分析了模拟多尺度物理问题对差分格式的要求,并指出,为了正确模拟小尺度物理量,要求采用高精度差分格式。文中利用文献[1]作者提出的耗散比拟思想,给出了构造具有耗散机制的高精度紧致差分格式的方法。并采用时间为三阶精度,空间为四阶精度的差分方法求解二维欧拉方程,给出了激波反射问题的计算结果。所得结果表明,本文给出的能模拟小尺度物理量的高精度差分方法也能正确模拟宏观尺度的激波。
    一类与马丢函数有关的强烈振荡积分的数值计算
    龙永兴
    1992, 9(S1): 506-510. 
    摘要 ( )   PDF(305KB) ( )  
    简要介绍了计算马丢函数积分的有关问题,计算实践表明,本文的计算方法不仅获得了满意的结果,而且节省大量的计算时间。
    某一类非线性反应-扩散方程组的奇摄动
    王国灿
    1992, 9(S1): 511-511. 
    摘要 ( )   PDF(81KB) ( )  
    本文考虑带泛函项的反应-扩散方程组的奇摄动。
    Navier-Stokes方程的线性协调元方法
    王申林
    1992, 9(S1): 511-511. 
    摘要 ( )   PDF(81KB) ( )  
    求解Navier-Stokes方程边值问题。
    复杂体系的对称系数计算
    龚新高, 郑庆祺, 潘卫
    1992, 9(S1): 512-514. 
    摘要 ( )   PDF(232KB) ( )  
    本文介绍了一种求原子集团对称系数的办法,它可广泛使用于实空间的量子化学计算。作为例子,简要给出了该方法应用于C60分子的一些DV-Xα计算结果。
    解抛物型方程的分组显式(AGE)方法同几种分数步长法的比较研究
    吴声昌, 郭革
    1992, 9(S1): 515-519. 
    摘要 ( )   PDF(326KB) ( )  
    本文讨论了解抛物型方程的分组显式(AGE)方法,给出了该法和几种分数步长法的实验模型的数值比较结果。
    守恒双曲方程的一类二阶精度TVB格式
    徐树荣, 龚大平
    1992, 9(S1): 519-519. 
    摘要 ( )   PDF(84KB) ( )  
    本文通过修改[1]中一类二阶精度TVD格式中的限制函数,得到一类二阶精度TVB格式,克服了TVD格式在极值点精度退化问题。数值实验表明该格式仍保留了原格式激波分辨率高的优点。
    计算衍射学中的Fourier方法
    何保利
    1992, 9(S1): 520-520. 
    摘要 ( )   PDF(79KB) ( )  
    计算衍射是计算物理学的最新分枝之一。本文主要论述计算衍射Fourier方法的的基本理论、数学物理模型及其在三维重构分析中的应用。在计算衍射Fourier方法中,计算衍射晶体点阵的每一组平面的衍射波可以分解为一个基波和若干谐波。基波给出一级Bragg反射,谐波给出次级Bragg反射。
    两相渗流问题的数值新方法及其理论分析
    梁栋
    1992, 9(S1): 521-525. 
    摘要 ( )   PDF(339KB) ( )  
    本文针对两相不混溶驱动和两相混溶驱动问题,提出了两类处理对流项的新方法:粘性分离算法和任意网格上的迎风算法,给出了理论分析,模型问题的试算得到满意的结果。
    变工况下离心三元叶轮内部二次流的数值分析
    席光, 王尚锦
    1992, 9(S1): 526-526. 
    摘要 ( )   PDF(78KB) ( )  
    研究离心叶轮内二次流对控制主流速度分布、减小流动损失有重要意义。目前国外,通过求解时均N-S方程,以相对气流角为指标研究叶轮内部二次流分布的工作已有报道,但目前国内外发表的文献中,研究变工况下离心叶轮内部二次流的工作尚很少见。本文通过求解三维、不同压时均N-S方程和K-ε双方程湍流模型,直接以垂直于主流方向的速度分量(即二次流)为对象,研究三种不同工况下闭式离心三元叶轮内部的二次流分布,考查变工况对二次流的影响。控制方程及数值求解方法见文献[1]。
    太阳热温室非稳态设计理论研究
    赵炳俊, 胡广芝, 张剑, 胡广耀, 刘五奎
    1992, 9(S1): 527-529. 
    摘要 ( )   PDF(228KB) ( )  
    本文用非稳态导热理论对太阳热温室进行动态模拟,建立数学模型,利用计算机进行设计计算,分析参数对温室性能的影响,以河南地区为例,得出了具有广泛实用价值的最佳设计值。
    半导体器件非线性泊松方程的一种新差分格式
    王国彬, 严荣良, 任迪远
    1992, 9(S1): 530-530. 
    摘要 ( )   PDF(83KB) ( )  
    本文对半导体器件非线性泊松方程▽2ψ=q/c(n-p-D),提出了一种新的差分格式。
    求解Euler方程的人工压缩加熵强迫方法
    金保侠
    1992, 9(S1): 530-530. 
    摘要 ( )   PDF(83KB) ( )  
    将文献[金保侠,处理接触间断的一种线性场修正方法,第五届全国流体力学数值方法讨论会议文集,重庆1991年,第36~38页]的思想进一步推广到非线性场的计算,并引进适当的熵强迫处理,以此进一步提高激波的分辨率及稀疏波波头和波尾的计算精度,且不出现非物理解。
    不同期合闸诱发高过电压的模拟
    赖定文
    1992, 9(S1): 531-532. 
    摘要 ( )   PDF(153KB) ( )  
    不同期合闸产生的过电压比同期合闸产生的约高出20~39%,对电力系统的危害更加严重。本文继续探索不同期合闸诱发高过电压的内因。
    非均匀介质区域二维场的边界元法
    宋斌, 傅君眉
    1992, 9(S1): 533-533. 
    摘要 ( )   PDF(96KB) ( )  
    电磁场边值问题往往是由一个偏微分方程和基本边界条件,自然边界条件构成的,在一些简单情况下,可用解析方法求解,但对于大多数问题,则需应用数值方法求解。本文采用边界元法对含有多种介质区域的二维静态场分布进行了研究。边界元法相对于区域性解法(如有限元法)具有结果精确、计算量小、处理问题灵活等特点。本文求解问题的过程为:将整个区域分成许多小区域,在每个小区域内介质是均匀的:对每个小区域取边界,进行元素分割,建立代数方程组;考虑到内部相邻两区域共同边界上的边界条件,把关于小区域的方程组联立起来构成系统的方程组并求解。文中在进行边界剖分时采用了二次元素,解的精度很高。
    已知爆速的爆轰装配
    刘邦弟
    1992, 9(S1): 533-533. 
    摘要 ( )   PDF(96KB) ( )  
    这是一种爆轰的数值模拟方法。把爆轰波视为没有厚度的数学平面,爆轰波以外的区域用差分格格式计算,爆轰波波阵面上的力学量,用冲击波关系式ρ(D-u)=ρ0D, P=ρ0Du, e=Q*+0.5P(1/ρ0-1/ρ), 确定。
    解奇异微分方程的一个新方法
    牟宗泽
    1992, 9(S1): 534-534. 
    摘要 ( )   PDF(88KB) ( )  
    对于奇异微分方程有奇异解的情形,目前尚无较好的数值方法求解,通常把奇点两侧分为二个区域分别求解,这样就需给出更多的定解条件。本工作给的方法不仅对解有奇异的情形有效,而且也适合解为非奇异的情形。
    弹性动脉血管内血液流场的有限元分析
    高应才, 王光辉
    1992, 9(S1): 534-534. 
    摘要 ( )   PDF(88KB) ( )  
    已建立的一组血液流动和管壁振动相耦合的非线性的带有奇异性的偏微分方程组及相应的定解条件如下:N-S方程、管壁振动方程组。
    轨道卫星自然辐射环境计算
    张殿辉, 陈雨生, 华鸣
    1992, 9(S1): 535-535. 
    摘要 ( )   PDF(102KB) ( )  
    近地空间轨道卫星的自然粒子辐射环境有三大来源:银河宇宙射线(GCR)、太阳宇宙射线(主要是太阳耀斑质子)和地磁捕获粒子。空间自然粒子辐射环境对卫星体特别是其电子系统的影响已成为危害卫星系统长寿命正常工作的重要因素。本文给出了一套预报任一卫星轨道上自然辐射环境参数的基本方法和一些典型轨道的计算结果。
    冲击波经多级扩散室的衰减
    郝保田, 田宙
    1992, 9(S1): 535-535. 
    摘要 ( )   PDF(102KB) ( )  
    本文通过数值模拟的方法,计算了单一和如图所示的多级扩散室对冲击波衰减的影响,得到了和试验实测一致的结果。
    烧蚀对电离边界层的影响
    曹文祥, 高铁锁
    1992, 9(S1): 536-536. 
    摘要 ( )   PDF(77KB) ( )  
    本文采用了化学非平衡边界层方程组来计算高速再入飞行器边界内流场参数和各组份浓度,壁面防热材料为含有微量碱金属的碳碳材料和泰氟隆亲电材料。壁面烧蚀参数计算采用烧蚀壁面和气体边界层耦合求解得到的传递系数关联公式。早期的传递系数关联公式没有考虑化学反应,误差很大,他们由空气-空气同质引射的边界层解或由简单薄膜理论解经过修正而得到的。
    金属凝固过程中温度场的有限元分析
    陈永进, 王洪纲, 程赫明
    1992, 9(S1): 537-538. 
    摘要 ( )   PDF(144KB) ( )  
    利用由金属物理学的基本概念导出的含一级相变非线性热传导方程,在微时段内线性化,由变物性线性化泛函,建立了金属凝固过程中温度场的有限元方程并对轴对称问题的金属凝固过程中的温度场作了数值模拟和测试。得出比较满意的结果。
    多元约束非线性规划的区间方法
    王海鹰, 张乃良, 刘蕴华
    1992, 9(S1): 539-541. 
    摘要 ( )   PDF(218KB) ( )  
    本文给出了求解约束非线性规划的区间方法。利用连续区间扩张,函数判别法等,删除一切不包含极值点的单元。在迭代过程中同时求得最优解和最优值。数值算例表明,这种方法是有效的。
    CAD在物理和生物化学中的振子方程的极限环上的应用
    石磊, 郭治安, 孙寿文
    1992, 9(S1): 541-541. 
    摘要 ( )   PDF(86KB) ( )  
    从形式上看,似乎这是由于方程所遵从的边界条件限定的,正象我们给物体运动限定了轨迹之后,不同的能量状态也要经过这个轨道,只是通过轨道的动能大小不同而已。
    广义GALERKIN法的校正算法
    张大凯
    1992, 9(S1): 542-544. 
    摘要 ( )   PDF(186KB) ( )  
    本文对广义Galerkin法给出了一类"低次元求解,高次元校正"的高精度算法-插值校正和误差校正算法。
    反演微分方程系数的有限分析-区间迭代法
    方保镕, 王如云
    1992, 9(S1): 544-544. 
    摘要 ( )   PDF(76KB) ( )  
    本文对一类微分方程给出一种反演系数的新算法。这个方法计算简单,归结为大范围最优化问题,且可估计出计算误差。
    A-光滑正则化算子的残差原则
    贺国强
    1992, 9(S1): 545-545. 
    摘要 ( )   PDF(73KB) ( )  
    X、Y为实Hilbert空间,AXY的非退化紧算子,则方程Ax=y(1)是不适定的。
    超几何级数余项估计
    陈宝振, 陈鹏
    1992, 9(S1): 546-546. 
    摘要 ( )   PDF(68KB) ( )  
    超几何级数F(α,β;γ;z)定义为下列收敛的无限和F(α,β;γ;z)=(Γ(γ))/(Γ(α)Γ(β(Γ(n+α)Γ(n+β))/(Γ(n+γ)) (zn)/(n!)(1)这里,Γ(α)是伽马函数:|z| < 1;γ不等于零和负整数。
    超越方程求根的Monte Carlo方法
    宫野
    1992, 9(S1): 547-548. 
    摘要 ( )   PDF(160KB) ( )  
    本文提供了两种求超越方程根的M-C方法,它不需要选择初始近似值,就能求出全部实根。所以,对那些使用其他方法无效的问题,而用此方法往往是有效的。
    射线跟踪方法在蒙特卡罗计算中的应用
    张立吾, 张玉琴
    1992, 9(S1): 549-550. 
    摘要 ( )   PDF(170KB) ( )  
    在分析了用蒙特卡罗方法解粒子输运问题时费机时的主要方面后,提出了一个提高计算速度的方法-射线跟踪法。这种方法在实际工作中已得到应用,取得了明显效果。
    用M-C方法和差分法混合编码研究电子束辐照效应
    周南, 丁升
    1992, 9(S1): 550-550. 
    摘要 ( )   PDF(91KB) ( )  
    本文利用Monte-Carlo方法与差分法相结合的方法对电子束辐照Al靶板的热-力学效应进行了数值模拟。
    无定形中靶位移原子深度分布的Monte Carlo模拟计算
    邵其, 潘正瑛, 霍裕昆
    1992, 9(S1): 551-554. 
    摘要 ( )   PDF(297KB) ( )  
    本文介绍一个模拟具有给定动能的离子在固体中传输及引起靶原子位移过程的Monte carlo程序,包括工作原理、物理模型和计算实例等。本程序采用两体碰撞近似研究α粒子轰击无定形靶的位移原子深度分布,并分析了位移原子深度分布与固体中能量沉积的关系。
    两种量子蒙特卡罗方法
    王忆平, 唐琳, 王玉清
    1992, 9(S1): 554-554. 
    摘要 ( )   PDF(80KB) ( )  
    对两种典型的量子蒙特卡罗方法(路径积分量子蒙特卡罗方法和扩散量子蒙特卡罗方法)进行了分析、比较和讨论。
    Monte Carlo方法计算晶体析出相的X射线发射强度
    何延才, 杨进, 陈家光
    1992, 9(S1): 555-559. 
    摘要 ( )   PDF(423KB) ( )  
    本文提出了用Mont Carlo模拟计算晶体中析出相X射线发射强度的理论和方法:严格计算晶体对析出相中电子散射的影响;用绝对法计算析出相的X射线深度分布函数;给出了析出相中X射线吸收及二次荧光计算公式。作为实例,对不同加速电压下Au、Cu基体中Cr-Fe-Ni合金析出相的X射线发射强度作了计算,并得出一些规律。本工作为析出相的电子探针定量显微分析(EPMA)研究奠定了基础。
    蒙特卡罗方法通用软件展望
    裴鹿成
    1992, 9(S1): 560-562. 
    摘要 ( )   PDF(229KB) ( )  
    这是一篇关于蒙特卡罗方法通用软件展望的概述性论文。内容包括:集团抽样方法;零方差技巧;自改善方法;向量蒙特卡罗方法。
    蒙特卡罗方法发展中的若干问题
    裴鹿成
    1992, 9(S1): 563-566. 
    摘要 ( )   PDF(332KB) ( )  
    这是一篇关于蒙特卡罗方法发展中若干问题的概述性论文。内容包括:深穿透问题;估计量无界问题;Metrpolis抽样方法的局限性和相关性问题;随机误差干扰问题和随机方程;通用软件的智能化和向量化问题。
    次级电子激发过程的Monte Carlo模拟计算
    韦亚一
    1992, 9(S1): 566-566. 
    摘要 ( )   PDF(89KB) ( )  
    次级电子发射是各种电子探测器和电子倍增器的工作基础。原初电子在次级发射体中激发次级电子的过程具有复杂性和随机性。过去的工作大多局限于定性描述和半定量估算,本工作用MC方法对原初电子在次级发射体内激发次级电子以及次级电子向表面运动直至逸出发射体的整个过程进行模拟:
    子母弹对目标毁伤概率的Monte-Carlo计算和分析
    黄启晋, 刘兴平
    1992, 9(S1): 567-568. 
    摘要 ( )   PDF(152KB) ( )  
    用二、三维Monte-Carlo方法计算子母弹对目标的毁伤概率,给出均匀抛撒时二维精确解析公式,计算结果表明和解析解是一致的。
    一个研究聚变中子辐照损伤的蒙特卡罗程序
    高兴华
    1992, 9(S1): 569-569. 
    摘要 ( )   PDF(107KB) ( )  
    自开展聚变堆研究以来,由于中子源强不够,直接用14MeV中子源研究聚变中子辐照损伤很难实现,人们寻找了各种代替手段,比如电子模拟、带电粒子模拟和计算机模拟等方法来研究聚变中子的辐照损伤问题,由于大型计算机的飞速发展,计算机模拟聚变中子辐照损伤的方法越来越受到人们的重视,本文介绍了用计算机模拟聚变中子辐照损伤的蒙特卡罗程序和研究方法,这个程序对其它堆型的中子辐射损伤同样也是适用的。
    非晶碳氢膜晶化凝聚的蒙特卡罗模拟
    刘启海, 杨功能, 刘志远, 何大韧
    1992, 9(S1): 569-569. 
    摘要 ( )   PDF(107KB) ( )  
    在材料科学领域中,关于非晶碳氢膜的研究,三十年来国内外一直处于活跃时期,非晶碳氢(a-CH)膜系类金刚石结构,有着许多特殊性能,电阻率高,红外及可见区透明,电介性能好,硬度大,抗蚀性强,还可掺杂。
    组分夸克模型中随Q2演化的强子的Parton分布函数
    钱路
    1992, 9(S1): 570-572. 
    摘要 ( )   PDF(221KB) ( )  
    用牛顿一拉夫逊方法求出了Q02时价子中的Parton分布函数,且采用把微分积分方程组形式的QCD演化方程离散为常微分方程组的方法,求出了任意Q2 > Q02时价子中的Paerton分布函数,最后通过卷积,得到了组分夸克模型中随Q2演化的强子中的Parton分布函数。
    用REDUCE语言计算非谐振子
    姜绍周
    1992, 9(S1): 573-577. 
    摘要 ( )   PDF(274KB) ( )  
    REDCE语言是用计算机推导解析公式的有力工具,本文用REDUCE语言计算了高次非谐项对非谐振子的效应。
    在BUU模型中引入集团效应的一种方案
    王志松, 霍裕昆
    1992, 9(S1): 578-578. 
    摘要 ( )   PDF(86KB) ( )  
    重离子核反应的Boltzmann-Uehling-Uhlenbeck(BUU)输运模型描写约化的单粒子相空间密度分布函数的平均演化行为。基于试验粒子方法[1],已经发展了众多版本的动力学计算程序[2],这些程序被应用于重离子核反应中粒子发射谱如硬光子、核子、Ⅱ介子等的计算,并与实验对比。结果表明,BUU模型和程序给出的能谱的高能端("硬尾")低于实验值,在这些计算方案的基础上,本文讨论一种能导致高能粒子发射的新机制。
    SNENM/RTSP: 一个用于压水堆堆芯瞬态分析的三维节块和子通道计算程序
    程平东, 邵天伟, 赵慧敏, 周全福, 申森
    1992, 9(S1): 579-582. 
    摘要 ( )   PDF(316KB) ( )  
    压水堆堆芯的运行瞬态和事故瞬态是伴随着复杂的热工水力过程的三维中子动力学问题。为了改进现代压水堆的安全设计,提高核电厂运行的机动性和经济性,对这类问题作直接的三维分析是必要的。然而,这种分析常常要花费很大的计算代价。因此,发展既有适当精度又能快速求解的三维模型和相应的计算机程序,就成为使这种计算分析具有工程实用价值的重要途径。SNENM/RTSP程序正是为了满足这种要求而研制的。
    混合堆放射性计算与FDKR程序研制
    冯开明, 阳彦鑫, 黄锦华
    1992, 9(S1): 582-582. 
    摘要 ( )   PDF(94KB) ( )  
    裂变产物、锕系元素、活化产物和氚是聚变-裂变混合堆内主要的放射性源。准确计算出各源项对放射性的贡献和危害对于混合堆设计的环境安全评价分析、结构材料的选择以及长寿命放射性核废料的处置是十分重要的。
    燃耗计算的一种有效方法
    钟文发, 罗经宇, 张瑞茵
    1992, 9(S1): 583-586. 
    摘要 ( )   PDF(304KB) ( )  
    本文介绍了反应堆物理计算中燃耗计算的一种方法,通过宏观截面随燃耗变化的参数表、插值方法及中子扩散方程的求解完成燃耗计算,计算结果与Shearon Harris堆的结果作了比较,它们基本上一致,该方法具有使用方便,计算时间较短,能满足一定精度要求等特点。目前该方法已用于实际的堆物理计算中,为工程设计提供了一种有力的工具。
    离散节块输运方法在二维(r,z)柱坐标系中的应用
    祝学华, 谢仲生, 吴宏春
    1992, 9(S1): 586-586. 
    摘要 ( )   PDF(80KB) ( )  
    离散节块输运方法(DNTM)在中子输运方程的数值解中获得成功的应用。本文讨论把DNTM方法推广到二维(r,z)柱坐标中的应用,并考虑到各向异性散射的普遍情况。
    大亚湾核电站低放废水低温废热排放和卷载效应的数值模拟
    张春粦, 黄祖湛, 孔令丰, 彭建伟
    1992, 9(S1): 587-592. 
    摘要 ( )   PDF(408KB) ( )  
    用ADI法计算了大亚湾核电站低温废水排放、低温废热排放在周围海域的放射性浓度场,温废场分布和吸水口的卷载效应。
    快速燃料管理程序包
    王英明
    1992, 9(S1): 592-592. 
    摘要 ( )   PDF(103KB) ( )  
    本程序包是由多功能数据处理、MUFT型54群快谱、THERMOS型30群栅元燃耗热谱、穿透几率组件均匀化、参数拟合和两维粗网结块燃料管理等六个程序组成。首先由数据处理程序根据堆芯栅元的几何、材料和结构参数计算出快谱和热谱程序所需的输入参数,再由这两个多群谱程序计算出各类栅元不同燃耗时刻的少(2~4)群截面参数,从而生成具有文字说明的标准格式栅元截面库。组件均匀化程序再根据组件的栅元分布特性,从栅元库中选出所需栅元的少群截面参数,作组件均匀化计算并生成格式上与栅元库完全相同的少群组件截面库。最后由两维粗网结块燃耗程序根据堆芯的组件布置,从组件库中选出所需组件的四群微观截面参数,解燃耗方程计算各燃耗步各元素的核密度并生成对应燃耗步的组件宏观截面,用它们作临界、燃耗和燃料管理计算。如果我们把这种燃耗计算称为“真燃耗”计算,为了节省机时,提高计算速度,可以用拟合程序从组件库中取出所需组件的随燃变化的宏观截面作拟合计算,得到各组件修正宏观截面随燃耗变化的拟合系数,这样在作两维燃耗计算时不用解燃方程,用拟合系数计算出各燃耗步的组件宏观截面,通常把用这种宏观截面进行的燃耗计算称为"伪燃耗"计算,伪燃耗比真燃耗计算节省机时约1/3。
    轻水堆物理分析方法的现状与发展
    章宗耀
    1992, 9(S1): 593-596. 
    摘要 ( )   PDF(386KB) ( )  
    本文从栅格物理和堆芯物理两大方面评述了轻水堆物理计算方法和程序的国内外现状,以及今后发展的主要动向。
    利用聚变-裂变混合堆转化放射性核废料的燃耗计算
    吴宜灿, 邱励俭
    1992, 9(S1): 596-596. 
    摘要 ( )   PDF(97KB) ( )  
    核废料的处置问题已成为核工业发展人们最关心的问题。"分离与转化"作为处置各种核设施尤其是裂变反应堆产生的高放射性核废料的一种方法具有永久处置的优点,同时能更加充分地利用铀钍资源。
    谱计算中的曲率反馈
    经荥清, 单文态, 罗经宇
    1992, 9(S1): 597-601. 
    摘要 ( )   PDF(344KB) ( )  
    本文研究了反应堆物理计算中的谱计算的曲率反馈问题。研究了一种实用的作谱计算曲率反馈的方法。通过计算讨论了曲率反馈对谱、少群常数的影响,曲率反馈迭代的收敛情况及谱区划分的影响。该方法已用于HTR-10 MW试验堆的物理设计中。
    He冷却混合堆第一壁的稳态、瞬态温度场计算
    王学人, 吕晓兰, 黄锦华
    1992, 9(S1): 602-606. 
    摘要 ( )   PDF(384KB) ( )  
    用有限差分热传导程序HEATING-5~1模拟和计算了混合堆He冷却第一壁的热工水力问题,介绍了求解稳态、瞬态热传导方程的数值方法。为了解决程序中的热平衡方程没有包括流体的流动项问题,先用解析法计算冷却剂在不同流动位置的温度分布,用作程序的输入参数。同时,还用热平衡方程带有强迫对流项的有限元热传导程序[2~4]计算了同一问题。其结果表明,用HEATING-5程序与AYER模拟冷却剂的强迫对流传热的计算结果符合得较好。
    穿透几率法在轻水堆燃料组件计算中的应用
    黄群英, 尹帮华, 谢仲生
    1992, 9(S1): 606-606. 
    摘要 ( )   PDF(92KB) ( )  
    由于轻水堆燃料组件结构的强烈非均匀性,应用扩散理论将带来很大的误差,而SN方法或蒙卡方法则耗时太多,不经济。碰撞几率法属于体积耦合型积分输运方法,各子区之间通过碰撞几率互相耦合,计算量仍然很大。这些年来发展的面流耦合的穿透几率法,各子区仅仅通过界面中子流及中子通量和相邻子区耦合,大大减少了计算时间,被认为是最经济有效的方法。
    放射性衰变测量区间中计数率的转换和校正
    李伟生, 李文新, 尹新民
    1992, 9(S1): 607-610. 
    摘要 ( )   PDF(338KB) ( )  
    提出了一个放射性衰变测量区间中计数率的转换公式,和带计数率校正的拟合多组分衰变曲线的分析方法及迭代程序。该程序已成功地应用于不同模拟衰变数据和有关实验数据的分析,并与通常使用的不带计数率校正的分析方法进行了比较,结果令人满意。
    放化工程屏蔽设计程序化的研究
    张洪猷
    1992, 9(S1): 610-610. 
    摘要 ( )   PDF(99KB) ( )  
    放化工程所处理的物料中含有多种放射性核素,其γ射线具有复杂的能谱,许多辐射源也不能等效成简单的点源,本课题就是为了适应放化工程中屏蔽设计的需要,改进了屏蔽计算方法和模式,将一般使用的能谱竞争线法发展成解析KT法,研制了相应的ESKT计算机程序。
    三温近似辐射传输计算的区域分解方法
    冯庭桂
    1992, 9(S1): 611-614. 
    摘要 ( )   PDF(253KB) ( )  
    本文讨论辐射传输计算的三温、单温区域分解方法,给出分解区域介面的差分处理方法和差分方程的整区域LU分解迭代方法,并讨论了迭代收敛问题。
    高分波碰撞强度的研究
    方泉玉, 蔡蔚, 李萍, 沈智军, 邹宇
    1992, 9(S1): 615-618. 
    摘要 ( )   PDF(293KB) ( )  
    本文通过研究Ne-Like Ge22+与电子碰撞激发的一些重要跃迁过程,对分波展开方法中的高分波作三种近似处理: 平面波近似,库仑Bethe近似以及下降等比级数近似。并且对这些近似进行了比较。从而解决分波收敛性的问题。
    Al10+离子的低能电子碰撞激发截面
    邹宇, 方泉玉
    1992, 9(S1): 619-620. 
    摘要 ( )   PDF(182KB) ( )  
    2态强耦合(2cc)和5态强耦合近似方法被用于计算Al10+离子的2s-2p,2s-3p和2s-3s跃迁的电子碰撞激发截面,分析了间接耦合的作用。与扭曲波近似结果比较可以看出,对于共振型2s-2p跃迁耦合效应不是很重要,但对2s-3s和2s-3p跃迁耦合效应是相当重要的。
    激光热烧蚀问题体能源体汽化模型的数值方法
    陶应学, 关吉利, 沈隆钧, 刘成海
    1992, 9(S1): 621-623. 
    摘要 ( )   PDF(220KB) ( )  
    本文对体汽化模型,根据等温平衡汽化原理,在相变区域以外采用Von Neumann计算方法;在相变区内建立特定的计算格式。该方法已应用在我们的总体程序LHAP-1DVG中,经过大量计算业已证明该计算方法是可行的。
    连续波氧碘化学激光器增益特性的模拟计算与分析
    陈丽吟, 楚泽湘
    1992, 9(S1): 624-628. 
    摘要 ( )   PDF(377KB) ( )  
    本文用一个简化模型从理论上系统地研究了CW OICL的增益系统及其与组份、压力、流速的关系,并得出了一些规律性的结果。研究的范围是:[O2*]/[I2]=10~200,[O2*]/[O2]=2~9,P=1~100 torr,U=5-100 m/s。详细探讨了流动守恒方程组对计算CWOICL的重要性。指出提高流速几乎不影响峰值增益,而有助于按与流速成比例地扩大可用增益区。
    三成分模型的气球模本征方程的数值计算
    隋国芳, 周海麟
    1992, 9(S1): 629-631. 
    摘要 ( )   PDF(218KB) ( )  
    用数值方法求解了三成分模型的气球模本征方程,研究了高速粒子对气球模不稳定性的作用。理论计算结果能够成功地解析等离子体物理实验中所观察到的物理现象。
    等离子体Cherenkov Maser的数值模拟
    王秩雄, 段耀勇, 陈雨生, 范如玉
    1992, 9(S1): 632-634. 
    摘要 ( )   PDF(244KB) ( )  
    用2(1)/(2)维、全电磁、相对论性的粒子模拟程序模拟了圆柱状相对论电子束(REB)与相同半径、内充满磁化等离子体的谐振腔相互作用,腔内激发准纵向的电子磁振荡,波长在荡,波长在微波波段,数值模拟结果和线性分析一致。
    利用映射方法研究带电粒子束在非线性薄透镜周期场中的运动
    郁庆长
    1992, 9(S1): 634-634. 
    摘要 ( )   PDF(92KB) ( )  
    利用映射方法可以讨论非线性薄透镜周期场中带电粒子束运动特性。设粒子在第i个透镜的r值为ri,在第i个透镜入射方的r'值为ri',此处r为径向坐标,r'=(∂r)/(∂z),z为轴向坐标,则有ri+1=ri+Lri+1',ri+1'=ri'-αri-βri3αβ为透镜参数,α > 0,L为透镜间距。
    欧姆加热等离子体的能量输运
    丁宁, 何启兵, 黄林, 邱孝明
    1992, 9(S1): 635-635. 
    摘要 ( )   PDF(108KB) ( )  
    尽管人们已经做了大量的欧姆加热等离子体的工作,但都没有明确给出安全因子q(r)的时间演化,也没有提出准确描述发生锯齿振荡q(0) < 1时,维持q(0)在1附近的与实验相符的模型。
    “双靶对接”类氖锗X光激光高增益实验结果的数值模拟
    张国平, 盛家田, 杨明伦, 张覃鑫, 邵云峰
    1992, 9(S1): 635-635. 
    摘要 ( )   PDF(108KB) ( )  
    90年4、5月,我们在神光装置上实现了"双靶对接"类氖锗X光激光的高增益,增益长度积达到15左右。为了分析实验结果和为下一轮的实验进行理论设计,我们编制了模拟X光激光在等离子体增益区中传播放大规律的XBY程序,和JB19程序、ALPHA程序一起组成一个完整的程序系列,来整体模拟实验室X光激光实验。
    激光产生X光的数值模拟研究
    张钧, 陈光南, 常铁强, 裴文兵, 尤希文, 隋成之, 徐少则, 古培俊, 章实, 张兴宏
    1992, 9(S1): 636-636. 
    摘要 ( )   PDF(102KB) ( )  
    本文简要地描述了激光X光转换涉及的物理问题和重要物理过程,概括地介绍了数值模拟程序的规模。激光聚变有两种途经,用激光直接驱动氘氚靶丸,使其达到高温高密度,必须实现理想的球面内爆,这就要求激乐靶表面要有很高的均匀性。但是由于激光的相干性达到上述要求是很困难的。通常称这种方式为激光直接驱动。如果用激光产生的非相干的X光辐照靶丸表面,就可达到极好的均匀性,能够实现理想的球面内爆,这被称为激光间接驱动。激光X光转换是间接驱动惯性约束聚变的重要课题和关键技术。数值模拟和地下核爆的小囊试验指出,实现点火和中等增益的热核聚变要求激光器的输出能量5-10MJ,功率~100TW。这将是本世纪末和下一世纪的目标。但是现有高功率激光器产生的X光在近期有多方面的应用。如高功率X光激光的泵浦源,在实验室中辐射流体力学现象的研究,以及高温高密度靶围内物质性质的研究等等。因此研究激光产生X光的物理规律不仅有远景目标,而且有重要的现实意义。
    强流电子束泵注入锁定XeCl准分子激光的理论模拟
    高占鹏, 杨国龙, 屠琴芬, 张建泉
    1992, 9(S1): 636-636. 
    摘要 ( )   PDF(102KB) ( )  
    本文模拟了在强流电子束泵浦条件下,Ne/Xe/HCl混合气体的反应动力学过程以及注入锁定非稳腔的激光生成、放大和输运过程,模拟结果表明,在激光生成阶段,注入种子光源后,明显地改善了激光的输出质量,使自发辐射部分受到了明显的抑制。得到了在不同注入强度、注入时间和放大倍数与激光输出的光谱特性、能量以及各种效率等参量的关系。
    强激光引起大气击穿的数值模拟
    吴乃清, 陈放, 王心正
    1992, 9(S1): 637-637. 
    摘要 ( )   PDF(94KB) ( )  
    脉冲高功率激光与大气非线性相互作用而产生的大气击穿效应,在激光武器系统研究中是一个很重的课题。本文从光致电场引起的多光子电离反向韧致辐射和串级电离过程等理论出发,考虑到大气主要成分氮和氧分子的分离、振动、电子的激发和电离,扩散和附着等机理及其反应速率,在此基础上,采用描述粒子守恒和能量守恒的有限量子波尔兹曼方程,再利用差分和数值积分相结合的方法,求解微分-积分方程,获得了满意的结果。
    强激光辐照铝靶的数值模拟研究
    王心正, 曹莉华
    1992, 9(S1): 637-637. 
    摘要 ( )   PDF(94KB) ( )  
    强激光与物质相互作用的重要特征是产生高温稠密的等离子体。其中非平衡辐射输运是一个重要的物理过程。目前国内外已有的激光打靶模拟程序多数仍然使用限流扩散近似计算辐射输运。最大熵变Eddington近似对具有复杂光学性质的辐射输运问题能给出更精确的结果。我们用该方法编制了一个非平衡辐射流体力学激光打靶模拟程序,得到了令人满意的结果。
    解相对论Dirac方程的低能电子弹性散射截面分波法
    谭震宇, 何延才
    1992, 9(S1): 638-638. 
    摘要 ( )   PDF(79KB) ( )  
    电子与固体相互作用是现代表面分析技术如SEM、EPMA及SAM等的物理基础。入射电子能量为几个keV或更低时,描述中能电子弹性散射的Rutherford截面不适用,需应用更为严格的量子力学分波法(PWM)。分波法计算电子弹性散射截面有解非相对论Schrodinger方程的pendry截面以及解相对论Dirac方程的Mott截面,但关于Mott截面的计算方法鲜见报导,且计算分析Mott截面、pendry截面及Rutherford截面之间的差异对描述低能电子在固体中散射物理过程十分重要。本文给出一个计算Mott截面的方法,由引入一变换使散射截面中相移计算由解复杂的Dirac方程转化为解一简单的一阶方程,应用Thomas-Fermi-Dirac(TFD)势函数,导出初值及相移计算式。
    一维地震勘探反问题的一种算法
    李志斌
    1992, 9(S1): 639-640. 
    摘要 ( )   PDF(142KB) ( )  
    对地震勘探中反演声阻抗的一维反问题提出了一种算法。先将反问题进行差分离散化为某一非线性方程组,然后对这一非线性方程组进行Newton线性化,再采用LMS迭代。最后给出一模拟算例,验证了本文算法的有效性。
    组件计算软件包HATRAP及其在秦山核电站临界计算中的应用
    沈炜, 尹邦华, 崔士钢, 程平东
    1992, 9(S1): 641-641. 
    摘要 ( )   PDF(103KB) ( )  
    随着粗网法和现代节块法的迅速发展,快速精确地确定燃料组件内中子通量分布以产生组件均匀化参数显得愈来愈为重要,这是因为粗网法和现代节块法是建立在假设堆芯粗网节块可以用等效均匀化参数描述的基础上的。对于全堆芯空间计算,由组件群常数的计算偏差引起的功率误差远远大于粗网、节块法方法本身所引起的误差。快速准确地确定组件群常数也是LWR物理设计及堆芯燃料管理的一个重要问题,本文介绍
    用MCNP程序计算核燃料废包壳缓发裂变中子形成的热中子通量
    高产峰, 张孝泽, 王瑞宏
    1992, 9(S1): 642-642. 
    摘要 ( )   PDF(68KB) ( )  
    核燃料废包壳中残存不少裂变物质,确定出它的含量,是核保障技术中的一个重要问题。实验测量的原理是:从加速器发出的14MeV中子,打入废包壳蓝中,在慢化过程中使其中的裂变核素产生出缓发裂变中子。通过探测到的这些缓发裂变中子,可以确定包壳中裂变物质的含量。为使探测有效,需要设计出使缓发裂变中子及其在探测器中形成的热中子通量最大的装置。为此,采用Monte Carlo方法来模拟各种实验装置安排下,缓发裂变中子所形成的热中子通量[1]
    Radon变换的计算及其地球物理应用
    李新祥
    1992, 9(S1): 643-646. 
    摘要 ( )   PDF(324KB) ( )  
    根据实际地震数字信号处理的要求,给出了Radon变换及其反变换的一种算法并具体地给出此变换在地球物理地震勘探数字处理中的应用原理和叠后去噪处理效果。
    恒星子午观测的模拟方法及其应用
    朱紫, 徐劲
    1992, 9(S1): 647-648. 
    摘要 ( )   PDF(158KB) ( )  
    本文通过恒星子午观测过程的分析,给出了输出信号所应满足的数学关系式。考虑到输出信号噪声成分的统计特性,得到恒星光子计数观测的计算机模拟方法。利用模拟方法,讨论和解决了两个子午观测中的实际问题。
    广义相对论对超新星星核坍缩的影响
    王贻仁, 姚进, 汪惟中
    1992, 9(S1): 649-653. 
    摘要 ( )   PDF(391KB) ( )  
    随着物态方程的发展,超新星星核可坍缩至核密度以上。广义相对论(简记GR)显示出重要的影响。GR效应使内星核坍缩至更高的密度,陷入更深的势阱,释放更大的能量。但GR的内星核较非广义相对论(简记NR)的内星核小。冲击波形成地点离星核边界较远,冲击波传输途径中光裂解损耗增加。中心密度达最大时,GR过程和NR过程两坍缩速度曲线相交点(静止质量坐标略大于NR声速点坐标)以内,GR坍缩速度分布高于NR的相应值。而于相交点之外,GR坍缩速度、密度分布皆低于NR相应值。这对于反弹冲击波的传播是有利的。
    椭圆波导的本征模序列
    张善杰, 沈耀春
    1992, 9(S1): 654-654. 
    摘要 ( )   PDF(107KB) ( )  
    椭圆波导在卫星通迅、雷达的天馈系统以及加速管的研究中具有广泛应用,其截止波长λc的计算是研究波导中本征模及共传播特性和进行波导设计的关键问题。
分享到: