In this paper, a third-order exponential time differencing scheme, named ETDRK3, was investigated for large time stepping in the computation of hypersonic non-equilibrium flow. The second-order Harten-TVD scheme was u...In this paper, a third-order exponential time differencing scheme, named ETDRK3, was investigated for large time stepping in the computation of hypersonic non-equilibrium flow. The second-order Harten-TVD scheme was used for the spatial discretization. The efficient implementation of the scheme with diagonalization of Jacobin matrix was established and carried out for the semi-cylindrical around flow. Current observations showed that the numerical results were in good agreement with those obtained by the classical explicit three-stage Runge-Kutta scheme (RK3) and implicit LU scheme. Efficiency assessments promised the effectiveness of the ETDRK3 scheme. The rationality of the application of this scheme was proved by its preferable accuracy and efficiency.展开更多
Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of d...Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of dry coal under gas adsorption equilibrium,gas flow and gas diffusion within wet coal under the generally non-equilibrium state are often ignored in the process of gas recovery.In this study,an improved apparent permeability model is proposed which accommodates the water and gas adsorption,stress dependence,water film thickness and gas flow regimes.In the process of modeling,the water adsorption is only affected by water content while the gas adsorption is time and water content dependent;based on poroelastic mechanics,the effective fracture aperture and effective pore radius are derived;and then the variation in water film thickness for different pore types under the effect of water content,stress and adsorption swelling are modeled;the flow regimes are considered based on Beskok’s model.Further,after validation with experimental data,the proposed model was applied to numerical simulations to investigate the evolution of permeability-related factors under the effect of different water contents.The gas flow in wet coal under the non-equilibrium state is explicitly revealed.展开更多
The high-temperature non-equilibrium effect is a novel and significant issue in the flows over a high Mach number(above Mach 8)air-breathing vehicle.Thus,this study attempts to investigate the high-temperature non-equ...The high-temperature non-equilibrium effect is a novel and significant issue in the flows over a high Mach number(above Mach 8)air-breathing vehicle.Thus,this study attempts to investigate the high-temperature non-equilibrium flows of a curved compression two-dimensional scramjet inlet at Mach 8 to 12 utilizing the two-dimensional non-equilibrium RANS calculations.Notably,the thermochemical non-equilibrium gas model can predict the actual high-temperature flows,and the numerical results of the other four thermochemical gas models are only used for comparative analysis.Firstly,the thermochemical non-equilibrium flow fields and work performance of the inlet at Mach 8 to 12 are analyzed.Then,the influences of high-temperature non-equilibrium effects on the starting characteristics of the inlet are investigated.The results reveal that a large separation bubble caused by the cowl shock/lower wall boundary layer interaction appears upstream of the shoulder,at Mach 8.The separation zone size is smaller,and its location is closer to the downstream area while the thermal process changes from frozen to non-equilibrium and then to equilibrium.With the increase of inflow Mach number,the thermochemical non-equilibrium effects in the whole inlet flow field gradually strengthen,so their influences on the overall work performance of the high Mach number inlet are more obvious.The vibrational relaxation or thermal non-equilibrium effects can yield more visible influences on the inlet performance than the chemical non-equilibrium reactions.The inlet in the thermochemical non-equilibrium flow can restart more easily than that in the thermochemical frozen flow.This work should provide a basis for the design and starting ability prediction of the high Mach number inlet in the wide operation range.展开更多
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.展开更多
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.展开更多
Measurement of bloodflow velocity is key to understanding physiology and pathology in vivo.While most measurements are performed at the middle of the blood vessel,little research has been done on characterizing the in...Measurement of bloodflow velocity is key to understanding physiology and pathology in vivo.While most measurements are performed at the middle of the blood vessel,little research has been done on characterizing the instantaneous bloodflow velocity distribution.This is mainly due to the lack of measurement technology with high spatial and temporal resolution.Here,we tackle this problem with our recently developed dual-wavelength line-scan third-harmonic generation(THG)imaging technology.Simultaneous acquisition of dual-wavelength THG line-scanning signals enables measurement of bloodflow velocities at two radially symmetric positions in both venules and arterioles in mouse brain in vivo.Our results clearly show that the instantaneous bloodflow velocity is not symmetric under general conditions.展开更多
The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a...The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.展开更多
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.展开更多
Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy ...Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.展开更多
文摘In this paper, a third-order exponential time differencing scheme, named ETDRK3, was investigated for large time stepping in the computation of hypersonic non-equilibrium flow. The second-order Harten-TVD scheme was used for the spatial discretization. The efficient implementation of the scheme with diagonalization of Jacobin matrix was established and carried out for the semi-cylindrical around flow. Current observations showed that the numerical results were in good agreement with those obtained by the classical explicit three-stage Runge-Kutta scheme (RK3) and implicit LU scheme. Efficiency assessments promised the effectiveness of the ETDRK3 scheme. The rationality of the application of this scheme was proved by its preferable accuracy and efficiency.
基金the National Natural Science Foundation of China(No.52079077)the Natural Science Foundation of Shandong Province(No.ZR2021QE069).
文摘Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of dry coal under gas adsorption equilibrium,gas flow and gas diffusion within wet coal under the generally non-equilibrium state are often ignored in the process of gas recovery.In this study,an improved apparent permeability model is proposed which accommodates the water and gas adsorption,stress dependence,water film thickness and gas flow regimes.In the process of modeling,the water adsorption is only affected by water content while the gas adsorption is time and water content dependent;based on poroelastic mechanics,the effective fracture aperture and effective pore radius are derived;and then the variation in water film thickness for different pore types under the effect of water content,stress and adsorption swelling are modeled;the flow regimes are considered based on Beskok’s model.Further,after validation with experimental data,the proposed model was applied to numerical simulations to investigate the evolution of permeability-related factors under the effect of different water contents.The gas flow in wet coal under the non-equilibrium state is explicitly revealed.
基金co-supported by the China Scholarship Council(No.202206840048)the Training Fund for Excellent Doctoral Candidates of Nanjing University of Science and Technology,China+4 种基金the Opening Foundation of State Key Laboratory of High Temperature Gas Dynamics,Institute of Mechanics,China(No.2021KF07)the National Key Laboratory Fund,China(No.2022-JCJQ-LB-020-01)the Foundation of Key Laboratory of Hypersonic Aerodynamic Force and Heat Technology,AVIC Aerodynamics Research Institute,China(No.XFX20220104)the China Postdoctoral Science Foundation(No.BX20200070)the Fundamental Research Foundation of the Central Universities,China(No.2022CDJXY012)。
文摘The high-temperature non-equilibrium effect is a novel and significant issue in the flows over a high Mach number(above Mach 8)air-breathing vehicle.Thus,this study attempts to investigate the high-temperature non-equilibrium flows of a curved compression two-dimensional scramjet inlet at Mach 8 to 12 utilizing the two-dimensional non-equilibrium RANS calculations.Notably,the thermochemical non-equilibrium gas model can predict the actual high-temperature flows,and the numerical results of the other four thermochemical gas models are only used for comparative analysis.Firstly,the thermochemical non-equilibrium flow fields and work performance of the inlet at Mach 8 to 12 are analyzed.Then,the influences of high-temperature non-equilibrium effects on the starting characteristics of the inlet are investigated.The results reveal that a large separation bubble caused by the cowl shock/lower wall boundary layer interaction appears upstream of the shoulder,at Mach 8.The separation zone size is smaller,and its location is closer to the downstream area while the thermal process changes from frozen to non-equilibrium and then to equilibrium.With the increase of inflow Mach number,the thermochemical non-equilibrium effects in the whole inlet flow field gradually strengthen,so their influences on the overall work performance of the high Mach number inlet are more obvious.The vibrational relaxation or thermal non-equilibrium effects can yield more visible influences on the inlet performance than the chemical non-equilibrium reactions.The inlet in the thermochemical non-equilibrium flow can restart more easily than that in the thermochemical frozen flow.This work should provide a basis for the design and starting ability prediction of the high Mach number inlet in the wide operation range.
基金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.
基金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.
基金funded by the National Natural Science Foundation of China(Grant/Award Numbers 62075135 and 61975126)the Science and Technology Innovation Commission of Shenzhen(Grant/Award Numbers JCYJ20190808174819083 and JCYJ20190808175201640)Shenzhen Science and Technology Planning Project(ZDSYS 20210623092006020).
文摘Measurement of bloodflow velocity is key to understanding physiology and pathology in vivo.While most measurements are performed at the middle of the blood vessel,little research has been done on characterizing the instantaneous bloodflow velocity distribution.This is mainly due to the lack of measurement technology with high spatial and temporal resolution.Here,we tackle this problem with our recently developed dual-wavelength line-scan third-harmonic generation(THG)imaging technology.Simultaneous acquisition of dual-wavelength THG line-scanning signals enables measurement of bloodflow velocities at two radially symmetric positions in both venules and arterioles in mouse brain in vivo.Our results clearly show that the instantaneous bloodflow velocity is not symmetric under general conditions.
基金Funding support from Heilongjiang"Open Competition"project(Grant No.DQYT2022-JS-758)is greatly acknowledgedFinancial support from the National Natural Science Foundation of China(Grant Nos.52304025 and 52174025)is acknowledged+1 种基金supports from Northeast Petroleum University and Guangdong Basic and Applied Basic Research Foundationsupport from the Heilongjiang Touyan Innovation Team Program.
文摘The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.
基金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.
基金financially supported by the Chang Jiang Scholar and Innovation Team Development Plan of China (IRT_15R29)the Basic Research Innovation Group Project of Gansu Province, China (21JR7RA347)the Natural Science Foundation of Gansu Province, China (20JR10RA231)。
文摘Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.
