管路系统气流两相瞬变流的矢通量分裂法

提出了一种新的求解气液两相瞬变流的方法－矢通量分裂法，将控制方程（质量方程、动量方程）中的矢通量按照矢量雅可比矩阵分裂为两个亚矢通量，采用迎风格式建立有限差分方程，该方法是稳定的并且可以自动准确地捕捉到两相瞬变流中的激波。     

管路系统气液两相瞬变流的矢通量分裂法

提出了一种新的求解气液两相瞬变流的方法——矢通量分裂法，将控制方程（质量方程、动量方程）中的矢通量按照矢量雅可比矩阵分裂为两个亚矢通量，采用迎风格式建立有限差分方程．该方法是稳定的并且可以自动准确地捕捉到两相瞬变流中的激波     

二维非定常Euler方程计算方法研究

目的研究二维非定常Euler方程计算方法。方法将修正的Steger-Warming矢通量分裂方法和NND格式相结合构成一种基于Steger-Warming矢通量分裂的差分方法。结果数值计算结果证明该算法有着良好的计算精度和计算效率。结论解决了Steger-Warming矢通量分裂方法在声速附近存在数值解振荡的问题。     

数值计算用于瞬态四波混频理论分析

采用差分方法数值求解瞬态四波混频耦合波方程，给出两种差分格式并进行数值计算，所得结果与实验基本相符。     

通风分裂算法在有限元非结构网格中的推广

将在结构网格中发展起来的通风分裂算法推广到有限体积非结构网格中,并成功地进行了定常跨声速二维欧拉方程的数值计算．此格式具有流通量矢量分裂格式的简单性及通量差分分裂格式的高阶精度．最后给出的ＮＡＣＡ００１２翼型绕流的结果与实验结果符合的较好．     

结合矢通量分裂的NND格式及其应用

将Van Leer矢通量分裂格式和NND格式相结合构成一种基于Van Leer矢通量分裂的差分方法,通过对激波管的数值分析证明该算法有着良好的计算精度和计算效率,且能自动捕捉激波.     

各向异性固体弹性介质中群速度的解析表达式

通过引入一个仅仅与介质弹性常数有关的矩阵M ,将Christoffel系数矩阵分解成M和一个与波传播方向有关的矢量 P两个部分的乘积 .由此推导出适用于计算各种各向异性固体弹性介质弹性波群速度的普适解析表达式 .群速度作为弹性波波矢的函数 .在特殊对称性方向上该解析表达式的结果与已有文献的结果一致 .利用Matlab软件数值求解Chritoffel方程 ,得到对应于任意方向上的相速度 .图 2 ,参 7.     

复杂射击规范下枪管温度场数值分析与试验研究

针对复杂射击规范下射击过程中的瞬态热冲击引发的枪管寿命问题,建立了枪管1维瞬态传热的数学模型.基于能量平衡法和交替隐式差分格式（ADI）,推导了内节点及边界点的有限差分方程,并以某5.8 mm步枪枪管为对象,在经典内弹道数值计算结果的基础上,编程求解了连续射击150发弹过程中枪管的温度场分布.最后,采用红外热像仪对枪管外壁温度进行了测试,数值仿真结果与实验结果吻合度较高,验证了所建模型及采用的算法在复杂射击规范下的有效性.      

适合渗透海床上的Boussinesq水波数学模型

为了数值研究波浪在渗透海床上的传播特性,基于近似到二阶完全非线性的高阶Boussinesq水波方程,通过引入考虑可渗海床影响的动量方程,将方程拓展到可适用渗透海床情况。对该方程常水深情况下进行了理论分析,并与解析解比较,讨论了方程的相速度及衰减率。在非交错网格下离散该方程,建立了预报—校正的有限差分数值模型。利用数值模型模拟了渗透地形上的波浪传播变形,并与相关实验结果进行了比较,实验结果与数值结果吻合较好,说明这种拓展Boussinesq水波方程至渗透海床情况是合适的。     

一类场矢量的并矢格林函数无旋和无散分解问题

将一类矢量偏微分方程相对应的并矢格林函数方程分解成无旋和无散部分,给出推导一类场矢量的并矢格林函数无旋和无散分散式的新方法。     

油水井瞬态流入动态关系曲线研究

基于质量守恒原理,建立了油水两相地层瞬态渗流的数学模型。通过定义油水两相拟压力函数,将渗流数学模型线性化,在获得数学模型解后,推导出在地层中为油水两相渗流的油水同产井的产能方程,从而建立了油水同产井瞬态流入动态关系曲线的理论,由此,可以作出油水同产井的瞬态流入动态关系曲线,为油水同产井的动态分析提供了理论基础。     

抽油机井杆管环空瞬态流场煤粉颗粒排出研究

对具有周期性运动边界的抽油机井杆管环空的瞬态固液两相流问题,提出采用动网格技术处理边界和流场内部网格节点运动的方法,保证较高的网格质量作为数值求解流场的前提;建立基于动网格技术的煤粉颗粒和水的固液两相流数学模型,采用基于非结构网格的SIMPLEC算法求解非稳态流场,得到煤粉颗粒瞬态的速度和浓度分布,进而求得杆管环空煤粉的质量流量及不同工况下煤粉颗粒的排出量,并通过排煤粉试验验证了数值计算的准确性。结果表明:抽油杆的往复运动对杆管环空流场产生的扰动作用,使得相同排量不同冲程冲次的抽油机井杆管环空的煤粉排出量不同;排量越小,煤粉粒径越大,煤粉排出量的差异就越大,最大相差达29%。     

三维飞船绕流问题的高精度迎风紧致差分格式

从迎风紧致逼近出发,给出了一个求解三维可压Navier-Stokes方程的一种高精度的数值方法,利用Steger-Warming的通量分裂技术将守恒型方程中的流通向量分裂成2部分,在此基础上据风向构造逼近于无粘项的三阶迎风紧致有限差分格式,对方程中的粘性部分采用通常的二阶差分逼近,所建立的差分格式被用来数值求解了三维飞船绕流问题。     

产水气井稳定点二项式产能方程及其应用

对于存在底水、边水、层间水的气藏,在开发过程中,易出现气水两相流动。准确评价气井产能对于合理开发气藏具有重要意义。从渗流力学的基本原理出发,结合质量守恒方程,引入水气比参数,建立气水两相拟稳定渗流的数学模型,推导出气水两相渗流的稳定点二项式产能方程。根据此方程,结合不稳定试井解释资料,以及一个稳定点的数据,能够确定气井开始产水时的无阻流量以及生产过程中不同时间段的无阻流量。作出对应的流入动态关系曲线图,只需在生产过程中测得一个稳定点资料即可,无需单独做产能试井,为准确认识产水气井的产能提供了理论基础。     

可动凝胶深部调驱流线模拟方法研究

可动凝胶深部调驱过程中流体的流动涉及复杂的物理化学过程,快速准确地模拟该过程是参数优化、方案设计的关键.在建立起可动凝胶深部调驱的油水两相多组分数学模型的基础上,推导出了考虑重力和毛管力的三维流线模型,将求解饱和度的三维问题转化为一维问题.采用算子分离技术将流线上的饱和度方程分解为典型的双曲型方程和抛物型方程,并采用分数步差分法求解饱和度方程,利用显示方法解决了隐式问题,提高了求解精度和速度,从而实现了可动凝胶深部调驱的快速准确模拟.     

可动凝胶调驱快速模拟方法研究

为了解决目前可动凝胶深部调驱数值模拟计算速度慢、数值弥散现象严重的问题，针对高含水期油田开采的特点，建立了可动凝胶调驱的油水两相八组分数学模型。采用流线方法将求解饱和度的三维问题转化为一维问题求解。流线上的饱和度方程属于典型的对流扩散方程，采用算子分离算法将该方程分解成一个双曲型方程和一个抛物型方程，并采用各自稳定性好的差分格式分别求解，提高了求解的速度和精度。设计了一个正韵律剖面模型。对影响可动凝胶调驱的注入参数和地层参数进行了敏感性分析。结果表明，随着聚合物浓度、聚胶比、注入体积、成胶时间的增加，采收率提高幅度均有所增加；纵向窜流程度增加，可动凝胶词驱效果变好。     

液氮垂直流动沸腾的双流体模型分析

采用双流体模型计算了液氮在垂直管内的上升流动沸腾过程,考察了壁面热通量和液体流量对流动及传热传质特征的影响.结果表明:垂直上升流动沸腾中重力压降占主导地位;根据截面液体温差的变化可判断沸腾模式的转变;壁面热通量与液相流量的相对大小决定了沸腾过程中的传热传质特征.     

气-液两相流界面迁移现象的数值模拟研究

借助于LevelSet函数,建立了气-液两相的统一控制方程组,并在交错网格中进行离散.用两种格式,即Superbee TVD格式和5阶WENO格式求解LevelSet函数的输运方程,用SIMPLER算法的思想对主流场控制方程的求解方法进行改进.数值实验结果表明,在求解LevelSet的控制方程时,5阶WENO方法比Superbee TVD格式的结果更准确;用改进的数值算法可成功实现对密度比大于1000/1的气-液两相流界面迁移问题的数值模拟.对几种典型大密度比气-液两相流问题的计算结果与实际问题的物理规律完全一致,验证了该方法的有效性和可靠性.     

直井气侵后气液两相参数分布数值模拟

钻井过程中发生气侵后,由于压力、温度、偏差因子等参数变化影响,气体沿井眼滑脱上升的同时会伴有一定膨胀,而加剧井底压差,导致气侵量不断加大直至井喷.气侵直至井喷过程中井筒内为瞬态气液两相流动,研究气侵后井底气液两相参数发展变化规律及井底压力变化特征对认识井喷发生发展规律有重要的指导意义.以甲烷-钻井液为循环介质,在接近井底实际工况情况下对气侵发生直至井喷过程进行了数值模拟,模拟计算了不同气侵量下气侵不同阶段的流型变化规律、井底压力变化规律、速度分布,并取得了初步认识.