结构网格方法对高升力构型的应用研究

采用"亚跨超CFD软件平台"(TRIP)数值模拟了梯形翼全展长与半展长高升力构型的复杂流场,并与试验做了对比分析。对应风洞试验是在NASA Langley 14×22英尺亚声速风洞(FST)和NASA Ames 12英尺增压风洞(PWT)中完成的。计算中采用一方程SA湍流模型和MUSCL-ROE格式,并综合运用对接/重叠/拼接网格技术,数值模拟了两种高升力构型的气动特性,给出了典型站位的压力分布,并对比研究了不同结构网格技术对此类高升力构型的计算差异。研究表明,与修正后的试验数据相比,数值模拟得到的气动力系数和典型剖面的压力分布均与试验结果吻合良好。     

DLR-F4翼身组合体的阻力计算

为了考察自行研发的CFD软件的计算能力和阻力计算精度，本文采用LU—SGS方法、MUSCL差分格式和Baldwin—Lomax代数湍流模型，数值模拟了AIAA阻力计算工作室提供的DLR—F4翼身组合体的绕流流场，综合分析了easel和case2的气动力的计算结果，并与NASA Christopher L．Rumsey采用CFL3D6．0和AFRL/VAAC Don W．Kinsey采用Cobalt60提供的两组计算结果以及AGARD提供的两种不同风洞的测力试验结果作了比较。计算结果表明，本文计算精度与国外CFD软件相当。为了提高激波，边界层干扰的模拟精度，今后要重点加强湍流模型的应用研究。     

高超声速飞行器钝舵缝隙流动数值模拟研究

目前,带缝隙钝舵的缝隙引起的流场结构和气动加热规律变化,还很不明确,需要研究缝隙诱导所形成旋涡的空间分布特征和旋涡运动对物面气动加热的影响规律。通过分析高超声速钝舵缝隙气动加热问题,基于无缝隙钝舵,建立一种带缝隙钝舵简化模型。使用有限体积方法求解可压缩Navier-Stokes方程,通量采用van Leer通量向量分裂方法计算。插值采用MUSCL方法,时间项采用LU-SGS隐式方法。结果表明:无缝隙钝舵流场结构相对简单,带缝隙钝舵流场结构同无缝隙钝舵相比要更为复杂,舵轴上游缝隙内会出现马蹄形涡串结构,相应地在缝隙的上下表面均会出现马蹄形高热流区;受缝隙诱导分离再附流动的影响,在舵轴迎风面以及舵体侧面后部均形成了局部高热流区。     

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

不同襟翼偏角梯形翼构型气动特性数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用二阶离散精度的单调迎风格式(MUSCL),结合剪切应力输运(SST)两方程湍流模型和γ-Reθ转捩模型,研究了梯形翼高升力构型襟翼偏角变化对气动特性的影响。主要目的是进一步确认Trisonic Platform(TRIP)软件模拟高升力梯形翼不同襟翼偏角引起微小气动特性变化的能力。首先,简要介绍了采用的计算方法;其次,介绍了两种襟翼偏角的梯形翼模型及风洞试验;最后,在网格收敛性研究的基础上,采用全湍流和转捩两种方式模拟了梯形翼构型不同襟翼偏角对气动特性的影响。与试验数据的对比结果表明,采用全湍流和转捩两种方式均可以较好地模拟不同襟翼偏角对气动特性的影响量,采用γ-Reθ转捩模拟方式可以提高梯形翼构型气动特性的模拟精度,失速迎角附近的气动特性模拟需要进一步研究。     

交互式棱柱网格生成中翘曲现象形成机制及消除算法

在基于径向基函数插值的交互式棱柱网格生成中,当物面网格单元空间法矢变化范围较大或变化剧烈时,棱柱网格生成会发生翘曲问题,降低棱柱网格质量。针对该问题,首先采用DLR-F6机翼外形算例对翘曲现象进行了介绍和初步分析,并结合平板外形棱柱网格生成,通过不同波长推进位移扰动模拟法矢变化对推进位移的影响,研究棱柱网格翘曲现象的形成机制。不同位移扰动下平板棱柱网格生成结果证实边界棱线网格点推进位移不光滑会通过空间插值放大,引起棱柱网格发生翘曲,且网格翘曲的表现形式与推进位移的分布有关。从而间接表明棱柱网格生成中翘曲现象是由法矢变化带来推进位移分布不光滑引起的。然后从翘曲现象的形成机制出发,基于拉普拉斯方法发展了推进位移光顺技术,并通过DLR-F6机翼算例进行了验证。测试结果表明,光顺技术能够完全消除棱柱网格的翘曲,获得光顺的高质量棱柱网格。最后通过DLR-F6翼身组合体外形棱柱网格生成,展示了添加光顺技术的棱柱网格生成算法的应用潜力。     

无尾布局支撑干扰数值模拟

采用TRIP2.0_SOLVER软件,开展了高、低速来流条件下,无尾飞翼布局支撑干扰的数值模拟技术研究.高速来流状态下支撑干扰的数值模拟,以尾支撑作为主支撑,腹支撑作为辅助支撑,完全模拟了实验中采用的两步法;低速来流状态下的支撑干扰的数值模拟,采用了有无尾支撑的方式得到支撑干扰量.笔者介绍了支撑干扰的典型高、低速数值模拟结果,并与高速实验结果做了初步的对比.     

多重拼接网格技术应用研究

利用多重拼接网格技术对三种不同类型的计算外形进行了内外流场的数值模拟,并将计算结果与单重对接和单重拼接的结果进行了对比,计算结果表明多重拼接网格技术在没有降低计算精度的情况下,大大简化了网格生成,同时可以把计算速度提高一个量级左右,取得了令人满意的效果。     

无尾布局支撑干扰数值模拟

采用TRIP2．0．SOLVER软件，开展了高、低速来流条件下，无尾飞翼布局支撑干扰的数值模拟技术研究。高速来流状态下支撑干扰的数值模拟，以尾支撑作为主支撑，腹支撑作为辅助支撑，完全模拟了实验中采用的两步法；低速来流状态下的支撑干扰的数值模拟，采用了有无尾支撑的方式得到支撑干扰量。介绍了支撑干扰的典型高、低速数值模拟结果，并与高速实验结果做了初步的对比。     

半展长襟翼梯形翼构型数值模拟技术研究

基于雷诺平均的 Navier-Stokes 方程和拼接结构网格技术,采用 MUSCL 格式和 SST 湍流模型,研究了网格密度对半展长襟翼梯形翼高升力构型的数值模拟结果的影响。相应的风洞试验是1998年在 NASA Ames 12英尺增压风洞(PWT)中完成的,试验结果包括了总体气动特性、压力分布。研究内容主要包括网格密度对收敛历程、气动力特性、压力分布和表面流线的影响,以及气动力特性随迎角的变化。研究表明,Ma=0.15,α=16.7°时,网格密度对收敛历程、典型站位压力分布和表面流态基本没有影响,气动力特性随网格密度单调变化;采用不同密度的网格,典型剖面的压力分布与试验结果吻合良好;与修正后的试验数据相比较,数值模拟得到的失速迎角前的气动力系数与试验结果吻合良好。