过刊目录

2003年, 第20卷, 第3期
刊出日期：2003-05-25 上一期    下一期
---
论文   

全全选选:

导出引用 EndNote Ris BibTeX

隐藏/显示图片

论文

Select

激光与靶相互作用体汽化模型的数值计算方法 陶应学, 关吉利, 陈发良 2003, 20(3): 189-192.  摘要 ( )   PDF(206KB) ( )   利用克劳修斯-克拉珀龙(Clausius Clapeyron)方程描述相变过程,利用von Neumann差分方法构造辐射流体动力学方程组的计算方法,数值求解了激光辐照金属材料的热效应、汽化过程和激光维持爆震波,并同作者前期的工作作了比较.结果说明,该计算方法对激光与靶材料相互作用体汽化模型的数值计算是行之有效的,且该汽化模型为相关问题的二维、三维问题的数值计算提供了简便有效的途径.

Select

三维非均匀不稳定渗流方程的自适应网格粗化算法 杨权一, 刘福平, 杨长春, 刘立峰, 卓磊 2003, 20(3): 193-198.  摘要 ( )   PDF(378KB) ( )   将渗透率自适应网格技术应用于三维非均匀不稳定渗流方程的网格粗化算法中,在渗透率或孔隙度变化异常区域自动采用精细网格,用直接解法求解渗透率或孔隙度变化异常区域的压强分布,在其它区域采用不均匀网格粗化的方法计算,即在流体流速大的区域采用精细网格.用该方法计算了三维非均匀不稳定渗流场的压降解,结果表明三维非均匀不稳定渗流方程的三维非均匀自适应网格粗化算法的解在渗透率或孔隙度异常区的压强分布规律与采用精细网格的解非常逼近,在其它区域压强分布规律与粗化算法的解非常逼近,计算速度比采用精细网格提高100多倍.

Select

耐腐蚀聚合药物基体上药物的扩散释放 S. Wasuwanich, N. Jinuntuya, P. Petpirom 2003, 20(3): 199-204.  摘要 ( )   PDF(371KB) ( )   主要对耐腐蚀聚合药物基体上药物的整个释放过程进行解析和数值研究.数学模型为基于Fick第二法则的边值问题,分为状态Ⅰ和Ⅱ.状态Ⅰ为具有移动性扩散波头的Stefan问题,而当此扩散波头消失时,即开始状态Ⅱ.该研究工作可作为数值研究更加复杂系统的基础之一.

Select

解三维扩散方程的积分内插法格式 陈光南, 李德元, 万正苏 2003, 20(3): 205-209.  摘要 ( )   PDF(248KB) ( )   研究三维扩散方程的数值模拟.在非正规六面体网格上,使用积分内插法建立扩散方程差分格式,涉及到27个相邻网格,适用于大变形网格上带间断系数的拟线性扩散方程的计算.叙述差分格式的建立,推导通量流和网格顶点温度的计算公式,给出了数值试验结果.

Select

向列相液晶沿面到垂面排列转变的分子场理论 向兰, 张志东 2003, 20(3): 210-214.  摘要 ( )   PDF(328KB) ( )   以摩擦基板间的液晶薄层为研究对象,用分子场理论研究了向列相液晶分子排列转变行为.分子质心固定在简单立方晶格的格点上.液晶由极性分子构成,与基板相接触的一层分子同时受到色散和极性两类表面作用.通过自洽的数值计算,获得3种相图,清楚地展示了摩擦基板间向列相出现的从高温沿面到低温垂面排列的转变;获得实现这类转变所需要的两类表面作用的参数范围.结果表明:基板的摩擦会改变基板表面色散作用,但不会影响基板表面极性作用;表面极性相互作用能引起基板间向列相液晶发生沿面到垂面排列转变.

Select

液晶胆甾相中的分子短程关联 张志东, 卢遵铭 2003, 20(3): 215-218.  摘要 ( )   PDF(208KB) ( )   向列相液晶的二粒子集团理论被推广应用于研究胆甾相二维模型.手征性分子固定在三维简单立方晶格的格点上,而分子取向限制在二维.理论结果表明,平衡态螺旋波矢依赖于温度的变化,且存在胆甾相到向列相相变.通过考虑分子间短程关联,二粒子集团理论的数值结果较平均场理论更接近Monte Carlo模拟结果.

Select

利用Ghost Fluid方法模拟激波与柱形界面相互作用 程军波, 唐维军, 李德元, 傅德薰 2003, 20(3): 219-225.  摘要 ( )   PDF(510KB) ( )   利用GhostFluid方法(后面简称Ghost方法)和γ-model方法,在同样的时空离散精度条件下,对激波与柱形界面相互作用的二维可压缩流场进行了直接模拟,并与实验结果相比较.从模拟结果看,在短时间内,Ghost方法和γ-model方法模拟的结果与实验结果基本相同,两种方法均正确地模拟出界面的位置、激波的强度和速度.但随着时间的发展,具有较大数值耗散的γ-model方法的计算结果与实验差别越来越大;而数值耗散较小的Ghost方法能较为正确地模拟界面的运动.

Select

开槽扰流片流动特性研究 张攀峰, 王晋军 2003, 20(3): 226-232.  摘要 ( )   PDF(421KB) ( )   用标准的k-ε湍流方程结合近壁流动的两层带状模型,在非结构网格上采用SIMPLE算法,对不同开孔度的开槽扰流片流动特性进行了数值模拟,并和测压实验结果进行了比较.结果表明开槽扰流片的流动结构取决于上下两部分流场的相互耦合作用,开孔度η=50%时两个流场互相没有干扰,开孔度η=20%,30%时两个流场相互耦合作用很强,扰流片后有较大的回流区以及较小的湍流动能.

Select

升力系数不变时跨音速机翼阻力优化设计 杨旭东, 乔志德, 朱兵 2003, 20(3): 233-238.  摘要 ( )   PDF(310KB) ( )   基于共轭方程的优化设计理论,应用三维欧拉方程进行了升力系数不变时跨音速机翼阻力优化设计研究,根据给定的目标函数推导了在物理空间上表述的共轭方程及边界条件,研究了共轭方程的数值求解方法及目标函数对设计变量的敏感性导数求解问题,发展了一种跨音速机翼阻力优化设计方法,应用该设计方法进行了跨音速机翼阻力优化设计研究,优化后机翼表面的激波强度减弱很多,有效减少了波阻.

Select

考虑扩散和吸附作用的聚合物驱替过程渗流数值模拟 侯健, 李振泉, 王玉斗, 陈月明 2003, 20(3): 239-244.  摘要 ( )   PDF(356KB) ( )   采用显式求解饱和度、隐式求解浓度的思路,对考虑扩散和吸附作用的聚合物驱替过程渗流模型进行了数值求解.饱和度方程求解应用了显式全变差递减(TVD)法;浓度方程求解过程中,空间项离散采用Crank Nicolson差分格式,时间项变量进行拟线性处理,保证了计算的稳定性.通过与解析解对比,验证了该方法的有效性.计算实例分析表明,扩散使聚合物在溶液中稀释,导致浓度传播分散;吸附使聚合物损耗,导致浓度传播滞后.同时,计算结果直观反映了聚合物驱重要的"油墙"形成机理.在段塞注入情况下,原油富集区在出口端的突破是介于聚合物浓度前缘突破和聚合物浓度峰值突破之间.

Select

二维小波变换及在青藏高原地形重构中的应用 戴新刚, 王国军, 汪萍 2003, 20(3): 245-254.  摘要 ( )   PDF(845KB) ( )   计算和检验了4种不同小波基底和不同阶小波基底对单脉冲和尖峰两个信号的分解及重构的性质,结果表明双正交小波的精度最高,Coifman小波的精度最低.通过设阈值对小波和快速Fourier变换进行截断后再重构,说明小波具有重构精度高、所用非零系数少以及误差被局域化等优点,这就削弱了"Gibbs现象",限制了它的影响范围.在此基础上,用二维快速Fourier和二维离散小波变换对实际青藏高原及其附近地形做了分解、重构及截断后重构等,说明二维小波变换能以较之FFT更少的非零系数和误差局域化等重构大地形,表明了二维小波变换在未来大气环流数值模式中的应用潜力.

Select

碳纳米管储氢性能的计算机模拟 程锦荣, 闫红, 陈宇, 张立波, 赵力, 黄德财, 唐瑞华 2003, 20(3): 255-258.  摘要 ( )   PDF(191KB) ( )   采用巨正则蒙特卡罗(GCMC Grand Canonical Monte Carlo)方法[11],系统地研究了锯齿(Zigzag)型单壁碳纳米管(SWNT-Single Walled Carbon Nanotubes)的储氢性能,得到了管径和管长与储氢能力的关系曲线,给出了氢在碳纳米管中的分布,并对碳纳米管储氢能力的表征提出了建设性意见.

Select

一层随机非球形粒子高阶散射矢量辐射传输迭代解的参数求逆 梁子长, 金亚秋 2003, 20(3): 259-263.  摘要 ( )   PDF(322KB) ( )   根据一层随机非球形粒子矢量辐射传输(VRT)方程的一阶与高阶Mueller矩阵解的差异,提出了在已知一些参数的条件下,用水平与垂直同极化和交叉极化后向散射测量反演非球形粒子复介电常数和单位面积粒子数的迭代方法,并用极化后向散射强度的模拟测量数据及机载合成孔径雷达(Air SAR)实测数据进行了反演试验,并对反演值与各参数的关系以及测量噪声对反演的影响进行了讨论.

Select

用遗传算法拟合Eu-151穆斯堡尔谱 孙雪峰, 林俊, 黄红波, Ali, 何云, 戴耀东, 夏元复 2003, 20(3): 264-266.  摘要 ( )   PDF(195KB) ( )   以Eu2O3(立方晶系)为例,论述遗传算法(GA,Genetic Algorithm)求解穆斯堡尔谱的过程,并将解谱结果与文献数据和Moss Win解谱结果进行比较.

Select

非线性薛定谔方程数值解法中傅立叶正逆变换选取的讨论 张晓光, 杨伯君, 俞重远 2003, 20(3): 267-272.  摘要 ( )   PDF(515KB) ( )   讨论了光波场正负频表示与正逆傅立叶变换形式的选取及非线性薛定谔方程的形式密切相关,指出了现在熟知的非线性薛定谔方程形式是选取负频表示的结果,因此正逆傅立叶变换形式的选取就不是任意的.在采用基于傅立叶变换数值方法解非线性薛定谔方程时,要注意所用编程语言中快速傅立叶变换正逆变换的形式.而在大多数编程语言中其快速傅立叶变换正逆变换形式是与非线性光学中应采用的傅立叶变换正逆形式恰好相反,因此时而出现一些错误的计算结果.分析了这些错误的原因,并给出了正确结果.

Select

非结构网格上的B-B单方程模型湍流计算 成娟, 黄明恪 2003, 20(3): 273-278.  摘要 ( )   PDF(403KB) ( )   研究如何在非结构网格上进行Navier Stokes(N-S)方程湍流计算.采用格心有限体积方法离散N-S方程.为了适应非结构网格,计算所用的湍流模型特别选用Baldwin Barth(B-B)单方程模型.此模型由一个单一的具有源项的对流扩散方程组成.为了能在非结构网格上求解B B单方程模型,提出一显式有限体积格式,并直接对带源项的格式进行稳定性分析,得到了相应的时间步长限制条件.最后以平板、RAE 2822翼型、多段翼型绕流等数值算例验证了计算方法的有效性.

Select

MEND程序中能谱的计算方法 蔡崇海 2003, 20(3): 279-282.  摘要 ( )   PDF(264KB) ( )   给出了在中能区(≤250MeV)多次粒子发射过程中计算各种出射轻粒子及各种反冲核之能谱的方法,指出它的计算效率远高于(两个数量级以上)传统方法,并给出了两个算例.