用多重网格法模拟三维流动

采用多重网格法模拟三维流动。控制方程在三维直角或圆柱坐标系中通过有限体积方法进行离散 ,SIMPLE算法对速度与压力进行耦合 ,亚松弛策略提高离散方程的收敛速度。在压力方程的求解过程中 ,界面速度应用由Rhie、Chow和Majumdar等人提出的动量方程插值法来剔除不合理的压力波动。三维流动算例表明采用多重网格法模拟三维流动可以得到良好的计算精度和收敛性。     

激波风洞气动力试验不确定度影响因素分析

参考北大西洋公约组织和AIAA推荐的风洞试验数据不确定度计算方法,结合激波风洞运行特点,确定激波风洞气动力试验的主要误差源,计算激波风洞13-2标模气动力测量结果的不确定度。采用改变单一变量的方法计算主要误差源对测量结果不确定度的影响程度,辨析对不确定度起主要作用的基本参数。计算结果表明：皮托压力和总压的测量结果对流场参数影响显著,皮托压力的测量结果比总压测量结果对流场参数与气动力测量结果影响更大;降低皮托压力和总压的偏离极限,有利于提高激波风洞气动力试验数据的质量。     

JF—10氢氧爆轰驱动激波风洞自由流的测量和诊断

针对ＪＦ－１０氢氧爆轰驱动风洞动行的高焓高速状态，采用吸收光谱技术、皮托压力瞬态测量技术、电离探针等的测试方法，测量了自由流中一氧化氮的含量和离子浓度、皮托压力及其沿喷管出口半径方向的分布。并根据以上各项测量结果判断了自由流的非平衡程度。     

注塑模熔体非牛顿流动的数值模拟

研究注塑模流动的数值模拟。用有限元-有限差分方法模拟非牛顿流体的压力场和温度场,用流体网格法（FAN）模拟熔体锋面的推进过程,这种方法具有通用性好和容易实现自动推进的优点。在486微机上对几个算例进行了模拟计算。比较了牛顿和非牛顿流体的流动特点和锋面推进自动处理的效果,得到了满意的结果。     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

基于DOE优化湍流模型的SERN分离点预测数值模拟与实验验证

在SERN(Simple expansion ramp nozzle)中,流动分离现象伴随着激波/边界层相互干扰、强剪切层、强激波等复杂的流动现象,导致通用的计算流体力学软件对SERN中分离点、分离区和再附点的计算精度不能满足工程需要。为了提高SERN中分离点的计算精度,基于Allamaprabhu等的研究成果,针对Menter的SST(Shear stress transport)模型,将修改湍流模型经验参数的方法运用到SERN分离流动的预测中,并利用试验设计(Design of experiment,DOE)优化方法得到了经验参数的最优组合,使计算和实验的无量纲分离点和无量纲分离点压力的误差分别降至1.99%,4.38%,压力分布均方根误差降至7.83%。     

基于指数插值的浸没边界法在可压缩流模拟中的应用研究

由于浸没边界法在处理复杂外形和动态问题时具有极强的灵活性和便利性,因而受到广泛的关注。在当前浸没边界法的研究和应用中,主要工作集中在关于中/低雷诺数和不可压流动的模拟。本文引入 Choi 等提出的指数插值法对目标点速度进行插值,并根据多项式插值和 Crocce-Busemann 关系式等确定目标点压力和温度,并借鉴直角网格法中的虚拟网格法计算物面内的流动变量,建立了适用于可压缩流动的浸没边界法。通过亚、跨声速段翼型绕流和超声速圆柱绕流算例对方法进行验证,结果显示在 k 值取值合适时模拟结果与实验数据较为吻合。通过对相关算例网格收敛性分析,证明了本文发展的浸没边界法具有较高的精度。文中各算例的模拟结果表明,此方法能准确捕捉亚、跨、超声速绕流问题中的各类流动结构,适用于可压缩流动的模拟。     

二等标准皮托静压管系数检定结果的不确定度评定

基于标准皮托静压管的工作原理,以标准皮托静压管系数检定结果和测量标准为研究对象,分析了二等标准皮托静压管系数检定时的测量不确定度来源,对其进行了测量不确定度评定。     

二维贴体曲线座标系的快速生成——多层网格法的应用

本文简要地叙述了多层网格法,并用它求解了一对非线性的二维贴体曲线座标系的生成方程。讨论了应用中的几个问题。本文将多层网格法与G-S方法以及带有最佳松弛因子的SOR方法作了比较,结果表明,多层网格法十分有效。     

多重网格法求解火焰稳定器钝头体管流问题

多重网格法是近年发展起来的一种软加速方法,用多重网格法可明显地加快数值计算的收敛速度,且网格越密其加速效果越明显。应用多重网格法求解有钝头体的火焰稳定器管内流动问题,计算结果与试验结果比较,吻合较好。     

部分膜弹流的膜厚与二次压力峰公式

用部分膜弹流模型来模拟线接触的混合润滑,提出浮动网格法,详细计算了二次压力峰的位置与大小,最小油膜厚度与中心油膜厚度。首次提出了线接触部分膜弹流的油膜厚度与二次压力峰计算公式。对浮动网格法,用 6种方法,即加、减、乘、除、对数与幂,进行了说明。理论分析的数学模型是基于Patir与 Cheng的平均流模型与微凸体接触模型。     

跨音速翼型计算的快速方法

用全速势方程求解绕翼型的跨音速流动时常用的计算方法有:SLOR方法、BGKJ方法(快速直接/SLOR混合方法)、隐式近似因子分解(AF)方法和多重网格法等。本文用其中的两种方法:BGKJ方法和AF2方法,以NACA 0012翼型为例进行了数值实验,比较了这两种方法的收敛速度。计算时,BGKJ方法用现有的计算程序,AF2方法则用最近自己编制的程序。数值实验的结果表明,BGKJ方法较AF2方法收敛速度快。     

静气弹中非线性气动力建模方法与分析

大展弦比低雷诺数气动布局容易较早出现气流分离,会带来明显的非线性气动力问题。针对此类布局提出了一种建立基于Kriging插值的非线性压力系数分布模型的方法。从Navier-Stokes方程计算的不同状态下飞机的压力系数中提取不同坐标的系数。利用Kriging插值函数建立CFD压力系数对迎角导数的响应面,将插值结果代入偶极子网格法(Double-Lattice Method,DLM)修正其线性方法。利用无限板样条(IPS)方法进行气动结构耦合,实现了有限元结构的非线性气弹响应分析。算例结果验证了方法对于静气弹分析的有效性,同时能准确地反映弹性带来的气动效率的降低和非线性力矩特征。     

绕任意刚性翼型,机翼非定常流动的跨声速欧拉方程解

本文在有限体积方法的基础上引入变系数隐式残值光顺及固接网格法求解了绕二维翼型及三维机翼非定常流动的跨声速欧拉方程,计算结果与同类结果及实验结果吻合。     

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

低超声速来流中机翼跨声速非定常气动力数值解法

本文给出一种低超声速来流中机翼跨声速非定常流的计算方法。该方法使用跨声速小扰动理论,正确地计入了沿流向平面流动特征线的压力波的运动边界条件。流动方程是使用时间精确交替方向隐式(ADI)有限差分解法求解。试算了M=1.1的矩形翼和M=1.1、1.34的F5机翼的定常和非定常气动力,并与国外的风洞试验和计算结果作了比较,符合很好。表明本方法是有效的。     

一种高超声速二元混压式进气道的研究

 针对飞行马赫数为6.00的二元进气道模型开展了高焓脉冲风洞试验研究,分析了进气道在不设置反压和设置反压两种情况下的激波结构、内通道皮托压分布及隔离段出口的性能,并结合数值仿真分析了通道内的流场特性。研究结果表明:在无反压情况下,进气道内通道激波反射明显,靠近下壁面的皮托压值均低于其他测点,在隔离段出口截面,靠近侧壁皮托压有所降低;在一定反压条件下,结尾激波系上传至隔离段内,结尾激波位置不对称;堵塞度为62%的反压条件下,结尾激波系位于喉道位置,隔离段出口截面下半部分已经是亚声速流动;在来流马赫数Ma=6.07,迎角α=4.5°无反压情况下,隔离段出口总压恢复系数为0.477,平均马赫数为2.72,增压比为44,流量系数为0.81,表明进气道性能良好。     

重合网格法在复合材料分层问题中的应用

使用重合网格法(S-FEM)分析了混合模态下CFRP层合板分层问题。首先将虚裂纹闭合技术集成到S-FEM程序中,求解了复合材料分层尖端处的应变能释放率。然后把重合网格法数值结果与基于梁理论的计算结果及实验数据进行了比较。结果表明,对于单方向及多方向CFRP复合材料层合板,重合网格法都能够得到更精确的数值,是一种有效的计算复合材料裂纹尖端应变能释放率的方法。     

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。