圆角化对方柱气动性能影响的流场机理

方形截面柱体的圆角化处理是常用的流动控制方法,但其流场作用机理尚未被澄清.采用大涡模拟方法,在雷诺数为2.2$\times$10$^{4}$时,考虑风攻角的影响,对均匀流作用下的标准方柱和圆角方柱的气动性能和流场特性进行了研究,定量分析了圆角化气动措施和风攻角变化对分离泡特性的影响规律,从流场角度澄清了圆角化气动措施对方柱气动性能的影响机理.研究表明:与标准方柱相比,圆角方柱的表面风压、气动力和涡脱强度呈整体下降的趋势,但圆角方柱的斯特劳哈尔数更高;圆角方柱的"分离泡流态'发生在更小的风攻角范围内,分离泡的出现会进一步造成方柱的尾流变窄,涡脱强度减弱;随着风攻角的增大,分离泡的长度会逐渐减小直至消失,分离泡的中心会逐渐向方柱前角(迎风向)和方柱壁面移动;与标准方柱相比,圆角方柱的气流发生初次分离的位置向下游移动,分离后的剪切层更贴近方柱,因而更易发生再附现象;方柱尾流宽度的减小和涡脱强度的减弱是导致圆角方柱气动力减小和斯特劳哈尔数增大的主要原因.      

流场非均匀性对非平面激波诱导的Richtmyer-Meshkov不稳定性影响的数值研究

基于Navier-Stokes方程组,采用可压缩多介质黏性流动和湍流大涡模拟程序MVFT (multi-viscousflow and turbulence),模拟了均匀流场与初始密度呈现高斯函数分布的非均匀流场中马赫数为1.25的非平面激波加载初始扰动air/SF6界面的Richtmyer-Meshkov (RM)不稳定性现象。数值模拟结果表明,初始流场非均匀性将会影响非平面激波诱导的RM不稳定性演化过程。反射激波加载前,非平面激波导致的界面扰动振幅随着流场非均匀性增强而增大;反射激波加载后,非均匀流场与均匀流场条件下的界面扰动振幅差异有所减小。进一步,定量分析流场中环量分布及脉动速度统计量揭示了前述规律的原因。此外,还与平面激波诱导的RM不稳定性进行了简单对比,发现由于非平面激波波阵面区域的涡量与激波冲击界面时产生的涡量的共同作用,使得非平面激波与平面激波诱导的界面失稳过程存在差异。     

基于颗粒离散元的沙纹演变大涡模拟研究

应用大涡模拟、"点球" 浸入边界法和基于"事件驱动" 模型和"弹簧-阻尼" 模型的颗粒离散元法, 数值模拟了明渠湍流中沙纹的演变. 通过对不同谢尔兹数下无沙纹床面的推移质输沙率进行计算, 并与经典输沙率公式进行对比, 验证了模型的精度和可靠性. 随后, 对明渠湍流中沙纹床面的演变过程进行了数值模拟, 计算了有沙纹床面的推移质输沙率、沙纹长度和高度、等效床面高度的最大值、最小值和平均值、沙纹形状阻力、体积流速随时间的变化曲线. 研究发现:初始平整的床面在较短的时间内(tU_b=h≈100) 发展出数条沙纹, 随后沙纹逐渐发展, 在tU_b=h为1 600~2 000 时, 沙纹发生合并. 在沙纹数量不变的条件下, 沙纹高度随时间近似呈线性增长, 而沙纹的长度的平均值却保持恒定. 随着沙纹高度的增大, 输沙率和体积流速逐渐降低, 沙纹形状阻力则逐渐增大;当沙纹发生合并时, 沙纹高度快速增加, 输沙率、体积流速和沙纹形状阻力也出现了大幅跳跃. 在同等的水流强度条件下, 有沙纹床面的输沙率小于平整床面的输沙率.      

应力历史对颗粒材料强度和变形特征影响的离散元模拟

主要关注了颗粒材料前期所受的应力历史对其后期宏观力学响应的影响。该应力历史由一段等比例加载应力路径以及卸载垂直方向应力至与水平方向围压相同的卸载段描述。具体工作为:基于PFC2D双轴压缩数值实验,调查了应力历史对颗粒样本的强度、变形特征、细观参量如组构的影响,得出颗粒样本的偏应力-应变、体积应变曲线及其发生破坏时的名义应变云图。数值结果表明:高低围压下样本分别发生剪切破坏和弥散破坏,高围压下弹性阶段的刚度受应力历史影响较大,而低围压下样本在刚进入塑性至应力峰值点阶段的弹塑性刚度变化较为明显。随着应力历史中等比例加载系数的增大,剪胀加快,变形局部化范围有所不同;另一方面,颗粒形状的不规则性会增强颗粒材料的各向异性,导致样本强度更高。     

全附体潜艇粘性流场的RANS模拟及场量和涡量的校验分析

采用相同拓扑结构和相近网格质量的4套网格和5种湍流模型,对全附体Suboff潜艇粘性流场进行RANS模拟,分析了网格密度、节点空间分布规律和湍流模型对计算精度的影响,详细校验了其力积分量、速度场量和涡量特征。结果表明:网格密度最大的G4网格(140万)计算精度最高,总阻力较实验值误差为0.723%,其采用SST湍流模型时最优。计算得到的压力系数和剪切应力系数分布均与实验值吻合很好;桨盘面速度等值线分布计算精度与文献相当,轴向相对速度0.9以上的计算半径稍大于实验值,其余半径与实验吻合较好;桨盘面上0.25倍半径处速度分量沿周向分布计算精度较文献高,轴向分量与实验值吻合较好,径向分量峰值稍小于实验值,但峰值所处周向位置与实验值一致。成功捕捉到了附体端面绕流诱导对旋涡、附体叶根截面下游处项链形涡对、尾翼端面尾缘上方附着涡蹄、附体马蹄涡系、尾翼截面通道流体挤压作用诱导涡以及桨盘面涡量汇集的潜艇涡量场特征,且围壳端面绕流诱导对旋涡沿流动方向持续稳定,不影响桨盘面涡量场,均与文献中由大涡模拟模拟得到的定性结论一致。研究表明,在网格密度较大、节点分布合理、网格质量较高、湍流模型选取适当和壁面函数使用有效的条件下,RANS模拟潜艇粘性流场的场量和涡量特征同样具有很高的计算精度,能够在工程应用中有力支撑新型艇型设计与性能分析。     

含液颗粒材料液固耦合分析的离散颗粒模型及特征线SPH法

对含液颗粒材料流固耦合分析建议了一个基于离散颗粒模型与特征线SPH法的显式拉格朗日-欧拉无网格方案。在已有的用以模拟固体颗粒集合体的离散颗粒模型[1]基础上,将颗粒间间隙内的流体模型化为连续介质,对其提出并推导了基于特征线的SPH法。数值例题显示了所建议方案在模拟颗粒材料与间隙流相互作用的能力和性能以及间隙流体对颗粒结构承载能力及变形的影响。     

流对波边界层湍流运动的影响

采用大涡模拟方法和Smagorinsky亚格子模型,求解三维Navier-Stokes方程,研究了波流边界层中的湍流特性．将大涡模拟结果与相应的直接数值模拟结果和实验数据进行比较,吻合较好．获得了不同波雷诺数,不同波流比情况下的大涡模拟数据库,并由此分析了波流边界层中各种湍流统计量,如速度廓线、剪应力、湍流强度等的变化规律．     

基于颗粒物质力学的铁粉末压制中摩擦特性对力链演化影响

结合颗粒物质力学理论,通过离散元法实现铁粉末压制过程模拟并通过压制方程进行验证,针对粉末体系中的力链演化问题,提出力链特征定量分析方式,进一步通过分析不同颗粒间摩擦系数、侧壁摩擦系数与颗粒运动状态转变的方式,探讨摩擦特性对力链量化特征的影响,从而建立摩擦行为与力链演化间的联系. 研究结果表明:随颗粒间摩擦系数增大,整体力链数目变少,力链方向系数、承载不均匀度及单位屈曲度均变大,而随侧壁摩擦系数增大,力链特征差异较小,则颗粒间摩擦系数较侧壁摩擦系数对力链特征演化具有更显著影响. 同时发现,颗粒接触状态的改变与力链特征演化间具有对应性. 研究成果将进一步拓展粉末压制中考虑摩擦行为及力链演化过程在内的粉体致密化行为理论.      

探测器对超音速刚性盘-缝-带型降落伞系统的影响

研究了流体初始马赫数为 2.0 时, 探测器的存在与否对刚性盘-缝-带型降落伞系统气动减速性能以及流场流体结构特性的影响. 对于非定常可压缩流体的数值模拟, 流场采用了三层块结构自适应网格加密技术, 配合混合形式的TCD (tuned center difference)和WENO (weighted essentially non-oscillatory)计算格式以及基于拉伸涡亚格子模型的大涡模拟方法来处理超音速流中的激波以及大尺度湍流旋涡结构等. 结果表明: 无探测器时, 降落伞系统的流场结构稳定, 扰动较小; 有探测器存在时, 探测器后端的湍流尾迹和伞衣内部逆向运动溢出的流体与伞衣前端的弓形激波周期性的相互作用, 使得激波位置发生前移、激波倾角变小, 伞衣内部流场难以达到平衡稳定状态. 这加剧了降落伞系统的气动阻力振荡脉动变化, 降低了降落伞系统气动阻力系数, 同时也使得降落伞系统流场尾迹结构更加复杂.      

舵几何特征对桨-舵系统尾流场演化的影响

发生在桨和舵之间的干扰会影响螺旋桨尾流的演化,导致尾流场中的湍流在下游增强,恶化船舶的振动和噪声性能,深入分析舵几何参数对桨-舵系统尾流场演化的影响能够为推进器尾流场的调节和减振降噪提供新思路.因此,从弦长、剖面和梯形舵入手分析不同的舵几何参数对螺旋桨尾流场演化特性的影响,使用大漩涡模拟方法模拟流场中的湍流结构,对不同舵弦长、剖面下的螺旋桨尾涡结构演化进行了分析,在舵弦长、剖面影响螺旋桨尾流场演化的研究的基础上分析了梯形舵对螺旋桨尾涡结构的影响,进一步分析了梯形舵影响下的螺旋桨尾流场中湍动能的分布.结果表明舵的弦长和剖面均会影响螺旋桨尾流场的演化,这种影响表现为更大的弦长和更厚的剖面会促进螺旋桨梢涡在舵压力面上的偏移,更薄的舵剖面会带来更强烈的螺旋桨毂涡偏移;涡管轮廓和舵表面脉动压力的对比均表明梯形舵会促进螺旋桨尾流场沿逆舵梯度方向偏移,从而导致螺旋桨的尾涡结构在舵两侧及下游呈现不对称分布,桨-舵系统下游的湍流结构与螺旋桨尾涡-舵碰撞过程、螺旋桨尾涡-舵随边涡干扰过程、螺旋桨梢涡-螺旋桨毂涡干扰有关,偏移更大的螺旋桨尾涡结构会在尾流场中更早地引起湍动能增强.