非定常风沙流中风速脉动对沙粒相瞬时浓度影响的实验研究

在大气边界层风洞内采用高速CMOS相机记录了时均风速线性上升的来流条件下运动沙粒的图像。根据二进制图像标记像素的连通性,从数字图像中提取相机拍摄区域内沙粒的数量N,由此计算得到风沙流的瞬时含沙量C。获取图像的同时,使用恒温式热线风速仪MiniCTA配合55R49热膜探头测量观测域上方主流位置和跃移层顶部附近的风速变化。实验结果表明,风沙流中的瞬时含沙量总是非稳定的,其变化频率至少在100Hz以上。风沙之间的相互作用存在一定的滞后效应,大气湍流脉动对含沙量的作用与沙粒粒径相关,粒径越小,风速脉动作用越大。当粒径≤160μm时,风速的微小变化即可引起含沙量的明显波动,二者高度相关。当粒径较大的时候,如300～500μm,只有较长时间较大幅度的风速脉动才能使含沙量有较大幅度的变化,研究时均风速变化对含沙量的影响更具有意义。     

水平槽道内气固两相湍流中颗粒行为的PIV实验研究

运用PIV技术对充分发展的水平槽道内稀疏湍流两相流中的颗粒行为进行了定量研究。实验中壁面雷诺数Rer=430,使用颗粒为直径60μm的聚乙烯微珠,测量了0．1％和0．5％两种质量荷载。研究发现,颗粒最大浓度出现在y＋≈10的位置上。对两种质量荷载,下壁面附近的颗粒体积分数最大值均远小于10^3量级,仍属于稀疏两相流。槽道下壁面附近,向上运动的颗粒概率大于向下运动的颗粒,两者概率的差异随着法向距离的增大而减小。在y’=10～30的范围内,颗粒相“喷射”事件（Q2）和“下扫”事件（Q4）概率分别显著地大于和小于当地流体对应值。颗粒相主要通过促进Q2且抑制Q4来实现对流体湍流的影响,颗粒的这一作用随着法向距离的增加而逐渐减弱。     

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率,对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序,该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想,验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合,说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动。为了提高计算规模及计算速度,作者对程序进行了并行化并针对微机机群系统进行了并行优化。实际计算表明本程序具有较高的并行效率。      

Phase coherence of foreshock MHD waves: wavelet analysis

The earth's foreshock is a region where particularly large amplitude MHD waves are commonly observed. They exhibit various waveforms, suggesting that nonlinear interaction between the waves is in progress. In a previous paper (Hada et al., 2003) we have introduced a method to quantitatively evaluate the strength of phase coherence among the waves from a given time series data. Here we further develop our method by applying wavelet filtering technique. From the analysis it was found that, although the turbulence is consisted of waves with a wide range of plasma rest frame frequencies, only those frequencies lower than the ion gyrofrequency are responsible for generating the phase coherence. This revised version was published online in August 2006 with corrections to the Cover Date.     

两方程湍流模型的可压缩性修正及其应用

本文在不可压两方程湍流模型的基础上进行可压缩性修正,给出了ｋ－ε、ｋ－ε两方程模型的可压缩性修正方程,并对亚声平板湍流流动和跨声喷管流动进行了数值模拟,得到了较好的结果。从而证实了两方程湍流模型经过可压缩性修正,能较好地模拟可压缩流动。     

数值模拟二维喷管激波/湍流附面层干扰流动

采用可压缩性修正两方程湍流模型,数值模拟了3种不同波前马赫数的跨声速二维喷管内激波/湍流附面层干扰流动,对流场中时均参数和脉动参数的计算结果与实验值进行了比较。结果表明可压缩性修正的两方程湍流模型准确地模拟了正激波/湍流附面层干扰流动的时均参数和脉动参数,无分离和有分离的激波/湍流附面层干扰流动的基本规律。      

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

带有双层壁扩压器的波瓣喷管混合流场的数值计算

采用三维贴体曲线坐标网格,边界网格加密且正交,在整个计算区域进行全场计算。为避免因波瓣造成的网格强烈的非正交引起解的发散,采用了Chen和Kim修正的k-ε湍流模型及同位网格SIMPLC计算方法,对带有双层壁扩压器的波瓣喷管混合流动进行了数值计算和分析。对波瓣及双层壁,采用大粘性的方法解决流固耦合。计算结果表明:在双层壁间有外界冷气流被引射进入,形成了壁面的冷却气流;自波瓣出口截面沿流向产生的环流速度场,强化了主次流的掺混,速度分布渐趋均匀。计算结果与实验数据符合良好,二者在离开波瓣52 mm的混合管内相对主流速度的最大误差为21.65 %。     

雷达天线罩廓形误差评定研究

 就弹载或飞行器上雷达天线罩的廓形测量及测量精度提高,探讨了应用鞍点规划的方法对天线罩廓形误差进行评定,给出了评定的线性鞍点规划模型,评定的最小条件,以及具体的操作处理方法,并介绍了评估例子,表明应用本文的评定方法获得了良好的结果。