预处理特征型边界条件及其回流处理技术

通过特征型预处理Euler方程推导出预处理特征型的边界条件处理方法,并发展了预处理的回流处理技术,增强计算的鲁棒性和精度。通过特征型预处理Euler方程的特征关系得到边界条件处理方法;在亚声速速入口和出口边界,根据其特征方程确定回流边界处理方法。通过数值模拟发现:带回流的特征型边界处理方法比Turkel方法收敛快1.5倍以上,而计算精度与SIMPLE方法相当;在出口带有回流的计算中,其压力误差比Fluent降低68.3%。由此可知带有同流处理的特征型边界处理方法具有较高的鲁棒性和精度,在实际的工程计算中有一定的应用价值。     

渐变边界与几种无反射边界的初步比较

简要介绍了通量渐变消失（gradually diminished flux,GDF）、完全匹配层（perfectly matched laye,PML）和牛顿阻尼边界条件,利用典型算例将这三种无反射边界进行了初步比较。     

边界条件对三维空腔流动振荡的影响

基于全隐式无分裂方法,数值求解三维非定常雷诺平均的N av ier-S tokes方程,紊流模型采用修正的B-L代数模型。运用此方法,本文数值模拟了三维开式空腔的超声速、跨声速、亚声速流动。研究了空腔口处剪切层的非定常运动现象,着重探讨了入口边界条件对三维空腔流自持性振荡的影响。数值计算结果表明,与层流速度型相比,入口条件为紊流速度型和特征边界型得到的结果能够吻合。     

二维多介质流动问题的界面处理方法

就二维可压缩多介质流动问题的数值模拟,给出了一种新的界面处理方法。通过在界面处构造Riemann问题,利用Riemann问题的解分别定义界面两边流体的边界条件,由于Riemann问题的解准确地描述了界面处流体的流动状态,因此得到了更加准确的界面边界条件。本文将由Riemann问题的解得到的界面速度外推到整个流场,重新定义速度场,避免了由于速度的大梯度变化而导致的Level—Set等值线相互交错,因而得到了更加精确的界面位置。利用该方法对水下激波与柱型气泡相互作用问题进行数值模拟,结果表明该方法能准确地捕捉各种物理现象。     

矩形变截面梁横向振动自振频率的传递函数渐近解法

分析具有固定宽度而高度成线性分布的楔形梁。引入状态变量将模形梁的自由振动控制方程及边界方程写成状态空间形式。定义小参数,并利用摄动方法,得到常系数微分方程,从而得到问题的摄动解。最后给出了一些数值算例,验证了方法的可行性。     

多层介质涂覆腔体电磁散射特性的IBC/FEM分析

将Leontovich阻抗边界条件(IBC)与矢量有限元法(EB-FEM)相结合对内壁涂覆多层介质开口腔体的电磁散射特性进行分析。将IBC和腔体口面边界积分方程分别转化为第三类边界条件形式，从而导出边值问题的相应泛函。表面输入阻抗通过广义反射系数得到，计入了介质内部电磁波的多次反射，提高了IBC／FEM混合法的精度。数值结果表明这种IBC／FEM方法在分析涂覆多层介质开口腔体电磁散射问题中的便捷性与精确性。     

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

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

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

矩形变截面梁横向振动自动频率的传递函数渐近解法

分析具有固定宽度而高度成线性分布的楔形梁。引入状态变量将楔形梁的自由振动控制方程及边界方程写成状态空间形式,定义小参数,并利用摄动方法,得到常系数微分方程,从而得到问题的摄动解,最后给出了一些数值算例,验证了方法的可行性。     

密实织物伞形状和应力计算研究

建立了充满稳定状态下密实织物伞形状和应力计算的理论模型。将充满的伞看成柔软的壳体，用力学原理建立伞绳（伞衣上径向加强带）和伞衣幅中线的力平衡方程。伞绳和伞衣处理成非线性弹性构件。考虑伞衣中线上的子午向应力，得到了一个双轴应力模型。结合伞充满状态下几何大变形的几何关系式，得到一组含６个微分方程的非线性微分方程组。从伞顶孔开始，利用伞顶处的边界条件，离散积分至伞底，校对伞底处的另一边界条件，如不满足再从伞顶开始进行迭代，直至满足伞底边处的边界条件。由此编成进行非线性迭代数值求解程序ＣＳＬＡＰ。选用一具典型的平面圆形伞进行了验证，结果很好     

基于嵌套网格和计入尾迹影响的旋翼流场数值模拟

给出了三维Euler方程数值模拟悬停状态旋翼流场的方法和模型。为充分考虑旋翼尾迹对流场的影响,采用了由动量理论导出的远场边界条件和由尾迹诱导速度代入滑移条件而构造的新翼面边界条件。为减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。针对嵌套网格中关键的贡献单元搜寻问题,本文采用一种伪贡献单元搜寻法。应用上述方法进行了算例计算,给出了有无尾迹修正时的旋翼桨叶表面压力分布和沿展向升力分布的计算结果,并与可得到的实验数据进行了对比。此外,从计算的等涡量线图上分析了旋翼尾迹流场畸变的发展趋势,表明了Euler方程法在旋翼涡尾迹捕捉方面的能力。     

双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

边界层中TS波及其高阶谐波的演化

基于抛物化稳定性方程，研究了边界层中TS波及其高阶谐波的线性和非线性演化问题。由局部法和Landau展开式导出初始条件，并计算了扰动幅值和速度型等的演化过程和特征，特别是非线性的重要作用。探讨了初始幅值、压力梯度、扰动频率对扰动演化的影响及其规律，这与边界层的稳定性和转捩研究紧密相关。算例结果与全Navier—Stokes方程的直接数值模拟结果一致。     

一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。     

翼型在自适应壁风洞中试验时的侧壁边界层三元效应研究

二元翼型试验时模型安装于侧壁转盘上会引起侧壁边界层三元效应，中国航空研究院西北工业大学和德国宇航院哥廷根流体力学研究所合作用计算和试验方法进行了研究。作者介绍了1995～1998年的研究，包含研究目的，方案，计算方法，在德国路德维希管风洞中的验证试验，结果分析和初步结论。     

三维弹性边界元分析中面力不连续和1/r积分奇异性问题的处理

本文将文[2]中的方法推广应用于解决三维弹性面积分中1/r的奇异性问题,从数值计算角度,较详细地讨论了消除1/r奇异性的原理。为了取得较高的数值精度及较少的运算时间,本文对高斯点数目的选择原则也作了分析与讨论,提出了一个实用的经验公式。针对三维弹性问题中常见的边界面力不连续问题,提出了一种简单的处理方法,即直接从边界离散化后的边界元方程出发,按力不连续点所处的积分单元分别进行处理。两个典型算例的数值试验结果表明,文中所用的方法是行之有效的。     

复合材料夹层圆柱壳的轴对称失稳问题

本文首先给出复合材料夹层圆柱壳大挠度问题的一般方程,然后在轴对称变形条件下进行简化,得到轴对称变形的非线性方程。在此基础上用摄动法,给出了复合材料夹层圆柱壳轴压轴对称屈曲方程,该方程与齐次边界条件一起构成分析轴对称失稳的特征值问题。通过详细分析得到了该特征值问题的解,无量纲的轴压临界值表达成无量纲的抗弯刚度参数,夹芯抗剪刚度参数和Badtoff几何参数的简单函数。通过这些无量纲参数,轴压临界值实际上是复合材料工程弹性常数,夹心横向剪切模量,各层纤维铺设角,铺层数,夹心相对于表层的厚度和圆柱壳参考中面几何尺寸的函数。     

NUMERICAL STUDY OF TWO INTERFACE METHODS FOR COMPRESSIBLE MULTI-FLUID FLOWS

INTRODUCTIONThe direct numerical simulation of com-pressible multi- fluid flows is an important re-search topic with various key applications.Therelatively dominant difficulty of simulation lies intreatment of the moving material interface. Ingeneral,there are two basic approaches to treatcontact discontinuities.The first,front trackingmethod[1,2 ] ,explicitly characterize the interfacemotion by the motion of individual points on theinterface. The other approach is front capturingmethod[3…