简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。     

k-kL两方程湍流模型的改进及验证研究

提高湍流数值模拟的准确性,从而明确湍流模型对数值模拟结果的影响具有重要的意义,应用K.S.Abdol-Hamid给出的尺度自适应k-kL两方程模型封闭RANS方程,并修改von Karman长度尺度的限制方法,通过平板、翼型、后台阶等流动的模拟,考察k-kL模型在湍流模拟中的准确性,及其反映主要流动特征的能力和网格收敛性,并对影响流动模拟准确性的因素进行讨论.结果表明:改进长度尺度限制之后的k-kL两方程模型无论是对附着流动还是分离流动都可以给出比较准确的结果.     

二维弯曲等截面管道中的激波串特性研究

利用Carroll的试验数据验证了数值方法的有效性之后,对二维等截面弯曲管道中激波串的特性进行了数值模拟试验,研究了管道弯曲对激波串的结构与特征长度、壁面沿程静压分布、出口截面马赫数与总压恢复、反压特性等的影响,研究中考虑了不同的进口马赫数和边界层厚度。结果表明,管道弯曲对流动的对称性有着明显影响,当马赫数较高时(如Ma₀=2.45)合适程度的管道弯曲有利改善直管道已有的流动不对称,使激波串长度缩短。管道弯曲能够有效抑制出口压力变化所导致的出口截面马赫数的大幅波动,考虑到低压比时(出口Ma_e>1)直、弯管道之间总压恢复系数存在明显差距,而当压比较高时(出口Ma_e<1)两者相当接近,因此亚燃发动机的超声速扩压器可适当使用大曲率以缩短管道的轴向长度。另外,鉴于弯曲管道与直管道内激波串长度之间的明显差异,已有的基于直管道的激波串长度经验公式不能很好地适用于弯曲管道。     

扩压器内跨音速湍流的数值模拟

采用Johnson-King非平衡代数雷诺应力湍流模型(J-K模型)和Baldwin-Lomax零方程湍流模型(B-L模型),数值模拟较强激波/边界层相互作用时扩压器内的分离流动。计算结果与实验值进行了比较,表明J-K模型比B-L代数湍流模型可较好地计算出分离流动的再附点位置,并且可更好地计算出激波强度和沿流程的压力分布,仅增加很少的计算量,并更易推广应用于三维湍流问题的数值模拟。      

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

高超声速绕钝体熵层基本流特性研究

基于N-S方程,针对超声速气体绕0°迎角钝锥流动采用三阶WENO(Weighed Essentially Non-Oscillatory)格式进行数值模拟,研究了熵层的基本流特性。研究给出了边界层及熵层外缘的确定方法,并分析了其适用性与局限性。数值计算结果显示基本流剖面存在"双广义拐点"现象,N-S方程与Euler方程数值结果对比发现,熵层区出现的广义拐点是由熵层区有旋运动引起的,而边界层的广义拐点由粘性作用导致。这是进行流动稳定性分析的基础。文章还分析了球头半径与马赫数对熵层特性的影响。结果表明,马赫数越大,熵层区的涡量越大,熵层影响的范围越大;雷诺数对无粘流的特征影响很小,熵层大小随球头半径几乎呈线性关系。     

应用一个新的具有转捩敏感性两方程低雷诺数湍流模型预测转捩流动

本文基于Launder-Sharma(LS)低雷诺数的 两方程湍流模型发展了一种新的具有转捩敏感性的低雷诺数湍流模型。针对LS模型的可实现性(realizability)问题和前缘滞止点湍动能预测过大的不足,新模型进行了改进。新模型只使用当地物理量,不需要求解壁面距离、Y+和边界层积分参数。新模型能够适用于广泛的流动,且容易应用到通用的CFD程序中。对具有详细数据的零压力梯度平板转捩边界层T3A实验的模拟结果显示,新模型能够预测转捩流动,并能对自由湍流变化给出合理的响应。 关键词：航空、航天推进系统;转捩模型;转捩流动;数值模拟     

盖板式预旋系统的压比和熵增特性

对预旋系统内的压力变化相关研究较少。基于理论分析、实验测量以及数值计算,对某盖板式预旋系统的压比及熵增特性进行研究。通过理论推导,对预旋系统内压比与无量纲温降的关系进行分析。在最高转速可达10000r/min的高转速实验台上,测量了转盘上的气流静压以及相对总温,进而获得压比及熵增特性。进行三维数值计算,将数值计算结果与实验结果进行了对比,并根据数值计算结果对预旋系统内的熵产分布以及各元件的熵增情况进行分析。结果表明:系统温降以及旋转马赫数大小决定了预旋系统的理想最大压比,而实际压比与理想压比的比值取决于系统内的熵增大小。采用数值计算以及实验测量所得结果对理论关系式进行了验证,最大偏差2.7%。旋转马赫数一定的条件下,随系统无量纲温降增大,系统压比逐渐减小。由于熵增影响,实测压比与理想压比最大相差约36%。预旋系统内的熵增主要发生在预旋腔静止壁面、接受孔前后、供给孔进口等气流旋转比发生剧烈变化的区域。预旋系统内主要元件的熵增随流量增大都呈逐渐增大的趋势,但接受孔处熵增最小值出现在喷嘴出口旋转比等于1左右时,流量过小或过大都会导致接受孔处熵增变大。     

固体火箭冲压发动机进气道三维数值模拟

基于RANS平均的N-S方程、标准k-ε双方程湍流模型,利用FLUENT软件对某型固体火箭冲压发动机进气道内外流场进行了三维数值模拟.研究了不同马赫数下进气道前体、入口、进气道出口流场和流动特性;考虑了弹体,缩比,网格划分对进气道流场出口参数的影响.模拟结果表明:随着飞行马赫数的增加,由前缘点产生的激波逐渐靠近弹体,进气道内部总压恢复系数降低,在进气道出口处压力比较高,出口流场基本均匀.      

一种新型的动态RANS/LES混合方法

为解决RANS/LES混合方法过渡依赖网格疏密的问题,使用熵函数fs作为交界面函数构造一种新型的动态RANS/LES混合方法。同时采用脱体涡模拟方法 (DES97),对平板边界层流动、超声速凹腔-压缩拐角流动以及超声速底部流动进行数值模拟,并与实验数据进行对比。结果表明,无论对于附体流动区域还是分离区域,fs都能够动态地指示出边界层的范围。构造的动态RANS/LES混合方法能够显著缓解DES97由于网格因素导致的雷诺应力模化不足的现象,同时保证分离区域涡结构的求解精度。另外,熵函数fs是从能量耗散的角度进行构造,并不依赖于特定的湍流变量,因此更具有普适性。     

熵,信息及其互补原理与守恒假说

本文从玻耳兹曼熵和香农熵的概念，推演出熵与信息的互补原理。并提出了熵量守恒与转化定律假说。     

复熵及其应用——Bayes-E决策分析法

 本文将熵的概念推广到复数域,并讨论了它的性质,进行了后验预分析,从而克服了Bayes风险决策方法对信息分析的不足,得到一个Bayes-E风险决策方法的新测度——传递熵。     

基于最大熵原理的顺磁性固体磁化研究

从最大熵原理出发,在经典和量子 2种情况下,研究了顺磁物质的磁化规律,得到了顺磁物质磁化强度的表达式,在高温弱磁情况下,磁化强度与外磁场成正比,与温度成反比,得到的结论与统计物理中的系综理论或最该然分布结论一致。该原理不仅适用于平衡态问题,而且也适用于非平衡定态问题。     

化学反应过程方式的判据及反应方向的判据

本文对化学反应过程的方式及方向判据问题进行了简要分析,并指出用△G~Q代替△G判断化学反应过程的粗略界限。     

DHGF与熵技术在风险评价中的应用

首先介绍DHGF原理与熵技术方法,同时针对在DHGF过程中主观因素,利用熵技术加以修正。在实例应用中,针对某商场的应急疏散系统,建立风险评价指标体系,应用DHGF和熵技术结合的方法进行风险评价,结果表明该方法确定的风险等级能够和保险公司承保的条件相符合,准确性较高。     

死亡是生命的必然归宿

本文根据物质系统演化的规律，平衡原理及熵理论，阐述了死亡是生命系统演化的必然归宿。     

基于融合信息熵距的转子裂纹-碰摩耦合故障诊断方法

针对转子系统的裂纹-碰摩耦合故障,提出了一种基于融合信息熵距的转子振动故障诊断方法.利用转子实验台模拟转子系统裂纹故障、碰摩故障及裂纹-碰摩耦合故障,并采集三种典型故障的振动加速度信号.利用时域的奇异谱熵、频域的功率谱熵、时-频域的小波能谱熵以及小波空间特征谱熵,计算融合信息熵距实现故障诊断.实例研究表明:这四种信息熵形成了综合评价转子振动状态的特征指标,多测点、多转速下的信息熵距曲线较好地区分了单一故障和耦合故障,有效地提高了转子振动故障诊断的准确性.测试信号与其对应的振动故障之间的信息熵距最小,信息熵距曲线位于坐标轴的最下方,达到了诊断单一故障和耦合故障的目的.      

转子振动故障的过程功率谱熵特征分析与定量诊断

针对传统旋转机械振动故障定性诊断的不足,提出了1种以信息熵理论为基础的转子故障特征分析与定量诊断方法。在转子试验台上模拟转子振动的4种典型故障,分别得到4个测点多转速下的振动过程故障数据;对这些故障数据进行分析和处理,提取反映其振动过程的故障特征——功率谱信息熵,建立能描述转子振动过程变化规律的多转速多测点下的故障信息熵矩阵,并对振动故障进行分析;通过对转子振动故障信号的实例计算和定量诊断分析,验证了该方法在转子振动故障分类和故障严重程度判断方面是可行的。     

基于相关信息熵的传感器故障检测研究

针对连续性过程系统中存在的大量强相关关系测点，根据测量数据间的关联特性，提出了一种应用相关信息熵进行传感器故障检测的方法。理论和实验分析证明，该方法具有简便、直观的传感器故障检测特点，较适合工程实际应用。     

熵损失下无失效数据的 Bayes估计

讨论了在熵损失下定时截尾数据的Bayes估计问题，将无失效数据看作逐次定时截尾试验的结果，进而得到无失效数据的Bayes估计，利用相关定理证明该估计是容许的。通过数值例子，进一步说明所得估计的合理性。