耦合飞行器运动的非定常粘性流场的数值模拟.pdfVIP

第十届全国计算流体力学舍议论文集 耦合飞行器运动的非定常粘性流场的数值模拟+ 袁先旭张涵信叶友达 国家计算流体力学实验室，北京100083 中国空气动力研究与发展中心，四川绵阳621000 捕要本文基于渭合通量分疆的思想，通过应用高斯一奏德尔迭代求解差分方程，构造了一种 时空二阶棒度的隐式选代NND算法，并特别讨论了时间精度．运动同格生成技术，几何守恒律的 应用以及动壁边界条件婷与非定常数值模拟密切相关的问题．作为检验算例，文中计算了二维气 动弹性腰柱层流脱涡绕流，并给出一个高超声速钝件动态绕流的结果。 荚t饲高斯一赛德尔选代方法，NND格式，混台通量分裂。 0弓I言 现代飞行器的设计日益追求超机动性和敏捷性．气动弹性干扰特性、特别是在跨声速范 围具有强烈非线性、非定常特征的气动弹性干扰特性在飞行器的设计中起着越来越重要的作 用，使得非定常流动的研究日益成为人们关注的重要课题之一。随着计算机能力的提高和计 算流体力学中的数值方法的很大发展，在气动设计中开展非定常弼方程的数值模拟也逐渐具 有了可行性．在气动弹性领域，一个具有拂战性的目标就是较精确地预测耦合流体一结构干 扰效应的非定常复杂流场。一般说来，数值模拟非定场粘性流场，采用时间高阶精度的显式 算法就会遇到小时间步长和过大计算量的问题，因此发展适于时间精确计算的陡式方法就有 较大意义．九十年代以来，国外相继发展了单时间步…和双时间步…隐式类牛顿迭代方法，并 逐渐被应用于数值分析颤振、直升机旋翼流场等非定常流动。 本文基于在定常流场计算中得到报好确证的算法一混合通量分裂简单隐式NND格式…，通 过采用高势卜奏德尔迭代求解非线性差分方程组，构造了一种时空二阶精度的隐式迭代nDlD算 法，使之适于计算气动弹性干扰问题。文中还讨论了时间精度及稳定性、运动网格生成技术、 几何守恒律的应用以及动壁边界条件等与非定常数值模拟密切相关的问题．作为检验算例， 文中计算了二维定态、弹性圆柱脱涡绕流和日本轨道再入试验飞行器作自由俯仰运动的高超 声速动态流场．并给出了气动特性参数的时间历程结果。 1数值方法 1．1控制方程 飞行器做非定常运动时，控制方程为描述气体运动的非定常NS方程和描述飞行器运动的 结构运动方程，对NS方程进行与时间相关的坐标变换后，即得到“1： ㈨… 西 a{a玎a( I af a打afJ C詈+筹+詈】 一烈j-iQ)+嚣+筹+嚣=土Re ’田家自然科学基金资助项目，批准号 —96— 塑Ot+K(S,CI，Cd，Cm)=o (2) ……… 。一7…Y程一(i堡量熏量神定义可参考文献㈨。方程(2)中，S为结构运动变量，具体形式根据 物理问题不同而改变a如果只考虑单自由度刚性俯仰运动，s具体形式如下： s=p，扫，，且口(，)=口一坼，口为攻角，∞，为平衡攻角。 1．2空间离散 无粘通矢量采用m格式离散，以流向为例，具体离敬如下： (黔__； _+；2萤j+屯+三minm呵哆；，哦÷]一严1 nm。《嗲。哦0 =三位也，)+托；蛆一五：；呱] 吆¨玑 ∥嚣 凹，盏 h验 0篙 巍 啦。瓣一；砒锁 恤嚣 哆蒜0簇