超声速湍流混合层实验图像的分形度量

以高时空分辨率的纳米平面激光散射（NPLS）实验技术为基础,在SML-1风洞中完成了超声速混合层的流动显示实验．相应的实验图像清晰地再现了层流、转捩及湍流区的流场结构,空间分辨率满足分形度量的要求．给出并比较了测量分形维数的两种常用方法,采用计盒维数法测量了超声速混合层转捩区和完全发展湍流区的分形维数．转捩区的分形维数随着湍流脉动的增强而增加．完全发展湍流区的分形维数并不会因为流场结构不同而有较大的变化,表面上杂乱无章的湍流界面具有基本相同的分形维数,体现了湍流流动的自相似性．     

Gao-Yong湍流模型对边界层转捩的适用性研究

准确预测边界层转捩对深入理解湍流边界层发展和飞行器翼型优化工业应用等具有重要意义。为提升Gao-Yong(G-Y)湍流模型对平板边界层发展及转捩特性预测的准确性与适用性，该文提出了采用G-Y漂移位矢Reynolds数(Re_T)表征当地湍流脉动强度的方法，建立实度系数(Cs)与Re_T间的函数关系，解决G-Y湍流模型中Cs取值间断、分布离散的问题。基于开源OpenFOAM软件开发了改进G-Y湍流模型的求解器，对经典平板边界层和转捩实验开展计算流体力学(computational fluid dynamics,CFD)数值模拟并验证。研究结果表明，相比k-ω等传统雷诺平均(Reynolds averaged Naviers-Stokes,RANS)模型，改进后的G-Y湍流模型不仅能准确模拟边界层流动特征，还能合理预测边界层发展过程及转捩位置，准确性较高。结果表明，经过改进的G-Y湍流模型可进一步用于对边界层转捩特性的研究。     

高升力多段机翼的转捩预测研究

本文采用Menter等发展的基于局部变量的γ-Reθt湍流转捩模式对平板和NLR7301翼型开展了转捩模式验证工作,进而对带有前、后缘襟翼的多段机翼流动转捩进行了数值模拟.研究结果表明,转捩对于多段机翼的气动特性具有较大影响,数值模拟可为多段机翼中襟翼偏角、间距等设计提供有益指导.     

用颗粒相双尺度二阶矩湍流模型模拟突扩两相流动

基于将颗粒脉动分成湍流引起的大尺度脉动和颗粒间碰撞产生的小尺度脉动的概念，应用单相流动多尺度湍流模型方法，从双流体模型出发，推导和封闭了颗粒的双尺度脉动雷诺应力方程、大尺度脉动能量传递率方程、小尺度脉动耗散率方程和两相脉动关联方程，用统一的湍流模型方法构造了包含颗粒两类脉动的新的颗粒双尺度二阶矩两相湍流模型．用该模型对突扩管内的两相流动进行了数值模拟，并和单尺度的二阶矩两相湍流模型进行了比较，计算结果和实验数据吻合较好，比单尺度的二阶矩两相湍流模型的计算结果有所改进．     

基于γ-Re_θ转捩模型的垂直轴风力机气动特性研究

在两方程的剪切应力输运SST k-w全湍流模型中添加γ-Re_θ转捩模型,以某两叶片垂直轴风力机为研究对象,数值计算得到叶片表面压力的分布规律,并与文献中的实验结果进行了对比.在此基础上,分析和对比了添加转捩模型后,风轮的转矩和推力系数、叶片表面极限流线分布和叶片表面摩擦阻力系数的变化规律.结果表明:相比SST k-w全湍流模型,添加γ-Re_θ转捩模型后的四方程Transition SST转捩模型预测得到的风轮的转矩系数和推力系数的均方根分别下降了6.13%和0.55%;计算时考虑γ-Re_θ转捩模型能够预测垂直轴风力机叶片表面边界层的转捩现象和更详细地捕捉叶片表面复杂的流动特征.     

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

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

翼面吸气混合层流控制的数值研究

本文采用Menter发展的γ-Reθ转捩/湍流模式预测平板和超临界RAE-2822翼型的转捩特性,验证了该转捩模式的有效性和可靠性,同时获得网格分布等规律.在此基础上开展某超临界翼型上表面及其前缘结冰时吸气层流控制后的转捩预测,分别获得干净和结冰外形下混合层流控制对转捩的影响规律.结果表明,合理设计的层流控制对干净翼型表面转捩推迟明显,能有效减阻;在结冰情况下,层流控制几乎失效.     

改进的k-w模型在亚音速漩涡流动中的应用

基于非结构网格,采用经典的wilcox k-ω模型和其改进的kω-Pω模型,建立了用于模拟大攻角旋涡流动的计算方法。以尖前缘的65?三角翼为例,模拟了旋涡的产生、发展、破裂过程,验证了wilcox k-ω模型和kω-Pω模型在典型的亚音速计算状态下对复杂涡系干扰的模拟能力。通过对多种计算的流场与气动力详细结果的比较分析,就两种湍流模型对大攻角复杂旋涡流动的预测能力和敏感性等进行了评估。结果表明：kω-Pω模型通过r值区分剪切层和涡核区域,从而对涡核区域的涡粘性进行修正,对最后的模拟结果有一定的修正作用,可以作为湍流模型修正的一个方向;RANS方法在预测涡破裂点位置和二次涡的强度及位置方面仍存在很大的缺陷。     

超声速边界层气动光学效应及控制研究

光学成像制导技术是精确制导武器研究的重要课题，气动光学效应的存在对制导精度产生严重的影响，使飞行器偏离目标位置，甚至脱靶。飞行器和来流之间的相互作用导致其周围流场结构非常复杂，激波、膨胀波以及湍流边界层的存在，使得流场的折射率在空间上分布不均匀，在时间上存在高频脉动。其中，边界层转捩过程引起的气动光学畸变是亟待解决的最大难题。边界层转捩过程是边界层由层流向湍流发展的过程，是一个多因素耦合影响的强非线性、复杂的流动现象，并且一直是湍流研究领域的热点问题。超声速边界层的转捩区具有非常强的非定常性、随机性，且脉动频率高，因此折射率分布极为复杂，光线透过后不仅会发生偏折，还会导致成像模糊和能量分散，严重影响成像制导的精度，通过光学校正的方法很难降低这种影响。本文试图通过延迟边界层转捩，来降低甚至消除超声速边界层转捩对光学传输性能的影响。以超声速旋成体飞行器光学窗口周围的边界层流场为研究对象，本文采用添加扰动片的流动控制方法对边界层气动光学效应变化规律进行了数值模拟研究。研究结果表明：随着攻角的增大，光程差分布由沿着对称面对称分布逐渐过渡为沿流向的变化，且逐渐平缓；施加了扰动片控制后，原本处于光...     

改进的K-ω模型在亚音速漩涡流动中的应用

基于非结构网格,采用经典的wilcox k-ω模型和其改进的kω-Pω模型,建立了用于模拟大攻角旋涡流动的计算方法.以尖前缘的65°三角翼为例,模拟了旋涡的产生、发展、破裂过程,验证了wilcox k-ω模型和kω-Pω模型在典型的亚音速计算状态下对复杂涡系干扰的模拟能力.通过对多种计算的流场与气动力详细结果的比较分析,就两种湍流模型对大攻角复杂旋涡流动的预测能力和敏感性等进行了评估.结果表明:kω-Pω模型通过r值区分剪切层和涡核区域,从而对涡核区域的涡黏性进行修正.对最后的模拟结果有一定的修正作用,可以作为湍流模型修正的一个方向.RANS方法在预测涡破裂点位置和二次涡的强度及位置方面仍存在很大的缺陷.     

考虑磁场影响的k-ε湍流模型推导及感度分析

针对电磁精炼及电磁连铸过程金属熔液流动现象,确定了磁场作用下导电流体湍流场中脉动电磁力的计算方法·结合Navier Stokes流体运动方程,推导了考虑磁场影响的k ε湍流模型·将修正后的k ε湍流模型应用到流动控制结晶器的流场计算中·结果表明,使用修正的k ε湍流模型计算的湍流动能分布区域缩小,平均流动动量略微变大·使用修正前与修正后的k ε湍流模型计算的速度分布差别较小,但更倾向于抑制湍流的产生·     

雷诺时均模拟的两种转捩模型比较

对两种基于RANS方程的转捩模型进行分析讨论,这两个模型都只需采用局部变量计算,其中SST-γ-Reθt模型基于实验经验公式,层流动能模型基于转捩现象。两个模型都被耦合到RANS低雷诺数湍流模型中,并通过不同工况下T3系列平板的数值计算对模型进行评估。结果表明转捩模型较原湍流模型需要消耗更多计算资源;相对层流动能模型,SST-γ-Reθt模型对y+更为敏感,在无压力平板算例中两个模型能够较为准确地预测转捩,只有在湍流度增大时,SST-γ-Reθt模型预测的精度才会下降;在有压力梯度情况下,两个模型预测转捩起始点都较实验值延后,当来流雷诺数较大时层流动能模型预测的转捩长度较实验值偏小;从平板上不同位置(层流区、转捩区和湍流区)湍动能的分布情况能够看出SST-γ-Reθt模型只是在数值上模拟转捩过程,并不考虑转捩内在规律,层流动能模型预测结果与实验值较为吻合。     

边界层人工转捩的大涡模拟及壁面脉动压力研究

为了验证边界层人工转捩技术实现船模边界层壁面脉动压力载荷相似的可行性,采用大涡模拟以及功率谱估计获得平板湍流边界层壁面脉动压力功率谱,并与试验值、半经验模型值进行了对比,验证了数值模拟的可靠性。采用大涡模拟获得了绊线上下游壁面脉动压力的均方值、自功率谱和波数频率谱;以脉动压力均方值、自功率谱特征为判据,对比了方形、锯齿形绊线的转捩效果。结果表明:绊线高度雷诺数大于吉宾斯雷诺数;两种绊线均可实现边界层转捩,且壁面脉动压力的功率谱特征、谱级相近;与同流速平板湍流边界层相比,即使在当地雷诺数较低时,绊线下游边界层壁面脉动压力自功率谱平台区仍然高出10 dB左右。设计的绊线可用于减弱船模边界层壁面脉动压力的尺度效应,并有助于实现壁面脉动压力载荷相似。     

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

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

由GF(pm)上的分组码构造p元域GF(p)上的分组码

Pless证明了三元（12,6,6）Golay码具有一种双层桔构,并据此给出了该码的快速硬判决译码算法．马建峰等人推广了Golay码的Pless结构,给出了由三元（n,k,d）线性分组码构造的三元（3n,n＋k,≥min（n,2d,6））线性分组码．本文证明了由任意域GF（p^m）上的（n,k,d）分组码构造p元域GF（p）上参数为（（m＋2）n,n＋mk,≥min{n,2d,2（m＋2）}）的分组码的可行性和方法．这种码具有很好的代数结构,可以快速译码．     

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

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

三峡水轮机尾水管涡带的CFD数值模拟

采用非定常雷诺应力湍流模式对三峡水轮机尾水管内的流动进行了数值模拟．分析了涡带的运动规律及其产生的水压力脉动的频率特性与相位特性和计算结果产生误差的原因,给出了改进的措施．研究结果表明,采用非定常雷诺应力湍流模式可以预测到尾水管内水压力脉动的频率特性,与现场试验的结果符合．     

使用PAFV湍流模型对叶栅流动的数值研究

为进一步发展更有效的湍流模型,提高对高负荷跨声速压气机流动的数值计算精度,以侧偏平均脉动速度和应变率张量为基础构造了新的侧偏平均脉动速度(partial average fluctuation velocity, PAFV)湍流模型。在使用Python计算机语言进行编程计算过程中,将流动控制方程在一般曲线坐标系下变换后,采用有限差分方法进行空间离散。方程对流项采用Steger-Warming通量向量分裂法,扩散项采用中心差分格式离散,时间导数项采用具有TVD(total variation diminishing)性质的二阶Runge-Kutta法离散。使用该湍流模型对德国航空航天研究试验院的L030-4叶栅亚声速和超声速流动进行了计算,在进行网格无关性研究基础上,获得了流场压力等值线云图、马赫数等值线云图,压力系数分布等结果。计算显示侧偏平均脉动速度在激波附近区域受到明显抑制,避免了激波附近区域出现过大的湍流黏性。数值计算结果与实验符合比较好,初步验证了PAFV湍流模型可以有效地模拟该类型的叶栅流动。     

非定常流动和转捩为湍流的数值模拟

欧洲流动、湍流和燃烧研究联合体（ERCOFTAC）于1990年举行了专题讨论会,来自12个国家38个科学和工业组织的90多位研究人员参加了这次会议,本书是该会议的论文集,1992年发行了第一版,本版于2008年出版。专题讨论会分为4个试验案例组：在S形导管中的边界层;水道中周期排列的圆柱;在自由流动湍流影响下边界层中的转捩;轴对称约束喷气的流动。最后一组分为两种情况：在稍微扩散的管道中的喷气和在开放空间中的喷气排放。这些试验案例反映科学界和工业界对于预测试验数据的各种模型和计算程序的兴趣。     

权分担三个集合的亚纯函数的唯一性

利用权分担集合的思想讨论了关于分担三个集合的亚纯函数的唯一性问题．证明了：设f与g是开平面上两个非常数亚纯函数,k：（i=1,2,3）为非负整数,n为不小于2的整数．若Ek1（{1,ω,ω^2,…,ω^n},f）=Ek1（{1,ω,ω^2,…,ω^n}g）Ek2（{0},f）=Ek2（{0},g）Ek3（{∞},f）=Ek3（{∞}g）且a,b,c,n满足（an-a-2）（bcn-b-f）〉2bcn,其中k1＋1=a,k2＋1=b,k3＋1=c,则f=tg（t^n=1）;或fg=s（s^n=1）,且0和∞为f与g的缺省值．