Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat ...Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.展开更多
The focus of this paper is on the measurement and calculation model of void fraction for the vertical upward co-current air-water slug flow in a circular tube of 15 mm inner diameter. High-speed photography and optica...The focus of this paper is on the measurement and calculation model of void fraction for the vertical upward co-current air-water slug flow in a circular tube of 15 mm inner diameter. High-speed photography and optical probes were utilized, with water superficial velocity ranging from 0.089 to 0.65 m·s^(-1)and gas superficial velocity ranging from 0.049 to 0.65 m·s^(-1). A new void fraction model based on the local parameters was proposed, disposing the slug flow as a combination of Taylor bubbles and liquid slugs. In the Taylor bubble region, correction factors of liquid film thickness Cδand nose shape CZ*were proposed to calculate aTB. In the liquid slug region, the radial void fraction distribution profiles were obtained to calculate aLS, by employing the image processing technique based on supervised machine learning. Results showed that the void fraction proportion in Taylor bubbles occupied crucial contribution to the overall void fraction. Multiple types of void fraction predictive correlations were assessed using the present data. The performance of the Schmidt model was optimal, while some models for slug flow performed not outstanding. Additionally, a predictive correlation was correlated between the central local void fraction and the cross-sectional averaged void fraction, as a straightforward form of the void fraction calculation model. The predictive correlation showed a good agreement with the present experimental data, as well as the data of Olerni et al., indicating that the new model was effective and applicable under the slug flow conditions.展开更多
The hydrodynamic study of the liquid film around Taylor bubbles in slug flow has great significance for understanding parallel flow and interaction between Taylor bubbles.The prediction models for liquid film thicknes...The hydrodynamic study of the liquid film around Taylor bubbles in slug flow has great significance for understanding parallel flow and interaction between Taylor bubbles.The prediction models for liquid film thickness mainly focus on stagnant flow,and some of them remain inaccurate performance.However,in the industrial process,the slug flow essentially is co-current flow.Therefore,in this paper,the liquid film thickness is studied by theoretical analysis and experimental methods under two conditions of stagnant and co-current flow.Firstly,under the condition of stagnant flow,the present work is based on Batchelor's theory,and modifies Batchelor's liquid film thickness model,which effectively improves its prediction accuracy.Under the condition of co-current flow,the prediction model of average liquid film thickness in slug flow is established by force and motion analysis.Taylor bubble length is introduced into the model as an important parameter.Dynamic experiments were carried out in the pipe with an inner diameter of 20 mm.The liquid film thickness,Taylor bubble velocity and length were measured by distributed ultrasonic sensor and intrusive cross-correlation conductivity sensor.Comparing the predicted value of the model with the measured results,the relative error is controlled within 10%.展开更多
Several applications such as liquid-liquid extraction in micro-fluidic devices are concerned with the flow of two immiscible liquid phases. The commonly observed flow regimes in these systems are slug-flow and stratif...Several applications such as liquid-liquid extraction in micro-fluidic devices are concerned with the flow of two immiscible liquid phases. The commonly observed flow regimes in these systems are slug-flow and stratified flow. The latter regime in micro-channels has the inherent advantage that separation of the two liquids at the exit is efficient. Recently extraction in a stratified counter-current flow has been studied experimentally and it has been shown to be more efficient than co-current flow. An analytical as well as a numerical method to determine the steady-state solution of the corresponding convection-diffusion equation for the two flow-fields is presented. It is shown that the counter-current process is superior to the co-current process for the same set of parameters and operating conditions. A simplified model is proposed to analyse the process when diffusion in the transverse direction is not rate limiting. Different approaches to determining mass transfer coefficient are compared. The concept of log mean temperature difference used in design of heat exchangers is extended to describe mass transfer in the system.展开更多
A detached eddy simulation(DES) and a k-ε-based Reynolds-averaged Navier–Stokes(RANS) calculation on the co-current spray drying chamber is presented. The DES used here is based on the Spalart–Allmaras(SA) turbulen...A detached eddy simulation(DES) and a k-ε-based Reynolds-averaged Navier–Stokes(RANS) calculation on the co-current spray drying chamber is presented. The DES used here is based on the Spalart–Allmaras(SA) turbulence model, whereas the standard k-ε(SKE) was considered here for comparison purposes. Predictions of the mean axial velocity, temperature and humidity profile have been evaluated and compared with experimental measurements. The effects of the turbulence model on the predictions of the mean axial velocity, temperature and the humidity profile are most noticeable in the(highly anisotropic) spraying region. The findings suggest that DES provide a more accurate prediction(with error less than 5%) of the flow field in a spray drying chamber compared with RANS-based k-ε models. The DES simulation also confirmed the presence of anisotropic turbulent flow in the spray dryer from the analysis of the velocity component fluctuations and turbulent structure as illustrated by the Q-criterion.展开更多
A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes...A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.展开更多
The flow field near a spur dike such as down flow and horseshoe vortex system(HVS)are susceptible to the topographic changes in the local scouring process,resulting in variation of the sediment transport with time.In ...The flow field near a spur dike such as down flow and horseshoe vortex system(HVS)are susceptible to the topographic changes in the local scouring process,resulting in variation of the sediment transport with time.In this study,large eddy simulations with fixed-bed at different scouring stages were conducted to investigate the changes in flow field.The results imply that the bed deformation leads to an increase in flow rate per unit area,which represent the capability of sediment transportation by water,in the scour hole.Moreover,the intensity of turbulent kinetic energy and bimodal motion near the sand bed induced by the HVS were also varied.However,the peak moments between the two sediment transport mechanisms were different.Hence,understanding the complex feedback mechanism between topography and flow field is essential for the local scour problem.展开更多
The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study are...The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study area which is extending along Karakorum Highway(KKH) from Besham to Chilas. Intense seismicity, deep gorges, steep terrain and extreme climatic events trigger multiple mountain hazards along the KKH, among which debris flow is recognized as the most destructive geohazard. This study aims to prepare a field-based debris flow inventory map at a regional scale along a 200 km stretch from Besham to Chilas. A total of 117 debris flows were identified in the field, and subsequently, a point-based debris-flow inventory and catchment delineation were performed through Arc GIS analysis. Regional scale debris flow susceptibility and propagation maps were prepared using Weighted Overlay Method(WOM) and Flow-R technique sequentially. Predisposing factors include slope, slope aspect, elevation, Topographic Roughness Index(TRI), Topographic Wetness Index(TWI), stream buffer, distance to faults, lithology rainfall, curvature, and collapsed material layer. The dataset was randomly divided into training data(75%) and validation data(25%). Results were validated through the Receiver Operator Characteristics(ROC) curve. Results show that Area Under the Curve(AUC) using WOM model is 79.2%. Flow-R propagation of debris flow shows that the 13.15%, 22.94%, and 63.91% areas are very high, high, and low susceptible to debris flow respectively. The propagation predicated by Flow-R validates the naturally occurring debris flow propagation as observed in the field surveys. The output of this research will provide valuable input to the decision makers for the site selection, designing of the prevention system, and for the protection of current infrastructure.展开更多
Brazilian pre-salt reservoirs are renowned for their intricate pore networks and vuggy nature,posing significant challenges in modeling and simulating fluid flow within these carbonate reservoirs.Despite possessing ex...Brazilian pre-salt reservoirs are renowned for their intricate pore networks and vuggy nature,posing significant challenges in modeling and simulating fluid flow within these carbonate reservoirs.Despite possessing excellent petrophysical properties,such as high porosity and permeability,these reservoirs typically exhibit a notably low recovery factor,sometimes falling below 10%.Previous research has indicated that various enhanced oil recovery(EOR)methods,such as water alternating gas(WAG),can substantially augment the recovery factor in pre-salt reservoirs,resulting in improvements of up to 20%.Nevertheless,the fluid flow mechanism within Brazilian carbonate reservoirs,characterized by complex pore geometry,remains unclear.Our study examines the behavior of fluid flow in a similar heterogeneous porous material,utilizing a plug sample obtained from a vugular segment of a Brazilian stromatolite outcrop,known to share analogies with certain pre-salt reservoirs.We conducted single-phase and multi-phase core flooding experiments,complemented by medical-CT scanning,to generate flow streamlines and evaluate the efficiency of water flooding.Subsequently,micro-CT scanning of the core sample was performed,and two cross-sections from horizontal and vertical plates were constructed.These cross-sections were then employed as geometries in a numerical simulator,enabling us to investigate the impact of pore geometry on fluid flow.Analysis of the pore-scale modeling and experimental data unveiled that the presence of dead-end pores and vugs results in a significant portion of the fluid remaining stagnant within these regions.Consequently,the injected fluid exhibits channeling-like behavior,leading to rapid breakthrough and low areal swept efficiency.Additionally,the numerical simulation results demonstrated that,irrespective of the size of the dead-end regions,the pressure variation within the dead-end vugs and pores is negligible.Despite the stromatolite's favorable petrophysical properties,including relatively high porosity and permeability,as well as the presence of interconnected large vugs,the recovery factor during water flooding remained low due to early breakthrough.These findings align with field data obtained from pre-salt reservoirs,providing an explanation for the observed low recovery factor during water flooding in such reservoirs.展开更多
基金supported by the National Natural Science Foundation of China(52304067,62273213)the Natural Science Foundation of Shandong Province of China(ZR2021QE073)+1 种基金the Natural Science Foundation of Shandong Province for Innovation and Development Joint Funds(ZR2022LZH001)the China Postdoctoral Science Foundation(2023M732111)。
文摘Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.
基金supported by National Key Research and Development Program of China ( 2018YFE011061)。
文摘The focus of this paper is on the measurement and calculation model of void fraction for the vertical upward co-current air-water slug flow in a circular tube of 15 mm inner diameter. High-speed photography and optical probes were utilized, with water superficial velocity ranging from 0.089 to 0.65 m·s^(-1)and gas superficial velocity ranging from 0.049 to 0.65 m·s^(-1). A new void fraction model based on the local parameters was proposed, disposing the slug flow as a combination of Taylor bubbles and liquid slugs. In the Taylor bubble region, correction factors of liquid film thickness Cδand nose shape CZ*were proposed to calculate aTB. In the liquid slug region, the radial void fraction distribution profiles were obtained to calculate aLS, by employing the image processing technique based on supervised machine learning. Results showed that the void fraction proportion in Taylor bubbles occupied crucial contribution to the overall void fraction. Multiple types of void fraction predictive correlations were assessed using the present data. The performance of the Schmidt model was optimal, while some models for slug flow performed not outstanding. Additionally, a predictive correlation was correlated between the central local void fraction and the cross-sectional averaged void fraction, as a straightforward form of the void fraction calculation model. The predictive correlation showed a good agreement with the present experimental data, as well as the data of Olerni et al., indicating that the new model was effective and applicable under the slug flow conditions.
基金supported by National Natural Science Foundation of China(42074142,51527805)。
文摘The hydrodynamic study of the liquid film around Taylor bubbles in slug flow has great significance for understanding parallel flow and interaction between Taylor bubbles.The prediction models for liquid film thickness mainly focus on stagnant flow,and some of them remain inaccurate performance.However,in the industrial process,the slug flow essentially is co-current flow.Therefore,in this paper,the liquid film thickness is studied by theoretical analysis and experimental methods under two conditions of stagnant and co-current flow.Firstly,under the condition of stagnant flow,the present work is based on Batchelor's theory,and modifies Batchelor's liquid film thickness model,which effectively improves its prediction accuracy.Under the condition of co-current flow,the prediction model of average liquid film thickness in slug flow is established by force and motion analysis.Taylor bubble length is introduced into the model as an important parameter.Dynamic experiments were carried out in the pipe with an inner diameter of 20 mm.The liquid film thickness,Taylor bubble velocity and length were measured by distributed ultrasonic sensor and intrusive cross-correlation conductivity sensor.Comparing the predicted value of the model with the measured results,the relative error is controlled within 10%.
文摘Several applications such as liquid-liquid extraction in micro-fluidic devices are concerned with the flow of two immiscible liquid phases. The commonly observed flow regimes in these systems are slug-flow and stratified flow. The latter regime in micro-channels has the inherent advantage that separation of the two liquids at the exit is efficient. Recently extraction in a stratified counter-current flow has been studied experimentally and it has been shown to be more efficient than co-current flow. An analytical as well as a numerical method to determine the steady-state solution of the corresponding convection-diffusion equation for the two flow-fields is presented. It is shown that the counter-current process is superior to the co-current process for the same set of parameters and operating conditions. A simplified model is proposed to analyse the process when diffusion in the transverse direction is not rate limiting. Different approaches to determining mass transfer coefficient are compared. The concept of log mean temperature difference used in design of heat exchangers is extended to describe mass transfer in the system.
基金Supported by the Ministry of Education Malaysia through RACE(RDU121308)and FRGS(RDU130136)
文摘A detached eddy simulation(DES) and a k-ε-based Reynolds-averaged Navier–Stokes(RANS) calculation on the co-current spray drying chamber is presented. The DES used here is based on the Spalart–Allmaras(SA) turbulence model, whereas the standard k-ε(SKE) was considered here for comparison purposes. Predictions of the mean axial velocity, temperature and humidity profile have been evaluated and compared with experimental measurements. The effects of the turbulence model on the predictions of the mean axial velocity, temperature and the humidity profile are most noticeable in the(highly anisotropic) spraying region. The findings suggest that DES provide a more accurate prediction(with error less than 5%) of the flow field in a spray drying chamber compared with RANS-based k-ε models. The DES simulation also confirmed the presence of anisotropic turbulent flow in the spray dryer from the analysis of the velocity component fluctuations and turbulent structure as illustrated by the Q-criterion.
基金This study was supported by the National Natural Science Foundation of China(U22B2075,52274056,51974356).
文摘A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.
基金supported by Shenzhen Science and Technology Program(Grant No.JCYJ20220818102012024)Hong Kong Research Grants Council(Grant Nos.T21–602/16-R and RGC R5037–18)。
文摘The flow field near a spur dike such as down flow and horseshoe vortex system(HVS)are susceptible to the topographic changes in the local scouring process,resulting in variation of the sediment transport with time.In this study,large eddy simulations with fixed-bed at different scouring stages were conducted to investigate the changes in flow field.The results imply that the bed deformation leads to an increase in flow rate per unit area,which represent the capability of sediment transportation by water,in the scour hole.Moreover,the intensity of turbulent kinetic energy and bimodal motion near the sand bed induced by the HVS were also varied.However,the peak moments between the two sediment transport mechanisms were different.Hence,understanding the complex feedback mechanism between topography and flow field is essential for the local scour problem.
基金financially supported by the Higher Education Commission of Pakistan (HEC) grant under National Research Program for Universities (NRPU) with No: (20-14681/NRPU/R&D/HEC/20212021)。
文摘The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study area which is extending along Karakorum Highway(KKH) from Besham to Chilas. Intense seismicity, deep gorges, steep terrain and extreme climatic events trigger multiple mountain hazards along the KKH, among which debris flow is recognized as the most destructive geohazard. This study aims to prepare a field-based debris flow inventory map at a regional scale along a 200 km stretch from Besham to Chilas. A total of 117 debris flows were identified in the field, and subsequently, a point-based debris-flow inventory and catchment delineation were performed through Arc GIS analysis. Regional scale debris flow susceptibility and propagation maps were prepared using Weighted Overlay Method(WOM) and Flow-R technique sequentially. Predisposing factors include slope, slope aspect, elevation, Topographic Roughness Index(TRI), Topographic Wetness Index(TWI), stream buffer, distance to faults, lithology rainfall, curvature, and collapsed material layer. The dataset was randomly divided into training data(75%) and validation data(25%). Results were validated through the Receiver Operator Characteristics(ROC) curve. Results show that Area Under the Curve(AUC) using WOM model is 79.2%. Flow-R propagation of debris flow shows that the 13.15%, 22.94%, and 63.91% areas are very high, high, and low susceptible to debris flow respectively. The propagation predicated by Flow-R validates the naturally occurring debris flow propagation as observed in the field surveys. The output of this research will provide valuable input to the decision makers for the site selection, designing of the prevention system, and for the protection of current infrastructure.
基金the support of EPIC-Energy Production Innovation Center,hosted by the University of Campinas(UNICAMP)sponsored by FAPESP-Sao Paulo Research Foundation(2017/15736e3 process).
文摘Brazilian pre-salt reservoirs are renowned for their intricate pore networks and vuggy nature,posing significant challenges in modeling and simulating fluid flow within these carbonate reservoirs.Despite possessing excellent petrophysical properties,such as high porosity and permeability,these reservoirs typically exhibit a notably low recovery factor,sometimes falling below 10%.Previous research has indicated that various enhanced oil recovery(EOR)methods,such as water alternating gas(WAG),can substantially augment the recovery factor in pre-salt reservoirs,resulting in improvements of up to 20%.Nevertheless,the fluid flow mechanism within Brazilian carbonate reservoirs,characterized by complex pore geometry,remains unclear.Our study examines the behavior of fluid flow in a similar heterogeneous porous material,utilizing a plug sample obtained from a vugular segment of a Brazilian stromatolite outcrop,known to share analogies with certain pre-salt reservoirs.We conducted single-phase and multi-phase core flooding experiments,complemented by medical-CT scanning,to generate flow streamlines and evaluate the efficiency of water flooding.Subsequently,micro-CT scanning of the core sample was performed,and two cross-sections from horizontal and vertical plates were constructed.These cross-sections were then employed as geometries in a numerical simulator,enabling us to investigate the impact of pore geometry on fluid flow.Analysis of the pore-scale modeling and experimental data unveiled that the presence of dead-end pores and vugs results in a significant portion of the fluid remaining stagnant within these regions.Consequently,the injected fluid exhibits channeling-like behavior,leading to rapid breakthrough and low areal swept efficiency.Additionally,the numerical simulation results demonstrated that,irrespective of the size of the dead-end regions,the pressure variation within the dead-end vugs and pores is negligible.Despite the stromatolite's favorable petrophysical properties,including relatively high porosity and permeability,as well as the presence of interconnected large vugs,the recovery factor during water flooding remained low due to early breakthrough.These findings align with field data obtained from pre-salt reservoirs,providing an explanation for the observed low recovery factor during water flooding in such reservoirs.
