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

基于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时,地层压力保持较好,且试验区开发效果最好.     

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

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

叉流式新风换热器数值模拟

换热器是新风系统的核心部件,以某品牌新风换热器为模型,对该换热器进行数值模型,分析换热过程中换热和流动特性,以及流速、流道高度对换热性能的影响。建立不同流道高度时新风侧温差、压差、努塞尔数、阻力系数、综合传热性能因子随速度的变化关系。发现在换热器内存在"换热死区",随速度增加"换热死区"范围增大,在流速0.5m/s～1.6m/s、流道高度3mm～5mm范围内,增加流速可提高换热能力、减小阻力系数,但压降同样增加;相同速度时,提高流道高度不仅提高了换热性能,还减小了压降。     

基于蒙特卡罗法的水下激光光幕探测性能研究

为了提高水下激光光幕的探测性能, 根据水下光束传播规律, 构建了水下激光光幕探测模型, 基于蒙特卡罗模拟方法对水下激光光幕探测性能进行研究。基于水下激光光幕探测模型, 利用MATLAB软件进行仿真, 分析海水衰减系数、初始功率及传输距离对水下激光光幕传输的影响。仿真结果表明:海水的衰减系数越小, 水下激光光幕传输率受到的影响越小。海水参数和传输距离一定时, 随着初始功率的增加, 只会影响到达探测端的最终功率, 但对传输率影响不大。当海水衰减系数一定, 传输距离为1 m时, 其传输率约为15%且变化稳定; 当传输距离增加到30 m时, 传输率在5%以下。