不可压缩流动的TG有限元模拟

采用三阶Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ（ＴＧ）有限元法离散Ｎ－Ｓ方程，并对离散格式中由高阶时间导数转化而来的附加项进行了缩并运算修正，使其起到了人工粘性的作用。由于该法隐含了流线迎风的耗散作用，可较好地应用于较高Ｒｅｙｎｏｌｄｓ数流动的计算。同时采用压力校正法求解不可压缩流场的各原始变量，并讨论了压力校正Ｐｏｉｓｓｏｎ方程的边界条件；最后给出了二维方腔拖带流的计算结果，证明了算法的正确性。     

可压缩N-S方程的特征有限元方法数值模拟

讨论了二维或三维可压缩N-S方程的特征有限元方法数值模拟,严格的理论分析表明这种方法是稳定的,并且具有最优阶误差估计.     

计算可压缩流体流动的LB多速度模型

建立了用于计算可压缩流体流动的三速度六角形ＬａｔｔｉｃｅＢｏｌｔｚｍａｎｎ模型，此模型要求粒子碰撞过程不仅满足质量和动量守恒，而且满足能量守恒．确定了此模型平衡态的形式，导出了此模型对应的宏观质量、动量和能量方程．此模型可用于计算具有温度变化的空气动力学问题．给出了一维压力波传播问题的计算结果．     

可压缩混合层流动的三维大涡模拟

为了加深对复杂流动混合现象的认识,为化工设备设计提供理论基础,用大涡模拟方法对可压缩气体混合层进行了三维数值模拟.计算结果显示了混合层流动的三维瞬态发展过程以及流动过程中涡的产生、生长、破裂以及涡列的形成与各涡的合并过程,揭示了可压混合层三维演变过程中同样具有自相似性,计算结果与实验结果相符.通过对混合层初始状态条件下不同扰动和不同初始速度比的计算表明:随机扰动量对混合层的失稳演变过程影响最大,而大混合层速度比同样可使混合层失稳加剧并提前.     

基于自适应笛卡尔网格的可压缩粘性流动数值模拟

基于自适应笛卡尔网格方法求解Navier-Stokes方程,网格以四叉树数据结构存储,固壁边界条件通过一种虚拟单元体方法引入。在几何外形确定的前提下自动完成网格生成、加密和流场的求解任务。对含有激波的NACA0012翼型的超音速绕流工况和含有回流区的双NACA0012翼型亚音速绕流工况进行了数值模拟,并与现有的非等距笛卡尔网格解和非结构网格解进行了对比。结果表明：基于自适应笛卡尔网格能够准确模拟可压缩粘性流动,同非等距笛卡尔网格相比,自适应技术的使用显著降低了网格量,但是同非结构网格相比,现有的自适应笛卡尔网格技术在边界层的分辨上效率较低,有待进一步发展。     

LU-AUSM~ 混合格式在三维可压缩流动数值分析中的应用

将一维AUSM 格式扩展到三维任意曲线坐标 .为了获得高阶精度 ,采用了三阶MUSCL格式 .将AUSM 格式与LU SGS格式结合 ,对凸包通道跨音速无粘流动、平面叶栅跨音速和超音速无粘流动 ,以及喷管超音速粘性流动进行了数值模拟 .计算结果与文献计算结果和试验数据相符很好 .     

用于模拟弱可压缩Navier-Stokes方程的格子气模型

给出压缩因子的形式．通过引入鬼场，消除压缩因子，得到了平衡态分布函数的一般表达式及可压缩的等温的Navier－Stokes方程，同时给出了自映射条件所要求的平衡态分布函数的范围．     

叶栅可压缩流动反问题的一种数值解法

本文提出一种解叶栅绕流反问题的近似方法.利用假想气体近似,将可压缩流动平面的反问题,变换为虚拟的不可压流动平面上的反问题.采用边界元方法造代求解出叶片型线,然后再变换到原物理平面.本文还对亚音速、无激波跨音速叶栅的三个例子作了数值计算.     

低雷诺数下翼型粘性绕流主动流动控制的数值模拟

通过数值求解Navier-Stokes方程模拟零质量射流对低雷诺数翼型绕流的控制作用.空间离散采用中心有限体积格式,时间推进为双时间推进方法.计算结果表明,采用零质量射流进行主动流动控制能有效抑制大攻角下大尺度的流动分离,改善翼型的气动特性 ,起到显著的增升减阻作用.     

解不可压缩流动N-S方程的隐式SMAC方法

该文基于SMAC(SimplifiedMarkerandCell)方法推导出了的一种直接求解不可压缩N-S方程的隐式数值方法。求解的基本方程是任意曲线坐标系下以逆变速度为变量的N-S方程和椭圆型的压力Poisson方程。压力Poisson方程用TschebyscheffSLOR方法交替方向迭代求解。N-S方程数值离散时对流项采用了Chakravaythy-OsherTVD格式。用该方法计算后台阶流场的结果与经典的实验结果相当吻合,表明该方法是可靠的,在合适的边界条件下求解不但是稳定的,而且能有效抑制流网扭曲大的地方产生较大的非物理振荡误差。     

水轮机蜗壳和固定导叶内部流动的数值模拟

为探求蜗壳内部的流动规律 ,优化蜗壳设计 ,对三峡混流式水轮机组模型蜗壳和固定导叶内部的流动进行了紊流解和 Euler解的全三维数值模拟。文中基于 Navier-Stokes方程的标准 k-ε紊流模型、Euler方程模型 ,采用贴体座标和交错网格系统 ,用 SIMPL EC算法 ,对蜗壳和固定导叶内部流动的联合计算结果进行对比分析 ,并与实验值进行了比较。还对流动的不均匀性进行分析 ,提出了改进和优化蜗壳和固定导叶设计的较为可信的解决方案     

多个转子的管道中黏性不可压缩流体的流动分析

研究了一类流固耦合问题,即带有多个转子的管道中不可压缩流体的流动问题。为了处理流体刚体之间的边界困难,将转子近似为黏性足够大的流体,从而把流固耦合问题转化为流体问题并得到近似解,然后通过黏性和运动方程确定出刚体的位置,从而将刚体与流体区域分开,最终得到该问题全局弱解的存在性。     

LU—AUSM^＋混合格式在三维可压缩流动数值分析中的应用

将一维ＡＵＳＭ＾＋格式扩展到三维任意曲线坐标。为了获得高阶精度,采用了三阶ＭＵＳＣＩ格式。将ＡＵＳＭ＾＋格式与ＬＵ－ＳＧＳ格式结合,对凸包通道跨音速无粘流动、平面叶栅跨音速和超音速无粘流动,以及喷管超音速粘性流动进行了数值模拟。计算结果与文献计算结果和试验数据相符很好。     

用分步有限元法求解三维不可压缩流动

将计算二维不可压缩流动的分步有限元格式扩展到三维情况 ,由于该格式没有引入新的高阶空间导数项 ,适用于多维空间的非线性问题。实际求解了三维空腔流 ,在 Re数较低的情况下得到的流态不随时间而变化 ;在高 Re数(Re=40 0 0 )的情况下 ,流态十分复杂并且是非定常的 ,在水平断面具有多个非定常旋涡 ,在垂直于流动方向的立面断面上可以模拟到 TGL涡。本格式的迎风效应是 Taylor展开式的高阶精度项 ,没有人工粘性引入 ,得到的结果可靠性好     

用扰动方程研究可压缩边界层中扰动的演化

本文对传统的可压缩流动中的通量分裂方法进行了修正，使其可以用于扰动方程的计算。利用这一方法，对来流马赫数Ma为4．5、雷诺数Re为10^5的可压缩平板边界层中扰动的演化进行了数值模拟。对于给定的基本流，在计算域入口加入T-S波扰动，研究扰动的空间发展演化。对于小幅值的扰动，计算得到的扰动幅值演化和扰动法向分布与线性稳定性理论的结果符合的很好。对于有限幅值的扰动，结果表明，各次谐波无论幅值还是扰动分布都与完整的Navier-Stocks方程计算的结果符合的很好，但由扰动方程得到的平均流修正比N-S方程得到的大一些，这也许是由于在N-S方程计算时，基本流与零阶谐波一起计算而产生的误差所至。     

可压缩气体合成喷射流数值模拟

为了更深入地研究合成喷射流性能与激励器结构参数的关系,在考虑气体可压缩性质的条件下,对合成喷射流进行了理论分析,得出了合成喷射流性能与激励器结构参数的函数关系式,并采用动边界条件模拟实际压电膜的振动,对合成喷外部射流场和激励器腔体内部射流场进行了联合数值模拟.计算结果与已有的实验结果及理论分析结果吻合较好,验证了该数值模拟方法的可行性.     

扰流片定常和非定常分离流动的数值模拟

扰流片的定常绕流和它作垂直方向振荡的非定常绕流,其流动特性有着本质的区别。本文利用有限差分方法求解二维不可压缩的定常和非定常Navier-Stokes方程。系统地研究了这两类流动问题。考虑了Reynolds数、扰流片高度和宽度、来流速度分布等因素对流动的影响。对于非定常绕流,还研究了振荡频率的影响。并对对流项采用一阶迎风、混合和幂函数等不同差分格式做了计算和比较。