渗透壁面竖直平板Blasius流熵产特性

对渗透壁面竖直平板Blasius流层流边界层流动熵产及传热熵产与总熵产的比值(Bejan数)进行了深入研究,将耦合的动量与能量偏微分方程组通过相似变换转换为非线性常微分方程组,用Runge-Kutta法求解获得了无量纲流函数与无量纲温度的相似解,进而研究了无量纲参数对熵产及Bejan数的影响.结果表明:吸入时最大熵产位于壁面处,吸入速度越大最大熵产越大,熵产随相似变量的增加单调下降;喷注速度越大壁面处熵产越小;近壁面处熵产随变黏度参数增大而增大,远离壁面处随喷注/吸入参数或变黏度参数的增大而减小;远离壁面处Bejan数随喷注/吸入参数或变黏度参数的增大而增大,近壁面处与此相反;熵产和Bejan数随着辐射参数和Reynolds数增大而减小,随着Biot数增大而增大;Biot数在0~2范围内对传热熵产的影响很剧烈,Biot数较小时总熵产由流动熵产控制,Biot数较大时总熵产由传热熵产控制;熵产随滑移参数的减小而增大,Bejan数随滑移参数的增大而增大;Grashof数对熵产的影响较大,在较大的Grashof数条件下熵产随相似变量的变化较剧烈;滑移参数越大,流动熵产越小,Bejan数越大.      

喷注壁面楔状流层流边界层速度与温度相似解

建立了壁面有喷注的楔状流层流边界层冷却数学模型,求解经相似变换得到的描述无量纲流函数和无量纲温度的常微分方程,获得了楔状流层流边界层无量纲速度和温度的相似解,给出了考虑壁面喷注边界条件的形式简洁、物理意义明确的无量纲流函数拟合式;用Runge-Kutta法求解无量纲速度与温度的常微分方程,获得了壁面有喷注的楔状流层流边界层无量纲速度和温度随楔形角、吹风参数、冷却介质温度的变化规律.通过计算0°,18°和36°楔形角的楔状流层流边界冷却速度与温度分布得到以下结论:楔状流速度边界层和温度边界层厚度都随楔形角的增大而变小,随着吹风参数的增大而增大;当冷却介质温度越低、吹风参数越小、楔形角越大时,靠近壁面处楔状流层流边界层内温度梯度越大.      

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

喷注壁面外掠平板层流边界层速度和温度的HPM-Padé解析解

运用HPM-Padé(同伦摄动-帕德逼近)法导出了当喷注流体速度与距平板前缘距离的平方根成反比时,不可压缩常物性流体外掠平板层流边界层内无量纲速度和无量纲温度解析表达式,无量纲流函数1阶导数的HPM-Padé解析解与4阶龙格库塔法的数值解结果一致.利用HPM-Padé解析解研究了喷注系数和流体Prandtl数对速度和温度分布的影响.结果表明:Prandtl数越大,温度边界层越薄,壁面温度梯度越大;喷注系数越大,速度边界层厚度和温度边界层厚度越大,壁面速度梯度和壁面温度梯度越小,壁面喷注对平板有隔热作用;当喷注系数为0.619时,壁面的速度梯度和温度梯度为零,高温来流向壁面的传热被喷注流体完全阻隔.     

振荡管流换热的解析解及强化换热特性

经求解随时间正弦变化的非定常压力梯度下,管内流体振荡时的运动方程和能量方程,得到了振荡管流换热的速度场、温度场以及当量导温系数的解析解.通过分析可知,影响管内振荡流轴向换热的无量纲参数主要有:无量纲振荡频率,无量纲振幅和流体普朗特数.其中,无量纲振荡频率直接影响流动的速度分布,且随着无量纲振荡频率增加,强化换热效果增强.随着无量纲振幅的增大,管内流动有效导热面积增加,等温面法向温度梯度也增加,使得强化换热效果增强.强化换热效果同样也随流体普朗特数的增加而增强.      

旋转状态边界层速度场及温度场特性实验

针对涡轮转子叶片内冷技术,使用TR-PIV(time resolved particle image velocimetry)技术与热线技术同步原位测量了壁面加热条件下旋转通道内边界层速度场和温度场特性。结果显示:旋转数大于0.48时前缘面附近出现了回流现象,并从受力分析的角度给出了解释;回流区一般出现于流场下游、较大密度比、较高旋转数下,可以利用回流区的影响达到增强前缘面换热的目的;得到了旋转条件下无量纲温度型、温度脉动量和努塞尔数的变化规律,可以看出湍流边界层内部的温度场分布在旋转效应的影响下产生了强烈的不对称性,与静止条件下的标准规律相比会产生一定的偏差。     

格林函数法求解含移动介质轴向导热的稳态温度场

用格林函数法求解了考虑移动介质轴向导热影响的稳态温度场,推导了第1类和第2类非齐次边界条件下求解稳态温度场的数学表达式,利用特征值和特征函数获得了格林函数解,证实格林函数法可以解决非齐次边界条件下的稳态传热问题.以平行板间和圆管内均匀速度移动介质的稳态传热为例,获得了在半无限长壁面和在有限长壁面有均匀热流密度两种情况下所对应的温度场的无量纲解析解,分析了移动介质轴向导热与Pe和x/H之间的关系.结果表明:在所分析的两种热流密度边界条件下,Pe越大无量纲温度越小;在半无限大区域有均匀热流密度条件下,无壁面热流密度区域温度变化越剧烈,有均匀壁面热流密度的区域无量纲温度随x/H的变化几乎不受Pe大小的影响;在有限区域有均匀热流密度条件下,与半无限大区域均匀加热的情况有显著差异,在全计算域内温度变化越剧烈.     

先进旋涡燃烧室速度场与温度场协同数值研究

为了探讨先进旋涡燃烧室流动传热特性,基于场协同原理,对不同来流速度、来流温度、壁面温度及燃气当量比下燃烧室的速度场、温度场及其场协同角分布进行了数值模拟。结果表明,协同角较小的区域主要分布在后钝体后侧、凹腔内部以及进气通道横向中心截面上。旋涡区可以强化换热。随着来流速度及来流温度的增大,场平均协同角呈递减趋势;随着壁温的提高,场平均协同角增大;当量比小于1.0时,场平均协同角随着当量比的增大而增大,而达到1.0之后变化不明显。对于速度场与温度场,燃烧室中心截面的场协同性能最好,且体平均协同角大于面平均协同角。     

齿轮瞬态温度场的仿真分析

以直齿圆柱齿轮为研究对象,基于能量守恒定律和傅立叶定律推导了齿轮瞬态温度场的导热微分方程,根据定解条件确定了齿轮各个界面边界条件,运用有限元方法和传热学理论建立直齿圆柱齿轮模型,加载边界条件,并对其瞬态温度场进行仿真,得到了不同周期的温度场分布和节点温度变化曲线,系统地分析了其温度场随时间的变化.结果表明:温度随着啮合周期的增多而增高;在啮合阶段节点温度有一急剧升高,在退出啮合后进入非啮合阶段,温度逐渐降低;啮合阶段温升大于非啮合阶段温度的下降,该节点温度总体趋势升高.分析结果符合实际,为齿轮的热分析奠定了坚实的基础.     

多喷嘴构型的轴对称氢氧加热器设计与仿真

为了使氢氧燃烧加热器满足自由射流试验台工作需要并获得均匀的出口气流参数,采用同轴剪切式7个喷嘴轴对称构型喷注器,利用CFD仿真软件对其进行了三维反应流场计算,燃烧模型采用氢氧单步反应模型,获得了设计工况下的参数.计算结果表明:燃烧效率随着中心喷嘴与外围喷嘴距离L与喷注面板半径R之比(L/R)的增大先上升后下降;喷注面板的温度随着L/R的增大而降低,最终维持在600K左右;加热器出口的氧气摩尔分数以及总温的均匀性基本不随着L/R变化而变化;出口主流区的马赫数在6左右满足设计要求.在各个喷嘴的影响区域大致相等时,加热器综合性能良好.氢氧速度比越大,完全燃烧所需区域越短,喷注面板温度越高.与单喷嘴、19个喷嘴的加热器比较发现7个喷嘴的构型较为合理.      

动能BGK算法在近连续流模拟中的应用

首先比较详细地介绍了动能BGK算法的基本思想,然后应用该算法,数值模拟了克努森数(Knudsen number)为0.0001～1.56,马赫数为1.96绕圆柱的二维粘性流动,研究了滑移边界条件和无滑移边界条件对圆柱阻力的影响,通过与实验结果[5]进行比较,验证了文献[1]关于由连续流到自由分子流的这种流动模型的划分.计算结果表明:采用滑移边界条件,直到克努森数为0.3,计算得到的阻力与实验都符合得很好.因此,采用BGK算法加滑移边界条件,能够模拟由连续流到近连续区的流动.     

机身边界层控制研究

对带壁面抽吸的三维可压缩机身边界层控制进行了计算研究。采用了非正交曲线坐标系边界层方程,通过修改矢量流函数关系式,将壁面法向速度转变成为控制方程的一个可调参数,从而消除了边界条件与方程未知参量之间的非线性耦合;所发展的横向分段推进数值技术,可用于计算横流速度在横向多次反号的复杂边界层流动;设计不同的壁面抽吸系数分布,以研究边界层控制对多种边界层特性的影响。结果表明:壁面抽吸对边界层的控制效果是显著的,合理的抽吸分布能有效地减小机身阻力。     

基于数值模拟的流场附面层边缘识别方法

附面层边缘通常取在速度达到主流速度0.99倍的位置,而复杂流场中主流流动往往并不均匀,给附面层边缘的准确识别造成了困难。为解决此问题,提出了用"参考主流"代替实际主流识别附面层边缘的方法:通过零剪切力滑移壁面边界条件下数值模拟得到不受附面层干扰的参考主流,在根据附面层定义确定附面层边缘时以该参考主流中的速度代替实际的主流速度。通过斜楔压缩和弯曲压缩两个超声速压缩流场对该识别方法进行了验证,所得到的斜楔压缩出口截面上附面层厚度与采用实际主流速度判断得到的厚度相对误差仅4.1%。根据该方法的识别结果对弯曲压缩型面设计进行附面层修正后,弯曲激波高度与无黏设计值之间的误差从修正前的2.0%降低至0.3%,压缩面末端压力的相对误差从修正前的6.6%降低至2.3%。该方法避免了指定主流速度的主观性,识别结果较为准确。     

扩散界面浸入边界法结合壁面模型在湍流模拟中的应用

提出了一种模拟高雷诺数湍流的扩散界面浸入边界法（IBM）。该方法采用壁面模型来减少壁面附近的网格量。为了实施壁面模型和边界条件，在物面外部设置了2个辅助层（一系列拉格朗日点）：外侧的参考层用于实施壁面模型来得到壁面剪应力，内侧的强制层用于实施离壁的边界条件。在实施壁面模型时，将动量方程沿物面法向积分，从而把物面切向的速度修正量与由壁面模型得到的壁面剪应力联系起来，而速度的法向分量则按二次曲线分布来近似重构。在实施边界条件时，为了严格满足无穿透条件，应用了基于隐式速度修正的IBM。最后，利用绕平板湍流和绕NACA0012翼型湍流来对方法的可行性进行了验证。     

吸附式压气机叶栅气动性能计算模拟研究

为考察附面层吸附技术在压气机静子势流区叶型上的应用,采用流场数值计算方法对吸气叶栅流场进行模拟.结果表明:(1)对于高亚声速压气机叶栅,采用吸力面附面层吸除可提高叶栅的扩压度,但不一定能减小流动损失.(2)对于中亚声压气机叶栅,采用吸力面附面层吸除不仅可提高叶栅的扩压度而且能减小流动损失.(3)如果叶片吸力面靠后缘处有流动分离,吸气位置在分离区的上游较远处可抑制分离,若在分离区附近可能不利于抑制流动分离.      

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。