油藏水平井数值模拟研究中一种新的求解方法

提出了底水油藏水平井油水两相渗流问题的物理模型和数学模型，并建立了差分方程，然后利用ＩＭＰＥＳ求解方法的思路，并对其不足之处进行了改进，从而形成了一种新的求解方法．这种方法是采用上游网格饱和度插值计算相邻网格间的水量交换，用本阶段时间中点的饱和度对压力方程进行再迭代．保持了ＩＭＰＥＳ方法简单方便、占计算机内存少、计算速度快的优点，而精度明显提高，结果更趋合理．实际计算表明，新方法效果良好，可广泛用于底水、边水或人工注水油藏油水两相渗流问题的数值模型研究，对水平井和垂直井均适用．     

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

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

以压力为变量解可压流动的非交错网格SIMPLE方法研究

提出一处以压力为基本求解变量的非交错网格ＳＩＭＰＬＥ方法，用来计算可压缩流动，本方法的特点在于将适用球可压缩流动计算的ＳＩＭＰＬＥ方法加以扩展，使之也适用于可压缩流动的计算，因而有较大的适用范围，对亚音速、跨音速及超音速等３个较典型的算例进行了数值计算，计算结果与实验及文献的结果符合较好。     

混流式风机流道内三维湍流流场的研究

应用高雷诺数ｋ－ε湍流模型进行了混流式风机叶轮机三维流动的计算研究，使用ＳＩＭＰＬＥ型算法求解了全三维雷诺平均质量加权Ｎ－Ｓ方程对流式采用了高阶精度的上流国权近ＴＶＤ类型的ＭＵＳＣＬ差分格式，采用非正交贴体坐标和同位网格系，速度的协变物理分量被选作为独立求解变量，为了消除非物理的压力和速度振荡，对界面上的逆变速度分量采用Ｒｈｉｅ，Ｃｈｏｗ动量插值方法，离散方程以能够使用三对角矩阵算法的形式书写，并     

固体火箭发动机燃烧室层流流场的数值模拟

该文发展了一种在任意曲线坐标系上求解层流Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值方法，该算法以ＳＩＭＰＬＥ为基础，采用了非交错网格，因而对原始算法中的压力修正方程进行了改进。用准定常方法数值模拟了固体火箭发动机燃烧室内的二维轴对称流场，计算结果能够反映流场内的旋涡与各参数的分布。计算表明，压力与速度等参数的分布明显受旋涡存在的影响，比传统的一维流场复杂得多。     

管式反应器中交混流场的数值模拟

对管式反应器中的交叉混合流场进行了三维数值模拟,采用分块结构化网格技术以及网络拼接技术实现了适用于复杂形体的分块结构化网格的计算程序,采用修正的ｋ－ε双方程湍流模型描述反应器的湍流流动,用类似ＳＩＭＰＬＥ法对控制方程进行求解,经冷态模拟试验检验,结果令人满意     

用Richardson外推法分析流动传热层流问题的不确定度

用ＳＩＭＰＬＥＣ算法计算了流动传热问题中２个有基准解的层流问题。对于２阶精度中心差分格式的２套网格数值解,用Ｒｉｃｈａｒｄｓｏｎ外推法可以得到４阶精度的解、其工作量比直接求解４阶差散方程的工作量减少许多;对于计算中的数值不确定度,用该方法分析了离散误差,并用逐次加密网格的方法探讨了外推法的合理性。     

高压离心压缩机轴向推力计算

详细分析了离心压缩机轴向力的组成,建立了离心叶轮外侧间隙内泄漏气体流动的计算模型,采用了低Ｒｅ数ｋε模型,用ＳＩＭＰＬＥＣ方法求解控制方程．对１台离心压缩机的轴向力进行了计算,结果令人满意．     

一种改进的SIMPLE算法

本文对传统的ＳＩＭＰＬＥ算法进行了改进，提出了一种新的协调型计算方法，并通过两个具体算例的计算与ＳＩＭＰＬＥ算法和ＳＩＭＰＬＥＣ算法进行了比较。结果表明，本算法有更好的收敛性。     

摩托车化油器腔体内部两相流动的数值模拟

针对摩托车化油器的工作原理,应用ｋ－ε两方程湍流模来描述气相场,并在欧拉坐标系中用非交错网格的ＳＩＭＰＬＥ算法求解;液相场在拉格朗日坐标系下用四阶Ｒｕｎｇｅ－Ｋｕｔｔａ算法求解,气－液两相的耦合采用ＰＳＩＣ法;最后对节气门全开、半开及经济工况下的流场进行数值模拟,得出了流场变化规律。     

计算流体力学与传热中的多重网格块修正算法

在流体力学的数值模拟计算中经常会遇到大型代数方程组的求解,而求解代数方程组的代价随着方程个数的增加而迅速增大。 主要将块修正多重网格算法引入ＳＩＭＰＬＥ通用程序中,并以流体力学中五个典型的实际问题为例,分别采用四种不同的迭代方法进行计算比较,展示了块修正多重网格法加速收敛速度的特性。     

轴流压气机转子内三维湍流流场计算

应用高雷诺数ｋ－ｅ湍流模型进行了轴流压气机叶轮内三维流动的计算．使用ＳＩＭＰＬＥ型算法求解了全三维雷诺平均质量加权Ｎ－Ｓ方程．对流项采用了高阶精度的上游加权近似和ＴＶＤ类型的ＭＵＳＣＬ差分格式．采用非正交适体坐标和同位网格系，速度的协变物理分量被选作为独立求解变量．为了消除非物理的压力和速度振荡，对界面上的逆变速度分量采用Ｒｈｉｅ和Ｃｈｏｗ动量插值方法．求解中考虑了周期边界条件和运动壁面边界条件．对两个轴流叶轮进行了计算，计算结果表明全三维湍流计算能有效地预测轴流叶轮对道内及无叶间隙处的分离，以便改进气动设计、几何设计以提高轴流压气机的性能．     

多重网格法及在HTR-10堆芯动态仿真

研究了ＨＴＲ－１０堆芯的动态仿真,为反应堆控制系统的分析和设计服务。由于堆芯模型方程是典型的非定常可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,要求快速准确的求解。求解算法采用了结合多重网格的ＳＩＭＰＬＥＣ算法,并且讨论了多重网格法的基本思想、算法步骤和计算流程。通过对求解域中的不规则区域的处理、合理的进行粗细网格的划分和校正、设置出口缓冲区等措施,克服了在ＨＴＲ－１０堆芯动态仿真中运用多重网格法的困难,并且所得的结果与单网格算法进行了比较。仿真结果表明：达到同样精度,采用多重网格方法,求解过程比单网格算法的迭代次数减少了约５／６,计算时间减少了约一半,为快速进行堆芯动态仿真提供了条件。     

NS方程的非结构化网格方法及其差分格式

采用同位非结构化网格上的有限体积法,对ＮａｖｉｅｒＳｔｏｋｅｓ 方程的ＳＩＭＰＬＥ 算法及差分格式进行了研究．提出了适用于非结构化网格的不用求解单元顶点变量值的二阶混合差分格式．该格式的优点在于：１） 减少了计算工作量;２） 避免了普通混合差分格式因使用简单的一阶迎风差分格式所引起的网格界面方向相关性的问题．最后,采用三角形网格,利用提出的方法及差分格式,对方腔内的驱动层流及绕翼型湍流进行了数值计算,计算结果与基准解或实验值的符合程度优于普通混合差分格式     

油水分离旋流器内部流场的数值模拟

本用数值方法模拟了油水分离旋流器内部的流场分布。油水混合物中的水被作为流体相，油被作为颗粒相。流体相的控制方程为连续性方程、Ｎ－Ｓ方程，其中结合了修正的Ｐｒａｎｄｔｌ混合长湍流粘性模型，颗粒相的控制方程是浓度扩散方程，流体相与颗粒相之间的相对运动由代数滑移模型确定，全部控制方程用Ｐａｔａｎｋａｒ－Ｓｐａｌｄｉｎｇ的ＳＩＭＰＬＥ算法求解。计算程序取名ＣＹＣＬＯＮ，用通用形式写出，用户只要输入旋流器     

三峡转轮内部流动的全三维欧拉解

水力机械内部流动非常复杂,是强三维流动,只有用三维数值方法计算才能获得比较满意的结果。本文对水轮机转轮内高度复杂的流动进行了三维无粘数值计算。计算建立在三维Ｅｕｌｅｒ方程基础上,速度场和压力场的求解方法为速度、压力修正算法即ＳＩＭＰＬＥＣ算法,并且对固壁条件的给出提出了一种新的方法。最后对一真实的模型转轮进行了计算     

粗网节块内功率的重构

粗网节块程序如ＮＧＦＭ、ＳＩＭＵＬＡＴＥ等只能给出平均功率。而反应堆设计中更关心的是堆芯内燃料元件的功率分布。本文讨论了获得堆内元件功率分布的重构方法。双二次多项式用于快群中子通量展开，与双曲函数项结合的双二次多项式用于热群中子通量密度展开。研究了２１项和１３项展开法，并与ＣＩＴＡＴＩＯＮ细网格计算进行比较。结果表明，本文中研制的２１项方法程序ＲＥＣＯＮ比ＳＩＭＵＬＡＴＥ－３ＮＥＭＯ程序中所用的１３项方法精确得多，在外围区域更明显。     

修正QUICK格式及其在弹状流中的应用

提出了一种适用于非均分网格的对流项差分格式优化ＱＵＩＣＫ格式,对数值传热学经典问题进行了验证性计算,得到较为理想的结果．应用这种优化ＱＵＩＣＫ格式结合ＳＩＭＰＬＥＣ算法,计算了垂直上升管中弹状流泰勒气泡周围的流场．在此基础上求解修正动量积分方程,从而修正气泡的形状,以满足对气泡内部压力恒定的要求,并选出合理的气泡上升速度．     

单井三维三相黑油模型的建立及自动历史拟合的实现

建立了一个柱坐标系下的三维三相单井黑油模型，通过隐式压力隐式饱和度（ＩＭＰＩＭＳ）方法对模型进行求解。获得了稳定的压力和饱和度计算结果。在此基础上，建立了一个带约束条件的非线性目标函数，运用最优化方法使目标函数达到最小，实现自动历史拟合，从而反求得各地层参数值。另外还对大港油田三口测试井的测试资料进行了自动历史拟合，求得了各自的地层参数值，同时将本文求得的地层参数值与美国ＳＳＩ软件解释结果及电测解     

复杂区域流动换热问题的一种新的网格处理方法

在阶梯形网格法的基础上发展了一种二维复杂区域流动与换热问题的新的数值方法。该法采和了同位网格下的ＳＩＭＰＬＥＮ算法,并借鉴了非结构化网格技术来进行复杂计算区域的网格划分及组织计算。用这种方法计算了方腔内的自然对流及圆内开缝正方形和热环体的自然对流,取得了很好的效果。