燃气源复杂通道三维冷流数值模拟

采用两方程RNGκ-ε湍流模型，通过非结构网格上的SIMPLE算法，对某远程火箭弹控制系统燃气源复杂通道内流场进行了冷流数值模拟，获得了通道内各气流参量的三维分布，了解了燃气源内各处特别是控制气流流经的迷宫形多层导热套环结构内的气流流动特性，计算结果与实验结果基本吻合。     

近壁湍流相干结构的数值模拟

提出了近壁湍流对称单个相干结构的理论模型,采用紧致有限差分和Fourier谱展开相结合的数值方法研究它的演化问题,得到雷诺应力分布及高剪切层的形成过程与实验结果一致。     

埋入式进气道流场的雷诺应力测量和频谱分析

本文给出了埋入式进气道内湍流流场的雷诺应力测量和频谱分析结果。在研究中利用“Ｘ”双斜丝探头测量了低速下埋入式进气道出口截面上六个雷诺应力分量沿径向和周向的分布,继而把两个单丝热线探头分别置于进气道的进、出口,通过ＣＦ－９２０频谱仪分析两个动态信号的自功率谱和互相关性质。结果表明埋入式进气道出口近壁区气流的雷诺应力呈明显的各向异性,其上游的气流脉动是按气流平均速率向下游传输的。     

应用FFT高精度快速计算瞬态响应

为了提高瞬态响应的计算速度和精度,构造了基于快速傅立叶变换的响应计算格式,并研究了这种格式的本质。理论分析和数值仿真表明：ａ所提的计算格式能够提高速度８ 倍;ｂ基于抛物线型求积公式的算法精度最高;ｃ只要频域采样点数 Ｍ 超过时域采样点数２ 倍,计算精度就与 Ｍ无关。     

模型的转捩流动计算分析简

 为了在黏性流动数值模拟中实现边界层转捩的自动预测,将γ-Re_θt转捩模型引入到三维非结构混合网格的雷诺平均Navier-Stokes方程求解程序（HUNS3D）。该转捩模型由两个依赖当地变量定义的关于间歇因子和当地化转捩起始动量厚度雷诺数的输运方程组成,其数值求解算法与流场求解程序中湍流模型的求解方法相同。为了考察和验证HUNS3D程序中γ-Re_θt转捩模型对航空工程中的常见附面层自由转捩问题的预测精度,对低速平板流动、Aerospatial-A翼型、NLR 7301超临界翼型和NASA Trap wing 高升力构型等典型外形的自由转捩流动进行了计算,并将计算结果与相关试验结果进行了对比分析。算例结果表明：γ-Re_θt转捩模型对于转捩位置具有很好的敏感性,能比较准确地预测自然转捩和分离转捩,可以有效提高HUNS3D程序对实际流动的模拟能力和预测精度。     

不同湍流模型对插板式进气畸变的数值模拟

为得出一种较为准确的湍流模型,应用三种湍流模型对插板式进气畸变进行数值模拟,通过与畸变实验实测数据的对比,研究不同湍流模型在插板式进气畸变计算中的预测能力。不同湍流模型计算出的出口流场畸变不同,重整化群k-ε模型得出高压区范围较大,S-A模型过渡区域范围大。结果表明:插板高度为畸变发生器内径的48%时,低压区范围角与实验值较为吻合,插板高度为畸变发生器内径的20%时,低压区范围角偏小。畸变指数的计算值以可实现化k-ε模型最为准确。     

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

结构误差对旋转稳定弹丸气动特性影响的数值模拟

基于剪切应力传输(SST)k-ω湍流模型,对考虑结构误差情况下的弹箭模型进行数值模拟计算。首先采用CFD和工程经验公式相结合的方法,得到了不同马赫数、不同转速情况下无结构误差模型的多种气动特性参数。通过与实验数据对比,阻力、法向力、俯仰力矩、压心位置误差在10%以内,验证了该方法的可行性和准确性。然后,建立了考虑结构误差即质量分布不对称、弹体不同轴和无结构误差模型,并分别进行了气动特性模拟计算。结果表明,质量偏心对滚转阻尼力矩系数、马格努斯力矩系数影响很大;弹体不同轴对法向力系数导数、俯仰力矩系数导数、马格努斯力矩系数导数和压心位置都产生很大影响。结构误差使模型的多种气动特性参数产生很大变化,将影响其飞行弹道和稳定性。     

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

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