振荡三角翼非定常气动特性的数值模拟

以欧拉方程为基础，采用有限体积方法，建立了用于进行非定常流动模拟研究的计算方法。文中引用了变系数的残值光顺方法，在保证原求解方程二阶精度的前提下，使得计算效率提高了近１０倍。利用所建立的方法对绕尖锐前缘三角翼作大振幅振荡的非定常旋涡流动进行了模拟计算，收到了良好的计算效果。     

不可压粘性流动N—S耦合方程的有限元解法—1理论

本文对不可压粘性流动非定常N－S方程组在耦合条件下应用有限元方法直接求解，并对在高雷诺数下如何提高计算的精度，稳定性以及收敛速度等方面进行了讨论，从而建立起不可压粘性流动N－S耦合方程的有限元解法。     

幂律流体平面渠道层流进口段流动阻力的分析

将Ｂ．Ｂ方法推广到幂律流体平面渠道间的流动，首先导出幂律流体平面渠道层流边界的能量积分方程，通过对能量积分方程的求解，对幂律流体平面渠道的进口段长度，以及进口段效应压力损失和流量修正系数进行了分析和计算，并且讨论了流动指数ｎ对进口段效应诸系数的影响，所提供的压力损失系数计算值，在ｎ＝０．７５时和在实验台所做的实验值进行了比较，二者符合得较好，因此验证了结果的可靠性。     

一种二维不可压粘性流动的求解方法

本文采用边界元法求解二维不可压粘性流动的Ω-ψ方程。通过适当的数学推导,导出了关于流函数ψ。涡量ω的边界积分方程,并把它们在区域边界上离散,求解。二维扩张管道的计算结果表明,这是一种较优越的求解二维粘性流动的方法。     

用改进的强隐式方法求解三维不可压粘性流动问题

提出了一种用改进的强隐式方法求解三维不可压缩粘性流动的交错网格有限体积法。该方法在ＳＩＭＰＥＲ算法的基础上，对Ｎ－Ｓ方程的对流项采用ＱＵＩＣＫ格式离散，利用ＭＳＩ方法求解代数方程，数值计算表明本方法具有计算稳定，收敛速度快等特点。     

微型直升机旋翼的不可压Euler方程数值计算

采用拟压缩方法求解低速不可压流动具有效率高、稳定性好的优点。通过求解拟压缩性修正后的Euler方程成功地模拟了悬停时微型直升机旋翼的绕流。计算网格采用O—H拓扑形式的结构网格，空间离散采用中心有限体积法，时间推进为五步Runge—Kutta方法。采用了当地时间步长、压强阻尼及隐式残值平均方法来加速收敛。应用此方法对低速条件下旋翼进行了数值模拟，计算结果与实验数据吻合得较好，经气动分析，证明拟压缩性方法运用于微型直升机旋翼完全可行。     

液体射流与空气交叉流动喷雾的数值模拟

采用气相的时均雷诺方程,在拉格朗日框架下求解离散相的运动方程,对水与空气交叉流动的喷雾流场进行了数值模拟研究.湍流模型为标准κ-ε模型,通过在气相方程中添加源项的方法,对两相间存在的质量、动量和能量交换过程进行了模拟;还添加了描述液滴的破碎及其碰撞聚合过程的计算模型;参考实验中的实测入口风速大小,添加用户自定义函数定义气相入口边界风速.计算时,交替求解气相和离散相方程直到两相都达到收敛来实现两相间的双向耦合.并将计算网格分割为两部分,在拥有两个处理器的工作站上并行运行.模拟结果与实验结果的对比显示,喷雾轨迹和液滴粒径大小与实验结果符合较好,也较合理地模拟出了液滴的空间分布情况.     

积分方程法计算三维磁场

本文研究了用积分方程法计算三维磁场的基本原理，列出了第二类Ｆｒｅｄｈｏｌｍ积分方程根据这些方程，采用离散方式在整个求解场域进行求解，另外还推导了耦合系数的解析式，最后应用本文方法对大型电力变压器大电流引线产生的漏磁场进行了计算，计算结果是令满意的。     

复杂边界水流高分辨率数值模拟

将基于流向量分裂技术的两步有限体积分量方向TVD格式应用于求解浅水方程的控制方程,建立了高分辨率的计算模型.实际算例表明,该模型能较好地计算浅水流动问题,为浅水流动数值模拟提供了精度高、稳定性好的数值计算模型.     

有限元方法在注塑流动压力场求解中的应用

:在分析注塑充填流动压力方程的基础上 ,采用有限元方法求解充填过程中非定常、具有移动边界的压力场 .提出的总体刚度阵合成方法与熔体前沿推进技术的结合可为一般的具有可变计算区域的有限元数值计算提供参考     

基于欧拉方程的一种机翼气动弹性计算方法

利用一种双时间方法求解三维非定常欧拉方程，采用无限插值理论生成O－H型代数网格，考虑了机翼变形时的网格生成问题，通过与气动力方程的联立求解，在时间域内用二阶龙－库塔方法求解机翼弹性运动方程。计算结果表明，本计算方法具有较高的计算效率，所计算的颤振临界速度与风洞实验一致。     

三维叶栅流道中不可压势流的边界元计算程序

本文采用边界元法求解三维无粘不可压流动问题。通过详细的数学推导,得到了三维边界元法在区域边界上数值离散后的代数方程,并以此方程编制了计算机程序。三维管道的计算结果表明,本文的尝试是成功的。用边界元法求解叶轮机内三维流动问题是一种现实可行的新方法。     

SIMPLE算法中压力松驰因子的自动生成

ＳＩＭＰＬＥ算法数值求解流体流动及对流换热问题的难点之一就是压力松驰因子的选取。基于使离散的动量方程残差的平方和最小的原理，建立了压力松驰因子的自动生成方法，程序编制时采用了三种具体的处理方法。对三个典型的算例计算表明，第三种方法最可靠、收敛速度快。     

SIMPLE算法中压力松驰因子的自动生成

ＳＩＭＰＬＥ算法数值求解流体流动及对流换热问题的难点之一就是压力松驰因子的选取。基于使离散的动量方程残差的平方和最小的原理，建立了压力松驰因子的自动生成方法，程序编制时采用了三种具体的处理方法。对三个典型的算法计算表明，第三种方法最可靠、收敛速度快。     

粘性不可压缩流体平面流动的有限差分解法

本文使用上风差分格式求解了二维非定常粘性不可压缩流体非守恒的Navier-Stokes方程,讨论了相应的差分方程的精度、收敛性和人工粘性.发现当雷诺数增大时,使用非守恒型方程比守恒型方程易于收敛、更为适宜.本文还对壁面拐角上点的涡量计算提出了一种新的方法.以曲折流道内粘性流体的流动为例进行了数值计算,结果表明与实验结果具有较好的一致性.     

非线性双曲型Euler方程的数值分析

采用LU-AUSMLW算法,对非线性双曲型Euler方程进行求解。该算法综合了LU-SGS方法与AUSMPW格式的优点。为了验证LU－AUSMLW算法的有效性,在曲线坐标系中,对亚音速、跨音速和超音速流动进行了数值模拟。数值试验的计算结果与文献结果相符很好。     

壁面凸体层流绕流流场的数值计算

采用贴体坐标变换方法将直角坐标系中不规则的流动区域变换成规则的计算区域．对求解原始变量（ｕｋ，ｐ）的ＳＩＭＰＬＥＲ方法进行了改进，应用有限分析法给出了由正常网格求解原始变量的选代方法．然后计算了台阶绕流和任意凸体绕流的流速、压力场．计算结果表明该方法简单有效，压力、速度收敛效果好，计算精度较高．可以用于计算各种流动区域的流场．     

一阶迎风差分格式求解非线性对流扩散方程的精度

采用一阶迎风格式分别对一维线性对流扩散方程和非线性对流扩散方程进行了求解,检验了一阶迎风格式用于求解一维线性对流扩散方程和一维非线性对流扩散方程的适用性.多个计算算例的结果表明:一阶迎风差分格式用于求解线性对流扩散方程的结果不甚理想,但用于求解非线性对流扩散方程时能获得相当精度.工程计算中,该格式可用于求解水流运动方程,但不宜用于求解被水流输移的物质对流扩散方程.     

超燃冲压发动机一维跨声速流动奇异初值问题

超燃冲压发动机一维跨声速流动方程在临界声速点存在奇异问题,现有的基于L’Hosp ital法则的求解方法在进行跨声速流动计算时可能出现较大的数值误差.从数值稳定性和误差特性两个方面对该算法进行了分析,结果表明问题主要在于方程的奇异初值误差,需要针对这一问题提出改进算法.     

弹性地基梁的数值计算方法

文章在半无限性体假设下的弹性地基梁道本方程的基础上，提出了一种便于计算机计算的解法，该解法将弹性地基反力进行分段性插值，利用积分变换，将弹性地基梁基本方程的求解变换成ｎ＋２阶的线性方程组的求解，该解满足所有基本方程。