初始压入位置对Ni基单晶合金纳米压痕影响研究

针对Ni基单晶合金建立初始压入γ 相的γ /γ' 模型和初始压入γ'相的γ'/γ 模型, 采用分子动力学方法模拟金刚石压头压入两种模型的纳米压痕过程, 计算两种模型[001]晶向硬度. 采用中心对称参数分析两种模型(001)相界面错配位错对纳米压痕过程的影响. 结果显示: 弛豫后, 两种模型(001)相界面错配位错形式不同, 其中γ'/γ 模型(001)相界面错配位错以面角位错形式存在; 压入深度在0.930 nm 之前, 两种模型(001)相界面错配位错变化不大, 压入载荷-压入深度及硬度-压入深度曲线较符合; 压入深度在0.930 nm之后, γ'/γ 模型(001)相界面错配位错长大很多, 导致相同压入深度时γ'/γ 模型比γ /γ'模型压入载荷和硬度计算结果小; 压入深度在2.055 nm之后, γ /γ'模型(001)相界面错配位错对γ 相中位错进入γ'相有阻碍作用, 但仍有部分位错越过(001) 相界面进入γ' 相中, γ'/γ 模型(001)相界面处面角位错对γ' 相中位错进入γ 相有更明显的阻碍作用, 几乎无位错越过(001) 相界面进入γ 相中, 面角位错的强化作用更明显, 所以γ'/γ 模型比γ /γ'模型压入载荷上升速度快.     

界面追踪方法中的激波限制器研究

针对多介质流体界面追踪(Front Tracking)方法,通过在界面处构造Riemann问题,研究激波限制器应用过程中的若干基本问题.针对气-气和气-水界面问题,通过比较平均守恒误差、L₁误差和激波强度随参数的变化情况,给出激波限制器参数可以在0.3附近选取.同时,通过理论分析和数值算例发现,当有激波接近界面时,选择激波波后状态作为界面处Riemann问题的初始状态,数值模拟结果较满意.     

基于压缩感知理论的非线性γ谱分析方法

γ谱分析是一种重要的放射性核素定量分析方法. 弱峰的检测和重叠峰的分解是γ 谱分析中的难点. 为了解决这两个问题, 基于压缩感知理论, 提出了一种新的γ 谱分析方法. 这一方法从谱仪对γ 谱调制的物理机理出发, 通过数学建模, 将γ 谱分析转化为一个以真实γ 谱为解的求逆问题, 并在压缩感知理论框架下, 运用γ 谱特征峰的稀疏性, 进行逆问题的求解, 直接获得γ 谱的估计结果. 数值模拟结果和蒙特卡洛模拟结果表明: 该方法能在降低统计涨落的同时, 有效减小谱仪调制带来的能谱展宽, 从而提高γ 谱分析精度. 关键词： γ谱分析')" href="#">γ谱分析 压缩感知 非线性 逆问题     

激波与SF6球形气泡相互作用的数值研究

本文基于大涡模拟方法, 采用高阶精度格式对平面入射激波以及不同反射距离条件下的反射激波与SF₆重气泡相互作用过程进行了三维数值模拟. 数值结果清晰地显示了SF₆重气泡在激波作用下诱导Richtmyer-Meshkov不稳定性过程, 揭示了入射激波以及反射激波在气泡界面聚焦诱导射流的过程, 详细分析了不同反射距离条件下反射激波与SF₆重气泡作用过程及流场结构.     

可压缩Kelvin-Helmholtz不稳定性低耗散锐界面方法直接模拟

采用低耗散WENO(weighted essential non-oscillatory)格式及锐界面方法模拟可压缩Kelvin-Helmholtz不稳定性问题.由于物质界面被描述成一种接触间断, 该方法可精确求解切向速度间断.基于优化模板对原始光滑指标进行正规化后, 得到一种低耗散WENO格式.修正后的方法显著降低了普通流动区域的过衰减问题, 保持了良好的激波捕捉性能, 并可获得与混合格式相当的求解精度.不同于以往求解单一流体或易混界面时, 通过初始设定有限宽度的剪切层或快速数值耗散以抑制高波数模态, 该方法允许高波数扰动的发展.计算结果表明, 高波数扰动展现出与以往理想Kelvin-Helmholtz不稳定性问题数值模拟或线化理论结果不同的特征, 但与有限厚度的剪切层结果相符.      

一类新型激波捕捉格式的耗散性与稳定性分析

高超声速流动是高复杂性的可压缩黏性流动, 其中存在激波、剪切层、激波/激波干扰、激波/边界层干扰、旋涡与分离流动等复杂流场结构. 对其进行准确模拟需要使用低耗散、强鲁棒性的激波捕捉方法. 本文基于一类新型的通量项分裂方法, 提出了一种耗散低且鲁棒性好的激波捕捉格式K-CUSP-X. 对该格式的耗散性和激波稳定性进行了详细的理论分析, 得到了格式激波稳定的数值条件. 推论认为, 迎风格式激波稳定的充分条件为速度扰动量具有衰减性, 数值实验验证了该推论. 研究表明, 该格式与Toro提出的通量分裂格式K-CUSP-T相比, 在保证精确捕捉接触间断的同时, 又具有更好的稳定性, 在激波处不会产生“红玉”现象.     

运动激波和气泡串相互作用的初步数值模拟

通过对激波和流体界面相互作用诱导的大变形界面演化的数值模拟,验证Level set方法精确模拟多个流体界面的有效性.采用2阶迎风TVD求解欧拉方程得到流场解,采用5阶WENO求解Level set方程追踪多流体界面,采用GFM方法处理流体内界面.利用文[1]的计算结果校核本文程序.在此基础上,对运动激波和气泡串相互作用过程进行了初步数值模拟,得到了不同时刻运动激波和圆管内的两个气泡作用后的演化图象,包括压力和密度等值线分布.计算结果表明:针对推广后的多界面Level set方程,该方法仍可高质量地捕捉多个流体界面.     

多孔介质中非等温互溶流体驱替过程的格子Boltzmann研究

利用格子Boltzmann方法,在孔隙尺度上对多孔介质内非等温混溶驱替过程进行数值研究,定量分析温度粘性膨胀系数（β_T）以及路易斯数（Le）对界面不稳定性和驱替效率的影响.结果表明：随着β_T的增大,界面不稳定性增强,驱替效率降低.当β_T>0时,随着Le的增大,界面不稳定性减弱,驱替流体与被驱替流体之间的界面趋于平缓,指尖残余率减小,驱替效率增大.当β_T<0时,Le对于驱替效率的影响相反.     

随机扰动下三维流体界面不稳定性的并行计算

对三维流体界面不稳定性的数值模拟引进了新的数值计算方法,并在MPI并行计算环境下进行了数值模拟.利用LevelSet方法确定界面位置,零水平集对应界面位置.对应离散LevelSet方程和界面两侧的两套Euler方程,借助于Ghost网格方法来完成离散.对最后网格点上的两套状态量的辨认依赖于该点的LevelSet值的符号.并进行了数值计算.     

基于LevelSet方法捕捉界面间断的Ghost技术

应用Ghost技术的LevelSet方法捕捉运动界面,采用五阶精度的WENO格式,时间离散采用TVD Runge-Kutta格式,完成了水池中油滴的上升与合并算例及火焰面燃烧算例的数值模拟,并与商用软件Fluent的计算结果进行了比较,肯定了采用Ghost技术的LevelSet方法捕捉产生间断的运动界面的合理性.同时,本文程序的编制过程体现了LevelSet方法无需进行复杂的界面重构,易于编程的优越性.     

脉冲源法测量次临界系统裂变衰减特性的理论分析和数值模拟

针对点堆动力学理论解释脉冲源法测试原理时存在的问题, 基于无源中子输运方程分析次临界系统总中子数、泄漏γ射线计数率随时间的变化关系。理论分析表明: 脉冲中子源作用结束后(无源条件下), 在一定时间范围内, 泄漏γ射线计数率和总中子数近似成正比, 两者随时间变化服从近似指数衰减规律, 反映系统本身的裂变衰减特性, 可以由总中子数和γ射线计数率求解瞬发中子衰减时间常数。基于蒙特卡罗程序构造类Godiva裸铀球次临界系统, 模拟脉冲中子源作用下中子和γ射线输运过程, 计算总中子数、泄漏γ射线计数率及两者比值随时间的变化关系, 结果与理论分析一致; 利用脉冲源法由总中子数、泄漏γ射线计数率计算瞬发中子衰减时间常数α₀, 得到与α-k迭代一致的α₀。说明总中子数、泄漏γ射线计数率可以准确反映系统本身的裂变衰减特性。此外, 根据理论分析和模拟计算给出脉冲源法可用数据的时间范围, 分析泄漏γ射线计数率和总中子数比值的影响因素。     

色散介质电磁特性时域有限差分分析的Newmark方法

根据Debye模型、Drude模型和Lorentz模型3种常见色散介质模型频域极化率的特点,利用频域到时域的转换关系jω→?/?t,将极化矢量P与电场强度E的频域关系转换成时域内关于P的二阶微分方程,其对3种色散介质模型皆适用,具有统一的形式.然后采用相比于中心差分具有更高精度的Newmark两步算法(Newmark-β-γ法)求解该方程,进而得到E→P的递推公式,再结合本构关系得到D→E的时域递推式.实现了色散介质电磁场量的时域有限差分迭代计算.数值计算结果表明该方法是适用于3种色散介质模型的通用算法,并且相比于移位算子时域有限差分方法等以中心差分为基础的离散方案具有更高的计算精度.     

少周期脉冲在稠密V型三能级介质中传播时的粒子布居演化

利用由预估计校正-时域有限差分(PC-FDTD)法求解Maxwell-Bloch方程得到的数值解,研究少周期激光脉冲在稠密V型三能级原子介质中传播时各能级粒子数布居的演化特性.结果表明,两跃迁偶极矩的比值(γ)对稠密介质中粒子数布居反转出现的时间和振荡的次数具有显著的影响,γ=1时,考虑洛仑兹局域场修正(LFC)时,粒子数反转出现的较晚,近偶极-偶极(NDD)相互作用延缓了粒子数的反转;γ>1时,考虑LFC时,粒子数反转出现的较早,NDD相互作用加快了粒子数的反转;电偶极矩不同时,粒子数布居振荡的次数比相同时多.粒子数布居在稠密介质中的演化特性与稀疏介质中显著不同,在入界处,稀疏介质中粒子数布居可以实现长时间的反转;随着介质密度的增大,粒子数布居产生准周期性的振荡,但考虑LFC时,振荡的次数比不考虑LFC时少.     

γ脉冲射线束在大气中的输运计算方法

在γ脉冲射线束测量中,束流波形是由许多不同能量的γ粒子叠加而成的,在大气介质传输中,由于γ粒子与空气原子的相互作用,测点的γ能通量强度随时间变化的波形、同源点的γ强度波形相比,在较大测点距离的情况下就会产生畸变。讨论了γ射线在大气中输运的传输函数以及用褶积计算测点γ能通量强度波形的方法,给出了解析表达式,对典型的γ源脉冲波形输运不同距离,计算了大气输运对测量的影响。最后讨论了减小或避免波形畸变的方法。     

金属Fe与间隙H原子相互作用的密度泛函研究

采用基于密度泛函理论的第一性原理方法, 研究了不同摩尔比下H在α-Fe和γ-Fe晶格中的间隙占位情况, 计算了稳态晶体的总能量、结合能、溶解热、电子态密度、电荷差分密度和电荷布居, 分析了间隙H原子和Fe晶格之间的相互作用, 讨论了H溶解对α-Fe和γ -Fe晶体电子结构的影响. 结果表明: H溶解引起α-Fe和γ-Fe晶体点阵晶格畸变, 体积膨胀率随着溶氢量的增加而增加. 从能量角度分析发现, H优先占据α-Fe的四面体间隙位, 而在γ -Fe中优先 占据八面体间隙位. 态密度、电荷差分密度以及电荷布居分析发现, 间隙H原子与Fe晶格的相互作用仅由H的1s轨道电子和Fe的4s轨道电子所贡献, 二者作用力相对较弱, 这是造成H在Fe晶格中固溶度较低的主要原因之一. 关键词： 金属Fe 间隙H原子 第一性原理 溶解热     

特征γ谱蒙特卡罗模拟方法的改进

针对中子-γ耦合输运问题,采用蒙特卡罗统计估计和期望值技巧,分别计算了被探测物的总能谱、时间谱和每个元素的特征γ谱,并与MCNP程序结果进行了对比.     

EV共聚物序列中的两种结构因素

系统描述了EV共聚物序列中两种结构因素对共聚物主链中亚甲基和次甲基¹³C NMR化学位移的贡献,将SCS方法与半经验方法结合起来,阐明了常数项δ_LE和δ₀的物理意义,提出了以聚乙烯链反式构象化学位移δ_t为基准计算上述亚甲基和次甲基化学位移的公式:δ_c=δ_t+Σn_iS_i+ΣP_jγ_j并在此基础上,讨论了聚乙烯,聚丙烯的特殊情况和乙丙共聚物的链结构.     

欧拉圆盘不同能量耗散机理之间的关联

本文研究了欧拉圆盘运动过程中盘厚度以及盘面与水平面夹角α两因素对能量耗散的影响. 得出圆盘厚度与直径之比x对能量变化中各项因子的影响: x很小时, 质心在竖直方向上的动能变化和重力势能变化是系统能量耗散的主要因素; 当x>0.4142时, 圆盘绕与之平行的轴的转动动能变化成为主要因素, 并给出圆盘厚度可忽略的条件. 模拟了滚动摩擦、空气黏滞等不同能量耗散方式与x,α的关系, 导出各种耗散方式在圆盘运动的过程中的转变规律, 并指出x=0.1733, α>18°时能量耗散形式为纯滚动摩擦, 这修正了文献[26]结论.