压缩折角激波-湍流边界层干扰直接数值模拟

进行了来流Mach数2．9,24°压缩折角激波．湍流边界层干扰的直接数值模拟,在上游的平板添加扰动以激发边界层转捩到湍流．计算得到的统计结果与实验吻合,验证了结果的可靠性．分析了角部分离区附近湍能的生成、耗散及分配机制．结果显示角部区域激波与湍流边界层相互作用造成大量湍动能产生,而湍动能的主要耗散区仍在近壁．湍流输运项起到了主要的平衡机制,把湍动能由外层输运到近壁区．通过对激波后及壁面瞬时压力的分析,认为激波低频振荡并非上游扰动弓』起,而是由于分离泡本身不稳定振荡产生的．     

二维激波边界层干扰的数值分析

该文利用二维可压缩雷诺平均Navier-Stokes方程，采用k-ε湍流模式，运用半离散有限体积法，对捆绑火箭助推器和芯级之间的激波边界层干扰进行了二维数值模拟．数值结果表明，在激波边界层的干扰下，壁面边界层流动出现了分离和再附过程，压力和温度出现激烈变化．     

可压缩效应对平板湍流边界层影响的直接数值模拟

采用高精度紧致型差分格式求解三维可压缩Navier-Stokes方程,通过在局部平板上引入周期性吹吸气小扰动的方法,直接数值模拟来流马赫数M∞为2.25的空间发展的可压缩平板湍流边界层.所得流场统计特征与相关理论和试验符合较好.在此基础上研究可压缩效应对平均流动特征以及湍能的生成和耗散特征的影响,指出在M∞=2.25的情况下,Morkovin假设基本成立,但可压缩效应使得拟序结构有明显差别.与不可压相比上抛、下扫对不同区域的作用有所不同.在近壁区,压力-膨胀项和压力-速度相关项仍是不可忽略的量,这导致声的产生和可压缩效应对湍能产生抑制作用.平板上声压的分布表明,随流动向下游的发展,边界层转捩过程激发了由物面脉动压力所产生的偶极子声源.     

超音速翼型激波边界层干扰的数值模拟

通过求解可压缩雷诺平均N-S方程（RANS）和S-A湍流模型,数值模拟了菱形翼型在超音速情况下的流动,并对在超音速飞行中襟翼偏转产生的激波一边界层干扰现象进行了分析和比较。计算得到的压力系数与理论计算的结果比较吻合,力矩系数也符合理论分析的结果。     

求解激波／边界层相互干扰问题的有限体积方法

提出了一种高效的求解二维激波／边界层相互作用问题的算法．应用该方法分析了超音速绕压缩拐角的流动情况,与有关实验数据比较,结果令人满意．本算法采用Ｒｕｎｇｅ－Ｋｕｔａ法对时间进行离散,湍流采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型．本算法可以推广到其他激波／边界层相互干扰的问题．     

超音速平板边界层湍流的直接数值模拟及分析

用时间模式对来流Mach数为4.5的平板边界层湍流进行了直接数值模拟. 发现当湍流充分发展后, 其平均流剖面、湍流Mach数、各种量的脉动均方根值和Reynolds应力等沿平板法向的分布都具有相似性. 但压缩性效应已较强而必需考虑, 强Reynolds比拟不再有效, Morkovin假说不再成立. 从转捩完成至湍流充分发展之间有一过渡过程, 其间上述相似性不成立.     

可压缩尖锥边界层湍流的直接数值模拟

采用高精度差分方法对来流马赫数0．7,来流Reynolds数250000／Inch,锥角为20°的尖锥边界层的整个空间转捩过程进行了直接数值模拟．对流项采用了7阶迎风格式离散,黏性项采用6阶中心格式离散,时间推进为3阶Runge—Kutta方法．对转捩形成的充分发展湍流进行了统计分析,包括平均速度分布,近壁湍流强度和雷诺应力等统计数据与平板边界层理论及实验吻合很好,验证了结果的正确性．显示了近壁湍流的典型拟序结构——高、低速条带结构并根据流向速度的周向相关量确定了条带的间距,以当地壁面尺度度量的条带间距沿流向并没有显著变化．给出了柱坐标下的可压湍动能发展方程,并据此对近壁湍动能的生成、耗散和输运机制进行了分析．     

边界层流动中湍流斑形成机制的数值研究

边界层流动中层流向湍流的转捩是一个古老但仍没有解决的问题．虽然对转捩过程的探索已付出了艰苦的努力,但依然有一个重要的物理过程还没有弄清楚,即转捩过程中湍流斑形成的原因以及湍流斑的运动特征是什么,这些问题正有待于人们做进一步的研究．文中提出以壁面局部脉冲形式的小扰动来模拟边界层流动中湍流斑的初始扰动场,采用紧致有限差分的方法,数值模拟了边界层流动中湍流斑的生成和演化规律;结果显示,它们在好多方面反映出湍流斑的基本特征．     

Mach8的平板可压缩湍流边界层直接数值模拟及分析

对来流马赫数等于8,壁温等于10.03倍来流参考温度的平板可压缩湍流边界层做了直接数值模拟,计算涵盖了从层流到转捩以及最终充分发展湍流的全空间演化过程.对湍流的统计特征做了详细的分析,结果表明,在当前的计算工况下,湍流边界层核心区平均速度剖面仍然满足对数率,且卡门常数基本不变;可压缩效应明显增强,由于采用近似恢复温度的等温壁条件,使得近壁区温度较高,导致当地声速增大,使得湍流马赫数绝对值较低,造成内在压缩性效应不强,与经典强雷诺比拟相比,除在数量上产生一些偏差外,强雷诺比拟关系近似成立,且Morkovin假设依然有效;对扩展自相似性和标度率分析表明,对于平板可压缩湍流边界层而言,高马赫数流动使得其适用范围减小;压缩性效应对近壁湍动能,条带结构,涡等值面分布的影响得到分析.     

平板边界层湍流的数值分析

对于不可压缩平板边界层的时间模式的直接数值模拟(DNS)而言,在转捩后继续计算,就可以得到湍流直接数值模拟结果.分析计算结果发现,当湍流充分发展后,其平均流剖面、各种量的脉动均方根值和雷诺应力等沿平板法向的分布都具有相似性.不过,从转捩完成至湍流充分发展之间有一过渡过程,其间上述相似性不成立.这一结果为简化但又较准确地计算工程技术问题中的湍流提供了一种可能性.此外,还分析了其他一些湍流特征,如形状因子和位移排移厚度随时间的演化规律.为了观察到相干结构,须将展向的计算域缩小一半,以增加分辨率,研究结果表明,湍流边界层近壁区存在相干结构,其主要表现是准流向涡或涡对.     

顺流平板的湍流边界层

本文作者将自己建议过的圆管湍流运劝的普适速度分布,用于不可压缩流动中顺流平板湍流边界层的分析,通过 Karman 动量积分,得到一组简明、准确的参数分析公式;进一步完善了内层变量法,给出了阻力系数 C_D 的积分式,以及相应的功应力公式和阻力公式。另外对部分重要阻力公式作了分析比较。     

平板紊流边界层探讨

本文是笔者已发表的“管中紊流新探索”一文的继续，笔者用该文提出的关于脉动切应力的假设，求解平板紊流边界层问题，给出了计算边界层厚度、位移厚度、动量厚度和能量厚度的公式和表格，所得的绕平板的阻力系数，无论在紊流水力光滑区，或是在紊流水力粗糙区，全与现有实验成果非常符合，对从水力光滑到水力粗糙的过渡区，本文提供的算法是首创．     

溢流反弧段紊流边界层变化规律的数值模拟

在溢流反弧段流线同心圆假定的基础上,通过理论推导,建立了表征反弧段动水压强、壁面势流流速、壁面切应力的相应表达式;将描述边界层发展变化规律的连续方程、动量积分方程、动能积分方程联立,并运用列主元高斯消去法将其转化为一阶微分方程组,采用双边法求其数值解,从而模拟出反弧段水深、边界层厚度以及流速分布指数的沿程变化规律。     

轴对称空腔内圆盘流动的数值模拟

针对轴对称空腔内圆盘流动，采用２种紊流模型（即混合长模型、混合长模型与高Ｒｅ数ｋ－ε方程相结合的分区模型）进行了数值模拟．计算结果表明：当圆盘旋转雷诺数Ｒｅθ较低时，流场中存在一个类似于固体旋转的核心区；在高Ｒｅθ数、小间隙比时，核心区消失．与实验值相比，ｋ－ε／ＭＬ模型得到较好的预测结果     

中浓纸浆悬浮液湍流边界层的模拟

介绍纤维悬浮液的分类及管道中纸浆纤维悬浮液的流动状态 ,重点进行中浓纸浆悬浮液湍流边界层的模拟研究 ,详细分析了相关的模拟结果 ,得出如下结论 :当管子壁面距离 0～ 9× 10 - 4m时 ,横过水环厚度的速度分布是线性的 ;水环厚度很小 ,水环中的剪切作用是层流剪切 ,剪切力作用于纤维网络表面 .突出网络塞体表面 (0～ 9× 10 - 4m )的纤维甚少 .据此可认为网络塞体表面较为光滑 ;管道中远离壁面处浓度高 ,壁面处浓度低 ,浓度分布具有递减的趋势     

一种含回流的复杂湍流流动的数值模拟

Ｔ－型绕流问题是一种典型的含回流的复杂湍流流动问题，具有一定的学术与工程意义。该文用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组结合ｋ－ε湍流模型数值模拟了这种流动，并与实验数据进行了详细对比。结果表明，壁面压力和时均速度场的预测精度是令人满意的，湍流量预测值的变化趋势和实验值一致，但在数值上有差异。文中指出了造成数值预测值与实验值差异的主要原因。     

液-液旋流分离管湍流流场的数值模拟

采用雷诺应力模型(RSM)对一种旋流管内部的单相流动进行了湍流数值模拟。通过数值计算得到了旋流管内部流动的速度矢量分布情况,计算结果与经典文献资料吻合,证明了模型的正确,为进一步数值模拟打下了基础。     

三维加力燃烧室湍流燃烧的数值模拟

本文对涡扇发动机三维加力燃烧室内的气相湍流燃烧过程进行了数值模拟。湍流模型采用标准ｋ－ｅ模型,湍流燃烧采用涡旋破碎（ＥＢＵ）模型,数值方法采用ＳＩＭＰＬＥ算法。计算结果定性合理。