再入湍流尾迹及其对雷达散射的影响研究

就高超声速再入体尾迹等离子体场而言,为进行其亚密湍流雷达散射截面的理论分析,提出了计算湍流尾迹脉动等离子体场强的理论模型及求解方法,即在研究高超声速尾迹流动特征的基础上,推导、使用包括化学组份浓度脉动强度的k－ε－g湍流模型,用以封闭高超声速粘性尾迹湍流运动时均控制方程组,并用全隐式有限差分法求解,算例小钝锥体的飞行条件为零攻角、M_∞＝21．3,Re_（∞D）＝1．33×10~6;M_∞＝20．5;Re_（∞D）＝3．22×10~5,计算结果得到了合理的参数分布;考虑电磁波在上述等离子体介质中的传播,以单电子多次散射模型─—畸变波Born近似方法,计算了湍流尾迹脉动等离子体雷达散射截面,给出散射能（RCS）在不同极化状态下的分布,分析了散射背景场脉动湍流对电磁波在其中传播的影响．     

两个角区湍流场及其尾迹的实验研究

绕两个翼型－平面的角区流动及其尾迹的实验是在低湍流度风洞中完成的．在零攻角条件下，对翼型－平面的角区流场内诸参数，如翼型表面和平板面上的压力分布、绕翼型及尾迹区内的平均速度、脉动速度、湍动能、二阶关联量ｕ′ｖ′及ｕ′ｗ′进行了广泛的测量．通过对比，分析了这两种模型与平面所构成的角区及其尾迹区内的流动特性     

高超声速三维化学非平衡流动N－S方程数值解

详细地叙述了模拟三维高超声速电离空气化学非平衡粘性流全流场（前体和底部近尾迹）的数值方法，化学反应粘性流求解是基于带有化学源项的Ｎ－Ｓ方程的守恒形式，总的连续方程由组元守恒方程组所代替．对于高温电离空气流动，存在如下７个主要组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ和ｅ￣－．研究发展了求解完全耦合的，与时间相关的偏微分方程组的数值方法．利用一类新的迎风ＴＶＤ激波捕捉有限差分格式求解守恒形式的控制方程组，对高超声速层流有攻角钝锥体绕流流场，在非催化壁条件下、对包括底部近尾迹在内的全流场各参数，如组元浓度，电子密度和温度，以及物面压力分布和热流率分布得出了计算结果．对化学反应气体和完全气体模型的计算结果进行了比较，计算的流场电子密度与飞行实验结果符合很好．     

高超声速三维化学非平衡流动N－S方程数值解

详细地叙述了模拟三维高超声速电离空气化学非平衡粘性流全流场（前体和底部近尾迹）的数值方法，化学反应粘性流求解是基于带有化学源项的Ｎ－Ｓ方程的守恒形式，总的连续方程由组元守恒方程组所代替．对于高温电离空气流动，存在如下７个主要组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ和ｅ￣－．研究发展了求解完全耦合的，与时间相关的偏微分方程组的数值方法．利用一类新的迎风ＴＶＤ激波捕捉有限差分格式求解守恒形式的控制方程组，对高超声速层流有攻角钝锥体绕流流场，在非催化壁条件下、对包括底部近尾迹在内的全流场各参数，如组元浓度，电子密度和温度，以及物面压力分布和热流率分布得出了计算结果．对化学反应气体和完全气体模型的计算结果进行了比较，计算的流场电子密度与飞行实验结果符合很好．     

k-D-a-B 模型和 Richtmyer-Meshkov 不稳定性的数值模拟

给出一个适用于流体力学数值模拟的湍流混合模型。由于引入湍流质量通量和密度脉动关联量的演化方程，它比一般的ｋ－ε模型能更准确地描述界面不稳定性引起的可压缩湍流混合过程，并且可以不需要人为给定初始湍流场。我们用同一套参数模拟了不同Ａｔｗｏｏｄ数的激波管实验，数值结果与实验基本一致     

k-D-a-B 模型和 Richtmyer-Meshkov 不稳定性的数值模拟

给出一个适用于流体力学数值模拟的湍流混合模型。由于引入湍流质量通量和密度脉动关联量的演化方程，它比一般的ｋ－模型能更准确地描述界面不稳定性引起的可压缩湍流混合过程，并且可以不需要人为给定初始湍流场。我们用同一套参数模拟了不同Ａｔｗｏｏｄ数的激波管实验，数值结果与实验基本一致。     

FV/MC混合算法求解轴对称钝体后湍流流场

介绍一种有限容积／Monte Carlo结合求解湍流流场的相容的混合算法．有限容积法求解Reynolds平均的动量方程和能量方程，Monte Carlo方法求解模化的脉动速度—频率—标量联合的PDF方程．将该算法发展到无结构网格，探讨了在无结构网格中实现两种方法的耦合，包括颗粒定位，颗粒场和平均场之间数据交换等问题．并以二维轴对称钝体后湍流流场作为算例，比较了计算结果与实验结果．     

LES/RANS方法混合函数特性研究

将两方程k-ω SST湍流模型和Sagaut的混合尺度亚格子模型通过一个混合函数相结合, 构造出一种混合大涡/雷诺平均N-S方程模拟方法(hybird large eddy simulation/reynolds-averaged navier-stokes, Hybrid LES/RANS), 采用这种混合模拟方法结合5阶WENO格式对Ma=2.8平板湍流边界层进行了数值模拟, 并在计算区域上游入口处采用“回收/调节”方法生成湍流脉动边界条件, 通过考查RANS区域向LES区域的过渡参数及网格分辨率对这种混合模拟方法进行了评价. 计算结果表明: 该文采用的混合模拟方法可以捕捉到湍流边界层中的大尺度结构且入口边界层平均参数不会发生漂移, 混合函数应当将RANS区域和LES区域的过渡点设置在对数律层和尾迹律层的交界处, 而过渡应当迅速以获得正确的雷诺剪切应力分布, 在该文采用的模型及数值方法的条件下, 流向及展向的网格小至与Escudier混合长相当时, 能够获得可以接受的脉动速度的单点-二阶统计值.     

90°弯管内流动的理论模型及流动特性的数值研究

从三维不可压缩雷诺时均Navier-Stokes方程出发,对90°弯曲管道内湍流流动进行数值模拟。网格划分采用六面体网格,湍流模型为RNGk-ε模型,在近壁区采用两层壁面模型进行修正,流场的计算结果与实验数据吻合较好。在此基础上,本文数值研究了来流方向对流场结构和流动特性的影响。得出在弯管流场中发生了分离现象,且随着来流侧滑角的增大,分离区范围增大。此外,随着来流从同一侧滑角变换至同一攻角时,横截面的二次流图像中也从具有两个对称主涡变成只具有一个主涡的现象。     

高阶紧致格式求解二维粘性不可压缩复杂流场

提出了一种求解二维不可压缩复杂流场的高精度算法．控制方程为原始变量、压力Ｐｏｉｓｓｏｎ方程提法．在任意曲线坐标下，采用四阶紧致格式求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组，时间推进采用交替方向隐式（ＡＤＩ）格式，在非交错网格上用松弛法求解压力Ｐｏｉｓｓｏｎ方程．对于复杂的流场，采用了区域分解方法，并在每一时间步对各子域实施松弛迭代使之能精确地反映非定常流场．利用该算法计算了二维受驱空腔流动，弯管流动和垂直平板的突然起动问题．计算结果与实验结果和其他研究者的计算结果相比较吻合良好．对于平板起动流动，成功地模拟了流场中旋涡的生成以及Ｋａｒｍａｎ涡街的形成     

涡度拟能传输惯性区二维湍流标量场方差谱

本文应用非平衡统计力学封闭方法,给出涡度拟能传输惯性区二维湍流标量场方差谱的完整表达式,数值计算该表达式中的比例系数B。由于传输过程的非局部性,B依赖于表征波数变化范围的局部化因子,不再是普适常数。     

关于具有剪应力的普通湍流理论发展

本文叙述了周培源1945年发表具有剪应力的普通湍流理论以后,各种湍流模式理论的主要发展概况.      

Decomposition of Turbulent Velocity Fields in an SI Engine

In this study, the turbulence filter, the phase averaging and the proper orthogonal decomposition methods are used to decompose experimentally measured turbulent velocity fields in an SI engine. The radial and circumferential turbulent velocity fields were measured using hot wire anemometer under motored conditions at different engine configurations. The decomposed results of each technique are compared with each other. In addition, the obtained organized and turbulence motions and their energy spectra are examined. Finally, coherent structures of velocity fields and their activities are investigated. This revised version was published online in July 2006 with corrections to the Cover Date.     

走近湍流

简要回顾了20世纪湍流研究的进展,特别是近30多年湍流研究的新成就。对于湍流研究的若干热点问题提出了一些看法。     

剪切湍流的涡黏张量模式

对剪切湍流提出了涡黏系数为四阶张量的涡黏张量模式。引入近代数学中Ｍｏｏｒｅ－Ｐｅｎｒｏｓｅ广义逆矩阵的研究结果，给出了构造涡黏张量各分量的计算公式。用平面后台阶流动验算了剪切湍流的涡黏张量模式，比ＲＳＭ和ｋ－ε模式更接近实验结果。提出了剪切湍流涡黏张量模式的应用设想。     

在我国开展工程湍流研究的设想

本文包括三部分:１）提出开展工程湍流研究的背景;２）工程湍流研究的特征和内容;３）开展工程湍流研究的设想,其中包括:ａ）湍流理论的概念鉴别和在工程中的应用,ｂ）湍流模式理论的应用和完善,ｃ）湍流测量和ｄ）湍流控制．      

Properties of the hybrid RANS/LES filter

One of the most important and challenging topics in the Large Eddy Simulation of turbulent flows is the connection of the LES technique to the well known and largely used RANS approach where the Navier–Stokes equations are Reynolds averaged. The hybridation of LES and RANS is not only important for its possible practical use, (a rational use of the computational means in different zones), but also from a theoretical point of view, and one possible procedure consists of blending RANS and LES models in the transition zone. In this paper a new filtering technique based on blending filters which transitions smoothly between LES and RANS is proposed and the associated universal model for the subgrid scale stresses is derived. PACS 47.27.Eq     

Reynolds stress and vorticity in turbulent wall flows

Experimental measurements in a boundary layer and a large-eddy simulation of plane channel flow have been used to study the dynamics of vorticity and mass transport in the nearwall region. It was found that Reynolds stress generation occurs in the vicinity of quasi-streamwise vortices, and that smoke particles tend to be ejected from the wall near these vortical structures.     

剪切湍流大尺度相干结构的模式研究

发展了一种计算剪切湍流大尺度相干结构的新模式．该模式的基础是认为大尺度相干结构为湍流场中流体脉动能量增长最快的那部分，且包含大部分的湍流脉动能量．在此基础上。通过对湍流相干能量方程的推演。建立了描述大尺度相干结构的特征控制方程，并应用Ｃｈｅｂｙｓｈｅｖ多项式方法求得湍流相干能量的最大增长率在波数空间的分布，从而获得对应的大尺度相干结构．应用该模式研究了槽流和一自然对流中的大尺度相干结构，得到的近壁区流动结构与实验现象十分接近．