Please wait a minute...

过刊目录

    1999年, 第16卷, 第6期
    刊出日期:1999-11-25
    上一期    下一期

    论文   
    全选选: 隐藏/显示图片
    论文
    定常Stokes问题的边界积分方程的高精度求积方法与外推
    吕涛, 黄晋
    1999, 16(6): 561-567. 
    摘要 ( )   PDF(257KB) ( )  
    借助Side-Israeli 的求积公式,给出解Stokes 问题的边界积分方程的机械求积方法。此方法精度较高、计算量少,近似解的误差有奇数幂的渐近展开,这表明使用Richardson 外推能够改善近似解的精度阶。
    用神经网络分辨“膝”区原初质子成分
    冯存峰, 孔繁敏, 张学尧, 何瑁, 戴志强, 张乃健
    1999, 16(6): 568-572. 
    摘要 ( )   PDF(212KB) ( )  
    利用不同的相互作用模型,对羊八井空气簇射阵列和乳胶室联合实验进行了Monte Carlo模拟,并利用模拟数据研究了利用人工神经网络方法进行原初宇宙线成分分辨的可行性。分析表明,利用该方法可很好地挑选出由原初质子引起的簇射,但所得结果与作用模型有关,并对相应的系统误差进行了估计。
    对宇宙线广延大气簇射Čerenkov光的计算机模拟
    查敏, 朱清棋, 丁林恺
    1999, 16(6): 573-579. 
    摘要 ( )   PDF(295KB) ( )  
    结合高能宇宙线广延大气簇射研究的需要研制了能够模拟广延大气簇射的强子成分,电磁成分和Čerenkov 光的Monte Carlo 软件。模拟程序的建立不仅能够针对梁王山阵列的实验结果进行相关的物理分析,而且具有共通性,可提供其它广延大气簇射实验产生对照的模拟样本。
    中子裂变链统计涨落问题的数值计算方法
    徐乃新, 汤敏君
    1999, 16(6): 580-586. 
    摘要 ( )   PDF(256KB) ( )  
    研究了弱中子源驱动下,裂变系统中子裂变链统计涨落问题的数值计算方法。进行了数值计算建模、数值方法分析、数值计算检验、一类问题概率分布函数的统计涨落特征量的数值计算示范。特例数值检验表明:只要数值解方程组阶数(截断) N足够大,数值解满足归一(守恒) 律、指数增长律,并与精确解析解一致。对于非定常裂变系统中子裂变链统计涨落问题提出了一维等效模型下数值模拟的方法。
    有分散凸起热源矩形腔内熔化的数值模拟
    蒋宏利, 陈钟颀, 徐通模, 惠世恩
    1999, 16(6): 587-591. 
    摘要 ( )   PDF(207KB) ( )  
    利用焓法求解了有三个分散凸起热源的矩形空腔内的熔化问题,并把计算结果与前人的实验结果进行了比较。计算和实验结果都显示,矩形空腔内的分散凸起热源对工质的熔化界面形状有显著的影响。
    纤维层中纤维的排列与散射相函数
    李炳熙, 别如山, 陆惠林, 杨励丹
    1999, 16(6): 592-598. 
    摘要 ( )   PDF(301KB) ( )  
    针对纤维层中两种典型的纤维排列求得其相函数:纤维在空间随机排列和平行于边界面的平面上随机排列。被考虑的两种入射辐射为垂直于边界面的入射和边界表面的表面辐射。纤维层的相函数以纤维层坐标系中的天顶角的函数形式表示。对入射辐射形式,尺度参数,纤维的排列及复折射率对辐射传递的影响进行了研究。
    磁场中反应热对电弧螺旋不稳定性的影响
    宫野, 张鹏云
    1999, 16(6): 599-605. 
    摘要 ( )   PDF(248KB) ( )  
    采用抛物型电导率、抛物型透明辐射和与时间有关的线性微扰理论,考虑存在反应热情况下,给出了热势分布(温度分布)、稳定性条件以及不稳定性增长率等定量结果。得出反应热对电弧螺旋不稳定性的影响是降低稳定性的结论。
    用MONTE CARLO方法计算托卡马克装置的中性氢原子密度分布
    张先梅, 李有宜, 万宝年
    1999, 16(6): 606-609. 
    摘要 ( )   PDF(172KB) ( )  
    简要描述了计算托卡马克中性氢密度分布的Monte-Carlo 方法。计算了HT-6B、HT-6M装置的中性氢原子密度分布,与实验结果进行了比较。
    激光腔靶中等离子体流场的数值模拟
    段庆生, 常铁强, 张维岩, 王光裕, 王重基
    1999, 16(6): 610-615. 
    摘要 ( )   PDF(271KB) ( )  
    针对激光腔靶中等离子体流场的复杂性,提出了一种网格的自适应算法以及活动网格方法。应用这些方法在不同的腔靶上进行了一些数值计算。得到的结果表明,这些数值方法是有效的。
    非交错网格下三维准稳态激光重熔熔池数值模拟
    曾大文, 黄开金, 谢长生
    1999, 16(6): 616-623. 
    摘要 ( )   PDF(322KB) ( )  
    对三维准稳态激光重熔问题进行数值模拟,模型考虑了熔池内金属液体的流动和换热。在非交错网格下,采用修正的动量插值技术和SIMPLE算法求解动量方程。数值结果表明表面张力温度系数∂γ/∂T 对熔池内金属液体的流动方式有重要影响。当∂γ/∂T 为-3-5×10-6 时,熔池中只存在沿顺时针方向流动的主对流环;当∂γ/∂T为-3.5×10-4 时,主对流环诱导产生了二次对流环。当∂γ/∂T为+3.5×10-4 时,与负的表面温度张力系数情形相反,熔池中只存在沿逆时针方向流动的主对流环。通过灰铸铁激光相变硬化实验,实测了硬化层深度和宽度。计算结果和实验结果进行了比较。
    准旋-角动量标量算子本征值与j壳层耦合态的完全分类计算
    李晓梅, 陈健华
    1999, 16(6): 624-629. 
    摘要 ( )   PDF(227KB) ( )  
    应用二次量子化和不可约张量理论,由4 个产生-湮灭算子耦合成准旋-角动量标量算子Z(K),K为参数。计算表明:Z(K) 的本征值有很强的分类功能,能对j≤15/2 壳层准旋-角动量耦合态完全分类(用准旋和角动量只能对j≤7/2 壳层耦合态完全分类)。讨论Z(K) 的构造和它的本征值的计算方法,给出j≤15/2 壳层耦合态的完全分类。
    周期三对角线性方程组的分布式并行算法
    迟利华, 刘杰, 李晓梅
    1999, 16(6): 630-637. 
    摘要 ( )   PDF(281KB) ( )  
    提出一种求解严格对角占优周期三对角线性方程组的并行算法(简称PAA 算法),新算法计算复杂性为O(8n),通讯复杂性为O(1).目前求解此类方程组的最优并行算法的计算复杂性为O(17n),通讯复杂性为O(logP).在SGI Indy 工作站网络环境下的试算结果表明,加速比呈线性增加,并行效率达到90%.
    采用多重网格法进行电子结构的计算
    王一中
    1999, 16(6): 638-641. 
    摘要 ( )   PDF(177KB) ( )  
    采用多重网格法(MG)和对称性理论以及密度泛函理论中著名的Kohn-Sham 方程,开发了一种计算具有一定对称性的原子簇的电子结构的方法,并对H2 分子的基态能进行了实际的计算,并与用其他理论计算得到的结果进行了比较。
    Si11+离子能级和振子强度
    闫公敬, 周忠源, 潘守甫
    1999, 16(6): 642-645. 
    摘要 ( )   PDF(162KB) ( )  
    用组态相互作用程序CIV3 计算了Si11+离子能级和振子强度。在计算中,考虑了一个电子保持在1s 轨道上,其余两个电子各种布居组态的组态相互作用。这种方法大大地改善了振子强度的计算结果。
    双轴分子构成向列型液晶中向列序的计算
    张志东, 张德贤, 黄锡珉
    1999, 16(6): 646-650. 
    摘要 ( )   PDF(199KB) ( )  
    通过格胞理论研究双轴分子构成的向列型液晶。发展新的数值方法计算液晶中的向列序。不用任何函数展开,通过迭代方法求解平衡态方程,得到精确的取向分布函数。计算了序参数和单分子内能对温度的依赖关系。给出了有关物理量在向列相-各向同性相相变点的值。数值结果较分子场理论更接近于Monte Carlo 模拟的结果。
    改进的PV-FD方法与沿井孔传播声波通过层的精细描述
    马俊, 董和风, 王克协
    1999, 16(6): 651-655. 
    摘要 ( )   PDF(288KB) ( )  
    采用改进的声压-速度有限差分(PV-FD) 数值模拟方法,着重研究水平薄层及裂缝对声波测井中全波列的影响。数值模拟了井外出现多种非对称水平薄层情况下的阵列声波记录,并模拟计算和绘制了移动测井仪器记录的实测声波变密度图。通过(斯通利波限幅) 放大,各类波通过界面和层的传播都得到精细的描述。
    圆柱空腔内旋转流动中轴对称涡破裂现象的数值模拟
    刘应征, 陈汉平
    1999, 16(6): 656-660. 
    摘要 ( )   PDF(240KB) ( )  
    对于上端盖旋转圆柱空腔内旋转流动来说,当流动达到一定条件时,就会出现旋涡破裂现象。目前针对这种现象已经进行了大量的稳态和非稳态实验测量,获得了精确的实验数据。采用数值模拟方法,重现了涡破裂现象及其发生发展过程,同已经发表的相关实验数据进行了比较
    非线性Schrödinger方程的守恒数值格式
    张鲁明, 常谦顺
    1999, 16(6): 661-668. 
    摘要 ( )   PDF(288KB) ( )  
    对非线性Schrödinger 方程提出了一种新的守恒差分格式,并证明了该格式的收敛性与稳定性。通过数值计算,对非线性Schrödinger 方程中非线性项的离散进行了讨论,获得如下结论:在取适当的参数后,所提出的差分格式精度上好于作为该格式特例的文[7] 中的格式。
    二维非结构网格生成及Euler方程计算的方法研究
    叶正寅, 杨永年, 钟诚文, 张仲寅, 陈迎春
    1999, 16(6): 669-674. 
    摘要 ( )   PDF(301KB) ( )  
    用离翼型表面最小距离作为阵面推进法中的参数选择依据,生成二维问题的非结构网格。这种方法抛弃了传统的背景网格观念,直接提供网格生成过程中所需的背景信息。在求解Euler 方程时,用格心格式的有限体积法作空间离散,用四步Runge-Kutta 方法作时间推进,采用不同的加速收敛措施获得定常流动。提出了两种边界条件的构造办法,并讨论了不同边界条件对结果的影响。
    基于车身外形CAD数据的表面气动网格生成
    姜乐华, 谷正气, 张丕付, 黄天泽
    1999, 16(6): 675-681. 
    摘要 ( )   PDF(385KB) ( )  
    提出一种在非均匀参数空间上生成表面数值网格的方法,开发了相应的表面网格生成软件SGRID,并应用于汽车空气动力学计算机仿真中具有复杂几何外形的车身表面网格生成。
    具有固支边的强厚度层合板的一种新解法
    盛宏玉, 范家让
    1999, 16(6): 682-687. 
    摘要 ( )   PDF(220KB) ( )  
    在文[1~3]的基础上,通过引入边界位移函数,对具有固支边的强厚度层合板建立了在任意荷载作用下的状态方程,给出静力问题的解析解。算例表明,该法具有收敛快,边界应力连续等特点。
分享到: