基于黏性涡域法的二维钝体绕流模拟

针对传统计算流体力学方法需要划分网格且计算耗时长等问题，一种确定性涡方法(黏性涡域法)逐渐发展并日臻成熟。本文利用该方法计算模拟了二维钝体的流场、气动系数及斯特劳哈尔数，并与其他成熟CFD软件计算结果对比来验证模拟精度。黏性涡域法的模拟过程为生成涡量、涡量演变及气动计算。文章首先推导了该方法各过程的基本方程，然后基于Lua语言编制了圆柱和方柱绕流计算程序，最后利用Gnuplot实现流场可视化。利用黏性涡域法模拟不同雷诺数下圆柱绕流、不同无量纲频率和无量纲振幅的振动圆柱气动力时程、不同风向角下的方柱绕流，结果表明：雷诺数为2×10²～1×10⁵时，基于黏性涡域法的圆柱绕流气动计算结果与其他方法的计算结果基本一致；基于该方法的不同风向角下方柱绕流气动计算结果与其他文献结果趋势相同。     

一种改进的浸入运动边界算法

浸入运动边界是数值处理流固耦合作用的主要方法之一，其原始权重函数根据经验确定且形式固定，难以准确地定量描述复杂的流固现象。本文提出了一种新的权重函数对格子Boltzmann方法框架下的浸入运动边界进行修正，通过假定零固含率处的权重因子多阶导数为0，引入可调参数b增强权重函数对固含率的关联，实现了对固体边界渗透性能的改进。对静止圆柱绕流、Taylor-Couette流和振动圆柱绕流3个算例开展数值实验，结果表明：随着雷诺数的增加，原始权重函数预测的升阻力系数均高于文献结果，严重低估了边界处流体组分的作用；当新权重函数中的参数b为3时，改进的浸入运动边界效果最佳；新权重函数增强了固体边界的渗透率，提高了流场的预测精度，使平均误差降低了38.5%；最后，运用新权重函数成功地复现了振动圆柱绕流的锁频现象，表明改进的浸入运动边界在处理复杂的流固耦合问题方面具有广阔的应用前景。     

旋转FGMs层合圆柱壳行波模态频率分析

为推进功能梯度（FGMs）材料在发动机、导弹和火箭等领域的应用，旨在研究旋转FGMs层合圆柱壳的行波模态频率。采用Voigt模型和Sigmoid体积分数描述FGMs层合圆柱壳的材料属性，考虑科里奥利力、离心惯性力、环向初应力以及热内力推导了FGMs层合圆柱壳的能量表达式。采用切比雪夫正交多项式构造位移容许函数，建立了任意边界旋转FGMs层合圆柱壳的模态频率方程，并探讨了组分含量、夹层厚度、温度梯度和弹簧刚度系数等对FGMs层合圆柱壳模态频率的影响。结果表明:夹层厚度相比Sigmoid体积分数对模态频率的变化更为敏感；高旋转短薄壁圆柱壳相比长薄壁圆柱壳对边界条件和失稳现象的影响更为敏感；轴向弹簧相比其他弹簧对模态频率的影响更大。      

一种求解低速跨流域流固耦合问题的离散统一气体动理学格式

运动边界及流固耦合问题是低速及高速、连续流域及稀薄流域流动中常面临的多场耦合问题。本文在原始离散统一气体动理学格式基础上，采用任意拉格朗日-欧拉方法及动网格技术，并耦合计算结构动力学，实现适用于低速流动运动边界及流固耦合问题的计算方法。在方法验证方面，进行连续流域圆柱绕流强迫振动和自由振动数值模拟，计算结果与参考文献结果吻合良好，准确捕捉圆柱在锁频区和非锁频区的流动特性。在稀薄流域范围，将方法中的连续流D2Q9离散速度模型替换成跨流域离散速度模型，对比低速圆柱绕流强迫振动算例在两种离散速度模型下的数值结果。虽主要关注低速流动问题，但发展的算法仍具备模拟高超声速运动边界问题的能力，对此进行了Ma=5.0稀薄流振动圆柱研究。低速和高速运动边界模拟能力显示了方法在考虑稀薄气体效应流固耦合问题方面的潜在应用价值。     

亚临界雷诺数圆柱涡致振动的实验研究

近年来，数值模拟结果表明圆柱涡致振动现象可在低至20的亚临界雷诺数下发生，但目前尚无相关实验证明该种现象的存在。本文首先研发了适用于低雷诺数流固耦合实验的旋转水槽，而后在此基础上开展了亚临界雷诺数下圆柱涡致振动实验，分别研究了支撑频率和雷诺数对涡致振动响应的影响。发现圆柱涡致振动存在的最低雷诺数为23，与数值模拟结果接近，证实了亚临界雷诺数涡致振动的存在；另外，发生亚临界涡致振动时，在尾迹区观测到了频率与圆柱振动频率一致的卡门涡脱，说明弹性支撑特性会使得流动稳定性降低，证实了相关数值模拟结果的正确性。