一种分离流数值求解方法

本文以弯管内分离流场为例，根据分离区与未分离的湍流动能差别较大的特点，提出了在分离区和未分离区采用不同的湍流模型来模拟流动的湍流效应，利用Ｎ－Ｓ方程进行了数值模拟，结果表明：利用这一方法，可以加速迭代的收敛，节省机时。     

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果     

弹性变形金属的力学化学效应与腐蚀过程相关性的研究

１根据力学，热力学、电化学的原理，就弹性变形对金属热力势和电机电势的影响，从而加速金属腐蚀的问题，进行了详细的理论分析与探讨。     

适用于激波/边界层相互作用的线性涡粘性湍流模式

本文选择了两个激波 /边界层相互作用诱导分离的跨声速问题 (轴对称圆弧突起和二维管道突起 ) ,采用五个有代表性的湍流模式 (BL模式 ,JL k-ε模式 ,k-ω模式 ,SST模式和双尺度模式 ) ,通过将数值计算结果和实验结果进行比较 ,对有关的湍流模式进行了评估和分析     

双转子临界转速的简易分析方法及应用

采用纯弯曲及弯扭耦合两种传递矩阵法和子结构传递矩阵法计算具有畸形结构的双转子系统的临界转速，提出了用线性插值法确定双转子临界转速的简便方法和由不平衡响应系统位能峰值转速确定双转子系统计及阻尼影响的临界转速，本文分析方法有较大的工程应用价值。     

绕大攻角翼型低速非定常流动的SUPG有限元数值模拟

采用既有良好稳定性又有较高精度的ＳＵＰＧ有限元方法,通过数值求解涡量一流函数形式的不可压缩流ＮＳ方程组,模拟大攻角翼型在高Ｒｅ数下几种典型情况的非定常绕流流场,同时验证了ＳＵＰＧ有限元方法在求解对流占优的二维非定常对流的－扩散方程时的有铲性。     

用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

跨声压气机动静干涉效应的数值研究

采用数值方法对一轴流跨声压气机在设计点的非定常流场进行了模拟, 对级内动静干涉进行了深入分析, 研究了尾迹和位势作用等对动静叶表面气动负荷的影响.计算结果表明:在压气机中, 上游动叶的尾迹等对静叶通道内部流动, 叶片表面静压的波动, 以及边界层流动损失的发生、发展和输运产生明显的影响, 需要在设计中加以考虑.      

新型高低速流统一算法及其叶轮机应用

 对于可压缩/不可压缩混合流动或高低速混合流动问题,目前没有较好的解决方法。基于Roe格式及预处理方法,提出了全速度Roe格式。全速度Roe格式是将Roe格式与低速Roe格式通过一个与本地马赫数相关的函数联合而得到的。与标准Roe格式相比,对亚音速流动,全速度Roe格式中声速对数值耗散的影响随着马赫数变化而变化。与预处理Roe格式相比,全速度Roe格式简单易于实现,而且没有计算稳定性的限制,不需要采用全局截断处理,从而在理论上对高低速流动具有一致精度。事实上,全速度Roe格式意味着可以通过直接改变Roe格式的非线性特征值,解决可压不可压流动统一计算精度的问题。这是对统一算法的新观点。数值模拟表明,全速度Roe格式具有与预处理Roe格式一样的激波捕获能力、极低速流动准确模拟能力。     

钉载作用下多裂纹板的裂纹闭合接触问题

 基于摩擦接触问题的数学规划解法,采用各向异性体平面弹性理论中的复势方法,建立了含裂纹群有限大各向异性板,在钉载作用下裂纹闭合或局部闭合问题的有效分析方法。忽略钉与孔间的摩擦,假设钉载沿孔边呈余弦分布,通过在可能闭合的裂纹边界引入互补变量函数并将其展成Fourier级数形式,以Faber级数为工具,应用保角映射技术和最小二乘边界配点法,导出裂纹面摩擦接触的线性互补模型,并通过算例验证了方法的有效性。数值结果表明,由于采用级数解描述板应力场和位移场,本方法具有较高的计算精度和效率,便于研究裂纹闭合对应力强度因子等断裂参数的影响。