基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

非线性双曲型守恒律的一类应用通量限制器的二步二阶格式

本文提出了计算非线性双曲型守恒律弱解的一类新的二步二阶精度有限差分格式,这类高精度非线性差分格式的构造应用了三点单调差分格式的数值通量和通量限制器。文章证明,在拟CFL条件限制下,格式为总变差衰减差分格式(TVD格式),因此不发生常用的诸如二步二阶Lax-Wendroff格式、Maccomark格式和其他古典二阶格式计算不连续解时在其周围出现的剧烈振荡现象。 为了减少Roe推广技术引起的计算复杂性,文中我们基于Davis,S.F.的思想,使用了Rusanov数值通量,把标量差分格式推广到方程组,构造了求解双曲型守恒律方程组的新的预测-修正格式,较之其他一般的预测-修正格式而言,仅增加了很少的计算复杂性,初步的数值试验表明,本文提出的格式无振荡和具有极好的激波分辨率。     

线性方程组的并行算法及在Transputer系统上的实现

在“一种有效的多Transputer系统的并行算法——ABC法”一文的基础上,本文进一步研究将ABC法用于变带宽矩阵线性方程组的求解问题,对线性方程组的系数矩阵采用了逐行一维存储方式,提出了相应的并行Gauss消元法,给出了该算法的效率.分析结果表明,带宽越大方程阶数越高,这种算法的效率就越高。因此本算法适用于高阶的大带宽线性方程组的求解问题. 根据本文的算法,编制了线性方程组的并行求解程序,并分别在一个、二个和四个T414系统上做了若干算例,结果表明本文分析的结论是正确的。     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

NiTi形状记忆合金的拉伸实验与物理模型

根据NiTi丝的拉伸实验结果,基于滞后元方法,建立了拉伸状态下应力-应变曲线的唯象物理模型,并进行了数值模拟。研究结果表明(1)SMA丝在马氏体相变,特别是R相变时出现一系列物理、力学性质的异常变化,其应力-应变关系呈现高度非线性;(2)在NiTi丝断裂之前没有明显的屈服现象;(3)物理模型的数值模拟结果与实验结果吻合较好。     

基于响应面法的跨声速翼型气动优化设计

基于响应面方法进行了跨声速翼型的气动优化设计。流场计算采用雷诺平均N-S方程。响应面模型采用二次多项式来构造，试验点的选取满足D优化准则。设计结果表明该方法只需通过较少流场计算就能对翼型的单个或多个设计点、在多约束条件下进行气动优化设计，设计质量较高，有较大的工程应用价值。     

堆积粉尘本构方程的实验测试

确定堆积粉尘的本构方程对讨论有关现象是十分重要的。设计了一种测试方法 ,利用阴影摄影和带示踪粒子的X光脉冲摄影 ,记录激波与堆积粉尘相互作用的流场 ,再通过对该流场的理论分析 ,以确定堆积粉尘的本构方程。以超细淀粉为例 ,对所测的方程做了实验验证。     

不对称喷管欠膨胀超声速射流的数值模拟

采用改进的MUSCL格式解三维可压缩平均雷诺纳维尔 斯托克斯方程组,湍流模型为Spalart Allmaras代数模型,数值模拟了四种调整片不对称喷管在喷口压比为2.0、2.8、3.4和4.0条件下的欠膨胀超声速射流场,获得了射流场流谱和参数分布,并给出了单调整片不对称喷管出口压强比与射流偏转角和扩张角的分布曲线,计算结果和实验数据符合良好。     

基于细以方法的颗粒增强复合材料弹塑性分析

针对高体积含量的颗粒增强复合材料，提出了一种细观结构模型：将颗粒简化为同质、同尺寸的弹性圆球，两颗粒之间的连接基本简化为一弹塑性短圆柱体，并假设细观应力、应变和塑性区均为轴对称分布。其体和颗粒的变形行为分别简化为类似于弹性地基和弹性半空是，以此为基础，建立了反映一对颗粒法向和切向变形协调关系的两个积分方程。数值求解这两个方程，可得到一对颗粒间基体中细观应力的分布形式，从而建立颗粒对的细观弹塑性本构，采用平均化方法，进一步推导出材料宏观的应力－应变关系。本文利用该模型对一种金属基复合材料的拉伸实验进行了模拟，理论预测与实验结果吻合较。     

Euler方程的动力学通矢量分裂格式

文中基于气体分子动力学理论，研究和发展了一类新型逆风差分格式——动力学通矢量分裂格式（ＫＦＶＳ），并论证了格式的逆风性质和Ｌ１稳定性。最后给出了高精度格式的构造，并就一维激波管问题作了数值试验，还与精确解作了比较，结果令人满意。