不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

壁湍流调制雷诺应力和拟序结构应变率空间形态的层析PIV实验研究

采用高时间分辨率的层析PIV技术,测量了水洞壁湍流三维速度分量空间分布的时间序列,应用局部平均速度结构函数概念对流向速度信号进行多尺度分解,以流向速度分量的局部平均速度结构函数过零点作为特征量检测壁湍流中拟序结构猝发的喷射和扫掠过程,应用空间相位平均技术提取拟序结构猝发的喷射和扫掠过程各速度分量、涡量分量、拟序结构速度应变率分量以及调制雷诺应力分量的空间相位平均拓扑形态。为了引入平衡态涡粘模型假设模拟调制雷诺应力,研究了拟序结构猝发过程中调制雷诺应力分量和速度变形率分量的空间分布形态,发现两者之间的空间相位分布不一致。由于存在时空相位不同步性,说明需要考虑大尺度拟序涡结构引起动量传递的时空弛豫效应。应用经典的线性平衡态下的Boussinesq涡粘模型不能准确地描述壁湍流拟序结构动量传递非平衡现象的物理机理。对于壁湍流拟序结构动力学方程中调制雷诺应力的模拟,应采用包含时空相位信息的复涡粘张量模型。由于雷诺应力与速度变形率的时空不同步性,对非平衡非局部湍流场的数值模拟提出了一个挑战性的问题,建议采用包含相位信息的复涡粘张量模型来模拟雷诺应力张量,从而更加符合雷诺应力演化的物理机理,这一模型有可能成为一个很有发展前景的封闭模型,从而更加准确地预测工业领域中广泛存在的非平衡湍流。     

应用修正湍流模型对旋转通道内换热的模拟

针对旋转光滑矩形通道分别应用针对旋转状态修正的 k-ε、标准k-ω以及提出的针对旋转状态修正的k-ω湍流模型进行流动和换热的数值模拟,通过与实验结论的对比,讨论了采用不同湍流模型对计算结果的影响.计算工况为旋转数Ro=0.24,流体进口雷诺数Re=25000.计算结果表明:采用所提出的针对旋转状态修正的k-ω湍流模型的计算结果要比采用针对旋转状态修正的k-ε以及标准k-ω湍流模型的计算结果更接近实验结论.     

空气龄的数值解法

提出了用ＣＦＤ法求解空气龄的具体方法，其步骤为先求解Ｎ－Ｓ方程，后求解空气龄。为求解Ｎ－Ｓ方程，采用了改进型ＨＳＭＡＣ法，并给出了压力和速度同时修正的方法。为求解空气龄，推出了空气龄输送方程式和给出了空气龄和速度成分的边界条件。计算机模拟验证了求解方法的有效性并由空气龄求出了通风效率。     

激光通信中的分集检测技术

在大气弱湍流模型下,应用Monte Carlo方法,通过产生符合光信号的对数幅度联合概率分布的数据,近似模拟计算了分集接收机的逐符号最大似然检测算法的似然比.研究了减少Monte Carlo方法中仿真点的个数对逐符号最大似然检测的分集接收机检测性能估计的影响.随后提出一种基于训练序列的分集检测方法.该方法在未知衰落分布相关信息的情况下,利用每帧内引入的一定长度的训练序列提供了某一较短时间间隔内的大气湍流信道的信息,基于与似然比表达式结构相似的算法获得改善的检测效果.仿真结果表明在各个分集接收机之间存在较强相关性的情况下,该算法仍可实现有效的检测.      

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

Bursty Bulk Flow Driven Turbulence in the Earth’s Plasma Sheet

During the last several years significant progress has been made in understanding MHD turbulence in the Earth’s plasma sheet. Due to the statistically transitory properties of fluctuations, finite size and boundary effects, however, issues of fundamental importance remain unresolved. Here we concentrate on such intrinsic features of plasma sheet turbulence as its origin and dynamical nature. In particular, we investigate bursty bulk flow driven multi-scale transfer of energy towards the dissipation scale, and provide evidence for the presence of non-linear interactions. We show that, in contrast with previous results, Alfvénic fluctuations together with 2D eddy interactions may appear as important constituents of turbulence in the plasma sheet.     

Alfvén waves in the context of solar-like star formation: accretion columns and disks

In this work we examine the damping of Alfvén waves as a source of plasma heating in disks and magnetic funnels of young solar like stars, the T Tauri stars. We apply four different damping mechanisms in this study: viscous-resistive, collisional, nonlinear and turbulent, exploring a wide range of wave frequencies, from 10⁻⁵Ω_i to 10⁻¹Ω_i (where Ω_i is the ion-cyclotron frequency). The results show that Alfvénic heating can increase the ionization rate of accretion disks and elevate the temperature of magnetic funnels of T Tauri stars opening possibilities to explain some observational features of these objects. This revised version was published online in August 2006 with corrections to the Cover Date.     

涡扇发动机进气总压畸变试验动态总压紊流度计算分析

通过对某型涡扇发动机进气总压畸变流场动态总压紊流度的计算，得到了在改变采样频率以及设计不同的数字滤波器对动态总压紊流度计算的影响结果，为分析和处理动态数据提供了一种方法。