Please wait a minute...

过刊目录

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

    论文   
    全选选: 隐藏/显示图片
    论文
    双参数波动方程反问题的数值解法
    张文飞, 李晓江
    1997, 14(2): 129-130. 
    摘要 ( )   PDF(254KB) ( )  
    提出在微分方程反问题的数值解法中可由一种反问题补充条件同时反求两个未知参数的观点和方法,并以一维波动方程为例推导了震源和岩性联合反演的详细算法。从概念上突破了传统的一种补充条件只能解一个未知数的反演理论的约束,解决了地震勘探中波动方程反问题的震源未知或测不准的矛盾,缩短了反演理论研究与工程实际应用间的距离。并给出相应的数值算例。
    三维旋转域内Poison方程的有限元线法分析
    陆坤, 袁驷
    1997, 14(2): 131-142. 
    摘要 ( )   PDF(281KB) ( )  
    应用半解析有限元线法(FEMOL)分析三维旋转域内的Poisson方程。通过环向采用Fourier级数展开,原问题被转化为一系列旋转面上的二维问题。构造了相应的FEMOL单元,并对由柱坐标系带来的r=0处的奇异性以及各类退化单元边、线作了细致处理。并给出数值算例。
    初至到时波前算法及在地震走时CT中的应用
    宋霖平, 刘浩吾
    1997, 14(2): 143-147. 
    摘要 ( )   PDF(356KB) ( )  
    基于惠更斯原理,给出了一个简练的初至到时算法,并应用于地震走时层析成像。结果表明,新方法计算效率较Dijkstra最短路径法高一个数量级。
    非正交网格上的九点格式在热传导问题中的改进
    段庆生, 袁国兴
    1997, 14(2): 148-154. 
    摘要 ( )   PDF(256KB) ( )  
    对文献[1]在二维Lagrange流体力学网格上构造的扩散方程九点差分格式做了进一步的讨论和改进,给出了一般形式边界条件的计算格式。数值试验的结果表明,这些改进提高了原格式数值结果的精度。
    三维SPH程序及其在斜高速碰撞问题的应用
    贝新源, 岳宗五
    1997, 14(2): 155-166. 
    摘要 ( )   PDF(447KB) ( )  
    叙述了平滑粒子流体动力学(简称SPH)的基本概念,重新推导了所用的基本方程。介绍了具体计算格式、算法以及策略等。对所编制的程序作了一定的检验。对SPH方法中的人为粘性和人为热流的作用作了数值考查。还提供了一个三维高速碰撞现象的算例。最后指出尚存在的若干问题。
    奇异摄动数值法及其在连铸结晶器内凝固传热分析中的应用
    赵西爱, 孙仁济, 曲英, 吕世意
    1997, 14(2): 167-170. 
    摘要 ( )   PDF(181KB) ( )  
    在奇异摄动数值法用于冶金过程中,提出了板坯连续铸锭在结晶器中凝固传热的一种数值模拟方法。把连续铸锭结晶器部分的传热模型归结为一边界层型奇异摄动问题,并用奇异摄动数值方法进行分析、计算,得出了结晶器内铸坯的温度分布、液相穴宽度、坯壳厚度等参数。
    数值模拟二维槽道内Oldroyd-B流体的流动
    王保国, 王栋, 武弘峨
    1997, 14(2): 171-178. 
    摘要 ( )   PDF(288KB) ( )  
    提出一种求解Oldroyd-B流体流动的数值方法。该方法采用有限体积离散,在求解连续方程时引进Chorin的人工可压缩性;求解本构方程时考虑了迎风处理的技巧;离散后的方程采用Gauss-Seidel线松弛迭代求解全流场。典型算例表明:该方法可用于非线性粘弹流体的数值计算,所得的二维算例结果符合非牛顿粘弹性流体的流动特征。
    非结构网格生成Bowyer-Watson方法的改进
    曾扬兵, 沈孟育, 王保国, 刘秋生
    1997, 14(2): 179-184. 
    摘要 ( )   PDF(283KB) ( )  
    完善了应用Delannay原理生成非结构网格的Bowyer-Watson方法。通过求解简化Pois-sion方程,灵活地控制网格的疏密;通过网格层次的区分,提高了网格的贴体性能;并通过分区,实现了适用于粘流计算的网格生成。
    辐射传输问题的P2扩散近似
    冯庭桂
    1997, 14(2): 185-189. 
    摘要 ( )   PDF(199KB) ( )  
    用渐近展开方法,给出和讨论辐射传输问题的P2近似扩散方程,并进行了数值模拟计算。结果表明,P2扩散方法优于普通扩散和限流扩散方法,其数值结果与输运方程结果十分接近。
    守恒方程的差分型格子气方法
    阎广武, 胡守信, 施卫平
    1997, 14(2): 190-194. 
    摘要 ( )   PDF(189KB) ( )  
    给出了守恒型模型方程的差分型格子气方法,通过定义广义的分布函数,将原方程的求解问题转化成关于分布函数的差分算法。该方法的意义在于建立格子气方法与差分法之间的联系,进而可以构造真正用于空气动力学的格子气模型,使格子气方法在可压缩流领域成为一种实用方法。分别构造了两步MacCormark格式和Warming-Beam格式的分布函数的差分方程。
    研究量子耗散配分函数的变分数值方法
    陆则慰, 李秀玲, 包景东
    1997, 14(2): 195-201. 
    摘要 ( )   PDF(262KB) ( )  
    提出用变分路径积分方法计算量子耗散系统的配分函数。将势能和耦合函数分别写成谐振子势(圆频率Ω2(x0))和线性函数(斜率μ(x0))与剩余项之和,其中所引入的两个试探函数Ω2(x0)和μ(x0)由有效势取极小值的条件所确定。则最终的配分函数为一个高斯密度函数下的平均值。并将本算法与蒙特卡罗方法进行了比较。
    非灰气体中辐射全交换面积的计算
    聂宇宏, 陈海耿, 杨泽宽
    1997, 14(2): 202-206. 
    摘要 ( )   PDF(189KB) ( )  
    将描述非灰气体辐射特性的指数宽带模型引入均匀离散射线法,推导了非灰气体参与介质条件下,辐射全交换面积的计算式,并计算了二维封闭体系的辐射全交换面积。结果表明,均匀离散射线法可用于模拟非灰气体中的辐射换热问题。
    超音速连续波COIL的二维半数值模拟
    张本爱, 陈云山, 李守先, 莫俊永, 刘广华
    1997, 14(2): 207-216. 
    摘要 ( )   PDF(426KB) ( )  
    给出了超音速连续波(CW)氧碘化学激光器(COIL)的理论模型,利用此模型对典型实验装置进行了二维半几何近似下的初步数值模拟研究,对所得结果进行了讨论,并与实验或文献结果作了比较。
    用于处理复杂对称结构电磁散射的扩展GIM技术
    朱峰, 任朗
    1997, 14(2): 217-220. 
    摘要 ( )   PDF(166KB) ( )  
    由于原广义镜像法只能够处理具有四面对称的电磁散射问题。为此,在广义镜像法原理的基础上,给出了处理一种具有特殊复杂对称结构的电磁散射问题的新方法,称之为扩展GIM方法。对比结果表明:利用扩展GIM,较之于原GIM,在对高于四面对称的散射体的求解过程中,能够节约更多的上机运算时间和运算所需的存储;还通过该方法对正六面柱的具体分析,给出了应用实例。
    200MW惯性约束混合堆燃耗计算
    刘忠兴, 伍钧, 吕荀, 黄文凯, 任捷, 曾文平
    1997, 14(2): 221-226. 
    摘要 ( )   PDF(282KB) ( )  
    建立不生产239Pu的惯性约束聚变裂变混合电站是可行的。我们使用BISON程序[4]计算了200MW电站主要特征量随时间的变化,计算结果表明裂变包层并不总是次临界的。讨论了既维持高的功率输出又能避免超临界的办法。
    求解极限环的谱展开法
    曾晓东, 齐翔林, 汪云九
    1997, 14(2): 227-232. 
    摘要 ( )   PDF(215KB) ( )  
    提出一种求解自治非线性常微分方程周期解的方法——谱展开法。它直接把解形式上写成Fourier级数,其频率和各谐波成份待定,而将原非线性常微分方程问题变换成求解相应的非线性代数方程组,还以van der Pol方程和神经元群平均场方程为例对该方法进行说明和讨论。
    一种高精度差分格式
    廖成, 赵愉深, 林为干
    1997, 14(2): 233-237. 
    摘要 ( )   PDF(263KB) ( )  
    在微分——Thompson变换结合时域有限差分(FDTD)技术计算复杂目标的电磁散射特性的方法中[1],差分格式的构造和选取与解的精度存在着密切关系。提出一种新的高精度差分格式,其数值实现进一步证实了该方法能精确模拟任意形状目标的电磁散射过程。
    三维寄生电阻电容直接边界元计算
    王泽毅, 周文明, 元彦宏
    1997, 14(2): 238-241. 
    摘要 ( )   PDF(199KB) ( )  
    应用直接边界元素法求解表征寄生电阻电容的混合边界条件拉普拉斯方程,处理角点法向电场间断是一个难点。作者曾提出多重法向导数方法处理2-D角点问题,取得良好效果[2,4]。本文详细介绍多重法向导数方法在3-D情形的推广,并给出数值例子。
    氮氧化硅与氧化硅玻璃结构,性能与量子化学计算研究
    闵新民, 邓志平, 赵修建
    1997, 14(2): 242-246. 
    摘要 ( )   PDF(161KB) ( )  
    用自洽场离散变分Xα(SCF-Xα-DV)量子化学计算方法研究了氮氧化硅及氧化硅玻璃,讨论了氮原子取代氧原子之后在结构、性能与化学键等方面的变化规律。计算表明氮氧化硅与氧化硅玻璃之间在性能上的差异,主要不是由Si-N和Si-O键的伸缩力常数之间的差异所引起的,也不是由Si-N和Si-O键的离子键强度之间的差异所引起的。Si-N共价键比Si-O共价键强,是N-Si-N键弯曲力常数比O-Si-O键大的主要原因,也与氮氧化硅玻璃比氧化硅玻璃具有更好的化学稳定性和更高的扬氏模量的实验结果一致。
    网格间距及正交性控制在翼-身-尾组合体三维绕流中的研究
    张正科, 庄逢甘, 朱自强, 罗时钧
    1997, 14(2): 247-252. 
    摘要 ( )   PDF(376KB) ( )  
    用椭圆型方程生成网格。在二维发展两种Hilgenstock方法实现对网格与边界间距及与边界正交的控制,并将其思想推广到三维,生成了翼-身-尾组合体与公共交界面或边界面正交的三维贴体分块网格,用于Euler方程流场计算,效果良好。
    KENO Ⅳ交互绘图功能的配制
    王先龙, 高彦锋
    1997, 14(2): 253-256. 
    摘要 ( )   PDF(145KB) ( )  
    介绍了为KENOⅣ配制的交互绘图功能的有关内容,说明了用户界面、几何可视化(绘制二维剖面图和三维立体图)、Keff结果图、运行进程显示等功能。
分享到: