HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

圆柱形汇聚激波诱导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时,地层压力保持较好,且试验区开发效果最好.