圆柱形汇聚激波诱导Richtmyer-Meshkov不稳定的SPH模拟

汇聚激波诱导不同物质界面的Richtmyer-Meshkov(RM)不稳定现象在惯性约束核聚变领域有重要的学术意义和工程背景.基于网格离散的宏观流体力学方法由于数值扩散问题往往需要高阶精度算法才能准确追踪界面演化,且对大变形和破碎合并等复杂界面追踪也极为困难.光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法采用纯拉格朗日算法,可以有效克服上述难点.但经典SPH算法需采用人工黏性处理强间断,在激波间断处往往会出现严重的非物理振荡,对于涉及强冲击不稳定性问题,很难达到理想的模拟效果.本文采用基于HLL黎曼求解器的SPH算法,实现了对强激波和大密度比物质界面的有效分辨和追踪.一维数值校核证明了代码的可靠性、健壮性,并进一步模拟了二维圆柱形汇聚冲击波冲击四边形轻/重气界面诱导的RM不稳定性问题,与已有实验结果进行了对比,发现模拟结果与实验结果吻合.通过分析界面演化过程中的密度及压力变化,发现本文所采用的方法可准确地追踪激波与界面作用的复杂界面和波系演化规律.研究结果为进一步理解和解释汇聚冲击条件下的RM不稳定性机理奠定了基础.     

旋转射弹高速倾斜入水多相流场与弹道数值模拟

基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程，结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响，然后对不同入水角下倾斜入水过程进行分析，得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性，但会降低弹体侧向稳定性，使射弹受到的阻力系数、俯仰力矩系数变小；入水角越小，形成的空泡越不对称，由射弹运动状态的改变引起的空泡形态变化越明显，在超空泡航行阶段，弹体运动较稳定，不同角度下流体动力系数差别很小，当弹体下表面刺破空泡壁沾湿时，弹体运动状态发生较大变化，流体动力系数迅速增大，此时入水角度过小，弹体容易失稳；弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。     

食品口腔摩擦学研究进展

食品在口腔加工的过程被认为涉及流变学和摩擦学机制，而口腔摩擦则被认为是口腔质构感知产生的一个重要机理，用来描述舌头与口腔上颚、舌头与食品之间相对运动以及由此产生的触觉感受. 本文中综述了近年来食品口腔加工领域的摩擦学研究进展，包括摩擦学测量技术与方法、体内口腔摩擦和感官知觉之间的联系等. 相关研究表明，摩擦学测量为几种与质构相关的感知提供了理论依据，但获得两者之间的量化经验关系仍然很有挑战性；同时，由于摩擦学和感官感知受到多种物理化学性质的影响，因此必须将摩擦学和其他表征技术相结合用以解释两者的相互关系及机制，从而促进食品口腔加工摩擦学的持续和深入发展，并为高品质食物和饮料的研究和设计提供重要工具和理论支撑.      

基底纹理对镍/铜纳米双层膜刮擦行为影响的分子动力学模拟

采用分子动力学模拟研究纳米尺度下的镍/铜双层膜的表面刮擦响应. 考虑基底铜膜材料的不同表面纹理，建立分子动力学模拟模型，分析了纳米双层膜不同基底纹理下的刮擦响应，并且对比了两种典型基底纹理下不同刮头半径、不同刮擦深度对纳米双层膜刮擦响应的影响. 研究结果发现:不同的基底纹理下，位错缺陷程度的不同会导致刮头前方的切屑体积不同，存在切屑体积最大的基底纹理；针对两种典型基底纹理，在刮擦深度或刮头半径一定时，刮擦力随刮头半径或刮擦深度的增加而增大；当刮擦深度或刮头半径超过临界值时，表面具有特定齿槽纹理的基底具有一定的减摩作用.      

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明，模型方程在小到中等Reynolds数下都适用，并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加，展向不稳定性增长率则随着Marangoni数增加而增加，随着Reynolds数增加而减小，流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明，在等热流密度加热条件下，液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构，液膜与气体的换热也因溪流状结构的出现而加强；在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强，液膜会提前破裂；在随机扰动的基础上引入占优势地位的流向最不稳定扰动时，液膜的换热会增强，但不会提前破裂；在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时，换热会加强且液膜会提前破裂。     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

超前注水开发指标变化规律及跟踪调整研究

X区块为低渗透油藏,目前采用超前注水开发方式,在数值模拟历史拟合的基础上对试验区超前注水开发指标的变化规律进行研究,并对注采系统进行调整.通过相对渗透率曲线,得到无因次采液指数与含水率之间的变化规律,结果表明,含水率低于65%时,无因次采液指数呈下降趋势,可适当降低注水强度;当含水率达到65%以后,无因次采液指数开始上升,此时可通过放大生产压差对试验区进行提液生产;而含水率达90%以后,无因次采液指数迅速上升,此时应以稳油控水为主.在确定产液量的基础上,在下降期和上升期分别设计了四个方案进行注水强度调整,结果表明,在无因次采液指数下降期和上升期,当注水强度分别为1.65 m~3/d·m和2.10m~3/d·m时,地层压力保持较好,且试验区开发效果最好.     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.