高速流非平行边界层稳定性研究

导出了可压缩流的抛物化稳定性方程（PSE）。针对高速流动,特别是超声速和高超声速流动的非平行边界层稳定性问题进行了研究。引入高效的边界层变换、全流场高精度的差分格式及预估校正迭代和推进求解法来求解PSE方程,使得PSE方法中至关重要的正规化条件得到了满足,确保了数值计算的稳定。采用高马赫数下对稳定性起支配作用的第二模式,研究了高速流边界层稳定性的演变和特征,分析了流动的非平行性、压缩性,以及壁面冷却等因素对流动稳定性的影响．所得结果与相关实验数据吻合较好。     

角接触球轴承分析模型的数值求解

建立了内圈静止、外圈旋转的角接触球轴承分析模型,考虑了滚动体的离心力、陀螺力矩以及弹性流体动 力润滑作用.建立了非线性方程组未知量的约束条件,解决了Newton-Raphson法求解方程易出现不收敛的问 题.以两个角接触球轴承为算例,分析了轴承内、外圈分别旋转对应力的影响,以及转速和径向载荷对轴承工作 性能的影响.结果表明,轴承外圈转时接触应力随转速增大而显著增大,内圈转时则无明显变化;低转速下外圈 转的接触应力明显低于内圈转;较小的径向栽荷下外圈转的接触应力高于内圈转;轴承微区滑动速度受转速和径向栽荷的影响显著.通过Jones球轴承拟静力学分析程序计算,验证了模型和程序的有效性.     

平面Poisson问题的间接变量等价边界积分方程

使用变分法确立平面Laplace外问题的确切形式；在此基础上，为了避免使用无限区域上的Green公式(不成立)，通过构造紧支集上的无穷可微函数，获得满足非齐次控制微分方程的特解的重要结果。使用非解析开拓数学方法，导出满足控制微分方程的位势函数集的等价边界积分表示式。获得Poisson问题的间接变量等价边界积分方程；通过一些例子，详细讨论了传统的外问题和边界积分方程的形式，表明它们不具有普遍适用性。     

Euler方程无网格算法及布点技术

研究了Bitina的显式无网格算法。在该算法的基础上,给出了Euler方程无网格离散形式,运用RungeKutta显式时间推进格式推进求解。文中的耗散项采用了与非结构网格上类似的守恒型耗散算子,边界条件的处理借鉴了结构化网格处理技术。此外,本文还描述了一种区域离散布点方法,研究了点云生成的选点准则,并成功地数值模拟了二维翼型的典型绕流。     

一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。     

Karman型四边夹紧正交异性矩形薄板的中等大挠度

利用Galerkin方法分析了Von-Karman型四边夹紧正交各向异性矩形板。所设的位移函数为梁振型函数，它不仅能精确地满足边界条件，而且具有正交的特性，从而把复杂的非齐次非线性偏微分方程组化为一组非线性代数方程组，通过非线性方程组的线性化和可调节参数的修正迭代解法找出问题的解。实践证明，梁振型函数收敛很快，只须取出级数的前几项即可满足精度要求。最后求出了不同复合材料的挠度和应力值并同已有的结果进行了比较。     

基于多GPU的大型线性和非线性方程组的求解(英文)

在处理工程问题时,常常需要对线性或非线性方程组进行求解。对于实际应用中经常遇到的大型方程组进行求解则需要相当长的时间。使用图形处理器(GPU)代替传统的CPU,将多块GPU通过操作系统进行协调,并将PBi-CGstab方法和Inexact Newton方法进行适合多GPU并行的改造以此作为多GPU求解器的核心算法,加速求解大型线性和非线性方程组。本文的多GPU求解器在成倍扩展了单GPU求解器允许的计算规模的同时取得了令人满意的加速比。     

不同温度和应变速率下BTi6431S新型钛合金流变应力行为

通过在870～900°C温度范围内和应变速率为0.000 1～0.01/s下,对BTi6431S新型钛合金薄板进行了单向拉伸试验,研究了BTi6431S新型钛合金在高温条件下流动应力与变形温度和应变速率之间的关系,并利用改进的Hooke law和Grosman方程建立了BTi6431S钛合金的应力-应变本构模型。研究结果表明:BTi6431S钛合金在870～900°C温度下,流动应力随着温度的升高而降低,随着应变速率的增大而明显升高,在进入塑性变形区间,基本保持稳态流动。拟合得到的本构方程能够较准确地反映了BTi6431S钛合金在高温下的流动应力变化情况。     

不同温度和应变速率下7022铝合金流变应力行为

通过在293-773 K的温度范围内和应变速率为0.001-0.1 s-1下对7022铝合金薄板进行温拉伸试验,研究了7022铝合金温拉伸性能,以及该合金在升温条件下流变应力与变形温度和应变速率之间的关系.并利用改进了的Hooke law和Grosman方程建立了7022铝合金在温拉伸时应力-应变本构模型.研究结果表明:7022铝合金的流变应力随温度的升高而降低,随应变速率的升高而升高;温拉伸试样的延伸率随变形温度的升高而升高,随应变速率的增大而减小.     

基于N-S方程的尾迹面法翼型气动阻力计算

研究了计算翼型阻力改进的方法——尾迹面法及尾迹面积分的位置和相应的积分技术.在亚跨声速时分别采用了表面积分和尾迹积分求得给定翼型RAE2822的阻力.结果显示,2种方法具有相同的结果,表明尾迹面积分方法是有效的,与实验值吻合较好.尾迹面法中,尾迹面位置应处于离后缘相对弦长距离0.6~1.0之间.尾迹面积分方法中积分结果不依赖于物体的详细几何外形,可以预计对曲率变化大的三维复杂外形,该方法有更大的优势.