薄片分离过程流场特性的数值仿真研究

利用弹簧近似法和局部网格重构相结合的非结构动网格方法,耦合求解Navier-Stokes方程和六自由度刚体动力学方程,数值模拟了三维薄片非定常分离过程。数值仿真结果有效反映出薄片分离过程中运动速度、轨迹、受力变化规律及流场特性的发展与变化。根据数值仿真结果对薄片分离机理、多体气动干扰作用进行了分析。分离过程中单个薄片与主体段周围形成复杂流动,并伴随着尾迹涡的运动、发展等非定常过程。气动载荷的非线性变化是薄片分离的主要原因。     

非结构动网格技术及其在超声速飞行器头罩分离模拟中的应用

针对激波诱导物体运动的二维动力学模型,比较了采用流动定常、结合非定常弹道方程的所谓"准定常"计算方法和非结构动网格计算方法,通过流场和运动特性分析,指出发展有相对运动的非定常流动模拟方法的重要性;另外采用修正弹簧近似和局部重构的非结构动网格技术,耦合求解Euler方程及6DOF弹道方程,数值模拟高度H=25km、飞行马赫数Ma=6.0、迎角α=0°状态下火箭抛罩的动力学过程.研究表明,稠密大气层内头罩受到的气动力比惯性力大几个数量级,导致头罩的运动特性完全由气动力控制.分离过程中,头罩和火箭头部之间的气动干扰在时间和空间上都表现出强烈的非线性特征,头罩姿态变化剧烈.本文还对分离过程中部分位置上的流场进行定常模拟,比较了定常和非定常的气动力,所得结论对于分析地面风洞可控轨迹系统(CTS)试验具有指导意义.     

基于非结构动网格的非定常激波装配法

从ALE方程出发,通过修正体积的新思路构建离散几何守恒律,基于非结构动网格的信息传递方法得到网格重构以后新网格的流场参数,采用一阶格式捕捉法得到的流场确定装配法的初始激波位置,成功地把激波装配法应用到非定常流动的模拟。模拟超声速飞行环境下整流罩分离的非定常流固耦合过程,在马赫数6条件下装配法和捕捉法的计算结果完全一致,随着马赫数增加,二阶精度格式的捕捉法需要进行调整限制器参数、时间步长等人工干预,但是装配法直到马赫数20计算过程没有出现异常。该基于非结构动网格技术的激波装配方法,与传统的激波装配方法相比,能够用于复杂外形产生的非规则形状激波和非定常流动产生的运动激波,解决二阶精度有限体积法用于高马赫数飞行条件下多体分离过程的流动模拟时遇到的稳定性问题。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

带反馈控制的基本带翼外形的耦合气动/运动数值模拟研究

介绍了一种先进的、时间精度的多学科数值模拟平台的开发和应用,该平台可用于预测飞行器机动过程中的非定常气动特性和飞行力学特性.通过耦合求解非定常Navier-Stokes方程、刚体动力学方程及控制律,对基本带翼外形进行了虚拟飞行数值模拟研究.求解器采用基于结构网格的有限体积法,时间方法采用双时间步方法.使用反馈控制器获得舵偏控制指令,重叠网格方法用来模拟物体间的相对运动.模拟结果显示控制器的增益必须正确选择才能获得正确的响应.仿真结果表明本文发展的虚拟飞行数值模拟方法能够预测包含动边界的复杂非定常流动,具有较强的工程适用性.     

安装角度对内埋式导弹分离特性的影响

基于弹簧变形理论和局部重构的非结构动网格技术,采用有限体积法求解ALE描述下的三维非定常Euler方程,并耦合刚体六自由度运动方程,模拟了超声速条件下三种不同安装状态导弹从弹舱中弹射分离的过程,分析了导弹质心相对位移以及姿态变化情况,发现增大初始安装角度能够改善导弹在超声速环境下的分离特性。     

外挂物投放问题的二维非结构动态网格技术研究

以动态网格守恒型方程为控制方程,在二维非结构网格基础上发展了弹簧式光顺法和局部网格重构法两种动态网格变形技术。结合运动轨迹计算方程,将动态网格法应用于飞机外挂物投放问题中,成功模拟了导弹自由投放的整个过程,获得了详细的包括弹体运动速度、运动轨迹和受力情况等在内的丰富的流场信息。成功的应用不仅验证了动态网格法的可行性;而且也表明,动态网格法的进一步完善与发展将成为解决外挂物投放问题的有力工具。     

扑翼微飞行器飞行时的数值模拟与分析

采用数值方法求解三维非定常不可压缩的N-S方程,生成了相应的扑翼微飞行器的三维动网格;模拟了扑翼微飞行器在不同频率、不同攻角、不同振幅、不同飞行速度、不同的翅膀初始位置下的非稳态流场;计算了扑翼微飞行器在不同运动模态下的升阻力系数;分析了影响扑翼微飞行器空气动力学的特性的不同因素;为扑翼微飞行器的飞行机制和建造提供了理论依据。     

基于动网格的非定常粘性空气动力学计算

提出了一种将阵面推进法和半扭转弹簧法相结合的动网格生成方法。运用阵面推进法决定附面层动网格的运动规律,半扭转弹簧法生成附面层外的非结构动网格。使得算法能够同时适用于非定常粘性与非粘性空气动力学计算。文中对于非定常可压缩N-S方程,在非结构动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代,计算了二维NACA0012翼型俯仰简谐运动的低雷诺数非定常空气动力,计算结果与参考文献提供的结果吻合良好。此外,本文还对三维1/4椭圆薄板的低雷诺数非定常空气动力进行了计算。     

多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

高速运动弹丸诱导的瞬态流场数值研究

采用基于任意拉格朗日欧拉(ALE)方程的动网格及嵌入网格法,结合二阶精度Roe格式和沉浸边界法,对弹丸由高压气体驱动从静止状态加速至超声速,射出膛口到完全飞离初始流场的整个过程进行了数值模拟.根据数值结果,详细讨论了初始流场、火药燃气流场的形成与发展以及与弹丸的耦合和相互作用过程.      

Enhancement on parallel unstructured overset grid method for complex aerospace engineering applications

In the present study, an efficient overset grid method by means of parallel implicit hole-cutting is proposed for the sake of simulating unsteady flows in aerospace engineering involving multiple bodies in relative movement. In view of the degraded computational efficiency and robustness for conventional overset grid assembly, several innovative techniques are developed within the overset grid assembly process, viz., a bookkeeping alternative digital tree method to speed up the donor-cell search...     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

聚合多重网格法在旋翼前飞流场计算中的应用

基于非结构动态嵌套网格, 将聚合多重网格方法引入到旋翼的前飞流场计算中.从初始的整体非结构网格出发, 快速建立相互重叠的两个子域网格, 同时无需建立背景网格便可以确定初始洞边界和初始插值贡献单元, 且动态嵌套时无需更新洞边界, 并采用了一种无需背景网格便可快速查找新的贡献单元的方法.分区利用聚合法生成多重网格, 利用双时间方法分区运用三层网格V循环求解非定常Euler方程, 两区的流场参数在重叠区实现耦合, 对旋翼的前飞流场进行了计算.结果表明, 多重网格方法能够有效地加快旋翼前飞流场收敛.      

螺旋桨非定常滑流数值模拟

运用非结构动态嵌套网格技术,通过求解三维非定常Euler方程,对螺旋桨飞机的螺旋桨滑流进行了真实数值模拟.动态嵌套网格方法通过建立运动区子域网格来模拟螺旋桨的旋转运动,静止的机身包含在静止区网格中,有效解决了螺旋桨滑流模拟中的相对运动问题.采用双时间方法求解控制方程,在每一个物理时间步上运用经典四步Runge-Kutt...     

基于动态边界控制的优化方法在高升阻比气动外形优化中的初步应用

本文是利用动网格技术和有限体积法来研究基于动态边界控制的气动优化的方法.该优化方法是通过在求解非定常流动的过程中计算优化的目标函数随设计变量变化的梯度.而利用参数化的结构动网格技术能有效、快速地生成与优化过程出现的外形对应的网格并进行计算,具有较高的效率.本文以典型的高升阻比外形优化为例,用非定常NS方程为流动控制方程进行优化,取得了初步的结果.     

飞机着陆过程中非定常地面效应数值模拟方法

针对飞机着陆过程的大位移运动,提出适于多块结构动网格生成的改进谢别德-无限插值混合方法,改进谢别德方法准确计算各网格块的边界网格变形,无限插值方法快速计算内部网格点变形,结合局部重构保证动态网格的合理分布.混合方法有效解决了单一无限插值法存在的动网格交错问题.进而将二维非定常地面效应N-S方程计算方法发展至三维情形,应用于带增升装置的翼身组合体着陆过程中非定常地面效应计算分析.结果表明:模型的升力系数随离地高度的减小而减小,非定常和准定常计算方法的结果存在一定差异,而考虑了非定常效应的算法更为准确.     

基于动态非结构重叠网格法的直升机前飞非定常流场数值模拟研究

 发展了非结构重叠网格技术，提出适合于直升机旋翼运动的动态非结构重叠网格算法，该算法使得在采用非结构重叠网格模拟直升机流场时实现程序自动处理网格，不需任何人工干预。流动控制方程采用非定常可压缩Euler方程，空间离散采用Jameson有限体积法，时间推进采用双时间步长法，为加速收敛，在双时间步长内循环迭代中采用预处理技术。为验证本发展的预处理解算器以及动态非结构重叠网格的算法的可信性和有效性，开展一系列的数值实验，模拟了不同速度范围的NACA0012翼型绕流,Caradonna旋翼悬停和前飞流场以及Robin直升机机身-旋翼组合的前飞非定常流场，结果表明采用本文算法得到的结果是可信的，网格处理是自动高效的，可以应用于直升机全机流场模拟，为直升机设计提供可信指导数据。