机翼静气动弹性数值模拟研究

通过求解基于三维非结构网格的欧拉方程计算机翼气动特性,采用结构柔度影响系数方法进行机翼静弹性变形的计算,利用三维非结构弹性网格技术进行几何外形变形下的网格修正,将以上三部分进行耦合迭代.实现了机翼的静弹性特性的数值模拟.采用开发的数值分析程序,进行了某支线飞机机翼静弹性特性的计算分析,计算结果合理,表明所发展的数值模拟方法可用于机翼静弹性气动特性的分析研究.     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切.使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格.应用特征恢复技术恢复凹面特征并改善凸面处的网格质量.流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量.加入Spalart-Allmaras一方程湍流模型以获得湍流解.上下对称高斯-赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性.算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算.     

基于图像处理思想的自适应网格技术

为了提高激波和旋涡等复杂流动现象的捕捉精度和计算效率,基于图像处理的思想,研究了非结构网格的自适应加密和稀疏算法.将流场计算网格与图像像素类比,借鉴图像处理的方法,提出了一种流场自适应区域的分割算法.该方法能自动对自适应区进行分层,中间区域网格可以达到原来多次自适应的效果,而且提高了生成网格的质量.将算法嵌入到流场求解器中,根据流场求解结果对网格进行自适应,并使用算例进行了验证.     

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

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

固体火箭发动机药柱的有限元网格划分研究

讨论了混合网格划分功能在对固体火箭发动机药柱进行有限元分析中的应用,对同一推进剂药柱进行了两种不同形式的网格划分,并对计算时间进行了比较,提出了推进剂药柱网格划分应注意的一些细节方面。     

贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

含偏心孔复杂盘体三维有限元分析

一、复杂三维盘体网格划分的基本策略　　图 1为含偏心孔的压气机盘实体 ,略去其与盘强度分析无关的若干细节 (如榫槽和花边 )。图 1　偏心孔的压气机盘实体对于图 1所示的压气机盘 ,根据其周期对称性取 30°扇区为分析对象 ,含半个偏心孔。这是一个非常复杂的几何域 ,其分网时遇到的主要困难是形体不规则和三维各向尺寸差异过大。广义上说用自由网格生成技术是可行的 ,但是实际上是行不通的。主要原因有 :( 1 )为了使单元自身的几何性质保持良好 ,对四面体、五面体单元畸变值要大于 0 .4,六面体单元畸变值要大于 0 .7。这势必导致节点数剧…     

用POISSON方程生成绕二维翼型及其它形状物体的网格

本文是在文献[1]的基础上,发展了用 Poisson 方程生成绕二维翼形和其它形状物体的有限差分网格的方法,对于非齐次项的处理,利用自动选择施加于内、外边界上的两个重要因子来实现:第一个因子控制与边界相交的网格线上网格点之间的距离,第二个因子控制网格线与边界的交角;并对计算方法和控制参数进行了详细的讨论。用 FORTRAN 语言编制了实现此方法的计算程序,给出了一组算例。     

多部件模型在全尺寸小型双函道涡扇发动机气流数值模拟中的应用

利用高阶单调Ｇｏｄｕｎｏｖ显／隐格式求解全尺寸小涡扇发动机Ｓ２流面上带粘性力项的非定常Ｅｕｌｅｒ方程组,获得了令人满意的稳定收敛解。对多部件模型方法在双函道小涡扇发动机气流数值模拟应用中存在的问题提出了建议。