关于三维Navier—Stokes方程的粘性项计算

借助于张量分析和张量计算，在贴体曲线坐标系下本文讨论了不同求解变量导致了粘性项个数上的重大差异和不同大小的计算量，并提出了便于粘性计算的最佳形式。文中借助于有限体积离散技术，通过引进两个对称辅助矩阵[A]和[B]，使粘性项的计算量大大减少，这对完成三维粘性流的数值计算具有重要的指导意义。借助于上述方法，本文完成了某型真实进气道两种工况的三维粘性Navier-Stokes方程计算（即M∞=3.0,a=0°和设计工况M∞=2.65,a=0°），获得了满意的全场结果；对于M∞=2.65的设计工况。同实验数据作了比较，符合良好。由于本文的方法明显的减少了粘性项的计算量且节省了大量内存，以致于使三维流场的N-S求解能在普通微机上进行。     

Navier—Stokes方程粘性项的有限体积离散法研究

Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程粘性项的有限体积离散法研究封建湖西北工业大学应用数学系［摘要］：三维Ｎ－Ｓ方程是描述真实气体流动的非线性偏微分方程组，如何快速准确地计算其粘性项至关重要。这不仅因为粘性项的计算量非常巨大，而且粘性项对急剧变化的流动梯度的...     

粘性流场三维NS方程有限体积法数值解

提出了有限体积法中的隐式ＮＮＤ格式数值求解三维全ＮＳ方程的薄层三维ＮＳ方程的方法，对绕带座舱的再入飞行器高超声速粘性流动进行数值模拟。     

粘性三维流动计算的任意拉格朗日—欧拉法

推导了旋转叶轮内流三维粘性流动计算的ALE格式。以网格的虚拟运动计算对流作用;采用粘性子循环时步、拉格朗日子循环时步加速收敛;改进上风格式为部分施主单元格式来减小扩散同时实现良好的计算性能;由通量限制控制耗散截断误差;计算中引入短波长阻尼以消除不规则性;以亚网格(SGS)湍流模式模拟湍流运动。      

基于粘性笛卡尔网格的N-S方程计算

本文采用基于叉树数据结构的粘性笛卡尔网格求解N-S方程。为了使网格的分布更为合理,采用基于物体几何外形的方法对网格进行自适应加密。在物面处采用光顺、投影、分割加密的方法保证了物面处网格的正交性与贴体性,从而生成高质量粘性笛卡尔网格。将有限体积法应用于粘性笛卡尔网格求解N-S方程,采用了SA紊流模型。算例计算结果与实验数据符合良好。     

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

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

基于Navier Stokes数值计算的风洞收缩段设计

运用Navier—Stokes数值模拟对不同收缩曲线和设计方案的风洞收缩段和试验段的流动进行模拟，并对收缩段的分离特性、出口速度均匀性以及试验段的速度均匀性和附面层厚度作出评价。计算格式在空间上采用中心有限体积离散，在时间上采用多步Runge-Kutta时间步长格式进行积分。     

用非线性涡粘性模式计算三维湍流边界层

针对非线性涡粘性模式在求解三维湍流边界层流动时的不足，本文从压力－变形率关联项中的快速项出发，考虑流动非均匀影响，在显式代数应力模式的雷诺应力表达式中引入了反映雷诺应力“松驰”效应的速度二阶导数项，构造了一个新的非线性涡粘性模式。通过对典型算例－翼体角隔流动的计算结果表明，新模式能较好地再现出三维湍流边界层内雷诺切应力方向的发展滞后于速度梯度方向发展这一流动特性。     

粘性超声速绕体流动的有限谱法计算

首先推导出粘性可压缩流体动力学基本方程组在广义坐标下矢通量分裂的一般形式，并将有限谱ENO格式用于求解超声速流动问题。通过典型算例，验证了这一格式的精度和稳定性。然后进行了圆柱和机翼超声速绕流的数值模拟，得到了与实验和其他计算相符的结果。这表明有限谱法不单具有精度高，应用灵活等多方面的优点，而且突破了传统谱方法全域性的局限，可以成功地应用到一般形状流场可压缩流动的数值计算中。     

叶轮机械三维粘性动静叶干涉的数值模拟

本文发展了 1套模拟叶轮机械三维粘性动静叶干涉的数值方法。在成熟的空间离散方法基础上 ,提出了 1种高效的双时间步隐式时间推进方法 ,并构造了一种保证守恒的动静叶交界面处的数值通量计算方法。本文对某三维亚音速涡轮的动静叶干涉进行了数值模拟 ,结果与实验符合较好。     

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

超声速绕凹角流动的压力分裂形式简化N－S方程数值解

本文从压力分裂形式的简化Ｎ－Ｓ方程出发，用Ｒｕｂｉｎ格式对超声速粘性流动进行了数值求解，给出了超声速绕凹角附着流及带有小分离气泡流动的结果，与其它数值方法所得结果相比，符合较好。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

某型串列静子的三维气动数值优化设计

针对某型高负荷轴流压气机串列静子近轮毂处存在的分离,利用自行开发的压气机三维气动数值优化设计平台,选择复合形优化方法,通过修改串列静子近轮毂处的三维积叠方式对串列静子进行了优化设计.结果表明,优化设计取得了良好效果,在消除分离的同时,有效提高了压气机的整体性能.      

三级风扇的三维流场计算

本文运用三维粘性流体计算软件Fine／NUMECA对某三级风扇多个转速下的三维流场进行了计算．并对其第一级三个典型转速进行了单独计算。分析了此风扇在不同转速下的特性变化,重点关注了设计点和近喘振边界工作点的流场特性,并分析了在设计转速时发生喘振的原因。结果表明,此风扇设计为典型的跨声速压气机设计．具有较高的效率,但仍有改进的余地。     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

明渠水流三维数值模型研究

本文在具有无尺度铅垂坐标的直角坐标系中采用完全的三维雷诺方程组描述明渠中的水流运动,利用二维沿水深平均的连续方程确定自由表面形状。文中给出的隐式求解方法保持了经典SIMPLE算法的全部优点,并具有高的收敛速度。算例表明,计算结果与实测值十分一致。     

跨声速轴流压气机三维黏性流场的数值模拟

改进了动静叶级间的掺混面处理方法,提出了用参数分布修正方法来保证参数在掺混面上的周向不均匀性,并使用亚松弛的参数分布系数来解决程序在动静叶轴向间距较小时的计算稳定性和收敛性问题。采用经改进的NAPA软件对NASA37级压气机的设计点和在不同转速下的流场进行了数值模拟,计算结果与试验结果吻合较好。     

复杂地形低空三维风场数值仿真方法

设计了飞机超低空飞行的复杂山地地形,基于球绕流理论,建立了低空复杂地形及其上空三维过山气流的工程化模拟方法,并进行了算例仿真和分析.该方法可以提供低空三维风场模型,并能满足分析飞机低空飞行动力学特性的需要.