计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

适用于气动学问题的DRP类有限差分方法

研究了一类低耗散低色散的高阶精度有限差分方法,目的是直接计算非定常欧拉方程用于气动声学问题.采用的数值方法是空间四阶、时间三阶精度的色散关系保持(DRP:dispersion-relation-preserving)类有限差分格式,通过波动方程算例验证了格式模拟波动问题的能力;采用加变幅值的高波数人工耗散项来抑制高频数值振荡,使得该格式可推广用于含激波、初始间断或非线性波动的问题.通过含激波、初始间断条件的算例验证了所用变幅值人工耗散的有效性,最后作为初步应用计算了二维亚、超音速均匀流中简单点声源辐射问题,得到了很好的结果,表明此类方法在气动声学问题计算上深有潜力.     

低速大迎角三角翼旋涡流动的有粘与无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Euler方程和雷诺平均Nayierstokes方程，应用Beam－Warming近似因子分解格式及其对角化形式，采用BaldwinLomax代数湍流模型。计算模型为70°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Euler方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

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

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

低速大迎角三角翼旋涡流动的有粘怀无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Ｅｕｌｅｒ方程和雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，应用Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因子分解格式及其对角化形式，采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。计算模型为７０°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Ｅｕｌｅｒ方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

CFD技术在动态压力校准中的应用研究

采用计算流体力学(CFD)流场可视化的手段,用数值模拟的方法研究了动态压力校准装置.研究工作主要集中于激波管动态压力校准装置和正弦动态压力校准装置的数值模拟,控制方程采用雷诺平均的N-S方程描述,方程的封闭采用精度相对较高的高雷诺数湍流模型,用有限体积法对控制方程进行空间离散和时间离散.并通过实验验证了控制方程和计算方法的正确性.     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

双调和方程在航空离心泵叶型优化设计的应用

为了减少构建代理模型的计算量，提高优化效率，基于双调和方程模型和混合人工鱼群算法对航空离心泵叶型进行优化改型设计。在Matlab平台下分别采用5点4次Bezier曲线和线性函数控制叶型圆周角分布和积叠变化规律，通过软件UG和Fluent的联合批处理方法对15组设计结果进行内流场数值仿真。在曲面插值获得计算空间边界条件的基础上，采用中心差分格式对双调和偏微分方程进行数值求解，建立超曲面性能代理模型。以效率最高为目标函数基于人工鱼群算法对设计变量进行全局寻优。结果表明:基于双调和方程的超曲面代理模型能够保证试验值与预测值完全一致，而优化后的离心叶轮流动中的尾迹效应被明显减弱，水力效率比原型泵提高5.4%。      

跨声速流场计算的时间分裂有限体积法

本文用时间分裂有限体积法计算了翼型和平面叶栅的跨声速流场,计算与实验值的对比表明,数值结果是比较准确的。计算应用了滞止焓沿流线不变的关系来代替能量方程后,给出了差分方程的稳定条件。分析了跨声速流计算中,差分格式引进人工耗散的必要性。此外,本文还对加速计算收敛的方法及其效果作了探讨。     

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

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

计算机翼跨声速压力分布的格林函数方法

本文给出一种用格林定理将跨声速小扰动速势方程化成积分方程式,然后进行数值求解的方法,文中对场元强度进行松弛迭代,一般只需20次左右即达收敛。算例表明,该方法结果与差分法结果一致。     

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

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

非对称铺层复合材料板室温形状的非线性有限元分析

由于一般铺层复合材料板壳存在拉伸-弯曲、拉伸-扭转等耦合效应,在温度变化时将出现垂直于中面方向的翘曲变形,给工程使用和理论分析带来了极大的不便。但另一方面,利用非对称铺层复合材料板壳可以有效地减轻结构的重量,并可能提供利用平面模具制造曲面层合壳的新工艺。     

二维振荡机翼跨声速非定常绕流的变域变分有限元解

本文首先对文献（１－３）提出的二维振荡机翼含激波跨声速非定常绕流的变分原理进行了改造，使之能适用于时间推进法；构造了三维时空的可变节点有限元来捕获自由尾涡面，而跨声流中的激波采用人工密度法捕获。在远场边界上使用简化的无反射条件和新型非定常Ｋｕｔｔａ条件。用本文方法求解了作俯爷振荡的ＮＡＣＡ６４Ａ０１０翼型的跨声速绕流，计算结果令人满意。本文方法可以推广到三维机翼及二维和三维叶栅的同类气动力问题上去     

解跨音速升力翼型的有限元法

 在应用解全速位方程的最小压强积分有限元法求解绕升力翼型的跨音速流动时,将不可压流中求解绕升力翼型的耦合单位环量流动和无环量流动的解法推广到可压流中。为了确定环量,本文所用Kutta条件是:在后缘处,气流流向平行于后缘角二等分线。因有限元法对网格无正交性要求,因而可在椭圆变换前后进行剪切和延伸变换。这种网格生成法易于构成适用于复杂形状的有限元网格。通过计算并将其结果与文献中的数据比较,表明这种方法应用方便且有较快的计算速度和较高的计算精度。     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。     

基于库伦摩擦效应的箔片轴承极限承载力理论与试验

基于有限单元法建立了考虑库伦摩擦的波箔型径向气体箔片轴承的箔片结构模型,采用有限差分法和有限单元法耦合求解Reynolds方程和气膜厚度方程,通过求解轴颈达到极限偏心率时的轴承极限承载力,研究了箔片结构库伦摩擦效应对轴承极限承载力的影响规律,并搭建了轴承极限承载力测试试验台,利用温度法测量了两个具有不同轴承壳内表面粗糙度的波箔型径向气体箔片轴承的轴承极限承载力.通过对比分析仿真结果与试验结果表明:轴承壳圆柱孔内表面粗糙度为0.4μm的轴承在10000r/min和20000r/min下,轴承极限承载力分别为15.5N和42.3N;而表面粗糙度为1.6μm的轴承极限承载力为10.9N和29.6N,这是由于波纹箔片和轴承壳体之间的库伦摩擦力增大了波纹箔片的刚度,因此增大箔片结构摩擦因数使得轴承极限承载力降低,并且仿真结果变化趋势和试验结果变化趋势吻合.      

一种无条件稳定的时域推进法与有限元混合分析瞬态热传导问题

本文将一种新的无条件稳定的时域推进法与有限元结合,用于分析瞬态热传导问题。根据微-积分型热传导控制方程,对时间变量在小区闻内插值积分,变初-边值问题为一系列离散时刻的边界值问题,再应用有限元法求解之。由于推进中的每一时刻的解,都严格满足原给定的初始条件,就消除了累积误差的影响。此外,这种时域推进法对时间变量的数值积分也优于数值微分。因而,与现行的基于微分型控制方程,采用时-空有限元模型同时离散,再逐步求解的直接积分算法相比,计算精度可大大提高。不难预料,当求解较长时间后的瞬态值时,本算法的优越性会更加明显。     

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好