叶片通道内二次涡的数值模拟

采用三维Ｅｕｌｅｒ方程数值模拟叶片通道内流向二次涡,及其对叶片出口气流角的影响。算例表明：采用这种方法可以计算出流向二次涡;该涡造成的出口气流偏转角与实验值比较,在远离壁面处相差较小。此方法可用于叶栅通道内流向二次涡的初步分析     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

一种求解欧拉方程改进的通量分裂方法

本文以Steger-Warming提出求解Euler方程的通量分裂方法为基础,利用改进的特征值分裂方法,建立了一种连续的通量分裂方法。该方法克服了Steger-Warming通量分裂法存在的分裂后的通量在个别点导数不连续的问题,同时方法简单,便于推广。对于通量分裂后的Euler方程,采用迎风格式,建立了相应的二阶隐式有限体积方程。为了减少求解所得块矩阵方程组所需要的工作量,本文利用特征变换,在不改变原方程组离散化精度的前提下,将所得块矩阵方程组化为分离的代数方程组进行求解,从而大大减少了计算量。文中分别对一维管流,二维绕NACA-0012和NACA64A-10翼型在亚声速和超声速来流的跨声速流场进行了计算,并将所得结果与有关实验和计算结果进行了比较,均符合良好。 本文所建立的方法具有占机内存小,计算时间短的优点。在普通中小型计算机上亦可使用。     

管道中复杂马赫反射流场的高精度数值模拟

本文分析了一类高精度迎风TVD欧拉方程差分格式的精度,并用240×80个网格点进行马赫反射的计算,获得了目前实验和计算不易得到的两个激波三叉点的复杂波系结构。本文还研究了激波管上壁面对马赫反射的影响,并计算了轴对称扩张管道中激波的运动及激波与半球形障碍物的相互作用。计算结果分析指出,本文所采用的高精度迎风TVD格式与计算方法可以用来模拟具有复杂波系的工程问题研究。     

跨音速颤振计算方法研究

采用快捷的动态弹性变形技术生成三维块结构动态贴体网格。结构运动方程采用二级精度的龙格—库塔时间推进。气动力求解则用非定常欧拉方程的双时间有限体积推进,外时间为物理时间,与结构运动方程同步;对每个真实物理时间步采用五步Runge Kutta时间推进进行虚拟定常迭代。针对可压流颤振计算中存在的质量不相似问题,采用变质量、变刚度的方法计算出质量匹配点处的颤振速压值,并根据颤振速压随质量或刚度倍数的变化趋势,得出可用的飞机跨音速颤振速压。对一飞机简化外形的跨音速颤振特性计算得到了合理的计算结果。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

某型发动机涡轮盘销钉孔结构分析与寿命评估

利用有限元方法计算了某型发动机涡轮盘径向销钉孔的应力应变分布, 分析了不同销钉数目、孔销之间的不同摩擦系数和涡轮盘是否施加扭矩对销钉孔周边应力的影响, 发现了前期试验状况与实际使用的差异.利用拉伸应变能寿命预测公式进行了寿命评估, 与试验有较好的一致性.有限元计算及试验结果说明了涡轮盘加扭疲劳试验的必要性, 为某型发动机高压涡轮盘的延寿工作提供了依据.      

基于准定常的飞行器动导数的高效计算方法

结合工程设计经验和工程计算方法,针对飞行器动导数计算方法展开研究。基于求解准定常欧拉方程发展了一种适合于亚跨超声速状态下,任意气动外形气动导数的高效计算方法,并将计算结果与一般非定常方法的求解结果以及国际标准实验数据进行了对比,研究结果表明该方法结果可靠、效率较高。     

基于CFD的机翼突风响应计算

利用网格速度理论,计算机翼在锐边突风和1-cos突风下的响应,研究气动非线性和自由度耦合对翼尖加速度和升力系数的影响.采用中心格式有限体积法进行空间离散,并用双时间推进法求解非定常Euler方程.计算了刚性(沉浮)和弹性机翼在锐边突风下的加速度响应,以及刚性机翼(沉浮+俯仰)在1-cos突风下的升力响应过程,并分别与片...     

基于混合CFD方法的旋翼桨-涡干扰噪声计算

建立了一个基于Euler方程和自由尾迹分析的新的、高效的混合计算方法,用于旋翼桨-涡干扰(BVI)气动和噪声特性的研究.在该方法中,仅针对单片桨叶进行计算,通过求解Euler方程来捕捉流场中旋翼近尾迹的影响以及桨叶表面的非定常载荷,流场空间中旋翼远尾迹的位置由自由尾迹模型给出;提出了广义网格速度方法以将旋翼远尾迹的影响...