Dynamic behaviour of direct spring loaded pressure relief valves in gas service: II reduced order modelling

A previous study of gas-service direct-spring pressure relief valves connected to a tank via a straight pipe is continued by deriving a reduced-order model for predicting oscillatory instabilities such as valve flutter and chatter. The reduction process uses collocation to take into account a finite number N of acoustic pressure waves within the pipe, resulting in a set of 2N+3 ordinary differential equations. Following a novel non-dimensionalization, it is shown analytically that the model can exhibit, at experimentally realistic parameter values, instabilities associated with coupling between the valve and acoustic waves in the pipe. The thresholds for each instability are such that for a given flow rate, the first mode to go unstable as the inlet pipe length increases is the quarter-wave mode, then a three-quarter wave, a 5/4-wave etc. Thus the primary mode of instability should always be due to the quarter wave. In the limit of low flow rates, a simple approximate expression is found for the quarter-wave instability threshold in the form of inlet pipe length against mass flow rate. This threshold curve is found to agree well with simulation of the full model. For higher flow rates there is a need to include fluid convection, inlet pressure loss and pipe friction in order to get good agreement. The reduced model enables the dependence of the stability curve on key dimensionless physical parameters to be readily computed.     

Steady-state critical two-phase flashing flow with possible multiple choking phenomenon: Part 1: Physical modelling and numerical procedure

In this paper, the modelling of pure fluid steady-state adiabatic flashing flow through a pipe involving an abrupt enlargement is presented. The thermal non-equilibrium two-phase flow model DEM has been adjusted for subcooled inlet conditions close to saturation and/or for inlet two-phase mixture state. The multichoking flow phenomenon occurs when two basic criteria are simultaneously fulfilled. A general procedure taking into account this possible occurrence is developed on the basis of two iterative algorithms. The first algorithm is applied to the mass flow rate upstream from the enlargement and the second algorithm is based on the length of the pipe downstream from the enlargement. The proposed methodology involves the improved physical one-dimensional two-phase flow model DEM and the global non-equilibrium model through the abrupt enlargement.     

Benchmarking of two-phase flow through safety relief valves and pipes

This paper compares calculated results for two-phase flows through safety relief valves and pipes using the TPHEM, CCFLOW and RRERSP computer programs. These studies were conducted to locate errors in the programs as well as to further our understanding of how each program worked. For most low-to-moderate viscosity flow examples, and for the frozen flow examples, the program results agreed despite differences in the calculation methods. Current thinking is that the TPHEM computer program gets better results for high viscosity flows through safety relief valves (nozzles) and for pipe down flows. This is because this program successively iterates nozzle flows to achieve results consistent with the choke pressure, temperature and quality and successively iterates pipe flows to achieve results consistent with the pipe length. Safety relief valve flows from the RRERSP computer program are reduced by the inlet pipe non-recoverable pressure loss. This effect is significant for several cases.     

Numerical investigation of flow erosion and flow induced displacement of gas well relief line

High-pressure particle-laden gas flow should be discharged through relief line of gas well timely to ensure safe test and exploitation. Erosion and vibration usually take place on the bend in relief lines, bringing a potential safety hazard to field operation. The majority of this paper investigates the factors affecting the erosion of bend and displacement of relief line in the downstream of bend using the computational fluid dynamics (CFD) methodology. A three-dimensional elbow pipe is selected as computational domain in this investigation. The kinematics and trajectory of discrete solid particles and liquid droplets are described by discrete phase model (DPM) while the hydrodynamic characteristics of continuous phase are obtained based on Reynolds-Averaged-Navier-Stokes (RANS) equations. An empirical erosion model is employed to predict the erosion rate of bend, and a fluid–structure interaction (FSI) model is adopted to calculate the displacement of relief line. The effects of types of multiphase flow (such as gas–solid two-phase flow and gas–liquid–solid three-phase flow), inlet flow rate and pipe diameter on erosion and displacement are discussed. The results show that large displacement and severe erosion present with large inlet flow rate in minor diameter pipe. The increase in liquid droplet content has less effect on flow erosion than that by the same increase in sand particle content.     

油水两相流弯管处安全分析

针对油水两相流经过弯管时的流向改变会导致流体速度和压力发生突变,造成发生静电事故和腐蚀事故风险上升 的不利影响,提出了RSM模型和Mixture模型相结合的安全分析方法。该方法对不同入口速度和含水率的油水两相流进行 数值模拟,并用Origin软件拟合了进口最大允许流速与管径及含水率的经验关系。结果表明,在含水率和入口速度一定 时,随着管径的增加,弯管处的最大速度呈现逐渐减小的趋势:当管径和入口速度一定时,随着含水率的增加,弯管处 的最大速度也逐渐减小。最大压力出现在弯管外拱壁处,最小压力出现在弯管内侧拱壁处。在实际生产中,增加弯管下 游直管段内侧壁的壁厚,可有效防止空化腐蚀所造成的危害;通过含水率来确定安全流速,可有效降低静电事故的风险 。     

袋式除尘器喷吹管的气流均匀性研究

针对袋式除尘器喷吹气流不均匀性的问题,利用CFD数值模拟建立含有不同结构参数的喷吹管气流分布模型,通过正交试验方法,选取了喷吹管管径D,喷嘴孔径d,喷嘴个数N,喷嘴间距L及喷嘴长度H为主要因素,选取L25（56）正交表得到25组方案,用方差分析法找到了影响喷吹气流不均匀性的主要因素,并采用流量修正和重复迭代的方法对喷嘴孔径进行优化设计,获得较好的气流量分布,为喷吹管设计提供了参考;对修正后的喷吹管,得到不同入口流量下的喷吹管气流不均匀的变化情况。     

The dynamic response of pressure relief valves in vapor or gas service,part I: Mathematical model

In order to prevent unstable operation, or “chatter”, of a pressure relief the API guidelines recommend limiting the irreversible pressure loss in the inlet line to a pressure relief valve to no more than 3% of the valve set pressure. This criterion is based on steady-state operating conditions and a typical blow-down pressure for the valve of about 7% of the set pressure. However, the stability of the valve is also influenced by other factors such as the dynamic response of the valve disk to the unsteady pressures and forces exerted by the fluid on the disk. A model for the opening lift dynamic response of a pressure relief valve in gas/vapor service is presented here which accounts for all of these effects through a set of five coupled nonlinear algebraic/differential equations. These equations are solved by a numerical method that can be implemented on a spreadsheet to predict the position of the valve disk as a function of time for given valve characteristics, operating conditions, and installation parameters. The model incorporates the influence of the various parameters on the stable/unstable nature of the disk response. An example is presented for a typical valve that illustrates the various modes of stable and unstable dynamic response that can be predicted by the model under various conditions. Two additional papers will be forthcoming: Part II – Experimental Investigation and Part III – Analysis of Data and Comparison with Model Predictions.     

深水测试地面流程节流油嘴段温压场研究

为对深水高压气井测试地面流程中的节流油嘴段的温压场及水合物生成情况进行研究,采用数值模拟的方式对整个节流管路进行气体流动特性分析.并利用P-T图回归公式法,得出不同温压场条件下的水合物生成情况,进行水合物生成范围对比.结果表明:针阀出口后端气体速度随着管路入口压力和针阀直径的增大而增大,随管路出口压力的增大而减小,而气...     

球接管容器气体爆炸尺寸效应的试验研究与数值分析

为了解尺寸对球形容器连接管道甲烷-空气混合物爆炸的影响规律,利用Fluent软件,采用κ-ε湍流模型、涡耗散模型(简称EDC模型)、壁面热耗散、热辐射模型及SIMPLE算法,建立了球形容器连接管道内甲烷-空气混合物爆炸的数值模型,对容器与管道内甲烷-空气预混气体爆炸的尺寸效应进行了数值模拟。结果表明:随管道内径增大,球形容器内最大爆炸压力逐渐增大,管道末端最大爆炸压力变化无明显规律;而随管道长度增加,球形容器内最大爆炸压力逐渐减小;改变管道内径,较大体积球形容器内最大爆炸压力均大于较小体积球形容器内最大爆炸压力,最大爆炸压力上升速率的规律则相反,容器体积对管道末端最大爆炸压力的影响无明显规律。     

人工煤气管道萘颗粒沉积规律研究

针对萘在人工煤气管道中沉积会造成管道堵塞,影响管道的安全运行的这一问题,以昆明人工煤气管道为例,运用计算流体动力学软件Fluent,选用离散相模型和雷诺应力模型,对水平直管、水平弯管和三通管进行萘颗粒沉积的数值模拟,对于不同的管径、弯曲比、管径比,分别分析萘颗粒直径、入口速度、温度及压力对萘颗粒沉积的影响。研究结果表明:水平直管、水平弯管、三通管中的萘颗粒沉积率与颗粒粒径成正相关关系,而与气流入口速度、压力成负相关关系;萘颗粒在人工煤气管道中的沉积率主要受颗粒直径、气流入口速度的影响;萘颗粒的沉积率随着水平直管的管径增大而增大,随着水平弯管的弯曲比增大而增大,随着三通管的管径比增大而先增大后减小;可通过适当增大管内煤气输送速度、压力,降低温度来降低萘颗粒在人工煤气管道中的沉积速度,进而减少萘颗粒沉积的发生。     

中压天然气管道泄漏扩散模拟研究

建立了埋地中压天然气管道发生泄漏时时的数学模型,将土壤视为各向同性的多孔介质,采用FLUENT对天然气在土壤中的扩散规律及浓度分布进行模拟,分析不同时刻地表的危险区域范围,并对比了不同管道压力、泄漏孔径大小、泄漏位置等工况下危险半径随时间的变化。结果表明:管道压力越大,泄漏的体积流量越大,同一时间危险范围越大;相同的泄漏压力下,泄漏孔径对危险半径没有很大影响;不同泄漏孔位置,泄漏初期向上开口时危险半径最大,一段时间后向下开口危险半径最大。     

快速接头内壁冲蚀规律的数值模拟分析

为进一步改善快速接头受连续相作用易失效、内壁冲蚀磨损严重的问题,采用CFD软件中的Realizable k-ε湍流模型,对不同入口速度下的快速接头冲蚀磨损情况运用数值模拟方法分析,得到快速接头内部的流场分布及其冲蚀磨损情况.结果表明:快速接头的冲蚀主要集中在管径缩小段的大小管道交界处,且在管径突变处流体的压力、速度会发...     

基于稳态与非稳态的喷吹管内流场的对比研究

针对等孔径喷吹管,利用CFD软件对其内流场进行数值模拟,得到不同入口流量条件下的流场分布特性。对喷吹管进行非稳态流场模拟,并与稳态模拟结果分析比较,模拟结果误差控制在1%以内,验证了将喷吹管非稳态内流场简化为稳态流场分析具有一定的可行性。沿着喷吹管气流流动方向,喷嘴出口流量分布不均匀,整体有增大的趋势;沿喷吹管轴向气流渐趋均匀,即越靠近端面,喷吹偏心角越小,甚至为0。     

混合药物粉尘输运过程静电特性实验研

兽药制药粉尘加工工艺过程中,由于粉尘颗粒之间或颗粒与设备、管壁之间的碰撞、摩擦,导致装置内部静电荷量积聚,激发静电放电,粉尘燃烧或爆炸的事故频发。实验主要通过包括粉尘与管材摩擦的漏电电流测试和静电放电火花对粉尘云点燃敏感性测试两部分。结果表明:单一药物药粉的静电漏电电流随着管材管径的增大,管长的增长,静电漏电电流逐渐变大;随着倾斜角的增大,静电漏电电流先增大后降低;镀锌铁管的漏电电流大于PVC管,电荷逸散速度更快。单一兽药粉的粉尘云放电火花最小点燃能量随质量浓度的变化,呈现二次曲线的变化趋势。混合兽药粉与单一兽药粉的漏电电流和粉尘云放电火花最小点燃能量的测试结果的变化趋势是一致的。     

Self-ignition of tetrafluoroethylene induced by rapid valve opening in small diameter pipes

This work investigates the ignition of tetrafluoroethylene induced by the adiabatic compression that can arise by activating a high speed valve separating two portions of a pipeline with a high pressure difference. In the tests performed the high pressure zone contained tetrafluoroethylene at pressures between 15 and 30 bar. For the low pressure zone, experiments with nitrogen, air and tetrafluoroethylene were carried out. The pressure range in the low pressure zone was comprised between 0.05 and 1 bar. The pipe diameters analyzed were 15 and 20 mm. For the analyzed geometries, special conditions were required in order to reach reproducible ignitions, namely air at temperatures of at least 105 °C had to be present in the compression pipe. Furthermore, a minimum length of the compression pipe had to be used. The current work describes the experimental setup employed for the tests and discusses the achieved results. Numerical simulations were performed in order to clarify unexpected findings.     

基于COMSOL的采空区埋管抽采参数优化数值研究

针对回采工作面上隅角瓦斯浓度超限问题,提出了回采工作面采空区埋管抽采的方法。以保德煤矿81307工作面为研究对象,运用数值模拟软件COMSOL模拟采空区无抽采和不同抽采参数条件下工作面内瓦斯分布规律,研究埋管抽采参数对上隅角瓦斯浓度的影响规律,确定最佳的采空区埋管抽采参数。同时进行现场抽采参数优化试验,对瓦斯浓度进行监测,研究结果表明:合适的布置间距、抽采负压和抽采流量能够有效解决上隅角瓦斯超限问题,试验期间内,上隅角瓦斯体积分数最大为0.74%,进风流中瓦斯体积分数最大为0.2%,工作面风流中瓦斯体积分数最大为0.45%,回风流中瓦斯体积分数最大为0.5%,均没有超过安全标准。     

The dynamic response of pressure relief valves in vapor or gas service. Part II: Experimental investigation

An experimental testing program designed to evaluate the opening stability characteristics of direct acting pressure relief valves (PRV) in gas/vapor service is described. Three different valve sizes from each of three different manufacturers were tested at two different set pressures to determine their opening characteristics (disk lift vs. time). The valves were tested with several different lengths of inlet piping as well as with and without discharge piping to determine the conditions under which unstable operation (chatter) would occur. Part I of this program described a mathematical model for predicting the dynamic response of PRV's, and the data described in this test program were used to evaluate the accuracy of the model, as described in Part III of this study to follow.     

弯管与盲通管冲蚀磨损对比分析研究

为了研究石油管道流向急剧改变处的冲蚀磨损问题,采用DPM模型,通过改变入口流速、颗粒质量流速、颗粒直径,对90°弯管与盲通管的流场分布及冲蚀情况进行数值模拟。结果表明:弯管与盲通管最大冲蚀速率随入口流速的增大呈指数增长,随颗粒质量流速的增大呈线性增长;在50~100 μm粒径范围内,最大冲蚀速率随粒径的增加逐渐减小,在100~300 μm粒径范围内,随粒径的增加而增大;在入口流速、颗粒质量流速、颗粒直径相同的条件下,弯管最大冲蚀速率明显高于盲通管最大冲蚀速率,盲通管的耐蚀性更强;由于流体在盲通管产生涡旋现象,增加了颗粒能量的耗散,从而减小了进入下游管线颗粒的速率,使得颗粒更易积存于盲通段形成堆积层,减小了下游管段冲蚀速率。     

阻力平衡器在除尘管网中的应用

通过对在除尘系统中作为流体阻尼元件的阻力平衡器和调节阀门的对比分析,阐述了除尘管网中采用阻力平衡器的意义, 并对阻力平衡器的结构、工作原理和设计计算进行了介绍,最后强调了阻力平衡器使用时应注意的事项.     

基于检验测试策略的PFD模型建立与分析*

为制定合理检验测试策略,提高安全仪表系统（SIS）在低要求运行模式下的安全性,提出要求平均失效概率(PFDavg)通用计算模型,引入检验测试分布因子和共因失效修正因子,表征部分和完全检验测试对SIS安全性的影响。结果表明:该模型适用于所有同构koon架构系统,可应用于周期性、非周期性部分检验测试及共因失效影响较大的场景,可以为企业制定检验测试策略提供理论依据。