An experiment on critical flow of water was conducted in two nozzles of 1.41 mm in diameter and 4.35 mm in length with rounded-edge and sharp-edge respectively, covering the ranges of inlet pressure of 22.1-29.1 MPa a...An experiment on critical flow of water was conducted in two nozzles of 1.41 mm in diameter and 4.35 mm in length with rounded-edge and sharp-edge respectively, covering the ranges of inlet pressure of 22.1-29.1 MPa and inlet temperature of 38-474 ℃. More than 200 data points were obtained and the characteristics and parametric trends were investigated. In the region of near and beyond pseudo-critical temperature the thermal-equilibrium is dominant, and the flow rate can be estimated by the modified homogeneous equilibrium model. In the below pseudo-critical region the results exhibit scattered feature as a result of hysteresis effect in the onset of vaporization, characterizing a bifurcation behavior. This effect is more significant in the nozzle with sharp-edge, especially at higher pressure. For temperature well below the pseudo-critical point, the flow is not at critical condition and the flow rate can be represented by the Bernoulli equation reasonably.展开更多
Incipient sediment motion plays a key role in scouring and bed load transport. However, the incipient sediment motion in the vegetated open-channel flows has yet to be fully understood. This study aims to quantify the...Incipient sediment motion plays a key role in scouring and bed load transport. However, the incipient sediment motion in the vegetated open-channel flows has yet to be fully understood. This study aims to quantify the critical conditions of the sediment particle movement in the presence of emergent and submerged vegetation. A new formula of the critical flow velocity was proposed to predict the incipient sediment motion based on the force balance equation of a sediment particle and the assumption that the velocity distribution in the bed roughness boundary layer fits the logarithmic law. Analysis of the derived formula revealed that the critical flow velocity for incipient sediment motion decreases with the increase in vegetation density. The proposed formula agrees well with the experimental data in the literature, thereby implying that the critical flow velocity can effectively quantify the incipient sediment motion in the vegetated open channel flows.展开更多
The response of an adjustable critical-flow Venturi nozzle is investigated through a set indoor experiments aimed to determine the related critical flow rate,critical pressure ratio,and discharge coefficient.The effec...The response of an adjustable critical-flow Venturi nozzle is investigated through a set indoor experiments aimed to determine the related critical flow rate,critical pressure ratio,and discharge coefficient.The effect of a variation in the cone displacement and liquid content on the critical flow characteristics is examined in detail and it is shown that the former can be used to effectively adjust the critical flow rate.The critical pressure ratio of the considered nozzle is above 0.85,and the critical flow control deviation of the gas flow is within±3%.Liquid flow can reduce the gas critical mass flow rate accordingly,especially for the cases with larger liquid volume and lower inlet pressure.The set of results and conclusions provided are intended to support the optimization of steam injection techniques in the context of heavy oil recovery processes.展开更多
Formulated are simple models for the flow in liquid film, formed by impinging jet, and a two-phase downward flow in pipe. The models are based on simplified equations of mass, momentum and energy. The solutions of suc...Formulated are simple models for the flow in liquid film, formed by impinging jet, and a two-phase downward flow in pipe. The models are based on simplified equations of mass, momentum and energy. The solutions of such conservation equations may have regular points belonging to one integral curve only as well as turning points can be found amongst them, which refer to extreme values in the appropriate co-ordinate system. The solutions can also have singular points belonging to none or more than one integral curve. Both the turning and singular points have a clear physical meaning. They could be linked to critical flow conditions in the pipe flow or to the so-called hydraulic jump. Analogy existing between critical conditions in the pipe and the flow of liquid films formed by the liquid jet have been shown in the paper.展开更多
Critical erosionflow rate is the key factor restricting the injection/production capacity of an injection/production well.At present,it is commonly calculated according to API RP 14E standard and its calculation resul...Critical erosionflow rate is the key factor restricting the injection/production capacity of an injection/production well.At present,it is commonly calculated according to API RP 14E standard and its calculation result tends to be conservative.So far,however,there is no definite laboratory experiment orfield data that can prove that critical erosionflow rate can be increased on the basis of API RP 14E.To deal this end,the concept of critical erosionflow rate was proposed based on corrosion rate for thefirst time in this paper.Then,a laboratory equivalent simulation experiment under real injection and production conditions was carried out by comprehensively taking into account the factors influencing string erosion(including temperature,pressure,gas component,water content,solid particle content and string material)while introducing the wall shear stress.Accordingly,the critical erosion coefficient(C)under experimental working conditions was calculated.Finally,a C value chart for three kinds of strings that are commonly used onfield(N80,SM80S and S13Cr)was established.And the following research results are ob-tained.First,solid particle content,water content,CO2differential pressure and wall shear stress are the main erosion controlling factors.Second,solid particle content is the most significant factor that affects the erosion of N80,SM80S and S13Cr strings,and erosion of N80 and SM80S strings is more sensitive to wall shear stress and water content.Third,as for S13Cr string,the C value can be 100 when the solid particle content is lower than 250 mg/L,180 when thefluid contains liquid but no solid particles,and 275 when thefluid is gas phase.Fourth,as for N80 and SM80S strings,the C value can be in the range of 100e180 based on different water content and wall shear stress when thefluid contains liquid but no solid particles,and 275 when thefluid is in gas phase.Fifth,in view that thefluid produced from the injection/production wells of Hutubi gas storage has a water content of 0.0010‰without solid particles,S13Cr is adopted as string material and the C value is set at 180.It is shown in the laboratory erosion experiments that no erosion trace occurs on the string samples under injection and production conditions and the erosion rate is extremely low without point erosion.In conclusion,the C value chart established in this paper is reliable and can provide the guidance for the scientific and reasonable determination of critical erosionflow rate.展开更多
The capability of 1he smoothed-particle hydrodynamics (SPH) method to model supercritical flow in circular pipe bends is considered. The standard SPH method, which makes use of dynamic boundary particles (DBP), is...The capability of 1he smoothed-particle hydrodynamics (SPH) method to model supercritical flow in circular pipe bends is considered. The standard SPH method, which makes use of dynamic boundary particles (DBP), is supplemented with the original algorithm for the treatment of open boundaries. The method is assessed through a comparison with measured free-surface profiles in a pipe bend, and already proposed regression curves for eslimation of the flow-type in a pipe bend. The sensitivity of the model to different parameters is also evaluated. It is shown that an adequate choice of the artificial viscosity coefficient and the initial particle spacing can lead to correct presentation of the flow-type in a bend. Due to easiness of its implementation, the SPH method can he efficiently used in the design of circular conduits with supercritical flow in a bend, such as tunnel spillways, and bottom outlets of dams, or storm sewers.展开更多
Water sensitivity phenomenon occurs during saline aquifer freshening process in seawater intrusion area, and clay particles released in the phenomenon can damage the infiltration capacity of the aquifer. In order to f...Water sensitivity phenomenon occurs during saline aquifer freshening process in seawater intrusion area, and clay particles released in the phenomenon can damage the infiltration capacity of the aquifer. In order to find out the factors and mechanisms for clay particle release, laboratory column infiltration experiments simulating saline aquifer freshening process were designed to measure the critical conditions(critical flow velocity, critical salt concentration and critical ionic strength) and force analysis for clay particle according to DLVO electric double layer theory was employed to illustrate the mechanisms for particle release. The research results showed that critical flow velocity for clay particle release is influenced by salt concentration of injecting solution. When salt concentration of injecting solution is very high, clay particles are not released, indicating that there does not exist a critical flow velocity in this situation. As salt concentration of injecting solution decreases, particles start to be released. The critical salt concentration for clay particle release is 0.052 mol L-1 in our work, which was determined by a constant-flux experiment for stepwise displacement of high concentration Na Cl solution. The critical ionic strength for clay particle release decreases as Ca2+ molar content percentage of the mixed solution of Na Cl and Ca Cl2 increases following the first-order exponential decay equation y = 0.0391e-0.266 x + 0.0015.展开更多
For the Liuzhuang gas storage reservoir with edge and bottom water,excessively high bottomhole pressure of production wells can cause edge and bottom water breakthrough,which reduces the effective storage capacity;whe...For the Liuzhuang gas storage reservoir with edge and bottom water,excessively high bottomhole pressure of production wells can cause edge and bottom water breakthrough,which reduces the effective storage capacity;whereas excessively low bottomhole pressure makes it difficult to achieve the gas production targets.Therefore,determining the reasonable wellhead pressure of production wells is of great significance for the efficient operation of gas storage reservoir.In this study,a comparison of seven critical liquid-carrying flow rate models,including the droplet inversion model and the liquid film inversion model,was made.The Belfroid critical liquid-carrying model was identified as the most accurate for evaluating liquid loading in the Liuzhuang gas storage reservoir.On this basis and in combination with nodal analysis,the relationship curves between gas production rate and critical liquidcarrying flow rate under different reservoir permeabilities and skin factors were obtained,leading to the development of wellhead pressure charts indicating the risk of liquid loading in gas wells of the Liuzhuang gas storage reservoir.This study provides technical support for field operations in diagnosing liquid loading in production wells.展开更多
Serrated flows are known as repeated yielding of bulk metallic glasses(BMGs)during plastic deformation under different loading conditions,which are associated with the operation of shear banding.According to the sta...Serrated flows are known as repeated yielding of bulk metallic glasses(BMGs)during plastic deformation under different loading conditions,which are associated with the operation of shear banding.According to the statistics of some parameters,the shear avalanches can display a self-organized critical state,suggesting a large ductility of BMGs.The emergence of the self-organized criticality(SOC)behavior in different BMGs is due to the temperature,strain rate,and chemical compositions.The SOC behavior is accompanied with the following phenomena:the interactions occur in the shear bands;the incubation time is longer than the relaxation time;the time interval is lacking of typical time scale;and the spatial or temporal parameters should display apower-law distribution.展开更多
Based on isentropic flow and thermal equilibrium assumptions, a model was derived to calculate discharge flow rate, which unified the rules of room temperature water discharge, high temperature and high pressure water...Based on isentropic flow and thermal equilibrium assumptions, a model was derived to calculate discharge flow rate, which unified the rules of room temperature water discharge, high temperature and high pressure water discharge, two-phase critical flow, saturated steam and superheated steam critical flow, and gave a method to calculate critical condition. Because of the influence of friction, the entropy is increased in the actual discharge process, and the discharge flow rate in thermal equilibrium condition can be obtained by the original model multiplied by an appropriate correction coefficient. The model calculated results agreed well with the experiment data of long nozzle critical flow.展开更多
Debris flows,a common natural disaster,cause annual casualties and property damage downstream.In mountainous watersheds,short-duration heavy rainfall generates significant hydraulic conditions in debris flow channels....Debris flows,a common natural disaster,cause annual casualties and property damage downstream.In mountainous watersheds,short-duration heavy rainfall generates significant hydraulic conditions in debris flow channels.The specified conditions trigger the mobilization of loose sediment,leading to the amplification effect that can ultimately overwhelm downstream mitigation measures.This study considers the influence of the fluid properties of viscous debris flow slurries on the forces exerted on buried coarse particles,aiming to accurately reflect real-world conditions and offer theoretical support for the amplification effect.A novel relationship was established between the initiation of motion for buried coarse particles and flow conditions,thereby broadening the applicability of existing models.The reliability of the model is assessed from two perspectives:comparison with physical model experiments and comparison with numerical simulations.A sensitivity analysis on the fluid property parameters of the debris flow slurry revealed that slurry density has the greatest impact on the cross-sectional initiation velocity for particle initiation,followed by fluid yield stress,fluid consistency index,and fluid behavior index,with the latter having the least influence.展开更多
In several previous studies,it was reported that a supported pipe with small geometric imperfections would lose stability when the internal flow velocity became sufficiently high.Recently,however,it has become clear t...In several previous studies,it was reported that a supported pipe with small geometric imperfections would lose stability when the internal flow velocity became sufficiently high.Recently,however,it has become clear that this conclusion may be at best incomplete.A reevaluation of the problem is undertaken here by essentially considering the flow-induced static deformation of a pipe.With the aid of the absolute nodal coordinate formulation(ANCF)and the extended Lagrange equations for dynamical systems containing non-material volumes,the nonlinear governing equations of a pipe with three different geometric imperfections are introduced and formulated.Based on extensive numerical calculations,the static equilibrium configuration,the stability,and the nonlinear dynamics of the considered pipe system are determined and analyzed.The results show that for a supported pipe with the geometric imperfection of a half sinusoidal wave,the dynamical system could not lose stability even if the flow velocity reaches an extremely high value of 40.However,for a supported pipe with the geometric imperfection of one or one and a half sinusoidal waves,the first-mode buckling instability would take place at high flow velocity.Moreover,based on a further parametric analysis,the effects of the amplitude of the geometric imperfection and the aspect ratio of the pipe on the static deformation,the critical flow velocity for buckling instability,and the nonlinear responses of the supported pipes with geometric imperfections are analyzed.展开更多
The stability and dynamical behavior of flexible and articulated rigid pipes conveying fluid have attracted the attention of many researchers in the field of fluid-structure interactions.The system of an articulated p...The stability and dynamical behavior of flexible and articulated rigid pipes conveying fluid have attracted the attention of many researchers in the field of fluid-structure interactions.The system of an articulated pipe composed of a flexible pipe and a rigid pipe is a class of hybrid flexible-rigid dynamical problems involving flow-induced vibrations.This paper establishes the governing equations of motion of a hybrid flexible-rigid pipe system based on Hamilton's principle,with the rigid pipe being hinged to the lower end of a flexible cantilevered pipe via a rotational spring.The coupling equations of motion are discretized via a Galerkin's approach.The mathematical model is validated by comparing the eigenvalue branches of a degenerated system by choosing extreme values of the parameters of the hybrid pipe with previous results.In the theoretical analysis,the critical flow velocities are calculated as a function of the stiffness of the rotational spring,mass ratio and length ratio of the rigid and flexible pipes.The unstable modes are detected from the eigenvalue branches and compared with those of a flexible cantilevered pipe.Numerical results show that the critical flow velocity is greatly influenced by several structural parameters.It is found that a small stiffness of the rotational spring tends to predict higher-mode instability,whereas a large rotational spring stiffness would generate a second-mode instability in most cases.In several system parameter spaces,the hybrid pipe may experience a transference of unstable modes with the increase of flow velocity.It is also shown that the hybrid pipe system may lose stability first in the fourth mode in some cases.Some of the fresh results obtained for the hybrid pipe system are expected to be helpful in understanding and controlling the dynamical responses of hybrid flexible-rigid fluid-conveying pipes.展开更多
Cantilevered pipe conveying fluid may become unstable and flutter instability would occur when the velocity of the fluid flow in the pipe exceeds a critical value.In the present study,the theoretical model of a cantil...Cantilevered pipe conveying fluid may become unstable and flutter instability would occur when the velocity of the fluid flow in the pipe exceeds a critical value.In the present study,the theoretical model of a cantilevered fluid-conveying pipe attached by an inerter-based dynamic vibration absorber(IDVA)is proposed and the stability of this dynamical system is explored.Based on linear governing equations of the pipe and the IDVA,the effects of damping coefficient,weight,inerter,location and spring stiffness of the IDVAon the critical flow velocities of the pipe system is examined.It is shown that the stability of the pipe may be significantly affected by the IDVA.In many cases,the stability of the cantilevered pipe can be enhanced by designing the parameter values of the IDVA.By solving nonlinear governing equations of the dynamical system,the nonlinear oscillations of the pipe with IDVA for sufficiently high flow velocity beyond the critical value are determined,showing that the oscillation amplitudes of the pipe can also be suppressed to some extent with a suitable design of the IDVA.展开更多
The wave method is introduced to vibration analysis of the fluid-conveying carbon nanotube. The constitutive relation of carbon nanotube on micro-scale is founded using the non- local elastic theory. The governing equ...The wave method is introduced to vibration analysis of the fluid-conveying carbon nanotube. The constitutive relation of carbon nanotube on micro-scale is founded using the non- local elastic theory. The governing equation on micro-scale is obtained. And the first five orders of the natural frequency of the carbon nanotube conveying fluid with various speeds are calculated through the wave method. Besides, the critical flow velocity when the carbon nanotube loses stability is obtained. Meanwhile, a contrast is made between the result obtained through tile wave method and that in previous researches.展开更多
The transient critical flow experiment with sharp-edged tubes as the break geometries is conducted in high pressure convective circulation test loop of Xi'an Jiaotong University. The initial steady operation press...The transient critical flow experiment with sharp-edged tubes as the break geometries is conducted in high pressure convective circulation test loop of Xi'an Jiaotong University. The initial steady operation pressure is up to 22.OMPa. An empirical correlation was made to obtain the critical mass flow rates, the critical pressure ratio and the thermal nonequilibrium number were correlated as the functions of the tube length to tube diameter ratio L/D. The predicted critical mass now rate gets a higher accuracy for short tubes with L/D 12.展开更多
A theoretical model is developed for the vibration and stability of a vertical pipe subjected concurrently to two dependent axial flows. The external fluid, after exiting the outer annular region between the pipe and ...A theoretical model is developed for the vibration and stability of a vertical pipe subjected concurrently to two dependent axial flows. The external fluid, after exiting the outer annular region between the pipe and a rigid cylindrical channel, is conveyed upwards inside the pipe. This configuration thus resembles of a pipe that aspirating fluid. The equation of planar mo- tion is solved by means of the differential quadrature method (DQM). Calculations are conducted for a slender drill-string-like and a bench-top-size system, for different confinement conditions of the outer annular channel. It is shown that the vibrations of these two systems are closely related to the degree of confinement of the outer annular channel. For a drill-string-like system with narrow annuli, buckling instability may occur in the second and third modes. For a bench-top-size system, however, both buckling and flutter may occur in the lowest three modes. The form of instability depends on the annuli size.展开更多
When a horizontal well is drilled through the natural gas hydrate(NGH)strata with different abundance,the NGH-bearing solid debris tends to lead to the deposit and adhesion of solid particles in the lateral sections a...When a horizontal well is drilled through the natural gas hydrate(NGH)strata with different abundance,the NGH-bearing solid debris tends to lead to the deposit and adhesion of solid particles in the lateral sections and consequently cuttings carrying is not smooth.In this paper,the critical return rate model of drilling fluid in the case of cuttings rolling(normal drilling)and saltation(pump off for sand settling)under NGH adhesion in the laterals was established according to the force and migration laws of cuttings under the condition of multiphase flow in a horizontal well and the particle migration theory.Then,numerical simulation was conducted and the influential factors and migration laws of cuttings starting in the case of normal drilling and pump off were analyzed.And the following results were obtained.First,the critical starting flow rate decreases with the rise of NGH abundance and it is higher when the NGH adhesion is taken into account.Besides,the higher the NGH abundance,the greater the effect of NGH adhesion.Second,the critical starting flow rate increases with the rise of drilling cuttings particle size when the NGH abundance is less than 85%;and it decreases with the rise of drilling cuttings particle size when the NGH abundance is more than 85%.Third,the critical starting flow rate decreases with the rise of drilling fluid density and viscosity.And fourth,under the same conditions,the critical return rate for saltation is about 1.28 times that for rolling.It is suggested that the rolling model should be adopted for the normal drilling while the saltation model for the recycle after pump-off sand settling.The research results are of great significance to the optimization of NGH drilling parameters and the reduction of drilling safety risks.展开更多
China’s unconventional gas fields have a large number of low-productivity and low-efficiency wells, many of whichare located in remote and environmentally harsh mountainous areas. To address the long-term stable prod...China’s unconventional gas fields have a large number of low-productivity and low-efficiency wells, many of whichare located in remote and environmentally harsh mountainous areas. To address the long-term stable productionof these gas wells, plunger-lift technology plays an important role. In order to fully understand and accurately graspthe drainage and gas production mechanisms of plunger-lift, a mechanical model of plunger-liquid column uplift inthe plunger-lift process was established, focusing on conventional plunger-lift systems and representative wellboreconfigurations in the Linxing region. The operating casing pressure of the plunger-lift process and the calculationmethod for the maximum daily fluid production rate based on the work regime with the highest fluid recovery ratewere determined. For the first time, the critical flow rate method was proposed as a constraint for the maximumliquid-carrying capacity of the plunger-lift, and liquid-carrying capacity charts for conventional plunger-lift withdifferent casing sizes were developed. The results showed that for 23/8 casing plunger-lift, with a well depth ofshallower than 808 m, the maximum drainage rate was 33 m3/d;for 27/8 casing plunger-lift, with a well depth ofshallower than 742 m, the maximum drainage rate was 50.15 m3/d;for 31/2 casing plunger-lift, with a well depthof shallower than 560 m, the maximum drainage rate was 75.14 m3/d. This research provides a foundation for thescientific selection of plunger-lift technology and serves as a decision-making reference for developing reasonableplunger-lift work regimes.展开更多
This paper introduces how crustal thickening controls the growth of the Himalaya by summarizing the P-T-t evolution of the Himalayan metamorphic core.The Himalayan orogeny was divided into three stages.Stage 60–40 Ma...This paper introduces how crustal thickening controls the growth of the Himalaya by summarizing the P-T-t evolution of the Himalayan metamorphic core.The Himalayan orogeny was divided into three stages.Stage 60–40 Ma:The Himalayan crust thickened to~40 km through Barrovian-type metamorphism(15–25°C/km),and the Himalaya rose from30°C/km).The three sub-sheets in the Himalayan metamorphic core extruded southward sequentially through imbricate thrusts of the Eo-Himalayan thrust,High Himalayan thrust,and Main Central thrust,and the Himalaya rose to≥5,000 m.Stage 16–0 Ma:the mountain roots underwent localized delamination,causing asthenospheric upwelling and overprinting of the lower crust by ultra-high-temperature metamorphism(30–50°C/km),and the Himalaya reached the present elevation of~6,000 m.Underplating and imbricate thrusting dominated the Himalaya’growth and topographic rise,conforming to the critical taper wedge model.Localized delamination of mountain roots facilitated further topographic rise.Future Himalayan metamorphic studies should focus on extreme metamorphism and major collisional events,contact metamorphism and rare metal mineralization,metamorphic decarbonation and the carbon cycle in collisional belts.展开更多
摘要An experiment on critical flow of water was conducted in two nozzles of 1.41 mm in diameter and 4.35 mm in length with rounded-edge and sharp-edge respectively, covering the ranges of inlet pressure of 22.1-29.1 MPa and inlet temperature of 38-474 ℃. More than 200 data points were obtained and the characteristics and parametric trends were investigated. In the region of near and beyond pseudo-critical temperature the thermal-equilibrium is dominant, and the flow rate can be estimated by the modified homogeneous equilibrium model. In the below pseudo-critical region the results exhibit scattered feature as a result of hysteresis effect in the onset of vaporization, characterizing a bifurcation behavior. This effect is more significant in the nozzle with sharp-edge, especially at higher pressure. For temperature well below the pseudo-critical point, the flow is not at critical condition and the flow rate can be represented by the Bernoulli equation reasonably.
基金supported by the National Natural Science Foundation of China(Grant Nos.52020105006,11872285).
摘要Incipient sediment motion plays a key role in scouring and bed load transport. However, the incipient sediment motion in the vegetated open-channel flows has yet to be fully understood. This study aims to quantify the critical conditions of the sediment particle movement in the presence of emergent and submerged vegetation. A new formula of the critical flow velocity was proposed to predict the incipient sediment motion based on the force balance equation of a sediment particle and the assumption that the velocity distribution in the bed roughness boundary layer fits the logarithmic law. Analysis of the derived formula revealed that the critical flow velocity for incipient sediment motion decreases with the increase in vegetation density. The proposed formula agrees well with the experimental data in the literature, thereby implying that the critical flow velocity can effectively quantify the incipient sediment motion in the vegetated open channel flows.
基金The authors would like to acknowledge the support provided by the National Natural Science Foundation of China(No.62173049)the open fund of the Key Laboratory of Exploration Technologies for Oil and Gas Resources(Yangtze University),Ministry of Education(Grant K2021-17).
摘要The response of an adjustable critical-flow Venturi nozzle is investigated through a set indoor experiments aimed to determine the related critical flow rate,critical pressure ratio,and discharge coefficient.The effect of a variation in the cone displacement and liquid content on the critical flow characteristics is examined in detail and it is shown that the former can be used to effectively adjust the critical flow rate.The critical pressure ratio of the considered nozzle is above 0.85,and the critical flow control deviation of the gas flow is within±3%.Liquid flow can reduce the gas critical mass flow rate accordingly,especially for the cases with larger liquid volume and lower inlet pressure.The set of results and conclusions provided are intended to support the optimization of steam injection techniques in the context of heavy oil recovery processes.
摘要Formulated are simple models for the flow in liquid film, formed by impinging jet, and a two-phase downward flow in pipe. The models are based on simplified equations of mass, momentum and energy. The solutions of such conservation equations may have regular points belonging to one integral curve only as well as turning points can be found amongst them, which refer to extreme values in the appropriate co-ordinate system. The solutions can also have singular points belonging to none or more than one integral curve. Both the turning and singular points have a clear physical meaning. They could be linked to critical flow conditions in the pipe flow or to the so-called hydraulic jump. Analogy existing between critical conditions in the pipe and the flow of liquid films formed by the liquid jet have been shown in the paper.
基金supported by PetroChina Major Scientifie Research&Techno-logical Development Project"Research on Key Technologies of Injection and Production Engincering of Underground Gas Storage"(No:2015E-4004).
摘要Critical erosionflow rate is the key factor restricting the injection/production capacity of an injection/production well.At present,it is commonly calculated according to API RP 14E standard and its calculation result tends to be conservative.So far,however,there is no definite laboratory experiment orfield data that can prove that critical erosionflow rate can be increased on the basis of API RP 14E.To deal this end,the concept of critical erosionflow rate was proposed based on corrosion rate for thefirst time in this paper.Then,a laboratory equivalent simulation experiment under real injection and production conditions was carried out by comprehensively taking into account the factors influencing string erosion(including temperature,pressure,gas component,water content,solid particle content and string material)while introducing the wall shear stress.Accordingly,the critical erosion coefficient(C)under experimental working conditions was calculated.Finally,a C value chart for three kinds of strings that are commonly used onfield(N80,SM80S and S13Cr)was established.And the following research results are ob-tained.First,solid particle content,water content,CO2differential pressure and wall shear stress are the main erosion controlling factors.Second,solid particle content is the most significant factor that affects the erosion of N80,SM80S and S13Cr strings,and erosion of N80 and SM80S strings is more sensitive to wall shear stress and water content.Third,as for S13Cr string,the C value can be 100 when the solid particle content is lower than 250 mg/L,180 when thefluid contains liquid but no solid particles,and 275 when thefluid is gas phase.Fourth,as for N80 and SM80S strings,the C value can be in the range of 100e180 based on different water content and wall shear stress when thefluid contains liquid but no solid particles,and 275 when thefluid is in gas phase.Fifth,in view that thefluid produced from the injection/production wells of Hutubi gas storage has a water content of 0.0010‰without solid particles,S13Cr is adopted as string material and the C value is set at 180.It is shown in the laboratory erosion experiments that no erosion trace occurs on the string samples under injection and production conditions and the erosion rate is extremely low without point erosion.In conclusion,the C value chart established in this paper is reliable and can provide the guidance for the scientific and reasonable determination of critical erosionflow rate.
基金financial support through projects TR 37009,TR 37010,granted by the Serbian Ministry of Education,Science and Technological Development
摘要The capability of 1he smoothed-particle hydrodynamics (SPH) method to model supercritical flow in circular pipe bends is considered. The standard SPH method, which makes use of dynamic boundary particles (DBP), is supplemented with the original algorithm for the treatment of open boundaries. The method is assessed through a comparison with measured free-surface profiles in a pipe bend, and already proposed regression curves for eslimation of the flow-type in a pipe bend. The sensitivity of the model to different parameters is also evaluated. It is shown that an adequate choice of the artificial viscosity coefficient and the initial particle spacing can lead to correct presentation of the flow-type in a bend. Due to easiness of its implementation, the SPH method can he efficiently used in the design of circular conduits with supercritical flow in a bend, such as tunnel spillways, and bottom outlets of dams, or storm sewers.
基金supported by the National Natural Science Foundation of China (Grant No. 41172209)National Public Welfare Scientific Research Project (Grant No. 201301090)
摘要Water sensitivity phenomenon occurs during saline aquifer freshening process in seawater intrusion area, and clay particles released in the phenomenon can damage the infiltration capacity of the aquifer. In order to find out the factors and mechanisms for clay particle release, laboratory column infiltration experiments simulating saline aquifer freshening process were designed to measure the critical conditions(critical flow velocity, critical salt concentration and critical ionic strength) and force analysis for clay particle according to DLVO electric double layer theory was employed to illustrate the mechanisms for particle release. The research results showed that critical flow velocity for clay particle release is influenced by salt concentration of injecting solution. When salt concentration of injecting solution is very high, clay particles are not released, indicating that there does not exist a critical flow velocity in this situation. As salt concentration of injecting solution decreases, particles start to be released. The critical salt concentration for clay particle release is 0.052 mol L-1 in our work, which was determined by a constant-flux experiment for stepwise displacement of high concentration Na Cl solution. The critical ionic strength for clay particle release decreases as Ca2+ molar content percentage of the mixed solution of Na Cl and Ca Cl2 increases following the first-order exponential decay equation y = 0.0391e-0.266 x + 0.0015.
基金funded by the National Science and Technology Major Project(No.2024ZD1406605).
摘要For the Liuzhuang gas storage reservoir with edge and bottom water,excessively high bottomhole pressure of production wells can cause edge and bottom water breakthrough,which reduces the effective storage capacity;whereas excessively low bottomhole pressure makes it difficult to achieve the gas production targets.Therefore,determining the reasonable wellhead pressure of production wells is of great significance for the efficient operation of gas storage reservoir.In this study,a comparison of seven critical liquid-carrying flow rate models,including the droplet inversion model and the liquid film inversion model,was made.The Belfroid critical liquid-carrying model was identified as the most accurate for evaluating liquid loading in the Liuzhuang gas storage reservoir.On this basis and in combination with nodal analysis,the relationship curves between gas production rate and critical liquidcarrying flow rate under different reservoir permeabilities and skin factors were obtained,leading to the development of wellhead pressure charts indicating the risk of liquid loading in gas wells of the Liuzhuang gas storage reservoir.This study provides technical support for field operations in diagnosing liquid loading in production wells.
基金Item Sponsored by National Natural Science Foundation of China(51371122)Program for the Innovative Talents of Higher Learning Institutions of Shanxi of China(2013)The Youth Natural Science Foundation of Shanxi Province of China(2015021005)
摘要Serrated flows are known as repeated yielding of bulk metallic glasses(BMGs)during plastic deformation under different loading conditions,which are associated with the operation of shear banding.According to the statistics of some parameters,the shear avalanches can display a self-organized critical state,suggesting a large ductility of BMGs.The emergence of the self-organized criticality(SOC)behavior in different BMGs is due to the temperature,strain rate,and chemical compositions.The SOC behavior is accompanied with the following phenomena:the interactions occur in the shear bands;the incubation time is longer than the relaxation time;the time interval is lacking of typical time scale;and the spatial or temporal parameters should display apower-law distribution.
摘要Based on isentropic flow and thermal equilibrium assumptions, a model was derived to calculate discharge flow rate, which unified the rules of room temperature water discharge, high temperature and high pressure water discharge, two-phase critical flow, saturated steam and superheated steam critical flow, and gave a method to calculate critical condition. Because of the influence of friction, the entropy is increased in the actual discharge process, and the discharge flow rate in thermal equilibrium condition can be obtained by the original model multiplied by an appropriate correction coefficient. The model calculated results agreed well with the experiment data of long nozzle critical flow.
基金supported by the National Natural Science Foundation of China(Key Project)(grant No.42230702).
摘要Debris flows,a common natural disaster,cause annual casualties and property damage downstream.In mountainous watersheds,short-duration heavy rainfall generates significant hydraulic conditions in debris flow channels.The specified conditions trigger the mobilization of loose sediment,leading to the amplification effect that can ultimately overwhelm downstream mitigation measures.This study considers the influence of the fluid properties of viscous debris flow slurries on the forces exerted on buried coarse particles,aiming to accurately reflect real-world conditions and offer theoretical support for the amplification effect.A novel relationship was established between the initiation of motion for buried coarse particles and flow conditions,thereby broadening the applicability of existing models.The reliability of the model is assessed from two perspectives:comparison with physical model experiments and comparison with numerical simulations.A sensitivity analysis on the fluid property parameters of the debris flow slurry revealed that slurry density has the greatest impact on the cross-sectional initiation velocity for particle initiation,followed by fluid yield stress,fluid consistency index,and fluid behavior index,with the latter having the least influence.
基金supported by the National Natural Science Foundation of China(Nos.11972167,12072119)the Alexander von Humboldt Foundation。
摘要In several previous studies,it was reported that a supported pipe with small geometric imperfections would lose stability when the internal flow velocity became sufficiently high.Recently,however,it has become clear that this conclusion may be at best incomplete.A reevaluation of the problem is undertaken here by essentially considering the flow-induced static deformation of a pipe.With the aid of the absolute nodal coordinate formulation(ANCF)and the extended Lagrange equations for dynamical systems containing non-material volumes,the nonlinear governing equations of a pipe with three different geometric imperfections are introduced and formulated.Based on extensive numerical calculations,the static equilibrium configuration,the stability,and the nonlinear dynamics of the considered pipe system are determined and analyzed.The results show that for a supported pipe with the geometric imperfection of a half sinusoidal wave,the dynamical system could not lose stability even if the flow velocity reaches an extremely high value of 40.However,for a supported pipe with the geometric imperfection of one or one and a half sinusoidal waves,the first-mode buckling instability would take place at high flow velocity.Moreover,based on a further parametric analysis,the effects of the amplitude of the geometric imperfection and the aspect ratio of the pipe on the static deformation,the critical flow velocity for buckling instability,and the nonlinear responses of the supported pipes with geometric imperfections are analyzed.
基金supported by the National Natural Science Foundation of China(Grant Nos.11902112,11972167,and 12072119)Hubei Superior and Distinctive Discipline Group of"Mechatronics and Automobiles"(Grant No.XKQ2021042).
摘要The stability and dynamical behavior of flexible and articulated rigid pipes conveying fluid have attracted the attention of many researchers in the field of fluid-structure interactions.The system of an articulated pipe composed of a flexible pipe and a rigid pipe is a class of hybrid flexible-rigid dynamical problems involving flow-induced vibrations.This paper establishes the governing equations of motion of a hybrid flexible-rigid pipe system based on Hamilton's principle,with the rigid pipe being hinged to the lower end of a flexible cantilevered pipe via a rotational spring.The coupling equations of motion are discretized via a Galerkin's approach.The mathematical model is validated by comparing the eigenvalue branches of a degenerated system by choosing extreme values of the parameters of the hybrid pipe with previous results.In the theoretical analysis,the critical flow velocities are calculated as a function of the stiffness of the rotational spring,mass ratio and length ratio of the rigid and flexible pipes.The unstable modes are detected from the eigenvalue branches and compared with those of a flexible cantilevered pipe.Numerical results show that the critical flow velocity is greatly influenced by several structural parameters.It is found that a small stiffness of the rotational spring tends to predict higher-mode instability,whereas a large rotational spring stiffness would generate a second-mode instability in most cases.In several system parameter spaces,the hybrid pipe may experience a transference of unstable modes with the increase of flow velocity.It is also shown that the hybrid pipe system may lose stability first in the fourth mode in some cases.Some of the fresh results obtained for the hybrid pipe system are expected to be helpful in understanding and controlling the dynamical responses of hybrid flexible-rigid fluid-conveying pipes.
基金The authors gratefully acknowledge the support provided by the National Natural Science Foundation of China(Nos.11622216,11672115 and 11972167).
摘要Cantilevered pipe conveying fluid may become unstable and flutter instability would occur when the velocity of the fluid flow in the pipe exceeds a critical value.In the present study,the theoretical model of a cantilevered fluid-conveying pipe attached by an inerter-based dynamic vibration absorber(IDVA)is proposed and the stability of this dynamical system is explored.Based on linear governing equations of the pipe and the IDVA,the effects of damping coefficient,weight,inerter,location and spring stiffness of the IDVAon the critical flow velocities of the pipe system is examined.It is shown that the stability of the pipe may be significantly affected by the IDVA.In many cases,the stability of the cantilevered pipe can be enhanced by designing the parameter values of the IDVA.By solving nonlinear governing equations of the dynamical system,the nonlinear oscillations of the pipe with IDVA for sufficiently high flow velocity beyond the critical value are determined,showing that the oscillation amplitudes of the pipe can also be suppressed to some extent with a suitable design of the IDVA.
基金the support of a grant from Aeronautical Science Foundation of China(2010ZA53013 and 2011ZA53014)the open funds of Key Laboratory of Advanced Design and Intelligent Computing(Dalian University),Ministry of Education(ADIC2010007)Northwestern Polytechnical University Basic Research Fund(JC201114 andJC20110255)
摘要The wave method is introduced to vibration analysis of the fluid-conveying carbon nanotube. The constitutive relation of carbon nanotube on micro-scale is founded using the non- local elastic theory. The governing equation on micro-scale is obtained. And the first five orders of the natural frequency of the carbon nanotube conveying fluid with various speeds are calculated through the wave method. Besides, the critical flow velocity when the carbon nanotube loses stability is obtained. Meanwhile, a contrast is made between the result obtained through tile wave method and that in previous researches.
摘要The transient critical flow experiment with sharp-edged tubes as the break geometries is conducted in high pressure convective circulation test loop of Xi'an Jiaotong University. The initial steady operation pressure is up to 22.OMPa. An empirical correlation was made to obtain the critical mass flow rates, the critical pressure ratio and the thermal nonequilibrium number were correlated as the functions of the tube length to tube diameter ratio L/D. The predicted critical mass now rate gets a higher accuracy for short tubes with L/D 12.
基金supported by the National Natural Science Foundation of China (Nos. 10772071 and 10802031)theScientific Research Foundation of HUST (No. 2006Q003B).
摘要A theoretical model is developed for the vibration and stability of a vertical pipe subjected concurrently to two dependent axial flows. The external fluid, after exiting the outer annular region between the pipe and a rigid cylindrical channel, is conveyed upwards inside the pipe. This configuration thus resembles of a pipe that aspirating fluid. The equation of planar mo- tion is solved by means of the differential quadrature method (DQM). Calculations are conducted for a slender drill-string-like and a bench-top-size system, for different confinement conditions of the outer annular channel. It is shown that the vibrations of these two systems are closely related to the degree of confinement of the outer annular channel. For a drill-string-like system with narrow annuli, buckling instability may occur in the second and third modes. For a bench-top-size system, however, both buckling and flutter may occur in the lowest three modes. The form of instability depends on the annuli size.
基金Project supported by the National Major R&D Program“New technology for test of marine natural gas hydrates through solid fluidization”(No.2016YFC0304008)the Key Program of National Natural Science Foundation of China“Theory and key issues for measurement and control of managed pressure drilling”(No.51334003).
摘要When a horizontal well is drilled through the natural gas hydrate(NGH)strata with different abundance,the NGH-bearing solid debris tends to lead to the deposit and adhesion of solid particles in the lateral sections and consequently cuttings carrying is not smooth.In this paper,the critical return rate model of drilling fluid in the case of cuttings rolling(normal drilling)and saltation(pump off for sand settling)under NGH adhesion in the laterals was established according to the force and migration laws of cuttings under the condition of multiphase flow in a horizontal well and the particle migration theory.Then,numerical simulation was conducted and the influential factors and migration laws of cuttings starting in the case of normal drilling and pump off were analyzed.And the following results were obtained.First,the critical starting flow rate decreases with the rise of NGH abundance and it is higher when the NGH adhesion is taken into account.Besides,the higher the NGH abundance,the greater the effect of NGH adhesion.Second,the critical starting flow rate increases with the rise of drilling cuttings particle size when the NGH abundance is less than 85%;and it decreases with the rise of drilling cuttings particle size when the NGH abundance is more than 85%.Third,the critical starting flow rate decreases with the rise of drilling fluid density and viscosity.And fourth,under the same conditions,the critical return rate for saltation is about 1.28 times that for rolling.It is suggested that the rolling model should be adopted for the normal drilling while the saltation model for the recycle after pump-off sand settling.The research results are of great significance to the optimization of NGH drilling parameters and the reduction of drilling safety risks.
基金the Fundamental Research Funds for the Central Universities of China(No.20CX02308A)CNOOC Project(No.ZX2022ZCCYF3835).
摘要China’s unconventional gas fields have a large number of low-productivity and low-efficiency wells, many of whichare located in remote and environmentally harsh mountainous areas. To address the long-term stable productionof these gas wells, plunger-lift technology plays an important role. In order to fully understand and accurately graspthe drainage and gas production mechanisms of plunger-lift, a mechanical model of plunger-liquid column uplift inthe plunger-lift process was established, focusing on conventional plunger-lift systems and representative wellboreconfigurations in the Linxing region. The operating casing pressure of the plunger-lift process and the calculationmethod for the maximum daily fluid production rate based on the work regime with the highest fluid recovery ratewere determined. For the first time, the critical flow rate method was proposed as a constraint for the maximumliquid-carrying capacity of the plunger-lift, and liquid-carrying capacity charts for conventional plunger-lift withdifferent casing sizes were developed. The results showed that for 23/8 casing plunger-lift, with a well depth ofshallower than 808 m, the maximum drainage rate was 33 m3/d;for 27/8 casing plunger-lift, with a well depth ofshallower than 742 m, the maximum drainage rate was 50.15 m3/d;for 31/2 casing plunger-lift, with a well depthof shallower than 560 m, the maximum drainage rate was 75.14 m3/d. This research provides a foundation for thescientific selection of plunger-lift technology and serves as a decision-making reference for developing reasonableplunger-lift work regimes.
基金supported by the Second Tibetan Plateau Scientific Expedition and Research program(Grant No.2019QZKK0703)the National Natural Science Foundation of China(Grant Nos.41972065 and 41888101)+1 种基金the Youth Innovation Promotion Association of the Chinese Academy of Sciences(Grant No.2022065)the State Key Laboratory of Lithospheric Evolution(Grant No.E152510201).
摘要This paper introduces how crustal thickening controls the growth of the Himalaya by summarizing the P-T-t evolution of the Himalayan metamorphic core.The Himalayan orogeny was divided into three stages.Stage 60–40 Ma:The Himalayan crust thickened to~40 km through Barrovian-type metamorphism(15–25°C/km),and the Himalaya rose from30°C/km).The three sub-sheets in the Himalayan metamorphic core extruded southward sequentially through imbricate thrusts of the Eo-Himalayan thrust,High Himalayan thrust,and Main Central thrust,and the Himalaya rose to≥5,000 m.Stage 16–0 Ma:the mountain roots underwent localized delamination,causing asthenospheric upwelling and overprinting of the lower crust by ultra-high-temperature metamorphism(30–50°C/km),and the Himalaya reached the present elevation of~6,000 m.Underplating and imbricate thrusting dominated the Himalaya’growth and topographic rise,conforming to the critical taper wedge model.Localized delamination of mountain roots facilitated further topographic rise.Future Himalayan metamorphic studies should focus on extreme metamorphism and major collisional events,contact metamorphism and rare metal mineralization,metamorphic decarbonation and the carbon cycle in collisional belts.