高超声速模型尾迹电子密度二维分布反演方法

弹道靶利用二级轻气炮将模型加速到高超声速状态,模型在靶室内超高速飞行时形成等离子体尾迹.为实现高超声速模型尾迹电子密度径向二维分布诊断,利用七通道微波干涉仪测量系统获得了高超声速模型尾迹截面不同位置处平均电子密度.该系统采用一发七收的方式,实现平面波照射等离子体及平面波接收,天线波束可完全覆盖尾迹径向范围.多通道微波干涉仪数据处理过程常将等离子体视为分层介质,考虑到分层界面上折射效应的影响,本研究利用射线追踪方式建立电磁传播模型,结合测量数据建立目标函数,通过遗传算法优化来反演高超声速模型尾迹电子密度径向二维分布.该数据处理方法的电子密度反演结果与相同来流条件下的数值模拟结果对比吻合较好,初步验证了该方法的有效性.分析了分层模型对电子密度分布特性的影响,结果表明利用七层模型对尾迹建模效果最佳,且适用于不同厚度尾迹,最大化利用接收通道数,确保了计算精度.利用该方法实现弹道靶高超声速球模型尾迹电子密度二维分布诊断,并给出了给定实验状态下模型尾迹电子密度二维分布的一些规律.     

急性髓细胞白血病基因筛选模型的贝叶斯分析

基因表达数据蕴含着大量的生物信息,在生物基因信息研究中,筛选表达水平发生显著变化的差异基因是认识疾病形成机理和辅助靶点药物研究的关键问题.根据急性髓细胞白血病(AML)的基因表达数据,构造基因均值差序列,建立贝叶斯分层混合模型,并为模型的参数赋予具有基因生物特征的先验信息.采用马尔可夫链蒙特卡洛(MCMC)算法对模型参数进行估计,并筛选出急性髓细胞白血病差异表达基因.在实际数据分析中,从美国生物信息中心(NCBI)的高通量基因表达数据库中获取急性髓细胞白血病基因数据集,从经过非特异滤波预处理的14688个急性髓细胞白血病基因中筛选出711个差异表达基因,差异表达基因数仅占急性髓细胞白血病基因总数的4.84%,这一结果与基因差异表达的生物学原理相吻合.     

基于模型的小域估计方法及其拓展——由单变量推广到多变量情形

小域估计是当今抽样调查领域的研究热点,通过建立小域估计模型,将所有调查区域通过公共的参数连接起来,利用相关区域数据作为辅助信息,解决多层次抽样调查中某些域样本量不足的问题,提高参数的估计精度。本文首先在目标变量是单变量的情形下,研究了域层次模型和单元层次模型的参数估计方法,进而推广到目标变量是多变量的情形,考虑变量之间的相关性,研究了多变量Fay-Herriot模型参数的经验最优线性无偏预测及其均方误差矩阵的性质,最后通过一个实例,验证了参数估计具有良好效果。     

光学诊断专题序言

光学诊断技术主要是以激光技术、光谱技术、光电探测技术、数据图像处理技术等为基础的一种综合性测试诊断技术可以实现复杂流场温度、组分浓度、速度、流场结构等参量信息的高时空分辨精确测量而且对测量流场无扰动.近年来光学诊断技术已开始从实验室基础研究走向于各类风洞、发动机等重要设备的工程应用这对于深入研究高超声速流动和燃烧化学反应动力学过程如高焓非平衡流动、高超声速边界层转捩、超声速燃烧动力学、汽车/飞机/火箭/卫星等发动机燃烧不稳定性和污染控制等具有重要意义.     

核电结构土-结相互作用分析分区混合计算方法

土-结构相互作用分析是核电结构进行抗震设计和安全评估的重要环节.在核电结构的土-结相互作用分析中,阻尼和非线性是影响结构反应的重要因素. 若采用频域分析,可以方便考虑阻尼,但需通过等效线性化来考虑非线性,不适合于强震作用下的土体非线性.若采用时域分析的逐步积分方法,适合于考虑非线性,但材料阻尼一般采用瑞利阻尼模型,除了紧靠指定阻尼比的少数几个振型外,其他振型的反应将受到瑞利阻尼模型所确定的大阻尼所抑制,造成地震反应与真实情形有较大差异.若采用时域分析的模态叠加法,可合理计入阻尼效应,但模态叠加法不能考虑非线性.因此,如何合理考虑阻尼和非线性是核电结构土-结相互作用分析需要关注的问题.基于此,本文提出一种模态叠加和时步积分结合的土-结相互作用分区算法.其中,出于安全性考虑,地震作用下核电主体结构一般不允许进入非线性,因此结构可采用模态叠加方法,以便合理考虑结构阻尼;土体和基础采用显式时步积分法,可考虑土体非线性;通过人工边界条件考虑无限域的影响 (辐射阻尼).通过简单算例对该方法进行了验证,并用于CAP1400核电结构的土-结相互作用分析中,对比分析了采用模态阻尼和瑞利阻尼时核电结构和场地反应的差异,结果表明结构阻尼模型对场地的反应影响不大,但对结构反应影响明显,在实际工程中应合理选取阻尼模型.      

CL-20/TNT共晶炸药的分子动力学研究

运用分子动力学(MD)方法,选择凝聚态分子势能优化力场(COMPASS),对六硝基六氮杂异伍兹烷(ε-CL-20)、2,4,6-三硝基甲苯(TNT)晶体及其等摩尔比的CL-20/TNT混合炸药和共晶炸药进行不同温度下恒定粒子数等压等温(NPT)系综模拟研究.结果表明,CL-20/TNT共晶的内聚能密度(CED)和结合能随温度的升高逐渐减小;共晶的CED比混合炸药的大,结合能是混合炸药的2倍多,预示其稳定性明显增强.对相关函数和局部放大结构显示共晶中组分分子间作用主要来自TNT中H和CL-20中O以及CL-20中H和TNT中O之间形成的氢键.通过波动法求得的弹性力学性能结果表明,CL-20/TNT共晶的拉伸模量(E)、体积模量(K)和剪切模量(G)介于ε-CL-20和TNT晶体之间,且随温度的升高而下降,符合一般预期;但共晶炸药的柯西压(C12-C44,Cij弹性系数)、K/G和泊松比(ν)均比其组分炸药ε-CL-20和TNT高得多,预示该共晶具有异常高的延展性和弹性伸长,主要是二组分呈层状交替排列且之间存在较强相互作用所致.     

基于车身绕流的低雷诺数湍流模型改进研究

本文将汽车绕流模块化为各典型局部流动,通过常用湍流模型对各典型局部流动进行数值模拟,结果验证了湍流模型对转捩的捕捉能力是准确模拟汽车绕流的关键. 在分析汽车绕流分离及转捩机理的基础上,优化了稳态和瞬态求解方法,改进了湍流模型对转捩的预测能力,进而提高了湍流模型在汽车流场模拟上的精度. 针对汽车绕流的稳态问题,将流线曲率因子及 响应阈值引入 LRN $k$-$\varepsilon $ (low Reynolds number $k$-$\varepsilon $) 模型,获得了一种能够更准确预 测转捩的改进低雷诺数湍流模型 (modified LRN $k$-$\varepsilon $),改善了原模型对湍流耗散率的过强依赖性及全应力发展预测不足等问题;针对汽车绕流瞬态求解,通过分析 RANS/LES 混合湍流模型的构造思想及特点,引入约束大涡模拟方法,结合本文提出的改进的 LRN $k$-$\varepsilon $ 湍流模型,提出了一种能准确捕捉转捩现象 的转捩 LRN CLES 模型. 分别将改进的模型用于某实车外流场和风振噪声仿真中,通过 Ansys Fluent 求解器计算,并将计算结果与常用湍流模型的仿真结果、HD-2 风洞试验结果和实车道路实验结果进行对比,表明改进后的湍流模型能够更准确模拟复杂实车的稳态和瞬态特性,为汽车气动特性的研究提供了可靠理论依据及有效数值解决方法.      

钝头体中的广义雷诺比拟关系

钝头体壁面的摩阻和热流分布规律不同,平板流动中的雷诺比拟关系在钝头体壁面失效. 文章在前期高超声速广义雷诺比拟理论研究工作的基础上,利用数值仿真的方法对不同外形和来流参数条件下的钝头体广义雷诺比拟关系开展进一步研究. 通过建立钝头体绕流边界层的理论分析模型,得到了钝头体壁面雷诺比拟系数的线性分布预示公式. 采用数值求解 N-S 方程的方法,计算了圆柱和幂次体壁面的摩阻和热流以及二者之间的比拟系数. 通过与前期数值和理论结果对比,以及计算收敛性和网格无关性检验,对数值方法进行了验证. 通过在不同雷诺数 ($Re_\infty = 3.98\times 10^2 \sim 1.59\times 10^6$) 和马赫数 ($M_\infty = 3\sim 12$) 条件下的计算结果对比分析雷诺比拟系数的分布,总结了钝头体中广义雷诺比拟关系受外形和来流条件的影响,评估了广义雷诺比拟理论的适用性. 研究发现,在较高雷诺数条件下,离驻点较远的下游 ($\theta > 60^\circ$) 部位,雷诺比拟系数的分布不同程度地偏离理论预示的线性规律. 相比于圆柱外形,幂次体壁面的雷诺比拟系数分布的线性规律相对较好,其分布斜率略低于圆柱壁面的结果. 研究表明,如果针对实际外形和雷诺数进行适当修正,可以提高广义雷诺比拟关系的预示精度.      

混合/流体润滑状态下原位离子液体添加剂的摩擦学性能研究

本工作重点研究了混合/流体润滑状态下原位离子液体添加剂的摩擦学性能,选用聚乙二醇(PEG-400)作基础油,将双(三氟甲基磺酰)亚胺锂盐(Li TFSI)溶解在PEG中原位合成离子液体.利用微型牵引力试验机测量在室温、60和80℃以及不同滑滚比下摩擦系数随卷吸速度的变化,研究离子液体添加剂的有效性以及离子液体添加剂对PEG流变行为的影响.本研究中将为深入研究离子液体的润滑机理提供一种新的研究手段,对于指导设计新型离子液体润滑材料具有较为重要的意义.