解对流扩散方程的有限元逆风因子法

 本文提出了一种解对流扩散方程的有限元逆风因子法。这种格式的数值解的稳定性很强,计算收敛速度很快,所引入的虚假扩散较少。一系列的数值试验证实了上述结论。     

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

前缘对进气道性能影响的数值模拟

应用质量平均的N-S方程和B-L代数湍流模型，在不同压缩面前缘半径条件下，以高超声速进气道的流场进行了数值模拟，对比研究了进气道压前缘半径对流场结构及性能的影响，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散。     

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

解不可压N—S方程的PCIMG法

本文提出了一种解不可压N-S方程的预测、修正、迭代、多重网格法。这种算法程序简单,结果精确及收敛快速。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

超声速对流推力矢量喷管的数值模拟

利用二次流对流剪切激励主射流方式来产生射流对流推力矢量方法是一种效率很高的矢量化喷流新方法,又称为对流推力矢量法。本文运用三阶精度MUSCL TVD有限体积格式求解二维非定常、可压缩、雷诺平均N-S方程和Spalart-Allmaras湍流模型,计算了对流推力矢量喷管的超声速射流问题,获得几种不同二次流吸气压强下对狭缝和喷管肩部壁面压强分布、剪切层诱导旋涡的结构特征和主射流矢量角等结果,计算结果与实验数据对比取得良好的一致性,并对计算结果进行了详细分析和讨论。     

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用

在非结构网格下,针对非定常三维N-S方程组发展了一种双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间上将r=3的加权ENO格式与强紧致格式相结合去处理N-S方程中的对流项以及离散方程的右端项,并用四阶精度的紧致格式去计算N-S方程中的粘性项。典型的3个算例从不同侧面对格式进行了考核。计算表明:该算法具有高效率与高分辨率的特征,所得的计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下模拟非定常流动的物理过程。      

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

线性GMRES算法在求解全速位方程中的应用

从全速位方程出发,利用压强极小积分得到有限元方程组。采用适用于有限元法的块结构网格生成技术,给出求解有限元总体方程组的线性GMRES算法。该算法比线松弛迭代法的收敛速度快得多。对M6机翼跨音速有升力情况进行了计算,其结果与实验值符合较好。     

基于浸入边界法的低雷诺数流固耦合数值模拟(英文)

提出一种基于SIMPLE算法的非定常流固耦合计算方法。流体Navier-Stokes方程空间采用非结构化网格有限体积法离散,时间项采用了欧拉隐式方法。利用浸入式边界方法模拟静止或者运动固体区域,流固界面作用力通过流体体积(VOF)方法进行处理。从而可以用固定网格求解任意复杂区域中的流固耦合作用。本文模拟了低雷诺数静止及振荡圆柱绕流,所得结果与文献中贴体网格计算结果吻合,从而验证了本文方法的合理性和正确性。     

固体中脉冲超声波传播的有限差分模拟

脉冲超声波传播过程的模拟研究对实际超声检测的信号分析具有重要的借鉴意义。对弹性固体中的声场方程应用时域有限差分法,进行了二维声场特性研究。通过完全匹配层方法对截断边界的吸收处理,解决了无限自由空间的计算问题。模拟计算表明该方法可用于各向异性介质的分析。分析了液固耦合情况下的超声波传播特性,取得了与光弹法物理试验相一致的结果。通过与平底孔试块实际缺陷回波信号的对比,验证了脉冲超声波有限差分法模拟计算的有效性和正确性。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

多钉连接件钉传载荷计算的一个解析方法

对多钉连接件钉传载荷的计算问题提出了一个解析分析方法,推导了求解钉载的线性代数方程组并给出了若干算例。该方法使用方便,并具有工程分析所要求的精度。     

斜坡增强混合的数值研究

采用隐式方法求解三维可压Ｎ－Ｓ方程及组分方程,数值模拟了后掠、无后掠斜坡的混合增强流场。比较了两种斜坡的混合增强效果及产生的流向涡量的强度。计算结果表明后掠斜坡产生流向涡的强度明显高于无后掠斜坡,其混合增强效果也更好。     

亚临界情况下翼型的非线性压缩性修正

本文从有限基本解方法出发,介绍一种计算亚音速流中非线性压缩性效应的方法。采用逐步逼近的办法,可以把有限基本解方法推广到求解非线性的无粘可压流体,从而为非线性方程的求解提供了新的途径。 编制了FORTRAN语言程序。与精确解和其他理论值的比较表明,结果比传统的有限差分法好,尤其是在翼型头部。