基于混合RANS/LES方法与FW-H方程的气动声学计算研究

在气动噪声的数值计算过程中,非定常流动的求解精度对声学计算结果有着重要的影响.以机体噪声计算的标准算例双圆柱绕流气动噪声问题为研究对象,采用基于非线性k-ε湍流模型的限制数值尺度(LNS)方法对双圆柱绕流进行了数值模拟,将计算得到的气动特性和流动特征与相应的试验结果进行了对比分析.为了求解远场观测点处的气动噪声,在精确求解双圆柱绕流流动的基础上结合基于FW-H(Ffowcs Williams-Hawkings)方程的声类比方法进行数值计算,并通过圆柱体的展向相关性对计算结果进行了修正,将得到的最终结果与相应的声学试验结果进行了对比,两者吻合良好,表明该数值方法是准确、可靠的.     

冠状喷口抑制涡扇发动机喷流噪声试验和数值研究

 以揭示发动机喷流噪声声场分布和冠状喷口抑制喷流噪声为目的,对涡扇发动机排气噪声声场和喷流流场进行了试验和数值计算,对比分析不同结构的冠状锯齿对喷流噪声的综合抑制效果。结果表明:采用大涡模拟(LES)、FW-H (Ffowcs Williams-Hawkings)声学模型和傅里叶变换的方法,可以预测出喷流噪声的空间分布和声压级值;冠状锯齿产生的阵列流向涡是强化流体间混合、降低喷流速度峰值、抑制喷流噪声的根本所在;在本文研究参数范围内,冠状锯齿结构不仅将排气喷流噪声由原型排气系统的159.6 dB降低到154.3 dB,而且冠状锯齿给喷管带来的损失控制在1.7%以内。锯齿冠状结构对分开排气式大涵道比涡扇排气系统喷流噪声有抑制作用。     

壁湍流调制雷诺应力和拟序结构应变率空间形态的层析PIV实验研究

采用高时间分辨率的层析PIV技术,测量了水洞壁湍流三维速度分量空间分布的时间序列,应用局部平均速度结构函数概念对流向速度信号进行多尺度分解,以流向速度分量的局部平均速度结构函数过零点作为特征量检测壁湍流中拟序结构猝发的喷射和扫掠过程,应用空间相位平均技术提取拟序结构猝发的喷射和扫掠过程各速度分量、涡量分量、拟序结构速度应变率分量以及调制雷诺应力分量的空间相位平均拓扑形态。为了引入平衡态涡粘模型假设模拟调制雷诺应力,研究了拟序结构猝发过程中调制雷诺应力分量和速度变形率分量的空间分布形态,发现两者之间的空间相位分布不一致。由于存在时空相位不同步性,说明需要考虑大尺度拟序涡结构引起动量传递的时空弛豫效应。应用经典的线性平衡态下的Boussinesq涡粘模型不能准确地描述壁湍流拟序结构动量传递非平衡现象的物理机理。对于壁湍流拟序结构动力学方程中调制雷诺应力的模拟,应采用包含时空相位信息的复涡粘张量模型。由于雷诺应力与速度变形率的时空不同步性,对非平衡非局部湍流场的数值模拟提出了一个挑战性的问题,建议采用包含相位信息的复涡粘张量模型来模拟雷诺应力张量,从而更加符合雷诺应力演化的物理机理,这一模型有可能成为一个很有发展前景的封闭模型,从而更加准确地预测工业领域中广泛存在的非平衡湍流。     

不同结构的冷却管脉冲射流流动响应特性

探究了结构参数对脉冲射流经冷却管内流动响应规律.在脉冲频率为1Hz,雷诺数为5000,脉冲射流占空比为50％的气动参数下,对冷却管管径、射流孔孔径、孔间距的变化对脉冲射流流动响应特性带来的影响依次分析,对结构参数对脉冲射流在冷却管内流动响应特性产生影响的内在机制进行综合阐述.针对上述结构参数拟合出冷却管射流质量系数随冷...     

基于OpenFOAM的超临界压力下低温射流大涡模拟

为了深入认识低温射流在超临界压力下的流动特性,将真实流体热物性模型写入开源CFD程序OpenFOAM平台中,开发出了用于模拟超临界压力下流动过程的均相求解器。针对液氮跨临界、超临界射流进行大涡模拟,研究了环境条件变化对于低温射流流动特性的影响。结果表明：采用的真实流体模型能够在广泛的温度与压力范围内准确计算流体的热物理性质,在此基础上开发的均相求解器能够准确描述超临界条件下的低温射流,其流动特性主要表现为稠密液体与环境气体之间的剪切层不稳定以及湍流混合,此类变密度射流在湍流充分发展区域同常密度射流一样具有自相似特性;超临界射流相比于跨临界射流具有更高的混合效率以及更短的液核穿透长度;环境温度升高以及环境压力降低均使得低温射流与周边气体之间的密度梯度增加,对剪切层中不稳定波动的发展起到抑制作用,导致了射流液核长度的增加以及扩张角度的减小。     

弹黏塑性材料稳恒扩展II型裂纹尖端应力场

采用弹黏塑性力学模型,对II型动态扩展裂纹尖端应力场的对数奇异性,进行了分析计算。详细地分析了黏性系数α、马赫数Ma²对裂纹尖端的应力场影响。指出了对数奇异性区域存在的问题,解释了过渡区的成因,对过渡区尖端场解的形式和求解方法做了合理的推测。     

圆自由射流多尺度涡结构迁移特性的实验研究

在流场背景湍流度比较大或流场中涡结构的尺度范围比较宽的情况下，为了研究湍流场中不同尺度涡结构的迁移速度，用IFA300双通道恒温热线风速计和两个热线探针测量圆自由射流中心线上距离一定的空间两点同一时刻的流向速度分量，用子波分析提取了每一个探针输出的不同尺度的脉动速度信号，计算了两个探针输出的同一尺度脉动速度信号的互相关函数，利用互相关函数达到最大值对应的延迟时间研究了湍流多尺度涡结构的迁移特性。实验结果表明：不同尺度涡结构具有不同的迁移速度，小尺度涡的迁移速度较大，大尺度涡的迁移速度较小。其中，一部分尺度的涡结构是湍流中的主要结构，该尺度涡结构对湍流的贡献也是主要的，其迁移速度与实验测得的湍流平均速度较为接近。     

适用于激波/边界层相互作用的线性涡粘性湍流模式

本文选择了两个激波 /边界层相互作用诱导分离的跨声速问题 (轴对称圆弧突起和二维管道突起 ) ,采用五个有代表性的湍流模式 (BL模式 ,JL k-ε模式 ,k-ω模式 ,SST模式和双尺度模式 ) ,通过将数值计算结果和实验结果进行比较 ,对有关的湍流模式进行了评估和分析     

三角翼大迎角风洞试验支架干扰数值模拟研究

现代战争要求战斗机能够在大迎角（AOA）状态下进行过失速飞行,对飞机大迎角绕流流场的研究主要的方法有风洞试验和数值模拟。在大迎角风洞试验中,常用的是尾支撑方法,支架的存在会对模型的试验结果产生一定的影响,本文通过数值模拟来对这个影响进行研究。以开源计算流体力学软件OpenFOAM 2.3为平台,采用PIMPLE算法求解Navier-Stokes（N-S）方程,PIMPLE算法是SIMPLE（Semi-Implicit Method for Pressure-linked Equations）算法和PISO（Pressure Implicit with Splitting of Operator）算法的结合体;采用基于有限体积的空间离散方法和空间二阶精度的线性插值方法,时间离散采用后向差分方法,湍流模型采用SA-DDES（Spalart-Allmaras-Delayed Detached Eddy Simulation）模型。为了验证方法的可靠性,首先对0°、10°、30°、50°、70°以及90°迎角下的有支架三角翼绕流流场进行计算,并将计算结果与试验结果进行对比,两者吻合较好。在此基础上,数值模拟了无支架的三角翼绕流流场,对比有/无支架情况下数值模拟结果,得到支架对三角翼绕流流场、背风面压强分布和气动力的影响。计算结果表明：大迎角情况下,有支架与无支架时相比,支架的存在会影响三角翼附近的流场（但是不会改变涡系等流动结构）、改变翼表面压强分布,从而导致三角翼的法向力系数和俯仰力矩系数发生明显变化。     

连续旋转爆轰发动机冷流场的混合特性研究

连续旋转爆轰发动机（CRDE）中燃料和氧化剂的快速掺混是实现爆轰波成功起爆和稳定传播的重要前提,然而目前国际上关于这方面的研究还相对较少。本文采用大涡模拟（LES）方法,对非预混CRDE中燃料和氧化剂的混合过程及其主要机理开展深入研究。研究结果表明,非预混CRDE流场中存在欠膨胀特征、大尺度涡结构,以及回流区等复杂的流动现象,其中由于Kelvin-Helmholtz（K-H）不稳定性产生的大尺度湍流涡结构是促进氢/氧混合的主要机制。此外,本文还考察了氧气喷注位置对非预混CRDE的流场结构和混合特征的影响,发现氧气喷注位置会影响射流剪切层形态、涡尺度,以及回流区分布等,进而影响氢气和氧气射流的混合过程和混合程度。与其他进气位置相比,氧气在靠近内壁面喷注时更有利于氢/氧的快速掺混。