Please wait a minute...

过刊目录

    1997年, 第14卷, 第S1期
    刊出日期:1997-12-25
    上一期    下一期

    论文   
    全选选: 隐藏/显示图片
    论文
    珍惜过去,开拓进取,扎实工作,再图发展——中国核学会计算物理学会第三届理事会工作报告
    符鸿源
    1997, 14(S1): 388-390. 
    摘要 ( )   PDF(120KB) ( )  
    十五年前,作为与国防尖端事业发展密切相关的计算物理学会从军工系统迅速地扩展到民用系统,进入国民经济发展的许多重要领域。计算物理学会正是适应这种形势发展的需要,得到了国防科工委的大力支持,于1982年8月23日在北京成立,隶属中国核学会领导。
    计算物理学会活动大事记(续)(1992-1996)
    张锁春
    1997, 14(S1): 391-395. 
    摘要 ( )   PDF(192KB) ( )  
    [47]1992年1月中旬发出第三届年会的第二次征文通知,顺告理事和有关单位,推荐新一届理事候选人,限于四月底将名单报来,同时也与各个专业组和有关重点单位打招呼,请他们推荐合适的人选。
    石油射孔弹的数值模拟
    杭义洪, 陈贤林, 董振, 于志鲁, 范中波, 耿沪男
    1997, 14(S1): 404-406. 
    摘要 ( )   PDF(173KB) ( )  
    介绍了一种能有效模拟石油射孔弹作用全过程的数值与解析耦合的计算方法及相应软件,给出了计算结果及图象,与实验结果吻合较好。实际应用表明,该软件具有较高的效费比。
    计算石油地质中的一些数学问题
    吴声昌, 袁益让, 白东华
    1997, 14(S1): 407-409. 
    摘要 ( )   PDF(132KB) ( )  
    主要讨论了在计算石油地质中有关油气资源评估,电法测井分析和不稳定试井分析的一些数学问题以及求解的途径和方法。
    激光靶耦合精密计算中的一些物理与计算数学问题
    常铁强
    1997, 14(S1): 410-410,409. 
    摘要 ( )   PDF(97KB) ( )  
    1精密计算与Benchmark实验高强度激光和靶的耦合是一个非常复杂的问题,它包括激光的传播和吸收,激光烧蚀产生高温等离子体,靶物质的非平衡电离与激发,X光的发射和输运,激光驱动的高压冲击波的传播等。但是,这个问题是激光惯性约束聚变研究的基础,确切地了解有关规律和精确地数值模拟十分重要。
    计算物理的发展和艾滋病数学模型的提出
    张锁春
    1997, 14(S1): 411-412. 
    摘要 ( )   PDF(97KB) ( )  
    (一)乍看起来,计算物理与艾滋病数学模型毫无关系。此话还得先从计算物理的由来说起。
    可压缩平面混合层稳定性数值计算
    王强, 傅德薰, 马延文
    1997, 14(S1): 413-416. 
    摘要 ( )   PDF(192KB) ( )  
    通过高阶差分求解三维可压缩稳定性方程,数值研究了时间发展可压缩平面混合层的稳定性特征。以对流马赫数为压缩性参数,以Reynolds数为粘性参数分析了二者对计算流场稳定性的影响。比较了二维和三维,无粘和粘性扰动波等情况。部分计算同相关文献进行了对比,取得了一致的结果。
    电离辐射在半导体器件中引起的光电流的数值模拟
    黄流兴
    1997, 14(S1): 417-418,416. 
    摘要 ( )   PDF(132KB) ( )  
    电离辐射在半导体中引起光电流问题是半导体器件辐照效应研究中一个带普遍性的课题,首先建立描述半导体中光电流过程的普遍性的物理模型,而后模拟一个实际pn结器件的光电流响应,给出模拟结果。
    氘氚燃烧等离子体的热不稳定性计算
    龙永兴, 石秉仁, 牟宗泽
    1997, 14(S1): 419-420. 
    摘要 ( )   PDF(93KB) ( )  
    由于氘氚反应率与离子温度有很强的依赖关系,氘氚燃烧等离子体是热不稳定的。讨论了描述稳态氘氚燃烧等离子体功率平衡方程的非线性特征。利用打靶方法,可以求得稳态温度分布和关于温度微扰方程的本征值。可是,由上述方法得出的平衡解会与随时间变化方程在t→∞时的解不一致。这是一类未解决非线性数学问题的一个典型例子。
    X射线激光全息的理论研究
    李敬宏, 彭惠民, 张国平, 雷广玉
    1997, 14(S1): 421-423. 
    摘要 ( )   PDF(144KB) ( )  
    在对实验室X射线激光空间相干性理论研究的基础上,提出了X射线激光全息的一种可行的原理性光路设计,对无透镜傅立叶变换全息的记录和再现过程进行了理论模拟。
    瞬态热传导问题的一种单纯边界元算法
    姚寿广
    1997, 14(S1): 424-425,423. 
    摘要 ( )   PDF(135KB) ( )  
    给出了一种二维瞬态热传导问题的单纯边界元算法,它完全在边界上离散数值求解,通过数值算例验证表明该方法充分体现了边界元法的独特优点。
    脉冲强激光产生的激波在靶材中的传播过程
    蔡宗良, 张建泉, 刘峰, 强希文
    1997, 14(S1): 426-428. 
    摘要 ( )   PDF(137KB) ( )  
    采用一维流体动力学模型,以烧蚀压力作为已知的边界条件,分别以解析和数值计算方法,描述了短脉冲强激光辐照固态铝、铜和碳酚醛靶材时,烧蚀加载在靶材中产生的激波的增强、持续和衰减过程,给出了激波峰值压力和波形等特征参数的时空关系。
    三维圆柱几何格林函数节块法中子扩散计算
    胡永明, 赵险峰
    1997, 14(S1): 429-430. 
    摘要 ( )   PDF(100KB) ( )  
    发展了中子扩散计算三维圆柱几何格林函数节块法。首先通过横向积分将中子扩散方程化为三个互相耦合的一维偏通量方程。对于径向偏通量方程,将径向扩散微分算符分解为平板几何的扩散微分算符和一个修正项之和,将修正项移到方程右端作为修正源项,这样,三个方程都化为平板几何的一维方程形式。再借助平板几何第二类边界条件格林函数,对圆柱几何相应体源作积分,建立偏通量积分方程。对于修正源项,通过分部积分方法将偏通量导数项转化为对格林函数的求导。通过源迭代法求解方程。基准计算表明,该计算精度高、速度快,可成为三维圆柱几何堆芯设计和燃料管理计算的有效方法。
    关于波函数模方逐步归一化修正的注记
    周忠源, 吴承埙, 丁培柱
    1997, 14(S1): 431-432. 
    摘要 ( )   PDF(108KB) ( )  
    采用R-K法和逐步归一化修正计算了一个简单例子,结果表明,归一化修正虽能使波函数模方保持归一,但其它物理量如能量和态的布居几率仍不合理,甚至变得更坏。
    N2双原子系统经典轨迹的辛算法计算
    石爱民, 母英魁, 丁培柱
    1997, 14(S1): 433-434. 
    摘要 ( )   PDF(107KB) ( )  
    采用辛算法和R-K法计算了N2双原子系统的经典轨迹。结果表明,辛算法保持系统的辛结构和固有性质不变,与实验和理论一致;R-K法则与实验和理论不符。
    增益与吸收对X光激光传播放大的影响
    吴建周
    1997, 14(S1): 435-437. 
    摘要 ( )   PDF(145KB) ( )  
    在光沿直线传播的假设下,研究了X光激光在等离子体中的传播与放大。对类氢炭等离子体,除非光线不断通过强吸收区,增益曲线的吸收部分不会严重影响X光激光的输出能量。
    多水平方法在二维山丘地形大气污染物扩散的数值模式的应用
    周振中, 陈历舒
    1997, 14(S1): 438-439,437. 
    摘要 ( )   PDF(131KB) ( )  
    对海风热力内边界层模式稍作修改后,不仅可模拟海风流场,也可以模拟河谷中的风场。在模拟河谷地区风场的实例计算中,数值计算结果和实测值基本符合。将多水平方法中用得最多的PCG方法(预条件共轭梯度法)用来数值解由海风热力内边界层模式离散后得到的线代数方程组,节省了计算时间。
    超短脉冲简并四波混频时间分辨研究
    陈谋智, 余力, 黄文达
    1997, 14(S1): 440-442. 
    摘要 ( )   PDF(140KB) ( )  
    研究了简并四波混频(DFWM)的时间分辨曲线,数值模拟了由延迟时间引起的两种不同光栅的时间分辨曲线以及曲线的约化半宽度与约化弛豫时间tr/tp的变化关系。
    广义Abel方程的Gauss数值积分法
    母英魁, 王毅, 丁培柱
    1997, 14(S1): 443-444,442. 
    摘要 ( )   PDF(134KB) ( )  
    采用Gaus数值积分法将光学厚柱对称等离子体径向发射系数的广义Abel方程离散成系数是上三角矩阵的线性代数方程组,计算特别简单、省时间,具有2阶精度。
    二维不可压无粘流动问题的流线——混合元数值模拟
    陈焕祯, 张怀宇
    1997, 14(S1): 445-447. 
    摘要 ( )   PDF(165KB) ( )  
    基于流函数-涡度表达式,对不可压无粘流动问题采用沿流线的差分格式逼近涡度函数,采用混合元方法同时逼近流函数及速度场。理论分析与数值实验表明:格式对上述三个物理量具有同时逼近、高精度、高稳定性的良好性质。这是传统的数值方法难已达到的。
    外磁场作用下液晶中向列序的计算
    张志东, 邵喜斌, 黄锡珉
    1997, 14(S1): 448-449,447. 
    摘要 ( )   PDF(133KB) ( )  
    在格胞理论的基础上计算外磁场作用下液晶中的向列序。考虑分子短程关联后,序参数和内能对磁场强度的依赖关系与计算机模拟结果较好地符合。克服了平均场理论遇到的困难。
    解对称带状矩阵特征值问题的二分法及其改进
    罗晓广, 李晓梅, 陈健华
    1997, 14(S1): 450-452. 
    摘要 ( )   PDF(169KB) ( )  
    提出了解对称带状矩阵特征值问题的一种二分法。当仅需计算指定的部分特征值及其特征向量时,该方法尤其适合。进一步,我们还对二分法作改进。改进策略是:先用二分法计算若干步,得到特征值的近似值;然后从该近似值出发进行Rayleigh商迭代,直至其达到要求的精度为止。
    异步电动机动态温升的数值计算法
    李大农, 马志云
    1997, 14(S1): 453-454,452. 
    摘要 ( )   PDF(129KB) ( )  
    采用求热传导方程数值解的方法,计算封闭式异步电动机的动态温升.所得结果的准确程度与稳态温升计算法的准确程度相同。
    电容层析成象的数值模拟
    常鸿森, 保宗悌
    1997, 14(S1): 455-457. 
    摘要 ( )   PDF(151KB) ( )  
    描述了两相流成象电容层析系统的有限元模型和利用反投影法进行数值模拟的方法,对典型的层流状模型重建出流体成分分布图象。
    动力发动机内外流场数值模拟
    李杰, 鄂秦, 李凤蔚
    1997, 14(S1): 458-459,457. 
    摘要 ( )   PDF(178KB) ( )  
    采用分区网格技术,以Euler方程为主控方程,数值模拟了民用飞机的动力发动机内外流场。
    电热油藏的数学模拟研究
    关继腾, 王殿生, 黄柳宾
    1997, 14(S1): 460-462. 
    摘要 ( )   PDF(146KB) ( )  
    给出了电热油藏的r-z能量模型,导出了电热油藏时温度分布和产量计算的数学公式。由数值计算结果,讨论了油藏电导率对温度分布的影响,以及加热功率和加热方式对原油累积产量的影响。
    周期势场中定向扩散流的蒙特卡罗计算
    包景东
    1997, 14(S1): 463-464. 
    摘要 ( )   PDF(137KB) ( )  
    提出一种稳定和精确的坐标变换和闭合积分方法,用以求解描述周期势场中布朗粒子作定向运动的乘性朗之万方程。
    颗粒物质随机厚度的蒙特卡罗研究
    王仲奇
    1997, 14(S1): 465-466. 
    摘要 ( )   PDF(97KB) ( )  
    利用蒙特卡罗方法处理随机几何问题,是核工程中提出的实际问题,也是蒙特卡罗方法在几何处理方面遇到的新的理论问题。颗粒物质随机厚度分布的研究,是其中的重要部分。
    耦合J-J型方程中的DEVIL STAIRS
    周天寿
    1997, 14(S1): 467-468. 
    摘要 ( )   PDF(100KB) ( )  
    通过数值计算,首次观察到J-J型方程的一些台阶,而且在chaos状态下,有devilstairs的存在。
    代数多重网格法在排队理论中的应用
    常谦顺, 马书清, 雷光耀
    1997, 14(S1): 469-471. 
    摘要 ( )   PDF(131KB) ( )  
    提出一种校正AMG方法来求解排队模型。该方法对一般AMG方法的限制算子进行了改进,从而在求解过程中保留了问题的奇异性。计算结果表明改进后的方法加快了收敛速度,提高了解的精度。
    j壳层核子准旋、同位旋、角动量分类因子化重复度计算
    陈健华, 李承祖
    1997, 14(S1): 472-474. 
    摘要 ( )   PDF(129KB) ( )  
    据群链[SO(5)⊃U(1)DSU(2)T]ⓧ[Sp(2j+1)⊃O(3)],对j壳层核子按广义准旋[s,t]、同位旋T、角动量J和核子数n=2D+2j+1分类,分类重复度可因子化为SO(5)⊃U(1)DⓧSU(2)TSp(2j+1)⊃O(3)约化重复度之积。给出上述两个约化重复度的简明计算公式,并列出部分结果。
    一维凝聚生长的计算机模拟
    吴锋民, 朱启鹏, 吴自勤
    1997, 14(S1): 475-476. 
    摘要 ( )   PDF(115KB) ( )  
    在一维随机成核生长(RSNG)模型的基础上,考虑粒子的扩散和粒子数密度的影响,在最近邻及次近邻条件下模拟得到一些聚集生长的图形,并分析了其分形特征。
    电子碰撞激发的全相对论扭曲波理论和计算
    陈国新, 邱玉波
    1997, 14(S1): 477-479. 
    摘要 ( )   PDF(132KB) ( )  
    给出计算高荷离子的电子碰撞强度的全相对论扭曲波Born理论和有关高剥离离子电子碰撞激发截面的计算结果。并将计算结果和其他工作进行了比较。
    用多极理论计算复杂截面的电缆电导
    郑勤红, 解福瑶, 李明, 李景天, 欧家呜, 王瑞丽
    1997, 14(S1): 480-482. 
    摘要 ( )   PDF(141KB) ( )  
    将多极理论用于复杂截面的电缆电导计算,为电缆电导计算提供一种新的计算方法。实例计算结果表明,多极理论的计算精度比边界元法的计算精度高得多。
    小波分析在寄生参数提取中的应用
    古江春, 王泽毅, 洪先龙
    1997, 14(S1): 483-484,482. 
    摘要 ( )   PDF(132KB) ( )  
    将小波分析方法应用于求解三维(3-D)边界元素法形成的线性方程组。通过数据压缩将稠密矩阵化为稀疏矩阵,当求解具有相同系数矩阵和多个不同右端向量的方程组时,可以节省计算量。
    磁六极球差校正系统的计算与分析
    周庆, 陈尔纲, 胡宗高
    1997, 14(S1): 485-486. 
    摘要 ( )   PDF(95KB) ( )  
    限制电子显微镜分辨本领的主要因素是电子透镜的球差,对一种实际的校正磁透镜的磁六极球差校正系统进行了计算和分析。计算结果表明利用该系统可以消除三级球差。
    核爆多群γ-EMP中的Keller Box方法
    尹文华
    1997, 14(S1): 487-489. 
    摘要 ( )   PDF(140KB) ( )  
    主要叙述核爆弹体释放出的γ射线所形成的电磁脉冲以及解Maxwel方程组的KelerBox方法,并且给出了γ-EMP程序同美国SHARP程序的计算结果。
    数字化卫星云图自动识别的动态聚类算法
    杨越奎, 张方友, 刘玉玲
    1997, 14(S1): 490-492. 
    摘要 ( )   PDF(127KB) ( )  
    根据数字化气象卫星资料的特点,设计了对卫星云图进行自动分类识别的动态聚类算法。利用GMS-5的10.5~11.5μm和0.55~0.9μm两个通道的资料,采用迭代自组织的数据分析算法(ISODATA算法)对云图进行了分类试验。
    先进双发民用飞机跨音速绕流数值模拟
    鄂秦, 李杰, 李凤蔚, 刘积仑, 于建中
    1997, 14(S1): 493-495. 
    摘要 ( )   PDF(170KB) ( )  
    采用求解椭圆型和抛物型偏微分方程的方法及多块网格技术,生成机翼/机身/短舱/挂架复杂组合体外形的贴体及与边界正交的计算网格。应用有限体积法高效求解欧拉方程,模拟先进双发民用飞机的跨音速绕流,并计入机翼边界层的粘性效应,获得存在短舱、挂架干扰的翼身组合体气动特性。计算表明,短舱、挂架对超临界机翼的气动特性有显著的影响。计算结果与实验模型的机翼表面压力分布(测量值)符合良好
    火箭喷流波系的数值研究
    杨宏伟, 马大为, 张福祥
    1997, 14(S1): 496-498. 
    摘要 ( )   PDF(144KB) ( )  
    根据MUSCL格式的基本思想,对格式进行了部分修改,数值模拟了火箭喷流流场的波系结构。
    低能电子束曝光中的Monte-Carlo方法及沉积能分布
    谭震宇, 何延才
    1997, 14(S1): 499-501. 
    摘要 ( )   PDF(154KB) ( )  
    建立一个低能电子在多元多层介质中散射Monte-Carlo模拟计算方法。低能电子在多元介质中弹性散射用Mot截面描述并应用作者提出的平均散射截面方法进行模拟计算;给出多元介质中修正的Bethe方程计算低能电子非弹性散射能量损失。计算了低能电子在电子束曝光胶中能量沉积分布。
    可视化技术在量子力学中的应用
    王无敌, 陈健华, 李晓梅, 胡庆丰
    1997, 14(S1): 502-504. 
    摘要 ( )   PDF(122KB) ( )  
    科学计算可视化是近年来迅速发展的一个新兴领域。介绍了一个量子力学可视化系统VISQUM,讨论了一些量子力学的可视化手段和技术。
    线性复值随机跳跃——扩散方程隐Miltein格式的渐近稳定性
    刘小清, 李镇华, 吴声昌
    1997, 14(S1): 505-506. 
    摘要 ( )   PDF(110KB) ( )  
    给出了线性复值随机跳跃——扩散方程隐Miltein格式的渐近均方稳定的必要条件,并揭示了离散步长的稳定区域与格式隐性之间的关系。
    解薄圆板动力响应的DQ半解析方法
    彭建设, 张敬宇, 杨杰
    1997, 14(S1): 507-509. 
    摘要 ( )   PDF(131KB) ( )  
    直接从薄圆板的控制微分方程出发,在空间域采用DQ法(Diferential Quadrature Method)在时间域取级数,采用时域配点得到求位移场全部待定参数的线性方程组,解一次方程组即得全域的响应位移场。
    一类分解刚性大系统的组合方法
    陈丽容, 刘德贵
    1997, 14(S1): 510-511,509. 
    摘要 ( )   PDF(134KB) ( )  
    针对分解的刚性大系统提出了组合RK-Rosenbrock(CRKR)方法,该方法分别采用Rosenbrock和显式RK方法在不同的处理机上并行求解刚性和非刚性子系统。文中讨论了算法的构造、收敛性以及数值稳定性,并在微机和多处理机上进行了数值仿真试验。
    物料混合问题的一种改进辨识方法——加权最小二乘法
    张承慧, 张天德, 谭震宇
    1997, 14(S1): 512-514. 
    摘要 ( )   PDF(128KB) ( )  
    改进文[1]提出的物料混合参数辨识算法,使其适用于时变参数的实时辨识,提高对时变参数的动态跟踪速度。给出计算实例,并与文[1]方法的计算结果进行比较。
    一类全对称周期Jacobi矩阵逆特征值问题
    孔翠芳, 吕烔兴
    1997, 14(S1): 515-517. 
    摘要 ( )   PDF(121KB) ( )  
    提出了关于全对称周期Jacobi矩阵的一类逆特征值问题。导出了问题有唯一解的一个充分条件;在条件满足时给出了计算公式。
    一类广义非线性Schrödinger方程的孤立子解及其性质研究
    曹庆杰, 张天德, Price G W, Djidjeli K, Twizell E H
    1997, 14(S1): 518-520. 
    摘要 ( )   PDF(135KB) ( )  
    研究一类广义非线性Schrödinger方程的孤立子解及其性质,研究非线性参数变化(α→0及α→∞)时孤立子性态的变化规律,同时研究该系统的数值解法,得到了一类广义非线性Schrödinger方程差分格式的收敛性和稳定性的条件。
    用有效势方法研究一个CI相变模型
    许爱国, 王光瑞, 陈式刚
    1997, 14(S1): 521-523. 
    摘要 ( )   PDF(138KB) ( )  
    使用一个具有实际意义的双阱势作为原子间作用势、受扰余弦势势作为外势,引入一个新的FK类模型,用有效势方法研究其相图。发现了一系列新奇的、复杂的,但同时又有一定规律可循的相结构。
    改进准静态方法的研究
    许云林, 经荥清, 李君利
    1997, 14(S1): 524-526. 
    摘要 ( )   PDF(166KB) ( )  
    处理反应堆三维动力学问题时,在时间方面采用了改进准静态方法,并加入了综合法改进,它通过对中子通量的因式分解,达到放大时间步长的目的,减少了进行空间计算的次数;在必要的空间计算中,该方法又利用了格林函数节块法的原理求解形状函数的空间分布,使空间网格距可以达到有限差分法的近二十倍。
    高频放电CO2激光器电子能量分布及输运系数
    王新兵, 徐启阳, 刘儿兀, 谢明杰, 李再光
    1997, 14(S1): 527-529. 
    摘要 ( )   PDF(149KB) ( )  
    利用文献中给出的截面数据,引入高频放电激励CO2激光器的Boltzmann方程,对其进行数值求解,研究了频率对电子能量分布及输运系数的影响。
    简化流通量分裂算法研究
    马汉东, 吴礼义, 李素循
    1997, 14(S1): 528-529,531. 
    摘要 ( )   PDF(150KB) ( )  
    针对Steger-Warming流通矢量分裂格式和其修正格式所构成的混合格式进行了研究,引入压差放大因子,使其对于有粘激波流场分辨能力有所提高,同时采用一种特殊分裂方法以简化矩阵运算。利用该方法研究了三种超声速流场:(1)平面激波反射问题;(2)超声速边界层流动;(3)高超声速压缩拐角流动。
    出壳γ射线特征能谱作为弹头指纹的可行性计算分析
    伍钧, 刘成安, 任捷
    1997, 14(S1): 530-534. 
    摘要 ( )   PDF(142KB) ( )  
    通过对核弹头部件中核材料自发衰变γ射线输运过程的分析,根据不同类型的核弹头有不同的出壳γ射线能谱形状,研究泄漏γ射线特征能谱作为核弹头指纹的可行性。还通过对高能和低能γ射线穿透谱的比较,分析球壳的厚度。
    晶体C60取向相变的Monte-Carlo研究
    刘英, 聂向富
    1997, 14(S1): 535-536,534. 
    摘要 ( )   PDF(142KB) ( )  
    对晶体C60的取向相变、热滞及压力诱导相变现象进行了Monte-Carlo模拟,数值结果与实验基本一致。
    非线性有限元结构分析增量法的并行处理
    周树荃, 彭立东
    1997, 14(S1): 537-539. 
    摘要 ( )   PDF(129KB) ( )  
    给出了非线性有限元结构分析著名的增量法的并行处理算法,并在多处理机系统Chalenge-L上进行了数值试验。
    反应堆探测器孔道的计算分析
    钟文发, 周全
    1997, 14(S1): 540-541,539. 
    摘要 ( )   PDF(131KB) ( )  
    针对反应堆物理设计中探测器孔道的设置,提出了计算模型和计算方法,并作了实际计算。结果表明,探测器孔道的设置使安全壳外和混凝土壳外的中子通量和γ通量都会增大,对屏蔽是有影响的。在屏蔽设计中应考虑它的影响,以确保反应堆的安全。
    三元合金异质结(AlP)x(Si2)1-x/GaP和(GaP)x(Si2)1-x/GaP的价带带阶
    蔡淑惠, 郑金成, 王仁智, 郑永梅
    1997, 14(S1): 542-544. 
    摘要 ( )   PDF(137KB) ( )  
    采用基于LMTO-ASA的平均键能计算方法和原子集团展开方法,研究了两组晶格匹配三元合金异质结(AlP)x(Si2)1-x/GaP和(GaP)x(Si2)1-x/GaP的价带带阶ΔEv(x)值。研究表明,两组异质结的ΔEv(x)值随合金组分x的变化都是非线性的,且表现出非单调的关系。
    超新星铁星核区对流不稳定的数值模拟
    张锁春
    1997, 14(S1): 545-546. 
    摘要 ( )   PDF(100KB) ( )  
    通过数值模拟发现,超新星铁星核区中~(0.74~0.91)M区间是十分重要的对流不稳定区域。
    Ⅱ型超新星铁星核质量对瞬发爆炸的影响
    王贻仁, 谢佐恒, 汪惟中, 张锁春
    1997, 14(S1): 547-549. 
    摘要 ( )   PDF(134KB) ( )  
    提出了将反弹冲击波在Ⅱ型超新星(SNⅡ)外星核区传输中的光裂解损耗看成为一个可调参量,用构造SNⅡ"瞬发爆炸模型"的方法研究前身星铁星核质量对SNⅡ瞬发爆炸的影响。并且给出了可发生SNⅡ瞬发爆炸的铁星核质量的范围。对WW(1993)SNⅡ前身星模型系列研究结果表明:由于铁星核质量过大,它们都不大可能发生瞬发爆炸。
    耦合海气动力学的初步研究
    姜达雍, 季良达, 杨大升
    1997, 14(S1): 550-552. 
    摘要 ( )   PDF(146KB) ( )  
    依据耦合海气系统观测研究事实和耦合海气动力学原理,研制了多个海气耦合模式(CGCMs);并利用实时观测资料,对季和年际时间尺度由海气相互作用所引起的气候变率,进行了数值模拟研究。
    用ALE方法数值模拟气缸内三维气体流场的研究
    梁树青, 张兆群, 魏德福, 田东风
    1997, 14(S1): 553-555. 
    摘要 ( )   PDF(145KB) ( )  
    用数值模拟的方法研究气缸内气体流动的全貌以探求合理组织气流的途径,这是当前国内外研究的重要课题,在开发KIVA Ⅱ程序的基础上研究适合该问题的三维ALE方法,并利用此方法数值模拟了一模型问题。
    溃坝洪水长波与地面障碍物作用的数值模拟
    王如云, 方国洪
    1997, 14(S1): 556-557,555. 
    摘要 ( )   PDF(150KB) ( )  
    利用TVD差分格式设计思想,给出了一种数值求解浅水波方程的数值方法。经过数值解与精确解的比较,表明该方法能很好地模拟溃坝洪水波间断面的位置和形状。进一步对溃坝洪水长波与地面障碍物的作用进行数值模拟,得到了合理的结果。
    具有边梁加固的板弯曲问题的差分方法
    曲小钢
    1997, 14(S1): 558-560. 
    摘要 ( )   PDF(131KB) ( )  
    由极小势能原理出发,构造了具有边梁加固的板弯曲问题的差分格式,给出了数值结果,定量地研究了板与边梁相互作用的规律。
    规整星型高分子回转半径的几率分布
    王治虎, 罗孟波, 许健民
    1997, 14(S1): 561-561,560. 
    摘要 ( )   PDF(134KB) ( )  
    采用简立方格点上的普通无规行走(NRW)模型,用Monte-Carlo模拟,研究了规整星型高分子回转半径的几率分布,得到它与总链段数、官能度和根均方回转半径之间的关系。
    有限差分法在CrNi钢淬火中的应用
    程赫明, 郭荣鑫, 王洪纲
    1997, 14(S1): 562-565. 
    摘要 ( )   PDF(132KB) ( )  
    将显式有限差分法用于CrNi钢淬火过程,求解非线性热传导方程的逆问题,从而可以确定钢淬火过程中的表面换热系数与温度的关系。
    正态分布的抽样方法
    宫野
    1997, 14(S1): 566-568. 
    摘要 ( )   PDF(133KB) ( )  
    给出了若干实现正态分布的随机抽样方法,并对各种抽样方法作了比较,可供在实际问题中选用。
    有限元方法在二维六角形组件堆芯临界问题计算中的应用
    陈昌友
    1997, 14(S1): 569-571. 
    摘要 ( )   PDF(129KB) ( )  
    应用有限元方法对二维六角形组件堆芯的少群中子扩散临界计算进行了研究,并编制了二维少群中子扩散方程有限元计算程序2DFME。以快堆300MW(e)-LMFBR临界问题[2]为例,作了数值计算。且与国内外有限元程序[1~3]和有限差分程序CITATION的计算结果进行了分析比较。
    三维多级扩散室对冲击波传播规律的影响
    郝保田, 张海波
    1997, 14(S1): 572-573,571. 
    摘要 ( )   PDF(160KB) ( )  
    运用三维FLIC数值方法研究了爆炸冲击波在三维不对称多级扩散坑道中的运动过程和传播规律。
    准一维无序系统的电子结构
    徐慧
    1997, 14(S1): 574-576. 
    摘要 ( )   PDF(132KB) ( )  
    利用负本征值理论的态密度计算方法,研究了准一维双链无序系统的电子结构。并针对系统大小,对角和非对角无序程度等各种参数,探讨了电子的局域化、系统能量分布范围等问题。在引进无序参量后,对角无序主要引起电子局域态的增多,而非对角无序则使系统的能量分布范围发生变化。
    低泄漏堆芯多循环优化方法及其程序系统
    程平东, 沈炜, 钱惠富
    1997, 14(S1): 577-578. 
    摘要 ( )   PDF(98KB) ( )  
    为低泄漏堆芯燃料管理发展了一种多循环优化方法和程序系统。多循环优化过程分三步完成:第一步,用线性规划确定能使多循环总寿期最长的相继各循环功率份额最佳分布和循环长度最佳分配;第二步,以前一步结果为条件,用线性规划或直接搜索法对多循环中相继的各个单循环进行燃料组件的优化布置;第三步,用可变容差法(FTM)实现可燃毒物最佳配置。
    计入分离涡时飞机机翼振动的极限环现象
    叶正寅, 杨永年
    1997, 14(S1): 579-580,578. 
    摘要 ( )   PDF(152KB) ( )  
    通过同时求解机翼弹性运动方程和空气动力方程,对机翼的流固耦合问题进行了数值仿真。在空气动力计算中仔细地模拟了机翼在较大迎角时从机翼前缘和侧缘处卷起的分离涡。通过不同机翼、不同状态的仿真结果,发现如下结论:(1)机翼前缘或侧缘分离涡将使机翼颤振临界速度降低,且机翼的基本迎角越大,颤振临界速度越低;(2)机翼前缘后掠角减小,分离涡对颤振性质的不利影响有所减缓;(3)在一定条件下,会出现稳定的极限环。
    VO332+53爆发时气体流的结构
    宓国柱, 沈金英
    1997, 14(S1): 581-583. 
    摘要 ( )   PDF(130KB) ( )  
    运用气盘和缓慢加速星风的混合模型,计算了由Be星和中子星组成的双星系统的X射线光度的变化,与VO332+53的观测资料作了较好的拟合。
    Euler方程的二速度三熵级格子Boltzmann模拟
    阎广武, 胡守信, 施卫平
    1997, 14(S1): 584-586. 
    摘要 ( )   PDF(131KB) ( )  
    提出了一个新的用于模拟可压缩Euler方程的格子Boltzmann模型。数值试验表明它可以模拟具有激波和接触面的问题。计算结果与经典算法的结果可以比较。
    线性和拟线性抛物型方程的格点模型
    李元香, 康立山
    1997, 14(S1): 587-589. 
    摘要 ( )   PDF(134KB) ( )  
    格点法是最近发展起来的一类模拟流体动力学的新方法,其基本思想是重新建模,推广格点法求解一般的数学物理方程,作为例子,建立线性和拟线性抛物型方程的格点模型,并运用它们模拟了几个典型的问题。
    振荡机翼跨音速非定常粘性绕流的数值计算
    郭广利, 杨永年
    1997, 14(S1): 590-591,589. 
    摘要 ( )   PDF(139KB) ( )  
    以非定常N-S方程为主管方程,采用LU-NND混合差分格式,C和C-H型贴体运动网格,B-L双层代数紊流模型,求解绕振荡翼型和三维机翼的跨音速非定常粘性流场。分别计算了NACA 0012翼型和M6机翼作俯仰振荡时跨音速非定常粘性绕流流场。研究了非定常绕流的气动特性。部分计算结果和风洞实验值作了比较。
    从混沌信号中识别振荡分量的一种算法
    刘健文, 董佩明, 李耀东, 王炳仁, 金维明
    1997, 14(S1): 592-594. 
    摘要 ( )   PDF(143KB) ( )  
    阐述了奇异谱分析方法的数学原理和分离振荡分量的重要性质。取100天的窗口长度研究我国东南地区地面气温序列,发现明显存在周期为35~40天和20天左右的主要低频振荡,相应频率段约占总方差的90%。
    网格划分对四周切向煤粉炉流场计算的影响
    汪箭, 刘团结, 范维澄
    1997, 14(S1): 595-596,594. 
    摘要 ( )   PDF(148KB) ( )  
    通过选用不同的网格划分方法,对四周切向煤粉炉流场进行数值计算。若选用常规的网格划分方法计算,虚假扩散的影响非常严重,势必造成计算结果较大的偏差;同时,给出适当的网格划分,以及相应合理的计算结果。
    解地下水流的Monte-Carlo随机法与积分插值法
    崔秀, 吴飞
    1997, 14(S1): 597-599. 
    摘要 ( )   PDF(136KB) ( )  
    分别采用积分插值法与Monte-Carlo随机法对地下水流问题进行了理论研究,提出了Monte-Carlo随机差分法。对这两种方法的计算格式及边值条件在格式中的处理进行了推导、分析及论述;并且选取了一个具有解析解的二维承压稳定地下水流数学模型作为例子,进行了数值计算。
    修正Burridge-Knopoff模型的地震前兆场数值模拟
    肖家鑫, 林家逖, 郑熙铭
    1997, 14(S1): 600-602. 
    摘要 ( )   PDF(142KB) ( )  
    对Buridge-Knopof模型作了非线性化与非均匀化两项重要修正,并利用四阶Runge-Kuta法对地震前兆场的演化过程进行数值模拟,得到与实例相似的结果。
    Levinthal矛盾的一种解释
    魏广川, 林家逖, 肖家鑫
    1997, 14(S1): 603-604,602. 
    摘要 ( )   PDF(135KB) ( )  
    研究蛋白质残基间相互作用对折叠过程中Levinthal矛盾的影响。结果表明此矛盾由两个原因所致:一是没有考虑转变几率与能量差的关系,二是没有计入残基间相互作用能的影响。
    二维模具成形过程计算机仿真
    李坚石, 陈天健, 梅光前
    1997, 14(S1): 605-607. 
    摘要 ( )   PDF(132KB) ( )  
    论述了二维模具成形过程力学模型及其接触搜寻算法和动态显式算法,给出一个电器模具应用实例。
    一个含F旋的IBM-2软件包
    张进宇, 龙桂鲁, 孙洪洲, 韩其智
    1997, 14(S1): 608-610. 
    摘要 ( )   PDF(133KB) ( )  
    介绍一个新的IBM-2程序。它以F旋对称的U(5)极限的波函数做为基,并将Hamiltonian写成F旋张量的形式。除了通常的显含质子中子指标的Hamitonian的输入外,还有F旋张量的形式输入。并且能够计算B(E2)和B(M1)。
    构造定常对流扩散方程高精度紧致差分格式的新方法
    田振夫
    1997, 14(S1): 611-613. 
    摘要 ( )   PDF(129KB) ( )  
    以一维定常对流扩散方程的高精度差分格式为基础,提出了一种构造二维定常对流扩散方程高精度紧致差分格式的新方法,并给出数值例子。
    激光等离子体特征参量和辐射耦合效应研究
    张建泉, 刘峰, 陈荣华, 强希文
    1997, 14(S1): 614-616. 
    摘要 ( )   PDF(141KB) ( )  
    采用一维单流体双温度辐射流体力学模型,在激光产生的等离子体做等温膨胀的假定下,研究了脉冲强激光辐照固态靶材时,计及再辐射后汽化等离子体特征参量与靶材物质和激光参数之间的定标关系,计算了等离子体再辐射对激光-靶材耦合效应的影响。
    六角形燃料组件内精细核密度分布计算
    陈超, 朱广宁, 王侃
    1997, 14(S1): 617-618,616. 
    摘要 ( )   PDF(138KB) ( )  
    介绍一种可用于六角形燃料组件的精细计算方法,在理论模型上编制了精细计算程序JXJS;JXJS计算结果与用文[1]提供的方程计算出的结果进行了比较(61Δ/六角形组件);并对结果进行了误差分析。
    非对称广义特征值问题的并行连续同伦算法
    薛长峰, 周树荃
    1997, 14(S1): 619-621. 
    摘要 ( )   PDF(131KB) ( )  
    非对称矩阵广义特征值问题的并行计算是大规模工程计算中的基础问题之一。迄今为止,国内外关于该问题的研究进行得很少,提出的有效算法不多。首次将连续同伦算法应用于非对称广义特征值问题的并行求解,提出并行连续同伦算法。
    用改进的射线踪迹法分析含漫射表面的半透明体表观光谱发射特性
    夏新林, 余其铮, 谈和平
    1997, 14(S1): 622-623,621. 
    摘要 ( )   PDF(147KB) ( )  
    通过反向应用射线踪迹法并与漫反射表面的有效辐射强度相结合来分析含漫射表面半透明体的表观光谱发射特性。漫反射表面的有效辐射强度通过迭代离散区域的能量平衡方程并结合反向应用射线踪迹法进行求解。将该方法与半透明大平板表观发射率的解析解及圆盘状半透明体表观发射率的蒙特卡洛方法数值解进行了比较,并用该方法对圆盘状半透明体的表观发射进行了分析。
    纳米金属颗粒复合物中偶极相互作用的数字模拟
    姚端正, 熊贵光
    1997, 14(S1): 624-626. 
    摘要 ( )   PDF(136KB) ( )  
    应用耦合偶极方程(CDE),描述了纳米尺寸金属颗粒复合物的偶极相互作用,计算了在准静态近似下纳米尺寸的金、银及铜的复合物中,偶极相互作用的频率变量和衰变参量的波长相关特征,并对结果进行了分析。
    半圆型光纤液面计结构的数值理论研究
    李英骏, 滕爱萍, 张广才, 彭翰生
    1997, 14(S1): 627-628,626. 
    摘要 ( )   PDF(149KB) ( )  
    采用数值模拟方法,对光纤纤芯半径分别是3μm,5μm,10μm,15μm,光纤半径都是20μm的光纤进行计算机模拟,并得到了相对反射光强在光纤入射平面上的分布曲线。其结果表明,对半园型液面计中反射式半圆形探头最佳耦合结构的半径满足ρmax=R-0.01004a3+0.37429a2-3.38095a+16.85714(μm)
    用快速傅里叶变换(FFT)计算任意平面分层介质中电磁散射
    胡俊, 聂在平
    1997, 14(S1): 629-631. 
    摘要 ( )   PDF(138KB) ( )  
    在非均匀介质电磁散射、逆散射等领域中,经常需要快速、精确地计算汉克尔变换来完成高效率的正演计算。针对任意平面分层的非均匀介质情形,以偶极子激励为例,应用快速傅里叶变换(FFT)来完成汉克尔变换计算。
    五对角阵逆特征值问题的一个数值解法及其应用
    孙合明, 吕烔兴, 赵淳生
    1997, 14(S1): 632-634,631. 
    摘要 ( )   PDF(205KB) ( )  
    给出一个由谱数据构造实对称五对角阵的数值方法及其在梁的振动逆问题中的应用。
    弹性波散射问题的边界元解法
    陈国荣, 姜弘道
    1997, 14(S1): 635-637. 
    摘要 ( )   PDF(125KB) ( )  
    采用边界单元法求解弹性波对孔口的散射问题。给出了弹性波散射问题的边界积分方程,针对数值计算中出现的奇异积分,提出了一种改进的把动力基本解分解为正则部分和奇异部分分别计算的方法。最后讨论了P波和SV波对圆孔和椭圆孔的散射而引起的动应力集中。
    直流电弧炉偏弧和控弧的计算机仿真
    洪新, 朱守军, 谢树元, 刘俊江, 傅润生, 罗建江
    1997, 14(S1): 638-640. 
    摘要 ( )   PDF(150KB) ( )  
    根据电弧弧柱在微元电流相互作用下的受力以及空间载流导线在电弧区所造成电磁力动态平衡的原理,建立了研究直流电弧电磁偏弧,多底电极分电流控弧和顶电极倾动控弧的数学模型,并通过计算机仿真实验对载流导线的电磁影响和一些纠正和控制偏弧的方法进行了考察和研究。
    一种快速收敛发散三角级数的有效算法
    俞慧丹
    1997, 14(S1): 641-642,640. 
    摘要 ( )   PDF(134KB) ( )  
    针对Bensimon方程中出现的无穷发散三角级数求和问题,提出一种快速收敛发散三角级数算法,并举例说明其应用。
    变压器励磁涌流场路耦合数字仿真
    肖鸿杰
    1997, 14(S1): 643-645. 
    摘要 ( )   PDF(161KB) ( )  
    变压器励磁涌流的准确计算是一个十分困难的问题。提出用非线性时步有限元法来模拟变压器励磁涌流的动态过程。采用的数学模型将非线性时变电磁场有限元方程与电压源下的电路微分方程耦合在一起,并考虑铁心的磁饱和与磁滞特性。仿真结果与实验曲线基本吻合。
    电子束辐照冲量的数值计算与实验的对比
    丁升, 周南
    1997, 14(S1): 646-648. 
    摘要 ( )   PDF(132KB) ( )  
    利用Monte-Carlo方法模拟了电子在材料中的输运过程,得到了材料中电子的能量沉积剖面,同时结合弹塑性流体力学方程组,对"闪光二号"脉冲电子束发生装置产生的电子束与Ly-12铝相互作用的热-力学效应逐炮进行了数值计算,着重研究了喷射冲量及质量亏损,并与实验测试结果进行了对比。
    辛算法在近地小行星轨道演化数值研究中的应用
    刘林, 廖新浩, 季江徽
    1997, 14(S1): 649-651. 
    摘要 ( )   PDF(133KB) ( )  
    采用改进的显式辛算法对近地小行星的轨道演化进行数值研究,在力学模型中除考虑各大行星的引力摄动外,还考虑了后牛顿效应,而在算法上则着重探索辛算法在近地小行星轨道演化研究中的应用前景,特别是当这类小行星与某一大行星靠近时辛算法的有效性。
    计算恒星结构演化的Henyey方法的代数形式
    罗国权, 黄润乾
    1997, 14(S1): 652-654. 
    摘要 ( )   PDF(131KB) ( )  
    基于恒星结构演化模型的方程组和两点边界条件,用矢量函数表述它们的差分方程组,并用Newton-Raphson方法展开,获得了可求解的代数方程和矩阵方程。
    一类带组合边界的变物性耦合场问题
    阎贵卿
    1997, 14(S1): 655-656,654. 
    摘要 ( )   PDF(136KB) ( )  
    在正交各向同性体框架内,讨论在热冲击作用下,具耦合效应的变物性热弹性方程的求解问题。它是近几年来关于热冲击下常物性耦合场讨论的继续。
    消除模式“负水汽”算法的研究与试验
    葛孝贞, 李峰
    1997, 14(S1): 657-659. 
    摘要 ( )   PDF(155KB) ( )  
    以中尺度模式MM4为基础,将二种保证正定并具高精度的平流格式引入该模式计算水汽平流,模式在分别采用不同水汽平流格式下(其它模式条件完全相同),通过实例风场下对不同格式进行了理论测试和多个暴雨过程的实例模拟。
    云模式的航空应用探讨
    王炳仁, 金维明, 李耀东, 王洪芳, 刘健文, 董佩明
    1997, 14(S1): 660-662. 
    摘要 ( )   PDF(239KB) ( )  
    以中尺度数值模式为动力框架与显式方案的云模式耦合,对云场的宏观分布和云中水含量、云粒子谱分布进行模拟,进而探讨其在航空方面的应用。
    用改进的格子Boltzmann方法模拟浅水长波问题
    施卫平, 胡守信, 阎广武
    1997, 14(S1): 663-665. 
    摘要 ( )   PDF(131KB) ( )  
    以正六边形的格子Boltzmann方法(LBM)为基础,通过气体动力学的Euler方程与浅水长波方程的比拟,修改了原有的LBM方法的平衡态分布函数,使LBM方法所对应的流体力学方程与浅水长波方程一致。一维问题的计算结果与离散速度法得到的计算结果作了比较。
    恢复函数的三点迎风紧致格式构造方法及应用
    王保国, 郭延虎, 沈孟育
    1997, 14(S1): 666-668. 
    摘要 ( )   PDF(137KB) ( )  
    在结构网格下,从三点迎风紧致逼近出发提出了一种适合于有限体积离散的恢复函数生成办法,在光滑区它具有三阶精度,并且在捕捉激波时有较高的激波分辨率。文中完成了典型算例。
    铁铝酸盐水泥矿物结构、性能与量子化学计算研究
    闵新民, 邓志平, 蔡克峰, 李嘉宇
    1997, 14(S1): 669-670,668. 
    摘要 ( )   PDF(139KB) ( )  
    应用自洽场离散变分Xα(SCF-Xα-DV)量子化学分子轨道计算方法研究铁铝酸盐水泥矿物,讨论了结构、性能与化学键之间的关系。计算表明,Al净电荷高于Fe,Al-O共价键弱于Fe-O,同时,随Al含量的增加,Al净电荷增高,Al-O共价键减弱,最低空轨道能量降低,是铁铝酸盐水化活性随铝含量的增加而提高的主要原因。
    格子Boltzmann方法的质量扩散模型
    孙成海, 王保国, 沈孟育
    1997, 14(S1): 671-673. 
    摘要 ( )   PDF(154KB) ( )  
    建立了多组份流体、多速度格子Boltzmann模型。该模型满足流体质量扩散方程和运动方程。通过正弦波的衰减过程测量了两流体间的扩散系数。对圆形区域内的静止流体在具有均匀来流速度的另外一种流体中的对流扩散问题进行了数值模拟。
    全方位离子注入柱形几何碰撞鞘层时空演化数值解
    于炯, 宫野
    1997, 14(S1): 674-676. 
    摘要 ( )   PDF(129KB) ( )  
    对使用流体动力学碰撞模型研究全方位离子注入柱形靶等离子体鞘层时空演化的求解过程进行了详细的讨论,得到了数值计算结果。
    发汗控制方程注记
    杨学实
    1997, 14(S1): 677-679. 
    摘要 ( )   PDF(131KB) ( )  
    指出了三种类型的发汗控制方程描述的物理对象和某些参数的性质。给出了部分数值仿真结果。为可控性证明提供了直观概念。
    用蒙特卡洛方法计算红外光学系统的杂散光
    夏新林, 谈和平, 余其铮, 王平阳
    1997, 14(S1): 680-681,679. 
    摘要 ( )   PDF(138KB) ( )  
    将蒙特卡洛方法应用于红外光学系统的杂散光计算。对光学系统的几何结构和部件物理类型进行了标准化处理,建立了有关概率模型并采用能有效提高计算速度的光线跟踪方法。在此基础上编制了通用的杂散光计算软件并对某型号红外相机的杂散光进行了分析。
    脉冲强激光冲击力学效应的数值模拟
    陈荣华, 张建泉, 李邦固, 刘峰, 强希文
    1997, 14(S1): 682-683. 
    摘要 ( )   PDF(106KB) ( )  
    利用一维理想弹塑性流体动力学模型和有限差分方法,研究了激光辐照固态靶材时产生的激波的形成和衰减过程以及动态断裂。
    数值方法形式矢量场的递推求解
    廖新浩, 刘林
    1997, 14(S1): 683-683. 
    摘要 ( )   PDF(65KB) ( )  
    给出了动力系统数值方法对应的形式矢量场的递推求解方法,这为数值方法的深入研究提供了一种手段。利用这一求解方法对具有下列Hamilton函数H=H(p,z)的动力系统进行了讨论,针对辛格式或非辛格式,给出了相应的形式矢量场。
    高比度同位素192Ir和60Co辐照的计算
    周全, 钟文发, 徐小琳
    1997, 14(S1): 684-686. 
    摘要 ( )   PDF(127KB) ( )  
    利用反应堆物理的方法,对使用高通量堆辐照生产高比度同位素192Ir和60Co的方案进行了计算。从同位素产量和对堆物理参数的影响两方面分析了计算结果,提出了一个方案作为实际生产的参考。
    光电子全息的计算机模拟研究
    韩继红, 陆薇, 顾昌鑫
    1997, 14(S1): 686-686. 
    摘要 ( )   PDF(58KB) ( )  
    光电子全息作为来自材料本体的电子全息,近年成为一个活跃的研究领域。X射线被原子芯能级吸收后,激发产生的光电子波在传播过程中被周围原子散射,带有周围原子结构信息的散射波(物波)与未被散射的光电子波(参考波)同源产生,有极好的相干性,这两列波的电子衍射图即是光电子全息图。通过适当变换,可实现全息图的重现,重构出三维空间原子像,获得原子级分辨率的三维结构信息。光电子全息作为直接揭示原子结构的有效手段而得到重视和发展。采用合理的物理模型,对金属Cu原子、C(2×2)S/Ni(100)表面中S原子以及YBCO原胞中Cu2原子近邻的一些原子簇,进行了光电子全息产生过程的计算机模拟,同时作了多重散射和球面小修正。全息重现则采用数值重现的方法,获得了准确的三维晶格结构参量。
    DNA分子非线性局域激发的数值计算
    杨显俊, 王仲奇, 卓益忠
    1997, 14(S1): 687-689. 
    摘要 ( )   PDF(185KB) ( )  
    用一简化模型对DNA非线性局域激发的机制作了初步的研究,用五阶Runge-Kuta数值计算来模拟DNA分子中N个非线性耦合碱基对的大振幅运动。通过适当选取H量中各参数及步长,调节色散效应与非色散非线性效应的比份以利于动力学稳定,使结果能初步给出非线性局域激发的形成与稳定传播的物理图象。并能看出非线性效应对DNA分子的"呼吸"效应的重要影响。
    平方守恒半拉格朗日法在水汽方程求解中的应用
    季仲贞, 王斌
    1997, 14(S1): 689-689. 
    摘要 ( )   PDF(60KB) ( )  
    平方守恒半拉格朗日法在水汽方程求解中的应用季仲贞王斌(LASG,中国科学院大气物理研究所,北京100029)摘要平方守恒半拉格朗日法是在迎风差格式和浮动算法的思想基础上发展起来的一种经济算法。它不仅具有类似半拉格朗日法的特点,而且还能保持原系统的重要物理特征-能量守恒性,可将它用于具体的短期气候模拟计算。对于IAPL0IT42气候谱模式,我们采用平方守恒半拉格朗日法来求解水汽方程。由于平方守恒半拉格朗日法具有平方守恒性和正定性,可以克服计算不稳定现象和出现负水汽现象,因此,即使在时间步长增大的情况下,仍可以较客观的描述水汽的输送,因而可以改善降水的模拟。用具体的实例计算,取得良好的效果。
    谐振势阱中有弱相互作用的玻色爱因斯坦凝聚
    郑金成, 严子浚, 林仲金, 陈丽璇
    1997, 14(S1): 690-692. 
    摘要 ( )   PDF(134KB) ( )  
    应用数值计算的方法计算了谐振势阱中有弱相互作用的玻色气体凝聚的临界温度和基态占据率。计算结果表明,二者都随着散射长度的增大而减小;但与理想玻色气体相比仅差约0.4%。
    延时反馈控制Van der Pol振子的混沌
    吴锋民, 施建青, 李巧文
    1997, 14(S1): 692-962. 
    摘要 ( )   PDF(61KB) ( )  
    VanderPol振子是力学、物理学和电子学中著名的数学模型,已有不少工作表明VanderPol振子在一定的条件下会出现混沌行为。在实际的工程问题中,混沌振动常会导致动力系统的不规则运作、机械系统的疲劳断裂以及电路系统的反常噪声,故这种混沌行为往往应设法避免。在此用德国学者K.Pyrages提出的延时连续反馈控制方法,通过计算机模拟,对VanderPol振子的混沌行为进行控制,结果表明,只要选择合适的延迟时间τ及控制信号幅值K,可以有效地抑制系统的混沌运动而得到稳定周期(或准周期)振动结果。
    非均匀圆柱光散射的彩虹现象
    吴振森
    1997, 14(S1): 695-695. 
    摘要 ( )   PDF(61KB) ( )  
    与球形粒子类似,当光照射具有大粒子尺寸参数圆柱时,也会出现明显的彩虹现象。根据第一,二级彩虹角位置和干涉条纹数,可以测定圆柱体的折射率和尺寸参数,应用于粒度分析中。利用平面波谱展开和矢量柱谐函数研究非均匀圆柱对2D高斯波束的散射,导出分层圆柱电磁散射系数的稳定,收敛的迭代公式和算法。采用32兆内存、586微机,其计算范围之大,可以使分层数与粒子尺寸参数的乘积达5×106。还计算了不同水柱,酒精柱和光纤的光散射以及彩虹角分布。同时根据射线理论导出彩虹角反演圆柱折射率的计算公式。上述正、逆两种问题与实验比较均有很好的吻合。
    最小二乘原理求解多维瞬态导热反问题
    白博峰, 郭烈锦, 陈学俊
    1997, 14(S1): 696-698. 
    摘要 ( )   PDF(143KB) ( )  
    采用最小二乘法原理,研究了确定未知边界条件的多维瞬态导热反问题求解方法.
    一类半导体问题的特征差分方法的L([0,T];L(Ω))模估计
    刘蕴贤
    1997, 14(S1): 698-698. 
    摘要 ( )   PDF(61KB) ( )  
    采用特征差分方法研究了一类半导体问题的电子和空穴浓度方程(e)/(t)-()/(x)(De(e)/(x))-μeq(e)/(x)-(∂μe)(x)qe-αμee(p-e+G)=-R(e,p)#br#(p)/(t)-()/(x)(Dp(p)/(x))-μpq(p)/(x)-(∂μp)(x)qe-αμee(p-e+G)=-R(e,p)其中未知函数是电子和空穴浓度e,p。用特征差分方法对一维有界区域上的非齐次牛曼边界条件(e)/(x)|Ω=g(x,t),(p)/(x)Ω=r(x,t)及初始条件e(x,0)=e0(x),p(x,0)=p0(x)。分别给出区域内部和边界的计算方法,并且用最大模原理给出建立在线性差值上的特征差分方法的L-模估计,误差界为Kt+h)。
    一类常微分方程边值问题二点格式的算法研究
    谭志海, 吴雄华
    1997, 14(S1): 699-701. 
    摘要 ( )   PDF(123KB) ( )  
    根据文[1]的思想对一类奇异摄动问题给出了高精度的任意不等距二点差分格式。在边界条件处理上对文[1]的算法作了改进,使算法不等距,在奇异点附近网格可任意加密并重点研究了离散方程组的计算方法。数值试验表明,本算法在处理奇异摄动问题时更灵活,更有效。
    可分线性哈密顿系统显式辛格式的守恒量
    吴承埙, 周忠源, 丁培柱
    1997, 14(S1): 702-704. 
    摘要 ( )   PDF(145KB) ( )  
    证明了对可分线性哈密顿系统的每个二次形守恒量,1阶和2阶显式辛格式存在对应的格式守恒量;由此可说明,步长适当小时,1阶和2阶显式辛格式保持量子系统的波函数模方守恒,是直接求解时间相关Schrdinger方程(TDSE),以研究量子系统时间演化的合理和有效的数值方法。
    以锕系元素为燃料的先进型次临界增殖反应堆研究
    武玉
    1997, 14(S1): 705-707. 
    摘要 ( )   PDF(146KB) ( )  
    以锕系元素为燃料的次临界增殖反应堆,在D-T聚变中子源驱动下,可以更有效地利用有限的铀资源;反应堆keff < 1,且功率密度较低,有较好的安全性;反应堆对D-T中子源要求低。50MWt~100MWt的聚变功率即可满足1GWe反应堆要求,包层可生产足够的氚供堆芯聚变使用,此外每年还可以增殖1000kg左右可裂变元素供裂变堆使用,支持同等功率3~4个裂变堆;反应堆每年燃烧3~4个裂变堆产生的锕系元素。
    非线性多点边值问题的插值矩阵法
    牛忠荣, 王秀喜
    1997, 14(S1): 708-710. 
    摘要 ( )   PDF(136KB) ( )  
    建立了插值矩阵法的基本理论,用于解非线性混合阶常微分方程组多点边值问题。制作了常微分方程组求解器IVMMS,可以支持计算力学中的有限元线法。
    平面叶栅跨音绕流并行数值模拟
    谢叻, 戴嘉尊
    1997, 14(S1): 711-713. 
    摘要 ( )   PDF(132KB) ( )  
    研究平面叶栅跨音绕流的并行算法。采用文[1]中差分格式,设计适合平面叶栅跨音绕流的并行算法,并在江南一2并行机上实现。
    国内外“计算物理”发展动态
    张锁春
    1997, 14(S1): 714-721. 
    摘要 ( )   PDF(268KB) ( )  
    我想趁计算物理学会第四届年会召开的机会,向从事计算物理工作的同仁们介绍国内外计算物理发展的一些情况,尽管知道得不多不全,但我想这样做有助于交流,对促进计算物理在中国的发展有好处。本文分为国际和国内两大部分,以正式出版物为依据(绝大部分资料在北京图书馆或中国科学院图书馆能找到),分为杂志类、图书类、丛书类和会议论文集类来叙述,从中看出计算物理的发展动态。本人稍加评注,更多的留给读者自己去评说。
    控制系统中的泛函微分方程转换方法
    计国君, 宋文忠
    1997, 14(S1): 722-722. 
    摘要 ( )   PDF(50KB) ( )  
    含有时延的控制系统代表的模型为泛函微分方程(或差分方程),最贴切实际的控制系统为x(t)=Ax(t)+Bu(t-τ)或x(t)=f(t,x(t),u(t-τ)),如何利用无时延的控制系统的现有结论去研究时延控制系统,给出了几种转换方法。将无时延系统中建立的方法转换到时延控制系统中为控制系统基本转换方法(如等效控制法,关联预估法等);利用泛函微分方程数值积分法为泛函微分方程基本转换法(如渐近方法,矩量法等);将时延项视为附加扰动输入法为无时延转换法;将时延项以灵敏度系数调节为灵敏度系数法;考虑具齐次初始问题的系统为降解维数法;利用Lyapunov泛函设计控制器或变结构控制为Lyapunov泛函法。有关时延系统近30年的研究表明,必须建立一套新的概念和方法,这有待于从事此方向研究的学者们进一步探讨。
    星载微波遥感地表特征参数反演的算法研究
    金亚秋
    1997, 14(S1): 722-722. 
    摘要 ( )   PDF(50KB) ( )  
    星载遥感的进展为人类了解自身所处的地球环境在广大的空间尺度和时间尺度上,提供了十分丰富的信息。由于微波遥感能全天候工作,并能给出地球复盖圈的内部信息,在近20年里微波星载遥感取得了很大的进展,并成为直至下世纪初最主要的星载遥感和空对地观察的主要手段之一。但是如何从星载遥感数据中提取信息,反演地表特征性参数,仍是一个留待深入研究的应用基础问题。它包括了遥感观察的正向模拟仿真,以及反向求逆反演二个方面的问题。介绍了星载微波遥感数值模拟和参数反演的一些进展。特别介绍星载微波遥感数据标定和验证算法上的一些新方法,展示一些新结果。这些成果已应用于全球特别是中国区域上地表特征性参数的反演提取。
    动力天文中常用多步方法的比较
    赵长印
    1997, 14(S1): 722-722. 
    摘要 ( )   PDF(50KB) ( )  
    在动力天文数值研究中,最早被广泛采用的是Adams-Cowel(简称AC)多步积分方法。90年,对称(简称SY)多步方法问世。国外的卫星定轨软件中,应用了KSG多步积分器。文中从几个方面对这三种常用的多步方法作了比较,得到的结论是:尽管对称方法能量误差随时间变化的特性好,但在相同条件下,它的计算精度较差,包括系统能量的精度,且不适用于显含速度的微分方程;KSG方法无论从能量误差的累积特性,计算精度,还是从计算效率上,均优于AC方法。这些结果为动力天文中积分器的选择提供了依据。KSG方法的辛性质值得研究。
    关于色散方程Ut=aUxxx的显-隐差分格式
    张大凯, 田益民
    1997, 14(S1): 723-723. 
    摘要 ( )   PDF(58KB) ( )  
    对色散方程Ut=aUxxx的混合问题,构造了一类显-隐差分格式(1)和(2)其中r=/h3是网比,τh分别表示时间步长和空间步长。它们的截断误差阶均为o(τ+h),稳定性条件分别为o < r≤(1)/(2)和-(1)/(2)≤r < 0.设Ih={x|x=kh,0≤k≤2M+1}。当a>0时,每层按节点的空间坐标交替使用(1)之上下式从左向右计算,对a<0时则使用(2)式从右向左算。避免了通常的差分格式在使用时对附加边界值的依赖问题,并具有显式计算的优点。数值例子验证了分析结果。
    格子Boltzmann方法的外力模型
    陈炬桦, 李元香, 熊盛武
    1997, 14(S1): 723-723. 
    摘要 ( )   PDF(58KB) ( )  
    在《基于格子Boltzmann方法的流体力学问题仿真》一文中详细介绍了格子Boltzmann方法(LBM)的通用模型;模拟了一些典型流体流动现象,证明LBM模拟流体流动的简单有效性,适合基于区域分解的分布式算法;首次提出并实现了基于LBM的流体力学问题的自动建模算法。现设计一个九点正方形LBM的外力模型。
    RLW方程的特征有限元格式及显式并行算法
    谢树森, 唐诗
    1997, 14(S1): 724-724. 
    摘要 ( )   PDF(57KB) ( )  
    考虑如下RLW方程的初边值问题ut+g(x,u)·▽u-δΔut=f(x,t),(x,t)∈Ω×(0,T];u(x,t)=0,(x,t)∈∂Ω×(0,T];u(x,0)=u0(x),x∈Ω.(1)其中Ω⊂Rd(d≥1,整数)为有界区域,g=(g1,…,gd),T,δ为正常数。
    MBE过程中杂质扩散系数新的确定方法
    刘英, 程文芹
    1997, 14(S1): 724-724. 
    摘要 ( )   PDF(57KB) ( )  
    分子束外延(MBE)过程中杂质扩散系数的确定,可以归结为动边界杂质扩散问题的反问题#br#ut=Duxx x>v(t0-t), t0>0, t>0 (1)#br#u(x,0)=n0 x>vt0, n0=Const.(2)#br#u(v(t0-t),t)=N0 t>0,N0=Const.(3)#br#u(x,t0)=N(x) x>0(4)#br#我们采用粒子数守恒定律精确地确定样品的生长速率;利用有限差分法计算杂质的扩散系数。对于GaAs中掺Be样品的计算结果表明:该方法可以简化传统的工序;确定生长速率的方法较振荡法精确;提高了扩散系数的计算精度。
    蒙特卡罗求积中对偶变数法的一个注记
    罗强
    1997, 14(S1): 725-725. 
    摘要 ( )   PDF(75KB) ( )  
    应用Monte Carlo求积时,对偶变数法是降低方差技巧之一。然而在一般文献中仅限于[0,1]区间的单调函数。为此推广了任意区间非单调函数的情形。
    Burgers方程的一个新的组显格式
    曾文平
    1997, 14(S1): 725-725. 
    摘要 ( )   PDF(54KB) ( )  
    Burgers方程的一个新的组显格式曾文平(华侨大学数学系,泉州362011)摘要以求解Burgers方程的显式逆风格式为基础,构造了一组新的分组显式格式。讨论了方法的线性化稳定性。数值结果表明,此方法优于Evans的分组显式格式。
    混合自旋模型动力学相距的时间演化
    刘策军, 胡嘉桢
    1997, 14(S1): 725-725. 
    摘要 ( )   PDF(75KB) ( )  
    考察混合自旋模型两种组态按Metro polis动力学随时间的演化。计算给出了T=1.55,1.65,1.75,1.85,1.95,2.05,2.15,2.25八个温度下其在平方晶格上动力学相距随时间演化的关系,即In<D(t)>-In (t)曲线。结果显示出T >1.95与T <1.85两个不同的高、低温度区域内曲线的不同走向,从而揭示该模型相变的存在。相变温度近似值即为Tc≅1.85-1.95。
    ECR等离子体中电子运动的计算机模拟
    宋汝安, 程先安, 祁祥麟
    1997, 14(S1): 726-726. 
    摘要 ( )   PDF(52KB) ( )  
    电子回旋共振(ECR)等离子体受到沿磁场递减方向的推力而加速,电子在运动中获取高能。对矩形微波谐振腔TE101模高频电磁场和Z轴向递减的非均匀磁场中的ECR等离子体的电子运动,用电子计算机进行模拟计算。从算出的电子运动轨迹和运动曲线可以看到:ECR运动激发的等离子体最宽是微波半波长。电子在运动中受磁场递减方向推力的最大位置不在共振点,而在共振点之前,随共振点的位置而稍有不同,推力作用区也不超过微波半波长。电子通过腔体从电磁场吸取能量,当共振点在微波驻波波腹处时,获得能量最大。
    泰氟隆烧蚀计算及化学反应模型对比分析
    高铁锁, 曹文祥
    1997, 14(S1): 726-726. 
    摘要 ( )   PDF(52KB) ( )  
    采用化学非平衡边界层方程组与烧蚀壁面耦合求解的方法来确定烧蚀壁面条件和化学非平衡绕流,针对泰氟隆材料进行了详细计算,研究了化学反应模型、化学反应速率、壁面催化特性等对流场和电离特性的影响,并和无烧蚀情况进行了对比,得到一些有用的结果。计算中采用的化学模型为31个组元,33~35个化学反应,数值方法为全隐格式。计算结果与有关文献及飞行试验进行了比较,结果合理。
    稳态对流扩散方程的斜迎风格式
    芮洪兴
    1997, 14(S1): 726-726. 
    摘要 ( )   PDF(52KB) ( )  
    考虑稳态对流扩散方程的Dirichlet边值问题-▽(a(x)▽u)+b(x)·▽u=0,x=(x1,x2)∈Ω,(1)其中b·▽u=b1(u)/(x1)+b2(u)/(x2。当0≤a(x)<<|b(x)|时,对流项占主导地位。通常的数值方法求解会出现严重的数值振荡,为此人们提出了迎风格式的技巧。
    振动系统的增容问题
    黄贤通
    1997, 14(S1): 727-727. 
    摘要 ( )   PDF(48KB) ( )  
    提出弹簧-质量振动系统的增容问题:增加现有振动系统的容量,即增加弹簧和质量的个数,使系统具有预先给定的固有频率。讨论了两种增容方式,给出了解存在的充分必要条件,且给出了数值算法。
    求解扩散方程的高精度差分格式
    张天德, 李久平, 王玮
    1997, 14(S1): 727-727. 
    摘要 ( )   PDF(48KB) ( )  
    用待定系数法只给出扩散方程的一个高精度差分格式,证明了差分格式的收敛性和稳定性的条件。数值实验结果表明,差分格式具有较好的稳定性和很高的精度,因而计算结果与精确结果吻合较好。
    半导体器件数值模拟的特征有限元方法
    车中原
    1997, 14(S1): 727-727. 
    摘要 ( )   PDF(48KB) ( )  
    半导体器件数值模拟需要研究扩散模型的非线性偏微分方程组的初边值问题,袁益让教授曾研究牛曼边界和混合边值问题,提出有限元和混合元的计算方法并得到H1误差估计。在周期性条件的情况下,用有限元来逼近位势方程,用特征有限元格式来逼近浓度方程,得到最佳阶的H1模误差估计和最佳阶L2模误差估计。
    静电式电子能量分析器的计算机辅助设计
    周庆, 马琳
    1997, 14(S1): 728-728. 
    摘要 ( )   PDF(52KB) ( )  
    在扫描电镜中,具有一定能量的电子束入射到被测样品上会产生很多信息。充分利用这些信息对于了解材料的结构及性能具有重要意义。静电式电子能量分析器是用于半导体材料及微电子器件分层体视及内窥的一种扫描电镜上的新型检测器。在此介绍了电子能量分析器的设计原理及能量分辨率的计算方法,并用计算机对电子能量分析器进行了辅助设计。软件用VisualBasic编写,界面友好,操作方便,可同时显示数值计算结果、电极形状及电子轨迹模拟结果,具有快速、准确、直观的优点。还给出一个用该软件设计的实例结果。
    洪特定则对同科P2电子组态的L-S耦合的精细结构同样适合
    段萍
    1997, 14(S1): 728-728. 
    摘要 ( )   PDF(52KB) ( )  
    洪特定则是关于从同一电子组态形成诸能级的上下次序问题的定则。它只适用于L-S耦合洪特定则的内容如下:"从同一电子组态形成的能级中(1)那重数最高的亦即S值最大的能级位置最低(2)具有相同S值的能级中,那具有最大L值的能级位置最低。"而洪特定则对同一电子组态且具有相同L值和S值,但J值不同的能级次序问题没有给出。以同科P2电子为例通过对各量子数的定量计算的方法确定出能级最低态,从而定出能级的次序问题。
    复杂目标雷达特性的可视化建模
    刘铁军, 周勇
    1997, 14(S1): 728-728. 
    摘要 ( )   PDF(52KB) ( )  
    提出一种新的把计算机几何辅助设计(CAGD)、计算机图形学及高频电磁散射法结合起来计算目标雷达特性的方法。首先构造目标外形的非均匀有理Bezier曲面或非均匀有理B样条曲面(NURBS)模型。设计特定的Phong光照模型,得到目标的用Z-Buffer方法消隐后的彩色图象,其上每个象素点的RGB值等于目标表面对应点法矢的三个分量。根据目标面上每个点的六维信息(X,Y,Z;Nx,Ny,Nz),利用物理光学(PO)、物理绕射理论(PTD)及阻抗边界条件等高频电磁散射方法,计算目标的雷达散射截面(RCS)。根据一维高分辨(宽带)雷达照射目标的物理图象,提出一种计算一维距离图象的可视化方法,并研制成功"军用目标几何造型及其雷达特征的可视化模型"软件,计算了大量军用目标的雷达散射截面及一维距离像,与实验相比,吻合较好。
    一类具有几何对称性的非协调矩形板元
    石东洋, 黄艾香
    1997, 14(S1): 729-729. 
    摘要 ( )   PDF(55KB) ( )  
    利用扰动方法,通过摆脱传统的ACM元的位移模式,推出了一类非常规矩形板元,其不仅协调性好,列式简单对称,保持了单元的几何不变性,而且还充分显示了很好的数值效果。
    地下核爆炸围岩能量耦合数值研究
    韩学安, 陈万瑞, 李振东
    1997, 14(S1): 729-729. 
    摘要 ( )   PDF(55KB) ( )  
    采用弹塑性流体动力学模型的一维球对称Lagrange有限差分方法,数值研究了地下核爆炸围岩介质中能量耦合。对1kt爆炸当量,求解了0.5~10m不同爆室半径下不同区域中耦合能量的分布。计算结果表明:(1)当爆室较小时,爆室初压较高,在围岩中可以产生汽化、液化、破碎、破裂、弹性等不同的效应区。随着爆室增大,汽化等效应区会逐渐消失,当爆室增大到一定程度,围岩仅保持单一弹性区。(2)地下核爆炸的大部分能量耦合在近区,少部分能量以弹性能形式输运到远区。(3)爆室由小变大,空腔中存留能量率也由小逐渐增大,而围岩能量耦合率由大变小。当爆室半径增大到一定程度时,爆炸总能量的绝大部分存留于空腔中,耦合到围岩中的能量很少很少。(4)地下核爆炸中,弹性区能量耦合率较小。随着爆室增大,能量耦合率出现局部增强效应,最后,能量耦合率越来越小,甚至可能趋于零,即出现完全解耦爆炸。
    小波分式滤波器
    匡正
    1997, 14(S1): 730-730. 
    摘要 ( )   PDF(49KB) ( )  
    采用有理分式形式的滤波器来构造小波。这种分式滤波器具有如下优点:1.不仅可以构造出包括Daubechies小波在内的小波函数,而且能包括B样条小波;2.可构造出具对称性的小波函数;3.通过恰当的构造,小波函数正则阶比Daubechies小波有显著的提高。
    二阶非线性偏微分方程组的一个算法
    张流生
    1997, 14(S1): 730-730. 
    摘要 ( )   PDF(49KB) ( )  
    针对应用牛顿运动定律建立的一类二阶非线性偏微分方程组没有解析解提出了一个基于速度加速度定义易于掌握和应用的简便算法。
    带电粒子束传输的相图函数计算
    郁庆长
    1997, 14(S1): 731-731. 
    摘要 ( )   PDF(50KB) ( )  
    粒子束的相图是束粒子在相空间所占据的区域。相图函数描述相图边界的形状。当粒子束在非线性电磁场中运动时相图函数不断变化,这种变化属于保辛变换。利用Lie变换方法可以直接计算粒子束相图的变化而无须跟踪各个粒子的运动。此方法也可用于解逆问题,即由系统入口和出口处的相图函数求系统相应的变换。逆问题有无限多解,但可根据实际情况选取技术上易于实现的解。
    ECR微波等离子体离子输运过程的数值模拟
    宫野, 温晓军, 张鹏云, 邓新绿
    1997, 14(S1): 731-731. 
    摘要 ( )   PDF(50KB) ( )  
    建立了ECR微波等离子体源离子输运的平板和柱面模型,对离子历经的空间区域的输运过程进行了数值研究。在平板模型下,采用M-C方法模拟了离子离开放电室后历经中性区、鞘层区,最后被加负偏压的工件表面吸收的全过程。考虑了离子与中性粒子的电荷交换碰撞和弹性散射,统一处理了中性区和鞘层区电势的衔接。先前方法是采用经验公式,连接不光滑,理论上也不自洽。采用曲线拟合,电势自洽的方法成功地把中性区和鞘层区衔接起来,得到自洽电势分布和鞘层区不同位置处的速度分布、能量分布及角分布。在圆柱模型下,提出了二次电子发射形成"屏蔽区"的概念,讨论了磁场、气压以及二次电子发射对鞘势的影响。
    ZH型蒸烘机的热力设计
    赵炳俊, 胡广芝, 胡广耀, 赵辉, 刘五奎, 师松健, 蒿建丽
    1997, 14(S1): 731-731. 
    摘要 ( )   PDF(50KB) ( )  
    利用传热、传质及热辐原理对ZH型蒸烘机进行了热力设计计算。ZH型蒸烘机是用于国内、外油脂工业中大豆一次浸出制油工艺的新型蒸烘设备。根据该蒸烘机的结构、工作原理及被干燥物料的特性,运用计算机计算该机输入、输出的热量;间接蒸汽的热量和消耗量及加热面积等。从理论和实践上证明了该蒸烘机具有溶耗低,脱溶、脱水效果好,产品蛋白含量高,干燥能力强,热效率高等优点和具有多吸脱溶、微负压运行优良性能。在油脂、饲料等工业中得到了广泛应用。
    板式射流烘干机的设计研究
    胡广芝, 赵炳俊, 胡广耀, 刘五奎, 赵辉, 张剑, 蒿建丽
    1997, 14(S1): 732-732. 
    摘要 ( )   PDF(51KB) ( )  
    板式射流烘干机是利用固体燃料高效热风炉供热,采用多孔板式热风射流装置实施高温作业的新型、高效、节能烘干设备,利用空气热力学,流体力学和传热、传质原理对该烘干机和射流热板换热器进行了热力设计计算,并以被干燥物料高岭土为例,经测试和计算,得出的主要技术经济指标:处理量为920kg/h;脱水能力达6~20%;干燥强度为2.52kg/m2;干燥能力为6.28t·%/h;热效率为21%等,具有国内先进水平。解决了在工业生产中,对众多非热敏性在高温下实施干燥作业的难题。具有广泛的推广应用前景。
    非归-分布随机抽样方法比较
    程锦荣, 段香梅
    1997, 14(S1): 732-732. 
    摘要 ( )   PDF(51KB) ( )  
    从实际计算的角度,对非归一分布的四种具有代表性的随机抽样方法进行了评价。通过选择八种离散型和连续型的非归一分布,采用Metropolis方法、期望方法、排队方法和舍选方法分别对它们进行随机抽样,并根据抽样分布、收敛速度、拟合性检验结果、抽样效率、花费CPU时间、抽样中进入新状态的样本数目和所产生子样的随机性质等指标对上述四种方法作了仔细的比较。计算结果表明,期望方法和舍选方法明显地优于Metropolis方法和排队方法,且期望方法更适合于抽取大样本,舍选方法更适合于抽取小样本。同时还就非归一分布随机抽样方法的应用技巧问题作了若干讨论。
    700hPa高度场距平的长期预报试验
    周小刚
    1997, 14(S1): 732-732. 
    摘要 ( )   PDF(51KB) ( )  
    象长期预报研究的是大尺度分量和平均行为,本质上是概率性的。研究方法可归结为统计学和动力学方法。动力学方法考虑了大气内在的运动规律,但由于其复杂性尚难应用于业务,如目前大部分动力长期预报误差是来自于模式的"气候漂移",即系统误差。经验统计学方法、特别是具有动力学意义的统计方法有巨大潜力。采用多通道奇异谱分析和传统相似预报的组合模型,探讨了利用700 hPa前期高度场距平直接预报后期30天平均格点值的可行性。
    正态分布近似式的研究
    郑金成, 蔡淑惠, 王仁智, 郑永梅, 程灿东
    1997, 14(S1): 733-733. 
    摘要 ( )   PDF(54KB) ( )  
    正态分布是一种重要的概率分布,在许多情况下,简单的近似解析方程有助于问题的解决。通过将理论计算与数值模拟相结合的方法给出了一个较好的近似式。
    有限区气候模拟侧边界嵌套误差分析及边界条件的研究
    颜宏, 赵俊英
    1997, 14(S1): 733-733. 
    摘要 ( )   PDF(83KB) ( )  
    为实现有限计算机条件下尽量高分辨率的区域气候模拟可以采用嵌套的方法。嵌套误差从本质上讲是由于侧边界条件的不适定性造成的,即对目前区域气候模式中多采用的静力平衡假定下的原始方程组,理论上讲并不存在适定的侧边界条件,各种侧边界处理方案都只是消除或抑制计算波使模式计算稳定的手段,因此嵌套误差的存在是必然的。区域气候模拟除存在以上嵌套误差的共性外,还有其特殊性,主要是长时间中积分误差的积累性。该文目的旨在对区域气候模拟中嵌套误差的种类、来源及其时空分布和时间积累性特点作一比较系统的分析,并在此基础上提出减少误差的新思路、新方法。
    有限元超收敛的新发现
    陈传淼
    1997, 14(S1): 734-734. 
    摘要 ( )   PDF(55KB) ( )  
    我们发现有限元超收敛的一种新结构,即所有单元的角点不是超收敛点。记unL2投影为uh,e=uh-u,D是平均梯度。Ω的剖分均匀,Ωh={x∈Ω,dist(x,∂Ω)≥ch|lnh|}。
    带记忆热传导方程的质量集中方法
    朱广冲
    1997, 14(S1): 734-734. 
    摘要 ( )   PDF(55KB) ( )  
    质量集中方法是一种修正的Galerkin方法,在具有相同精度的前提下,有着计算格式简洁,计算量小的优点。其基本思想是,在抛物型偏微分方程的标准Galerkin格式中,对于时间导数项采用一个近似内积作逼近。这实际是一种数值积分作的近似,经过这样处理之后,在离散格式等价的代数方程组中,质量矩阵成为对角矩阵,因而求逆矩阵变得十分简单,实际计算时大大减少了计算量。
    涡流检测中的微扰模型
    段耀勇, 陈荣华, 盛剑霓, 陈雨生
    1997, 14(S1): 963-695. 
    摘要 ( )   PDF(146KB) ( )  
    根据低频电磁场在导体中的趋肤效应及基准模型中激励线圈相对较小的特点,将轴对称场与三维场相结合,建立了以裂纹形成的微扰场作为求解对象的微扰模型
分享到: