单裂缝中携砂液流动规律研究

裂缝中携砂液流动是一种固液两相流,携砂液的运移与支撑剂的铺置是水力压裂裂缝保持导流能力的关键.本文基于FLUENT流体计算软件,采用双流体模型,将颗粒看作拟流体,携砂液按照牛顿流体处理,分析了支撑剂体积分数α_s、阿基米德数Ar、颗粒雷诺数Re以及裂缝入口边界对流动规律的影响.研究结果表明:携砂液在裂缝中的流动过程中,发展成为支撑剂体积分数不同的四个区域,包括砂堤区、颗粒悬浮区、颗粒滚流区和无砂区;支撑剂的沉降程度随着支撑剂体积分数和阿基米德数的增加而增加,而随着雷诺数增加而降低;入口为网眼型时,进入裂缝后过流面积的增加导致流速突降,使得支撑剂更容易在入口处产生堆积,在同一入口流速下,较均匀入口的工况铺砂高度大.     

气液两相流压力波传播速度研究

将双流体模型用于绝热无相的管道气液两相流,依据小扰动线化分析原理,导出了压力波波数Ｋ方程通过对不同空隙率下肉体上压力波小随角频率变化的计算,研究了虚拟质量力和狭义相间阻力对压力波波速及其人色散性的影响。对泡状流和弹状流压力波波速的计算结果与前人的测量结果作了比较,两者符合良好。     

可侵蚀地表上方含尘大气运动的数值模型

在稀相气固两相流的双流体模型框架下，导出含尘大气运动的基本方程，其中的源项考虑相间双向耦合作用；通过对基本方程进行无量纲化，求得控制气相和弥散相动力学行为的相似准则，包括弗洛德数、罗斯贝数、颗粒雷诺数、颗粒惯性参数、颗粒剪切横移参数、颗粒旋转横移参数和沙尘质量载荷率等. 作为一个模型问题，研究可侵蚀地表上方充分发展湍流大气边界层流动对土壤颗粒的气动卷扬作用，其中计及气动阻力、萨弗曼升力和重力等因素的影响. 为了克服由于沙尘跃移运动轨迹交叉导致流动参数多值性所造成的困难，引进拉格朗日坐标下弥散相连续方程. 对两种不同风速和4种不同粒径条件下的沙尘运动进行数值模拟，给出相应情况中沙尘运动特性和密度分布剖面并讨论风速和粒径等参数的影响，还细致分析了含尘大气边界层中沙尘与气流之间的能量传递过程. 所得的方程、准则和方法可以为土壤风蚀和沙尘暴等自然灾害的预报提供理论基础.     

视密度加权的时平均统一二阶矩两相湍流模型

在气粒两相湍流的双流体模型中,颗粒相的视（表观）密度是有脉动的,在时平均的统一二阶矩（USM）模型中出现了和颗粒数密度或视密度脉动有关的项和方程,使模型方程比较复杂。实际上,用LDV或PDPA测量的流体（用小颗粒代表）和颗粒速度都是颗粒数加权平均的结果。因此,在视密度加权平均基础上推导两相湍流模型更为合理。通过推导和封闭了视密度加权平均的统一二阶矩模型（MUSM）方程组,改进了两相速度脉动关联的封闭,并引入了颗粒遇到的气体脉动速度及其输运方程。MUSM模型可以减少所用方程数,节省计算量。视密度加权平均的统一二阶矩两相湍流模型是一种对原有时间平均的统一二阶矩模型和改进和发展。     

平流式沉淀池液固两相流力学特性的立面二维数值模拟

采用双流体模型,选取RNG k-?紊流模型封闭双流体两相流时均方程,对平流式二次沉淀池液固两相流力学特性进行立面二维数值模拟。采用有限体积法求解双流体微分方程;紊动能、紊动能耗散率方程均采用Quick离散格式;耦合求解速度与压力时使用了压力隐式算子分裂PISO(Pressure-Implicit with Splitting of Operators)算法。通过模拟表明:水流经过挡板后,底部是急剧扩散的主流,流速大,上部是回流;流速沿流程逐渐减小,到中部近似均匀流;出口附近水流收缩流出,流速变大。进水口附近紊动动能较大,而后紊动动能沿流程逐渐减小,在出水口附近,由于流速的增大紊动动能也随之增大。将模拟值与实测值进行了对比分析发现,流速的最大误差为10%,出口悬浮物相对浓度的最大误差为15%。这表明该模型能够较好地模拟沉淀池中的水力特性分布规律,对沉淀池的设计具有一定的参考价值。     

明渠气二相流的双流体模型

从解决平均二相流基本方程的封闭问题出发,分析了水利工程中气水二相流的特点,介绍了据此提出的可用于水利工程稀疏气泡流计算的双流体模型,并对明渠气水二相流进行了算例验证。     

论气体混合物和两相流中的压强和热流

本文首先给出单流体模型和双流体模型中压强和热流等物理量的两种不同定义。然后通过引入流场中某面元dS的动置流密度M和能量流密度δ来讨论压强和热流等量的物理意义。指出,在描述混合物运动的守恒方程中以采用单流体模型定义的各量为宜,在描述组元运动的守恒方程中以采用双流体模型定义的各量为宜。本文还详细讨论了稠密混合气体中的碰撞压强和碰撞热流,说明它们的来源,它们在动量方程和能量方程中的表达形式,以及它们同通常的分压和分热流的本质差别。 气固两相流可视为一种稠密的混合气体流动。两相流中关于动置方程中是否应有“惯性耦合项”的长期争论,可以用分压的两种不同定义予以澄清。     

垂直上升圆管内等温泡状流的数值模拟

介绍了描述等温状态下垂直气泡流流动的总体平衡方法(population balanceapproach). 将平均气泡数密度控制方程引入到双流体模型中实现总体平衡. 介绍了描述气泡聚并与破裂机制的Yao和Morel理论模型以及Hibiki等人的实验方法. 利用商业计算流体力学软件ANSYSCFX10, 对圆管特定位置上的5个基本变量:气含率、气泡平均直径、相间表面积浓度以及气体和液体速度沿半径分布进行了模拟计算. 预测结果与实验数据的比较表明,预测结果和实测数据之间有很好的一致性.      

基于相变松弛法则的水下爆炸空化模型研究

水下爆炸过程中存在着大量的空化现象,空化的产生、演化及其溃灭过程对于水下冲击波传播、爆炸气泡运动以及水下结构物冲击损伤都会产生重要影响。本文基于多相可压缩流体理论模型,考虑空化发生过程中汽-液两相流体亚平衡状态下两相之间发生的热力学-化学平衡机制,分析汽-液两相介质之间的质量和热量交换,从而实现对相变过程的自动捕捉。该系统的控制方程采用分步法处理,首先利用二阶MUSCL-Hancock格式和HLLC黎曼求解器来求解齐次双曲型方程,再采用牛顿迭代法求解相变方程。数值测试结果表明,本文的计算模型对于空化相变过程具有较好的捕捉能力。最后将该模型应用到水下近水面爆炸空化的数值模拟当中,研究发现空泡的溃灭压力峰值约为冲击波压力峰值的15%,有效作用时间是冲击波载荷有效作用时间的2倍以上。本文的空化相变模型能够为水下爆炸空化现象的机理研究提供重要支撑。     

一种修正的低温流体空化流动计算模型

为了更准确地预测低温流体的空化流动特性, 基于Kubota空化模型, 对蒸发和凝结源项进行修正, 建立了一种考虑热力学效应的空化模型. 分别采用原始和修正的Kubota空化模型, 计算了绕对称回转体液氮的空化流动, 通过与实验结果的比较对修正的空化模型进行了评价. 结果表明, 与原Kubota空化模型比较, 修正的空化模型由于考虑了热力学效应, 计算获得的蒸发量减小, 凝结量增大, 空穴长度减小, 空穴界面形态呈模糊状态.计算结果与实验结果更加一致, 说明修正的空化模型能准确的描述低温流体空化过程的质量传输过程, 能够更准确模拟低温流体中的空化流动特性.      

Analytical considerations of flow boiling heat transfer in metal-foam filled tubes

Flow boiling in metal-foam filled tube was analytically investigated based on a modified microstructure model, an original boiling heat transfer model and fin analysis for metal foams. Microstructure model of metal foams was established, by which fiber diameter and surface area density were precisely predicted. The heat transfer model for flow boiling in metal foams was based on annular pattern, in which two phase fluid was composed by vapor region in the center of the tube and liquid region near the wall. However, it was assumed that nucleate boiling performed only in the liquid region. Fin analysis and heat transfer network for metal foams were integrated to obtain the convective heat transfer coefficient at interface. The analytical solution was verified by its good agreement with experimental data. The parametric study on heat transfer coefficient and boiling mechanism was also carried out.     

Turbulent mixed convection heat transfer to liquid sodium

The influences of buoyancy on turbulent heat transfer to a liquid metal flowing in a vertical pipe are considered. A theoretical model is presented which provides a criterion for the conditions under which such influences become significant and which predicts the impairment of heat transfer for upward flow and enhancement for downward flow. The variation with Peclet number of the maximum impairment of heat transfer and conditions under which it occurs are established. A generalization of the model leading to an equation for the entire mixed convection region is proposed. From this an equation for turbulent free convection to liquid metals is obtained.     

自由场中液氮单空泡动力学特性的实验研究

本文专门设计搭建了低温介质空泡演化实验测试平台, 对液氮单空泡非定常演化过程和动力学特性开展了实验研究. 实验中利用电火花瞬态放电激发液氮汽化形成单空泡, 通过高速摄影系统对单空泡的瞬态特征进行了精细化捕捉. 为了进一步揭示低温介质独特的物理性质以及强热力学效应对单空泡演化过程的影响机制, 对比分析了在相同环境压力下, 77.41 K液氮和298.36 K水单空泡的演化过程和动力学特性. 基于实验得到空泡半径与界面速度等定量数据, 阐明了液氮单空泡球形与非球形演化阶段的非定常特性. 研究结果表明: (1) 在相同输入电压下, 液氮单空泡的整体尺寸比常温水更小, 当输入电压为400 V时, 液氮空泡的最大半径约为常温水空泡的0.69倍; 同时, 液氮单空泡经历了膨胀阶段?收缩阶段?振荡阶段以及上升阶段的演化过程. (2)液氮空泡的收缩过程主要由相界面的热传导主导, 没有明显的塌陷现象, 收缩阶段液氮空泡的最小收缩半径约为常温水的5.5倍. (3)在液氮空泡振荡初期, 空泡相界面传热增强, Rayleigh-Taylor不稳定与热力学效应共同引起了空泡界面的表面粗化效应; 在整个振荡阶段, 空泡界面附近存在破碎的小泡. 当输入电压较高时, 空泡底部的小泡数量显著增多. (4)由于液氮空泡浮力系数较大, 液氮空泡在演化后期空泡整体向上迁移显著, 液氮空泡底部收缩更快产生凹陷, 促使空泡变为环状.      

Numerical simulation of the modulated flow pattern for vertical upflows by the phase separation concept

The multiphase heat transfer could be enhanced by creating thin liquid film on the wall. The phase separation concept is called due to the separated flow paths of liquid and gas over the tube cross section to yield thin liquid film. Our proposed heat transfer tube consists of an annular region close to the wall and a core region, interfaced by a suspending mesh cylinder in the tube. The heat transfer tube is a multiscale system with micron scale of mesh pores, miniature scale of annular region and macroscale of tube diameter and length. Great effort has been made to link from micron scale to macroscale. The Volume of Fluid (VOF) method simulates air/water two-phase flow for vertical upflow. The three-dimensional system was successfully converted to a two-dimensional one by using three equivalent criteria for mesh pores. The non-uniform base grid generation and dynamic grid adaption method capture the bubble interface. The numerical results successfully reproduce our experimental results. The numerical findings identify the following mechanisms for the enhanced heat transfer: (a) counter-current flow exists with upward flow in the annular region and downward flow in the core region; (b) void fractions are exact zero in the core region and higher in the annular region; (c) the liquid film thicknesses are decreased to 1/6–1/3 of those in the bare tube section; (d) the gas–liquid mixture travels much faster in the annular region than in the bare tube; (e) three-levels of liquid circulation exists: meter-scale bulk liquid circulation, moderate-scale liquid circulation around a single-elongated-ring-slug-bubble, and microliquid circulation following the ring-slug-bubble tails. These liquid circulations promote the fluid mixing over the whole tube length and within the radial direction. The modulated parameters of void fractions, velocities and liquid film thicknesses in the annular region and three-levels of liquid circulation are greatly beneficial for the multiphase heat transfer enhancement.     

Numerical modeling of cryogenic chilldown process in terrestrial gravity and microgravity

In this paper, a numerical model is developed to predict the cryogenic chilldown process of a vertical tube for both terrestrial and microgravity conditions. The flow field is covered by four distinct regions, which are single-phase vapor region, dispersed flow region, inverted annular flow region, and single-phase liquid region. Heat transfer mechanisms are dictated by the flow patters. A two-fluid model is employed to analyze the dispersed flow region and the inverted annular film boiling region. Gravity effect on the chilldown process is also investigated. Model results indicate that film boiling heat transfer decreases with decreasing gravity level for the bottom flooding condition. The model results show a good agreement with the experimental data.     

开区注氮采空区自燃温度场的数值模拟研究

为了研究开区注氮采空区遗煤自燃和被抑制的力学过程，用非均质多孔介质中的渗流连续性方程、气体弥散方程和综合传热方程的联立，建立了采空区注氮非定常数值模型．结合实例，用迎风格式有限元方法求解．给出了采空区的漏风流态、氮气流态，描绘了采空区自燃过程中氧、CO浓度和温度分布的变化过程．计算考虑了瓦斯涌出和工作面推进的影响．得到注氮条件下采空区高温区形状变窄；随着注氮量的提高，使自然发火期逐渐变长，直至不自燃．理论计算与实际情况吻合．     

Experimental study of jet structure and pressurisation upon liquid nitrogen injection into water

The rapid phase change and heat transfer obtained by direct contact heat exchange between a cryogen and water can generate high rates of pressurisation, which is of interest to a number of applications. A visualization study of liquid nitrogen injection into water is conducted in this work, with synchronized pressure and temperature measurement, to obtain insight into this complex phenomenon. High speed imaging reveals a four-stage evolution of liquid nitrogen jet structure upon injection into water, with a thick vapour blanket forming around a liquid nitrogen core and break-up brought on predominantly through impact with the vessel wall. Maximum pressurisation rate occurs in the third stage of injection due to a combination of heat and mass transfer. Pressurisation rates in excess of 350 bar/s are recorded and found to vary proportionally with injection pressure. The scenario of gaseous nitrogen injection is also investigated, and compared with liquid nitrogen injection. A clear advantage of liquid nitrogen injection is elucidated from the point of heat transfer and pressurisation, and implications for use in a cryogenic engine are discussed.     

Application of a lattice-gas numerical algorithm to modelling water transport in fractured porous media

A 2D lattice-gas numerical algorithm was used to simulate liquid infiltration into unsaturated porous media with a parallel sided central crack. The chosen model, the interacting liquid-gas model of Appert and Zaleski, has dynamic properties leading to a phase transition and liquid and gas phases can be simulated by the model. These two phases are used to simulate biphasic flow in porous media. Fourteen numerical experiments of liquid infiltration were carried out which differed in the morphology of the microporous matrix, in the aperture of the central crack and in the amount of liquid supplied. For the same microporous matrix, the infiltration dynamics in the dual media depended upon the ratio between the amount of liquid supplied and the crack aperture. Variations in water storage over time and liquid flow regimes within the cracks are discussed.     

Pure steam condensation model with laminar film in a vertical tube

A new physical model for calculating the liquid film thickness and condensation heat transfer coefficient in a vertical condenser tube is proposed by considering the effects of gravity, liquid viscosity, and vapor flow in the core region of the flow. To estimate the velocity profile in the liquid film, the liquid film was assumed to be in Couette flow forced by the interfacial velocity at the liquid–vapor interface. For simplifying the calculation procedures, the interfacial velocity was estimated by introducing an empirical power-law velocity profile. The resulting film thickness and heat transfer coefficient from the model were compared with the experimental data and the results obtained from the other condensation models. The results demonstrated that the proposed model described the liquid film thinning effect by the vapor shear flow and predicted the condensation heat transfer coefficient from experiments reasonably well.