叶轮机械流道中粘性可压缩流动的数值计算——(Ⅱ)旋转叶列中轴对称流动的计算

本文把第一部分中所提出的计算叶轮机械流道中粘性可压缩定常层流问题的数值方法应用到旋转坐标系中。首先导得旋转叶列中粘性可压缩流动所满足的基本方程组,然后通过在任意非正交曲线方向上的焓梯度方程、能量方程和熵方程之间的迭代计算,得到叶轮机械旋转叶列流道中粘性轴对称流动的数值解。     

计算二维粘性流动的流线迭代法

本文提出一种利用流线迭代法来计算任意形状流道中定常粘性层流流动的数值方法。通过任意非正交曲线网格中的压力梯度方程、能量方程和熵方程之间的迭代计算,可以得到整个流道中定常粘性可压缩(或不可压缩)流动的数值解。本文导出了二维(包括轴对称)流道中的基本方程,并详细地叙述了本方法的计算步骤。利用本方法对一些流道进行了数值计算,计算结果与其他巳知的数值解符合得很好。     

计算粘性流动的半隐式流线迭代法

引言对于任意形状流道中的定常粘性层流问题,作者曾在文献中提出了一种流线迭代法来进行数值计算。由于流线迭代法是在任意非正交网格中进行有很差分计算,它可以根据流动区域的特点来画出不等距的非正交网格,从而可以在保证精度的前提下减少计算工作量。正因为流线迭代具有方法简便,计算速度快等优点,这种方法可以广泛地应用于工程计算之中,作者曾用流线迭代法计算了叶轮机械中的粘性流动。(如文献)     

船用燃气轮机进气滤清器惯性级内的流场计算和实验验证

本文提出了一套求解船用燃气轮机进气滤清器流道流场的数值方法,成功地计算了流道内流场的状态,给出了各种不同型号流道的气动特性,对指导滤清器的设计有较大的现实意义, 在这套方法中,我们应用上风差分来逼近二维、非定常、粘性、不可压缩流体非守恒型的N-S方程,提出了一种可计算雷诺数高达上万的粘性流的差分方程,考察了这种差分方程的稳定性,收敛性、精度和人工粘性,本文还提出了处理一些边界拐点处壁涡的计算方法,实际算例表明,使用本文提出的差分方程和壁涡处理方法给出的计算结果和实验吻合良好。     

奇点法在叶轮机械任意回转面亚音速可压缩流动计算中的应用

本文在[3]的基础上,进一步将传统奇点法的应用范围从计算不可压缩流动推广到叶轮机械任意回转面亚音速可压缩流动的计算,并在数值处理方面做了有益的改进.算例和设计实例表明本方法计算可靠并具有实用性,为叶轮机械准三元流动计算中S_1和S_2两类流面间的反复迭代提供了一个有力的计算工具.     

二维定常不可压缩粘性流动N-S方程的数值流形方法

将流形方法应用于定常不可压缩粘性流动N-S方程的直接数值求解,建立基于Galerkin加权余量法的N-S方程数值流形格式,有限覆盖系统采用混合覆盖形式,即速度分量取1阶和压力取0阶多项式覆盖函数,非线性流形方程组采用直接线性化交替迭代方法和Nowton-Raphson迭代方法进行求解.将混合覆盖的四节点矩形流形单元用于阶梯流和方腔驱动流动的数值算例,以较少单元获得的数值解与经典数值解十分吻合.数值实验证明,流形方法是求解定常不可压缩粘性流动N-S方程有效的高精度数值方法.     

黏弹性流体扩展及收缩流动的格子Boltzmann模拟分析

本文应用格子Boltzmann 方法(LBM)并结合Oldroyd-B 模型,讨论了不可压缩的 Navier-Stokes 方程和平流扩散本构方程的解耦及各自求解方法,以及两类问题的边界处理格式,实现了黏弹性流体在二维1:3 扩展流道以及3:1 收缩流道中的流动的数值模拟．获得了不同雷诺数Re 和维森伯格数Wi 以及黏度vs 下流动的流线分布,计算给出了漩涡的涡心位置和大小,并分析了参数Re、Wi 和vs 对流动特点的影响．模拟结果表明本文所采用模型和边界处理方法具有良好的精度和稳定性．     

高精度迎风紧致差分格式与热流计算研究的注记

利用高精度差分格式求解了可压缩 N-S方程球头热流问题。分析了不同差分格式在对球头粘性绕流热流计算中存在的问题 ,并分析了相应的网格雷诺数。在利用高精度迎风紧致 [1 ] 格式求解粘性绕流热流问题时 ,采用 Steger-Warming[2 ]的通量分裂技术将守恒型方程中的流通向量分裂成两部分 ,在此基础上据风向构造逼近于无粘项的高精度迎风格式。对方程中的粘性部分采用中心差分格式。数值结果表明 :高精度差分格式能在较大的网格雷诺数下较好地计算球头驻点热流     

轴对称二维定常可压缩流体湍流边界层的计算

本文根据应力张量表示的动量方程,推导出了轴对称坐标系中无旋转定常轴对称流动可压缩湍流边界层的动量方程.可以看出,如果流道内流动要全部按粘性流体来处理,那末流道的曲率,以及曲率沿流动方向的变化率都不能太大.这一点在已有文献中都没有论证. 根据文中提出的方程组编制了一个计算程序并计算了一个例子.     

可压缩多介质粘性流体的数值计算

将考虑热传导和粘性情况下的Navier Stokes方程描述的物理过程分解成3个子过程进行数值计算,即把整个流量计算分解成无粘性流量、粘性流量和热流量3部分,采用多介质流体高精度parabolic piecewise method（PPM）方法、二阶空间中心差方法和两步Rung-Kutta时间推进方法相结合进行数值计算。给出了激波管中Riemann问题和二维、三维Richtmyer-Meshkov界面不稳定性的Navier Stokes方程和Euler方程对比计算结果,显示了粘性对界面不稳定性的影响。     

水泵叶片表面三维边界层计算

本文用有限差分法计算混流式可逆水力机械水泵工况叶片表面的三维边界层。水泵叶轮中主流区的三维势流由直接边界元法计算。对于叶片面附近的粘性流动。用三维半正交坐标系中的边界层方程表示。为了提高计算精度采用贴体坐标技术生成边界层区域的计算网格。并利用Cebeci等变换函数及Keller差分格式离散方程。用分块解法求解。计算叶轮叶片表面的压力分布与相应试验结果进行了对比。     

基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

一种用边界元方法求解二维不可压缩粘性流动的方法

本文用边界元方法求解了二维不可压缩粘性流动的涡量——速度方程,利用求解区域边界上的速度法向导数和速度值直接得到了涡量的边界条件,克服了利用涡量方程求解二维不可压缩粘性流动时涡量边界条件(主要是壁面边界条件)难以给定的困难,算例表明:这种方法比较简单且计算结果基本上是令人满意的。     

用迭代法模拟注塑成型的三维流动

Hele-Shaw模型是模拟注塑成型过程的常用模型,它的主要缺点是不能模拟一些重要的物理现象及引入实际制件中并不存在的“中面”概念,为了消除这些不足,本文开发了真三维的流动分析程序。建立了粘性、不可压缩的非牛顿流体流动的控制方程,为了避免同时求解耦合的压力场、速度场,本文引入拟稳定技术独立地求解这些方程,用迭代法求耦合方程的解。这种方法可以减少内存并提高数值方法的稳定性。算例表明数值结果与实验结果吻合较好,这种方法成功地模拟了注塑成型流动过程中的重要特征。     

用流线迭代法计算同时含有内流及外流场的粘性流动问题

本文提出一个通过坐标变换把外部流场问题转化为内部流场问题来处理的方法。从而可将计算内部流场粘性流动的流线迭代法推广应用到外部流场中。用这种方法,本文对同时含有内流场和外流场的粘性流动问题进行数值计算。最后给出了计算实例。     

粘性不可压缩流体定常流动的差分解法

本文研究了定常N-S型方程和压力泊松方程的耦合求解。提出了一种处理压力自洽边界条件的方法,结合文[7]中给出的自调差分格式,可以对一些较复杂的粘性不可压缩流进行数值求解。     

计算流体动力学的现状和展望

本文综述了计算流体动力学的进展,侧重物理问题(而不是方法)。首先,讨论了描述质量、动量和能量守恒的流体动力学控制方程。其次,在接下去的两节中叙述了无粘流和边界层流的近期计算工作。最后,结合近期文献中的几个例子阐述了数值求解Navier-Stokes方程的粘性流计算。     

剧变截面圆管内渗流的数值计算方法

对于剧变截面圆管的渗流问题写出不可压缩渗流的基本方程组,对直接求解原始变量(速度和压力)的数值计算方法作出改进。先由非主流方向的运动方程计算压力,后由主流方向的运动方程计算主流方向的速度分量,再由连续性方程计算非主流方向的速度分量。这样可以避免在一般的求解原始变量方法中由连续性方程计算压力时出现的困难和麻烦。根据本方法和剧变截面圆管的特点,采用半交错不等距非正交贴体混合网格系。本文详细写出差分方程和迭代计算公式,对剧变截面圆管内的渗流算例进行数值计算。本方法的优点是简单和实用,在工程上具有较大的应用价值。     

可压缩、粘性翼前缘分离涡研究

从N-S方程出发,推得用于描述飞行器分离涡的可压缩、粘性涡核方程,给出了对于亚、跨、超音速都适用的数值解法,它可以计算分离涡或尾迹涡特性。导出了可压缩分离涡的破裂条件,并用于计算破裂位置。得到和分析了压缩性等参数对分离涡运动影响的一系列感兴趣特性。     

时间推进方法在叶轮机械内部流场计算中的进展

时间推进方法目前广泛应用于叶轮机械内部流场的计算,这种方法可以划分成显式时间推进方法和隐式时间推进方法．本文简要介绍了时间推进方法用于叶轮机械内部流场计算的起源及这种方法的优点和缺点,重点记述了国内外叶轮机界应用显式时间推进方法和隐式时间推进方法计算叶轮机内部流场的最新进展,并简要介绍世界上一些著名的研究机构开发的时间推进方法计算程序所采用的数值计算方法