二维平板水漂运动数值模拟

探索平板水漂运动的流体力学现象机理对飞行器着水问题的研究有重要的参考价值。基于有限体积法和k-ε RNG湍流模型求解非定常雷诺平均Navier-Stokes（URANS）方程,采用流体体积分数（VOF）模型与速度入口造波法构造数值波浪水池,并结合整体动网格方法,完成二维平板在静水面与波浪水面的水漂运动数值模拟。在与试验值和理论值对比的基础上,讨论初始姿态角、投掷角度与投掷速度对水漂运动的影响。进一步研究不同波浪参数和波浪位置对平板水漂的影响,并从能量角度展开分析。结果表明：20°姿态角平板能以最小的投掷速度实现水漂运动;运动中平板能量的相对损失率受初始投掷速度的影响较小,主要受投掷角和姿态角作用,随投掷角或姿态角的增加而增大。波浪情况下,在上行波位置触水的平板能够获得更大的接触面积,而更长的触水时间发生在波谷处;因此,这两个位置触水的平板在水漂运动过程中的相对能量损失大,其数值变化比波峰位置大5%左右;在上行波位置触水时,平板的速度衰减与相对能量损失随着波高的增加而增大,而波峰位置则有相反的变化趋势。     

高速平板着水数值模拟

探索和揭示物体入水冲击的流体力学现象与机理对飞行器水上迫降问题的研究有重要的参考价值。对高速平板着水涉及到的复杂物理问题展开数值模拟,采用有限体积法求解非定常雷诺平均Navier-Stokes(URANS)方程和标准k-ε湍流模型,流体体积(VOF)模型捕捉水气交界面,整体动网格技术处理平板与水面的相对运动。在二维楔形体入水冲击的算例验证基础上,详细研究平板高速着水引起流体喷溅、射流、空气垫等现象和平板底面压力变化历程,结果表明:空气垫现象明显,俯仰角4°平板下表面出现规律的空气泡,10°时则不存在;平板下表面的水体沿壁面运动,当俯仰角为10°时,壁面水体的运动速度显著增加;在大俯仰角的情况下明显出现负压区。     

基于光滑粒子流体动力学的波浪中航行体入水数值模拟

航行体在波浪条件下高速入水的运动响应和载荷特性是其研发设计过程中需要重点考虑的问题,为了对该问题进行精准预测,采用无网格光滑粒子流体动力学（smoothed particle hydrodynamics, SPH）方法,提出了一种新型周期性波浪边界技术,利用四元数法计算物体六自由度运动,建立了波浪条件下入水模拟的数值水池。通过对静水中方块体垂直落水和航行体倾斜入水运动轨迹和冲击载荷的模拟,对比于实验参考结果,验证了数值模型的计算精度。随后,在SPH数值波浪水池中对航行体在不同波浪相位角下的高速入水过程开展研究,结果表明航行体弹道稳定性受到波浪相位角影响显著,0°相位角入水时弹道最为稳定。该新型SPH数值水池能够实现航行体波浪中入水过程的精确预报。     

基于边界造波法的波浪数值模拟

<正>海上风力机塔架的波浪力计算是一项比较复杂的工作,其与风力机安装点的波浪、水深、及塔架的几何形状等密切相关。研究近海风力机在波浪作用下的载荷分析模型,并通过适当的计算方法得到正确的波浪载荷,是进行海上风力机载荷分析的前提。传统的物理实验和理论分析一直是这类问题的主要研究手段,随着计算流体力学的发展,数值仿真越来越多地应用于海上风力机风波联合作用下载荷计算。波浪的数值模拟是分析波浪与塔架相互     

波浪情况下民机水上迫降性能数值分析

为了研究波浪情况下民机水上迫降性能,选用计算流体力学的有限体积法求解非定常不可压缩RANS方程,基于VOF法、整体动网格法、斯托克斯五阶波浪模型和自适应网格技术,以空客A320-200为研究对象,按照适航规章相关要求与建议,构建了民机水上迫降数值仿真模型。首先,对比分析了静水面和波浪水面下飞机水上迫降过程。结果表明：相同初始飞行参数下,波浪水面下最大水平过载为2.42 g,是静水面的1.09倍;最大垂向过载为4.82 g,是静水面的2.82倍。2种情况下,空气垫效应和潜水现象明显,碰撞初段,飞机均受到吸力作用;波浪情况下,跳跃现象出现,但并未对飞机的运动过程产生剧烈影响。其次,研究了波浪参数对水上迫降性能的影响。分析结果表明：同一波峰撞击阶段内,波高越高,则最大水平过载和最大垂向过载越大;波长越长,则最大垂向过载越小;波高越高,波长越长,则最大下沉速度越大。     

非结构聚合多重网格法流场数值模拟研究

 在非结构网格上应用聚合多重网格技术来加速 Euler方程的收敛过程。设计了一种高效率的网格聚合方法。采用四重网格,在每一层网格上采用有限体积法进行计算。通过对多段翼型绕流的数值求解验证了多重网格加速收敛的高效性。     

基于数值模拟的气动载荷计算方法研究

飞机气动载荷的获取一般要根据多项实验结果,经过飞行参数的仿真模拟,实验结果的处理等过程才能获得,而在飞机的初步设计阶段不可能进行大量的实验,但需要计算飞机的载荷和寿命.本文提出了一种利用数值模拟来确定基本气动载荷的方法,对于偏副翼的情况,用工程计算方法获得偏转引起的附加量,加到基本气动载荷上得到偏副翼时的气动载荷,该方法已用于某飞机型号设计.     

微射流作动器三维外流场数值模拟

采用有限体积法求解Faver平均N—s方程，对三维、粘性、非定常、不可压微作动器外流场进行数值模拟。计算结果揭示了微射流产生、发展与耗散的过程，和实验中观察到的现象比较吻合。旋涡对沿展向的不稳定性是涡环破裂、旋涡对的相干性消失的主要原因。二次涡的发现为以后研究微射流旋涡对的细微结构打下良好的基础。     

阵风响应及减缓的非定常数值模拟

基于网格速度和动网格方法,通过求解非定常N-S方程计算翼型副翼偏转时的阵风减缓响应,研究开环控制律参数和副翼缝隙对阵风减缓效果的影响。对于单个1-cos阵风,随着偏转角度的增大,减缓效果变大;由于粘性的影响,CFD计算得到的升力峰值减小量明显呈非线性的趋势,减缓效果低于准定常线性理论。对于连续sin阵风,通过调整副翼偏转的频率和最大偏转角,升力能达到很好的减缓效果,此时,副翼偏转频率应接近阵风频率,最大偏转角存在一个最优值。在阵风环境下,副翼间隙对升力系数影响不大,而力矩系数增大一倍,间隙会降低副翼偏转对升力的减缓效果。     

非结构多重网格流场数值模拟研究

本文在非结构网格上应用多重网格技术来加速Euler方程的收敛过程。此方法在一系列粗网格和细网格之间进行值的传递，并且粗网格和细网格之间结构上无关。在每一重网格上用有限体积法进行计算。采用本文的方法数值求解了NACA0012翼型的绕流。     

基于动态边界控制的优化方法在高升阻比气动外形优化中的初步应用

本文是利用动网格技术和有限体积法来研究基于动态边界控制的气动优化的方法.该优化方法是通过在求解非定常流动的过程中计算优化的目标函数随设计变量变化的梯度.而利用参数化的结构动网格技术能有效、快速地生成与优化过程出现的外形对应的网格并进行计算,具有较高的效率.本文以典型的高升阻比外形优化为例,用非定常NS方程为流动控制方程进行优化,取得了初步的结果.     

三维亚声速冲击射流流场的数值模拟

利用有限体积法对三维可压缩的N－S方程进行离散，对喷嘴的亚声速垂直冲击射流和斜冲击射流进行了数值模拟。网络的划分采用非结构性网络，湍流模型为RNGκ-ε湍流模型，在近壁处采用壁面函数进行修正，流场计算结果与实验值吻合一致。计算可得到垂直冲击流场结构的三个区域：自由射流、滞止区和壁面射流区域，并显示出垂直冲击点附近区域的卷吸回流、速度分布。冲击射流近壁处存在着回流现象，由于边界层流动和近壁回流引起距离固壁2D内的速度梯度变化较大。斜冲击在近壁处平行壁面速度的最大值，以及壁面上的总压最大值都偏离冲击点，并且相对应。     

基于HLLC格式的理想MHD方程数值模拟

针对理想MHD(magnetohydrodynamics)方程数值求解困难的问题,基于原始的HLLC (Harten-Lax-Van Leer Contact wave)近似黎曼解方法,发展出一种新的适用于理想MHD问题的通量计算格式。控制方程采用有限体积法离散,时间推进采用隐式的LU-SGS（lower-upper symmetric Gauss-Seidel）格式,并且采用双曲型散度清除技术来抑制磁场散度的累积。通过一维激波管问题的数值模拟表明,HLLC-MHD格式能准确的分辨并捕捉复杂的磁流体力学波系结构,保证数值计算的精度,最大的数值计算误差不超过10%;通过二维的Rotor问题的数值模拟研究表明,HLLC-MHD格式能够应用于多维理想MHD问题的数值模拟,并且能够准确捕捉多维情况下磁流场中的阿尔文波;通过对比有、无散度清除的Rotor问题计算结果表明,双曲型散度清除技术可以将磁场散度峰值从50降低到2,有效抑制高磁场散度区域的散度累积,同时也会将误差传播到低散度区域,并且引起边界处散度的累积,影响计算的稳定性。      

Numerical simulation of a high-speed propeller flowfield based on Euler equations

Ｔｈｅｐｒｅｄｉｃｔｉｏｎｏｆｆｌｏｗｆｅａｔｕｒｅｓｏｖｅｒｈｉｇｈ－ｓｐｅｅｄｐｒｏｐｅｌｅｒｉｓａｖｅｒｙｉｎｔｅｒｅｓｔｉｎｇａｎｄｉｍｐｏｒ－ｔａｎｔｐｒｏｂｌｅｍ,ｗｈｉｃｈｈａｓｃｏｎｓｔａｎｔｌｙｃｈａｌｅｎｇｅｄｍａｎｙａｅｒｏｄｙ...     

初始投放条件对内埋式导弹分离轨迹的影响

采用有限体积法对欧拉方程进行离散,耦合求解六自由度运动方程,运用非结构动网格技术计算了导弹在三种不同初始条件下从三维开式弹穴中投放后的分离轨迹和姿态变化情况,分析了各种现象产生的原因,给出了较为理想的内埋式导弹发射方式。     

超声速欠膨胀冲击射流的数值模拟

利用有限体积法对轴对称可压缩的N－S方程进行离散，对喷嘴的超声速垂直冲击射流进行了数值模拟，网格的划分采用非结构性网络，得到详细的波系结构，激波形状和位置与实验相吻俣，大压比下，壁面上冲击中心区域压力相对低，分布较平坦，压力最大点所在位置偏离中心冲击点，并且在近壁处有回流。     

大水滴撞击壁面的动态特性数值模拟

采用 VOF 方法建立了大水滴撞击壁面的计算模型,模拟了大水滴以不同直径、不同速度撞击光滑壁面的动态撞击过程和撞击特性。计算结果表明：在大韦伯数情况下,水滴撞击光滑壁面会在铺展过程中发生边缘水滴飞溅;在水滴撞击壁面的收缩阶段,随着水滴直径的减小和水滴速度的提高,会发生液膜缓慢收缩、边缘液环和液膜分离、中心处部分液膜和边缘液膜分离、液膜完全破裂等不同结果;当水滴直径和撞击速度增大时,同一时刻水滴的铺展半径、最大铺展半径、最大铺展系数增大;水滴在壁面上达到最大铺展系数所用的无量纲时间随水滴直径增大而增大,同一直径水滴不同初始速度下达到最大铺展系数所用的无量纲时间变化较小。