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.展开更多
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.展开更多
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.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
Non-equilibrium solidification structures of Cu55Ni45 and Cu55Ni43Co2 alloys were prepared by the molten glass purification cycle superheating method.The variation of the recalescence phenomenon with the degree of und...Non-equilibrium solidification structures of Cu55Ni45 and Cu55Ni43Co2 alloys were prepared by the molten glass purification cycle superheating method.The variation of the recalescence phenomenon with the degree of undercooling in the rapid solidification process was investigated using an infrared thermometer.The addition of the Co element affected the evolution of the recalescence phenomenon in Cu-Ni alloys.The images of the solid-liquid interface migration during the rapid solidification of supercooled melts were captured by using a high-speed camera.The solidification rate of Cu-Ni alloys,with the addition of Co elements,was explored.Finally,the grain refinement structure with low supercooling was characterised using electron backscatter diffraction(EBSD).The effect of Co on the microstructural evolution during nonequilibrium solidification of Cu-Ni alloys under conditions of small supercooling is investigated by comparing the microstructures of Cu55Ni45 and Cu55Ni43Co2 alloys.The experimental results show that the addition of a small amount of Co weakens the recalescence behaviour of the Cu55Ni45 alloy and significantly reduces the thermal strain in the rapid solidification phase.In the rapid solidification phase,the thermal strain is greatly reduced,and there is a significant increase in the characteristic undercooling degree.Furthermore,the addition of Co and the reduction of Cu not only result in a lower solidification rate of the alloy,but also contribute to the homogenisation of the grain size.展开更多
The effect of the gradient content of Co element on the solidification process of Cu-based alloy under deep under cooling conditions was explored.The non-equilibrium solidification structure of the under cooled alloy ...The effect of the gradient content of Co element on the solidification process of Cu-based alloy under deep under cooling conditions was explored.The non-equilibrium solidification structure of the under cooled alloy samples were analyzed.It is found that the rapidly solidified alloy has undergone twice grain refinement during the undercooling process.Characterization and significance of the maximum undercooling refinement structure of Cu60Ni35Co5 at T=253 K were analyzed.High-density defects were observed,such as dislocations,stacking faults networks,and twinning structures.The standard FCC diffraction pattern represents that it is still a single-phase structure.Based on the metallographic diagram,EBSD and TEM data analysis,it is illustrated that the occurrence of grain refinement under high undercooling is due to stress induced recrystallization.In addition,the laser cladding technology is used to coat Co-based alloy(Stellite12) coating on 304 stainless steel substrate;the microstructure of the coating cross-section was analyzed.It was found that the microstructure of the cross-section is presented as columnar crystals,planar crystals,and disordered growth direction,so that the coating has better hardness and wear resistance.By electrochemical corrosion of the substrate and coating,it can be seen that the Co and Cr elements present in the coating are more likely to form a dense passivation film,which improved the corrosion resistance of the coating.展开更多
Long-term permeability experiments have indicated that sorption-induced swelling can switch from internal to bulk depending on the evolutive sorption status.However,this sorption swelling switch mechanism has not been...Long-term permeability experiments have indicated that sorption-induced swelling can switch from internal to bulk depending on the evolutive sorption status.However,this sorption swelling switch mechanism has not been considered in current analytical permeability models.This study introduces a normalized sorption non-equilibrium index(SNEI)to characterize the sorption status,quantify the dynamical variations of matrix swelling accumulation and internal swelling partition,and formulate the sorption swelling switch model.The incorporation of this index into the extended total effective stress concept leads to an analytical transient coal permeability model.Model results show that the sorption swelling switch itself results in the permeability switch under stress-constrained conditions,while the confined bulk swelling suppresses the permeability recovery to the continuous reduction under displacement-constrained conditions.Model verifications show that current experimental observations correspond to the early stages of the transient process,and they could be extended to the whole process with these models.This study demonstrates the importance of the sorption swelling switch in determining permeability evolution using simple boundary conditions.It provides new insights into experimentally revealing the sorption swelling switch in the future,and underscores the requirement of a rigorous model for complex coupled processes in large-scale coal seams.展开更多
Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks th...Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.展开更多
文摘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.
基金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.
基金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.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
文摘Non-equilibrium solidification structures of Cu55Ni45 and Cu55Ni43Co2 alloys were prepared by the molten glass purification cycle superheating method.The variation of the recalescence phenomenon with the degree of undercooling in the rapid solidification process was investigated using an infrared thermometer.The addition of the Co element affected the evolution of the recalescence phenomenon in Cu-Ni alloys.The images of the solid-liquid interface migration during the rapid solidification of supercooled melts were captured by using a high-speed camera.The solidification rate of Cu-Ni alloys,with the addition of Co elements,was explored.Finally,the grain refinement structure with low supercooling was characterised using electron backscatter diffraction(EBSD).The effect of Co on the microstructural evolution during nonequilibrium solidification of Cu-Ni alloys under conditions of small supercooling is investigated by comparing the microstructures of Cu55Ni45 and Cu55Ni43Co2 alloys.The experimental results show that the addition of a small amount of Co weakens the recalescence behaviour of the Cu55Ni45 alloy and significantly reduces the thermal strain in the rapid solidification phase.In the rapid solidification phase,the thermal strain is greatly reduced,and there is a significant increase in the characteristic undercooling degree.Furthermore,the addition of Co and the reduction of Cu not only result in a lower solidification rate of the alloy,but also contribute to the homogenisation of the grain size.
基金Funded by the Basic Research Projects in Shanxi Province(No.202103021224183)。
文摘The effect of the gradient content of Co element on the solidification process of Cu-based alloy under deep under cooling conditions was explored.The non-equilibrium solidification structure of the under cooled alloy samples were analyzed.It is found that the rapidly solidified alloy has undergone twice grain refinement during the undercooling process.Characterization and significance of the maximum undercooling refinement structure of Cu60Ni35Co5 at T=253 K were analyzed.High-density defects were observed,such as dislocations,stacking faults networks,and twinning structures.The standard FCC diffraction pattern represents that it is still a single-phase structure.Based on the metallographic diagram,EBSD and TEM data analysis,it is illustrated that the occurrence of grain refinement under high undercooling is due to stress induced recrystallization.In addition,the laser cladding technology is used to coat Co-based alloy(Stellite12) coating on 304 stainless steel substrate;the microstructure of the coating cross-section was analyzed.It was found that the microstructure of the cross-section is presented as columnar crystals,planar crystals,and disordered growth direction,so that the coating has better hardness and wear resistance.By electrochemical corrosion of the substrate and coating,it can be seen that the Co and Cr elements present in the coating are more likely to form a dense passivation film,which improved the corrosion resistance of the coating.
基金supported by the Australian Research Council(Grant No.DP200101293)the National Natural Science Foundation of China(Grant No.42202286)the Zhejiang Collaborative Innovation Center for Prevention and Control of Mountain Geological Hazards(Grant No.PCMGH-2017-Z-02).
文摘Long-term permeability experiments have indicated that sorption-induced swelling can switch from internal to bulk depending on the evolutive sorption status.However,this sorption swelling switch mechanism has not been considered in current analytical permeability models.This study introduces a normalized sorption non-equilibrium index(SNEI)to characterize the sorption status,quantify the dynamical variations of matrix swelling accumulation and internal swelling partition,and formulate the sorption swelling switch model.The incorporation of this index into the extended total effective stress concept leads to an analytical transient coal permeability model.Model results show that the sorption swelling switch itself results in the permeability switch under stress-constrained conditions,while the confined bulk swelling suppresses the permeability recovery to the continuous reduction under displacement-constrained conditions.Model verifications show that current experimental observations correspond to the early stages of the transient process,and they could be extended to the whole process with these models.This study demonstrates the importance of the sorption swelling switch in determining permeability evolution using simple boundary conditions.It provides new insights into experimentally revealing the sorption swelling switch in the future,and underscores the requirement of a rigorous model for complex coupled processes in large-scale coal seams.
基金supported by National Natural Science Foundation(52204050)Sichuan Science and Technology Program(2021ZHCG0013,22ZDYF3009)。
文摘Bottom water coning is the main reason to reduce the recovery of horizontal bottom water reservoir. By water coning, we mean the oil-water interface changes from a horizontal state to a mound-shaped cone and breaks through to the wellbore. Autonomous inflow control device(AICD) is an important instrument maintain normal production after bottom water coning, however, the resistance increasing ability of the swirl type AICD is insufficient at present, which seriously affects the water control effect. Aiming this problem, this paper designs a multi-stage resistance-increasing and composite type AICD. The separation mechanism of oil-water two phases in this structure, the resistance form of oil-water single phase and the resistance-increasing principle of water phase are analyzed. Establishing the dual-phase multi-stage separation and resistance-increasing model, and verified by measuring the throttling pressure drop and oil-water volume fraction of the AICD, it is found that the composite type AICD has the effect of ICD and AICD at the same time, which can balance the production rate of each well section at the initial stage of production, delay the occurrence of bottom water coning. In the middle and later stages of production, water-blocking can be effectively increased to achieve water control and stable production.After structural sensitivity analysis, the influence law of various structural parameters on the water control performance of composite AICD was obtained. The simulation calculation results show that,compared with the existing swirl type AICD, composite AICD has higher sensitivity to moisture content,the water phase throttling pressure drop is increased by 4.5 times on average. The composite AICD is suitable for the entire stage of horizontal well production.