跨音松弛法在叶栅绕流中的应用

文章介绍跨音松弛法应用于叶栅绕流的定常与非定常流场计算。在定常流中,采用弦线法向小扰动假定;在非定常流中,采用线化非定常分量小扰动假定,从而简化了数值求解过程。由于松弛法大大地节省机时和内存,便于工程上推广应用。     

“跨音速非等熵全位势方法及粘流/无粘流相互作用的计算”通过鉴定

流体所朱自强教授为首的课题组自1986年历经5年研究的“跨音速非等熵全位势方法及粘流/无粘流相互作用的计算”课题于1991年底通过校级鉴定,1992年3月初通过部级鉴定。课题组发展了一套以非等熵全位势方法解和三维可压缩湍流边界层积分方程反方法解相耦合     

绕大钝头短体外形跨音速大迎角湍流流动计算

采用数值模拟求解N-S方程方法计算与分析了绕大钝头短体外形的跨音速大迎角时的湍流流动.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD(Total Variation Diminishing)性质,采用中心差分格式来离散粘性通量项.计算中采用了B-L湍流模型.算例定性地表明绕流的三维分离具有准周期性的非定常特性.     

守恒型有限元法及其在跨音速流中的应用

本文提出一种新的定常无粘跨音速流动的有限元计算方法——守恒型有限元素法。 这方法以积分形式的通量守恒方程作为问题的数学模型。从这一方程出发,推导了单一元素上的有限元类比方程,再将这些有限元类比方程进行最小二乘方迭加。最后用最优化的方法进行求解。由于所采用的优化方法可起到捕捉激波的作用,本方法没有引入人工粘性项。 最后本文给出了几个算例,并将计算结果与实验数据进行了比较。算例表明,二者吻合较好。     

跨音压气机叶栅流场的附面层修正

本文发展了二维压气机叶栅亚音、跨音流中粘性——无粘相互作用程序。无粘计算在弦向大扰动、弦法向小扰动的假设条件下用跨音松弛法求解;粘性附面层用Nash方法计算。考虑了叶栅尾迹的影响。最后,连同通道区和尾迹区一起,求得了包括冲波影响的叶栅的总损失系数。计算结果表明,在跨音情况下,粘性对冲波强度、位置以及全流场的速度分布都有不可忽视的影响。     

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

求解跨音位流问题的改进SIP方法

本文研究Sankar提出的,用全位势方程求解跨音流场的强隐式方法(Strongly Implicit Procedure)。针对强隐式方法在求解跨音流场中的问题作改进,即采用了新的人工压缩密度计算方法,计算绕孤立机翼、叶栅通道的跨音流场,降低了在跨音位流计算中激波前后的数值振荡,并提高了算法的稳定性,     

民用通风机气动设计中弯掠长叶片的数值模拟

采用数值模拟的方法,分别对弯掠叶片造型和直叶片造型两种情况下高性能民用通风机叶片内三维紊流流场进行了计算;计算程序中采用有限体积法,在任意非正交曲线坐标系下,用SIMPLE算法求解三维N-S方程和K-ε湍流方程.计算结果表明,弯掠叶片造型的民用通风机的流动参数和总性能参数更接近设计指标,而且采用弯掠叶片可以更有效地控制气流的流动,对减少二次流损失起了很大作用.      

复杂目标散射强度的三维分布特性

建立复杂目标几何外形的数学模型,对独立目标使用外法向判断可见面元,对复杂目标使用Z-buffer算法处理自动遮挡问题.通过求解电磁场边值问题求解目标表面的等效电磁流分布,并求解散射强度的三维分布特性.仿真计算结果被若干测试结果验证,证明该方法的精度较高、实用性强、可靠性好.算法可以直观地显示目标散射强度的三维分布图.      

尾流激振情况下叶片强迫响应瞬态分析方法

对叶片强迫响应问题,提出应用瞬态分析的方法,针对转子叶片在前排静子叶片尾流激振情况下的位移和应力进行预估.流场进行全场三维非定常求解,将三维非定常气动力引入转子叶片有限元结构计算中,对尾流激振情况下叶片强迫响应问题进行瞬态分析,比较了不同转速下叶片的强迫响应,得到响应位移及应力变化的分析结果.采用叶片强迫响应瞬态分析方法,从流场的求解到结构响应的求解均应用成熟的通用软件,为工程应用该方法进行尾流激振情况下叶片振动应力预估及阻尼的研究提供了一个可行的方法.      

带副翼的翼身组合体绕流的 Euler方程解

提出了一种求解带副翼偏转的翼身组合体绕流的Euler 方程计算方法.将对接分区网格与分区求解算法相结合,有效地求解了绕此外形的复杂流动.提出了一种满足通量守恒的内边界耦合条件.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种限制器(Limiter)函数以保证TVD性质.数值算例表明本文方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

Durbin法在阻尼梁动响应求解中的运用分析

运用铁木辛柯梁理论和K-V阻尼理论，研究了非比例阻尼梁在冲击载荷作用下的频域振动求解方法。推导采用了传统拉普拉斯正变换和基于Durbin公式的拉普拉斯反变换策略（统称拉普拉斯法），发展了阻尼梁系统的动力学方程解法。拉普拉斯法的推演同时涵盖了3种典型的梁边界条件，具有广泛的适用性。数值法的验证采用了特殊构造的比例阻尼点条件，并与基于模态叠加法的求解结果进行了对比分析，且数值算例充分考虑了数值参数和系统参数的影响。计算结果表明:在不同边界条件和受载状态下，拉普拉斯法与模态叠加法均能合理地计算出基本阻尼梁系统的动响应曲线，且两者的求解精度保持在同一量级；同时，捕捉到拉普拉斯法的求解精度会受到系统长细比等参数的影响。拉普拉斯法具有比传统实、复模态叠加法更易操作的特性，但其精度受到了算法固有参数和阶跃外载型式的影响，稳定性仍需进一步提高。      

高超音速化学非平衡流的隐式格式迭代解法

针对高超音速化学非平衡流计算中的隐式TVD格式给出了一种迭代解法.该方法是在原有的Harten-Yee格式基础上,在每一时间步对无粘部分进行迭代运算,目的是为了加速收敛.借助对模型方程的收敛性分析,证明该迭代算法是收敛的.二维数值实验表明该方法对于高马赫数时的化学非平衡流能明显提高隐式格式的计算效率,其计算结果是满意的.     

固体火箭发动机燃烧室三维流动数值计算

采用k-ε双方程湍流模型,以SIMPLE计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流N-S方程.在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构.计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡.      

固体火箭发动机燃烧室三维流动数值计算

采用ｋε双方程湍流模型,以ＳＩＭＰＬＥ 计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流ＮＳ 方程．在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构．计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡．     

飞行器结构缝隙电磁散射问题的研究

提出了求解飞行器表面结构缝隙电磁散射的方法,由场等效原理及场连续条件,运用矩量法解关于缝隙口面上等效磁流为未知数的方程,由此得到长直缝隙口面上的等效磁流的数学模型,基于上述理论,由座舱结构缝隙建立无了大导电曲屏面上的缝隙的数学模型,将曲屏面上缝隙给合理划分成若干个似的直缝隙,求解每一段直缝隙的等效磁流,由辐射积分方程求解该等效磁流的散射,由这些散射场的叠加得到曲屏面上的缝隙的散射场,最后给出飞行器     

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

三维多段机翼复杂流场的计算和分析

采用中心差分的有限体积方法和分区技术求解了N-S(Navier-Stokes)方程,在分区求解技术中采用了满足通量守恒的内边界耦合条件.分析讨论了三维襟翼的绕流流场,对"剪刀差"处的三维分离进行了探讨,初步揭示了绕流的主要特征及对附近物面压强的影响.同时,分析了三维襟翼前缘的分离流动.      

一种前飞状态直升机可达区域快速计算方法

直升机可达区域可用于生成空域安全态势,提高直升机在低空复杂环境中的飞行安全水平.为了满足求解可达区域的实时性和准确性要求,提出了一种快速计算方法.根据能量机动性求解了直升机控制域,结合运动学方程建立了可达区域数学模型.基于最优控制原理,实现了边界点快速计算;提出了基本可达区域的概念,将实时可达区域构造问题转化为少量边界点的计算问题.以AH-1G直升机为对象进行了二维和三维的算例分析,并与微分包含法计算结果进行了对比.算例结果表明,直升机控制域是时变的,需根据实时飞行状态进行计算;该算法能实时计算前飞状态直升机的可达区域,有较高的计算精度.