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

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

方形截面弯管内的二次流动及其分叉

[摘要]基于曲线坐标系内的流动模型及其数学描述,采用Galerkin方法对方形截面弯曲管道内充分发展的稳态不可压层流流动进行数值模拟,分析管道曲率参数,c及流动条件Re数对管内二次流动特性的影响;给出管内二次流动的分叉状态图,并精确地求得出现四涡解的分叉点．     

恒定有压扩散流的大涡模拟

采用大涡模拟（LES）方法对恒定有压扩散流进行了数值模拟，对大尺度漩涡直接求解Navier-Stokes方程，小尺度漩涡采用标准Smagorinsky亚格子模型模拟．针对不同雷诺数下扩散段内的水流进行数值模拟，得到了各雷诺数流动时的时均速度场，分析了不同流态下扩散段内恒定有压扩散流时均流动特性．通过对不同流态下扩散段内连续瞬态流场的分析，探讨了不同雷诺数流动时扩散段内瞬态流场结构特点以及分离区内漩涡的变化规律，指出在恒定流状态下扩散段内存在流动的局部不稳态现象（主流的非周期性摆动），并对此现象进行分析．计算结果和试验结果吻合较好．     

应用滑移界面方法求解动／静相干 叶排内非定常流动

基于滑移界面技术，发展了求解二维动／静相干叶排内非定常粘性流动的计算方法和软件，该方法克服了混合平面方法所引入的方法误差．压力与速度耦合采用非迭代PISO算法；计算区域的离散采用非交错斜交网格．计算了二维不可压缩Couette流动及无限长相对运动动／静叶排内的非定常流动，验证了用该方法和软件计算动／静相干叶排内非定常流动的有效性．     

松散颗粒介质水的渗透率与孔隙率实验研究

以松散颗粒介质渗透率实验出发点，考证多孔介质内流动滑移的同时，将此实验作为考察微孔道流动的方法提出，避开单一微槽道流动测量的困难．实验结果表明，大于0．1mm的孔道内，观测不到水的偏常规层流流动．     

应用控制容积有限元法分析注射充模过程中塑料熔体的流动与传热

将Ｂａｌｉｇａ和Ｐａｔａｎｋａｒ提出的控制容积有限元法推广用于分析注射充模过程中熔体的流动与传热，并叙述了其特点．应用此方法计算分析了包括合不同几何形状嵌件的模腔内充模过程的流动与传热问题．得到了一些有益的结论．这将有利于指导在复杂型腔充模过程流动分析时的简化处理．     

基于泵内纸浆悬浮液数值计算的纸浆泵设计

在由吸水室、前盖板、半开叶轮和涡壳组成的过流部件内，对纸浆悬浮液的流动进行了三维不可压两相湍流流动计算．计算采用双流体模型及κ-ε湍流模型，Phase Coupled SIMPLE算法进行压力速度耦合计算．根据计算结果，对两相流动的速度、压力及体积分数进行了对比和分析研究．发现在叶轮叶片数较少的情况下，叶轮流道内速度变化不均匀．并着重对前盖板与叶轮间隙内的流动进行了分析，揭示了半开式叶轮输送两相流介质的优点所在．该计算结果为建立少叶片半开式叶轮纸浆泵的设计方法提供了理论依据．根据新的设计方法设计出的纸浆泵，效率经检测高出国家标准17．5％．     

基于多孔介质模型的刷式密封泄漏流动特性研究

采用数值求解基于多孔介质模型的Reynolds-averaged Navier-Stokes方程技术，对刷式密封内泄漏流动特性进行了详细的数值研究．根据刷式密封的泄漏量试验数据，确定了刷丝束多孔介质的渗透率系数．利用所确定的渗透率系数，数值计算和分析了7种压比和5种径向间隙条件下的某轴端刷式密封的泄漏量和泄漏流动形态，并且与迷宫式密封进行了比较．结果表明：在相同的压比下，刷式密封的泄漏量小于迷宫式密封；压比影响刷式密封的泄漏量，但对泄漏流动形态的影响可以忽略；径向间隙不仅影响刷式密封的泄漏量，而且影响泄漏流动形态；在相同的径向间隙下，压比越大则泄漏量越大；在相同的压比条件下，径向间隙越小则泄漏量越小．研究工作可对刷式密封在汽轮机中的应用提供理论依据和技术参考．     

微小矩形通道内气体充分发展层流流动的数值模拟

基于贴壁层概念和已经求得的贴壁层内气体粘性系数变化规律 ,建立了数理模型 ,求解了矩形微小通道内气体流速分布完全发展的层流流动 ,得到了不同截面高宽比 ξ和不同 Knudsen数下的量纲一速度分布和流动阻力因数 ,并对结果进行了分析 ,为进一步研究不可压气流在微小矩形道内的层流换热提供了分析基础     

离心泵蜗壳内部清水流动测量

利用ＬＤＶ分别测量了最优和小流量工况下矩形断面蜗壳内３个断面内的清水平均流动．实 验发现,蜗壳内部液体绝对速度的圆周分速度比径向分速度大一个数量级;蜗壳内存在螺旋运动; 蜗壳断面内流体角动量不守恒;无论在最优还是在小流量工况,蜗壳内部流动都是扩压流动;在小 流量工况,扩压更为严重．     

τ方法在无粘轴对称流稳定性研究中的应用

用谱方法中的τ方法研究了不可压无粘轴对称流的空间稳定性问题．以旋拧型（ｓｗｉｒｌｉｎｇｆｌｏｗ）尾涡流和射流（ｊｅｔ）为例，分析了不稳定特性，并与已有的结果作了比较，分析了τ方法的优缺点．     

汽轮机GX-1型调节阀流动特性的试验与数值研究

以汽轮机常用的GX-1型调节阀为研究对象，通过在阀座喉部、阀座渐扩段等阀体内各关键部位设置测点，利用微小型高频动态压力传感器及其采集系统进行多工况和多方位的试验研究，同时利用有限差分方法进行三维可压缩湍流数值模拟．通过理论分析、试验研究和数值模拟深入研究了阀内复杂的流动特性．结果表明，中小压比时阀内气流会出现流型反复交变，阀喉部附近流场最不均匀且不稳定，至使阀门的主要损失集中在喉部附近，而且可能引起气流的脉动．     

摇荡物体周围粘性流的数值模拟

提出了一种新的求解非定常不可压Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值方法，在每个时间步长上采用人工可压性技术，通过对动量方程取散度，从而获得离散化的Ｐｏｉｓｓｏｎ方程，当关于虚拟时间的迭代过程达到收敛时，便得到了无散速率场，采用该方法模拟了一个纵摇椭圆柱体周围的非定常不可压性流动。     

离心泵压水室内部流动数值模拟

基于标准的κ-ε双方程紊流模型,采用SIMPLEC算法,在贴体坐标系下,通过求解三维Navier—Stokes方程,对一离心泵压水室内部流动进行了数值模拟和分析,得到了内部流动的主要特征、研究表明：按常规水力设计方法设计的压水室,在压水室径向断面内存在明显的二次流．速度矩也有明显的梯度．这些流动现象可能是导致压水室损失较大的主要原因．     

水轮机蜗壳内流动的数值研究

在分析研究现有的水轮机蜗壳设计理论的基础上，对蜗壳在不同工况下的内部流动情况进行了三维不可压定常粘性流动的数值模拟。基于N－S方程和标准的κ－ε模型，采用贴体坐标、非结构化网格和SIMPLE算法，对蜗壳内部流动进行了计算分析，其结果与实验值吻合，可作为蜗壳优化和固定导叶设计、改型的依据。     

基于RNG k-ε湍流模型钝体绕流的数值模拟

采用RNG k-ε湍流模型对钝体绕流流场进行了数值模拟．计算采用非交错网格的有限体积法求解二维不可压N-S方程．计算结果表明，RNG k-ε湍流模型可以成功地模拟绕钝体的不稳定、非定常和剧烈分离流动。     

求解不可压缩流动的一种新算法

在Jameson有限体积法的基础上发展了一种模拟不可压流动的数值方法。该算法用四步Runge-Kutta时间推进法直接求解非线性的动量方程,充分利用了良好的非线性收敛特性,并从连续方程出发推导了求解压力的方程,解决了压力与速度的耦合问题。最后用该方法在同位网格上,对二维平板边界层流、二维顶盖驱动方腔层流及二维扩压器管道湍流进行了数值模拟。计算结果与理论值、标准解及试验数据进行了比较分析,从而证实了该方法对模拟不可压流动的适用性,并显示出了编程简单、易于推广、计算效率高等优点。     

砂岩脆—韧性变形转化的特征

通过砂岩的三轴试验，发现：在常温下，随着围压增加，多孔砂岩由脆性破裂转化为韧性碎屑流动；围压和孔隙在此转化过程中起了决定性的作用，砂岩脆性破裂是一个自组织过程，产生了分形的破裂面和分形分布的微破裂；砂岩的韧性碎屑流动是砂粒内的微破裂的连续滑动；砂岩的韧性碎屑流动导致了孔隙的破裂。     

三维低雷诺数后台阶流动的Lattice Boltzmann Method模拟

Lattice Boltzmann Method (LBM) 是近年来兴起的新型流体计算方法,该方法已成功应用模拟各工程中的流体流动.本文利用该方法研究在低雷诺数的不同雷诺数下三维后台阶流动的计算,并与前人的研究结果进行了比较,计算表明:与传统的计算流体力学方法相比,该方法具有极易处理复杂的边界条件和卓越的并行计算性能.已证明该方法在模拟恒温、不可压和低雷诺数条件下的流体流动时,有非常好的计算精度.     

砂岩脆－韧性变形转化的特征

通过砂岩的三轴试验，发现：在常温下，随着围压增加，多孔砂岩由脆性破裂转化为韧性碎屑流动；围压和孔隙在此转化过程中起了决定性的作用；砂岩脆性破裂是一个自组织过程，产生了分形的破裂面和分形分布的徽破裂；砂岩的韧性碎屑流动是砂粒内的微破裂的连续滑动；砂岩的韧性碎屑流动导致了孔隙的破裂。