基于混合方法的二维圆柱绕流气动噪声数值预测

采用计算流体力学与声比拟理论相结合的混合方法,探讨了二维圆柱模型在低马赫数下的非定常绕流诱发的噪声算问题。首先采用雷诺时均Navier-Stokes方程计算了二维圆柱模型的非定常流场,与相同参数条件下的计算数据和实验数据进行对比,得到了可靠的流场数据;然后对流场数据进行处理,利用声比拟理论将流场数据转化为声源,通过求解柯尔(Curle)方程计算了远场噪声。将声场的计算结果与文献中的结果进行了对比,得到了较好的噪声计算结果。     

基于深度神经网络的含运动边界非定常流场预测方法研究

针对飞行器设计过程中对流固耦合系统快速预测的需求,探索基于数据驱动的非定常流场建模策略,缩短流场演化求解耗时,从而加快流固耦合系统模拟速度。流固耦合系统中流场演化求解部分等价于含运动边界的非定常流场演化。本文提出了一种基于神经网络的流场预测模型,来学习并预测运动边界附近的非定常流场演化过程。此神经网络可以基于当前流场状态及边界运动信息预测出下一个时刻的流场状态。通过不同振动频率及振幅下的运动圆柱绕流问题,测试了本文提出的神经网络模型的预测精度及泛化能力。神经网络的预测结果和计算流体力学仿真结果中流场结构吻合度较高,且通过对预测的流场数据中边界上的压力积分得到的气动力也具有较高的精度。测试结果展示了该神经网络模型具有良好的预测性能,因此该方法可以用于快速、准确获得运动边界周围非定常流场状态。     

基于比拟理论计算圆柱/翼型的气动噪声

采用计算流体力学(CFD)与"声比拟"相结合的方法,通过两个步骤模拟流动的声学远场。第一步,在包含所有声源的近场区域内,通过求解URANS方程获得控制面处的非定常流场参数;第二步,采用基于可穿透数据面的FW-H方法模拟声学远场。该方法与传统的半经验方法相比,计算精确,易于工程实现,并且可以计算非线性噪声。根据气动声学计算积分公式,构造了高阶的时间离散格式,采用高斯积分公式计算圆柱/翼型的气动噪声。在时域和频域上计算了观测点处的声压及声压级随时间步数的变化情况,其结果与国外实验结果对比取得了较好的一致性,此外研究了翼型尾流部分界面对计算精度的影响。     

类升力体外形俯仰阻尼特性数值研究

采用有限差分方法求解薄层近似的非定常Navier-Stokes方程，定常流场采用交替方向隐式分解的NND格式，非定常流场采用四步Runge-Kutta方法，在保证时、空二阶精度前提下引入变系数残值光顺技术提高非定常流场计算效率，复杂带翼外形的空间网格通过求解抛物化的椭圆型方程生成，最后在Etkin理论下给出球锥及类升力体外形的俯仰阻尼导数计算结果。     

耦合POD重构舰面流场的直升机舰面起降数值模拟

为解决舰面非定常流场数据量过大的问题,采用本征正交分解(POD)方法对舰面流场进行重构,发展了一种耦合POD重构流场的直升机舰面起降数值模拟方法。首先采用计算流体力学(CFD)方法计算舰面非定常流场,获得离散数据样本;然后提取流场的POD模态,并截取能够捕捉到原流场主要特征的少量模态对原流场进行重构;最后建立耦合重构舰面流场的直升机高阶飞行动力学模型。以直升机返航进场为例进行数值模拟,并将计算得到的操纵量和飞行状态与飞行试验结果进行对比。结果表明:使用POD方法重构后的舰面流场数据约为原始样本数据的8.5%,且重构流场与原始流场吻合良好,POD方法能够解决舰面非定常流场数据量过大的问题。与飞行试验数据的对比表明,本文方法捕捉到了舰面非定常流场对直升机的影响,可用于直升机舰面起降研究。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

压气机叶片流固耦合数值计算

采用流固耦合方法对压气机单排跨声叶片进行气动弹性数值模拟.首先对压气机流场进行三维定常求解,并与试验结果进行比对.在定常流场分析的基础上,进行流固耦合计算,采用非定常方法分别求解流体域Navier-Stokes(N-S)方程和固体域结构动力学方程.耦合过程中,流体域和固体域在物理边界上进行能量传递.计算结果表明:所采用的流固耦合数值方法用于判断叶轮机械叶片气弹稳定性是可行的,可以得到叶片的瞬态响应,从而判断叶片是否发生颤振.      

高超声速滚转阻尼导数数值模拟

采用非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程描述物体简谐振动流场,并在Ｅｔｋｉｎ理论下给出绕定轴转动时滚转阻尼导数的计算公式,定常流场的计算采用空间二阶精度的交替方向隐式分解的ＮＮＤ格式,非定常流场的计算采用时、空二阶精度的Ｒｕｎｇｅ－Ｋｕｔｔａ多步格式,采用代数和方法生成物体静、动网格。最后给出高超声速下钝体外形滚转阻尼导数的计算结果,以及滚转力矩系数随瞬时振幅的变化曲线。     

基于流场/声爆耦合伴随方程的超声速公务机声爆优化

基于自主研发的大规模并行结构化网格CFD求解器PMB3D以及并行化伴随方程求解器PADJ3D，开展了流场/声爆伴随方程的求解研究。首先采用标准算例，对内部CFD代码PMB3D软件和声爆预测代码进行了声爆计算可信度验证，以及声爆强度对近场声压梯度的校核。针对并行环境下多块对接网格的近场声压提取操作的复杂性，提出了"包围盒"的方法实现并行环境下近场声压装配单元编号、网格块编号以及对应的进程编号确定，基于声爆计算坐标将并行传递的数据进行一维排序，为声爆预测、伴随方程以及梯度求解提供输入条件。通过线性插值雅克比矩阵实现均匀坐标系梯度信息向非均匀坐标转换，并进一步根据结构化网格特征提出了插值原则，简化了近场声压转换雅克比矩阵的变分。通过装配单元记录，实现声爆强度对流场守恒变量的变分结果向各个进程装配，将装配结果作为流场伴随方程的右端项实现流场声爆耦合伴随方程的求解。此外，对小型超声速公务机开展了声爆优化，对比分析了设计前后的声压及其频谱特性。     

基于本征正交分解和代理模型的流场预测方法

 为了实现流场的快速求解,基于本征正交分解(POD)与代理模型提出了一种全新的流场预测方法。其基本原理为:首先采用本征正交分解将一定数量的样本流场分解为同等数量的基模态流场;然后用少数包含了绝大部分样本流场特征的基模态流场拟合所有的样本流场;最后用代理模型建立起决定样本流场的输入参数与拟合系数之间的近似函数关系。针对几何外形不同的二维翼型定常流场预测结果表明:在亚声速情况下,预测误差收敛的模态数量不超过20个,继续增加模态的使用数量不能明显提高预测精度;在跨声速情况下,预测误差收敛的模态数量为26个,当使用的模态数量达到前10个时,继续增加模态的使用数量能提高绝大部分流场区域的预测精度,但同时会在激波附近引入"噪声"激波特征而降低该区域的局部预测精度。在这两种情况下,预测流场所需时间均不到高精度计算流体力学(CFD)方法的1/200。     

超临界翼型尾缘噪声影响因素研究

使用大涡模拟(LES)和FW-H相结合的方法对超临界翼型尾缘噪声进行了研究。针对超临界翼型后缘的不同设计参数(后缘厚度、后缘角和几何形状)对尾缘噪声的影响,使用ANSYS FLUENT的LES湍流模型计算声源,采用FW-H积分方法求解远场噪声总声压级。首先完成二维非定常流圆柱绕流的流场及其噪声的验证计算,计算结果与实验值符合,证明了求解器设置和网格生成的合理性。然后基于此正确的求解器设置和网格生成,使用LES/FWH对比了典型超临界翼型的不同后缘设计参数对远场总声压级的影响,所得结论对低噪声超临界翼型优化设计有参考意义,同时为进一步的噪声控制优化设计提供了基础。     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

高速列车绕流的数值模拟

给出了用低速Ｎ－Ｓ方程求解高速列车绕流的数值计算方法，采用显隐式结合的方法来加快收敛速度，基于该方法编制了计算软件系统Ｔｒａｉｎ ３Ｄ，用该软件系统对高速列车外形的外流场进行了数值模拟，获得了该外形的空气动力特性和流场结构。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

基于运动嵌套网格的前飞旋翼绕流N-S方程数值计算

 通过求解 Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间上采用中心平均的有限体积法进行离散,时间方向应用含隐式子迭代的双时间法推进求解。为了验证程序的正确性,数值计算了一有升力悬停流场,旋翼桨叶表面压力分布的计算值与实验值吻合很好。     

细长体简谐振动绕流的数值模拟

用数值模拟方法研究了椭圆锥细长体简谐振动非定常过程的流场特性。数值模拟的出发方程是可压缩流的全Ｎ－Ｓ方程，计算格式为ＴＶＤ格式。数值模拟中，研究了细长简谐振动绕流的非定常流场特性，不同振动缩减频率和压缩性效应的影响，揭示了非定常流场的特性和流动机理。     

基于非线性谐波法评估压气机转静干涉诱发的强迫响应

针对压气机转静干涉诱发的强迫响应问题,对某型风扇试验件进口可变弯度导向叶片进行流固耦合数值模拟。基于NUMECA的非线性谐波法,求解在后排转子势流作用下可调导叶的非定常流场,从而获得谐波压力,并对其作用在前排可调导叶相应模态下的强迫响应进行了计算,与试验验证结果吻合。结果表明:基于非线性谐波法的流场计算以及强迫响应分析,为工程应用建立了转静干涉诱发强迫振动响应的数值评估方法。     

高超声速圆球模型飞行流场的数值模拟和实验验证

对高超声速圆球模型飞行场进行数值模拟，分别采用空气空气气体模型，平衡气体模型以及热化学非平衡11组元气体模型求解非定常轴对称N－S方程组，使用有限差分时间相关法捕捉激波，得到了定常流场的解，差分方程稳式部分采用了LU－SGS方法以避免矩阵运算，对化学反应和振动能量源项采用预处理矩阵以解决刚性问题，由计算结果处理得到的阴影图和干涉条纹图与再入物理弹道靶实验照片进行了对比分析，验证了实验中圆球飞行流场大部分区域接近于平衡状态。     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数