Low-speed flow experiments in which ultra-fine copper tubes are used to simulate micro-fractures in carbonate strata are conducted to analyze the variations of gas flow state in fractures of different fracture heights...Low-speed flow experiments in which ultra-fine copper tubes are used to simulate micro-fractures in carbonate strata are conducted to analyze the variations of gas flow state in fractures of different fracture heights,determine flow state transition limit and transition interval,and establish the calculation method of flow state transition limit.The results show that the ideal Hagen-Poiseuille flow is the main form of gas flow in large fractures.Due to the decrease of fracture height,the gas flow in the fracture changes from Hagen-Poiseuille flow with ideal smooth seam surface to non-Hagen-Poiseuille flow,and the critical point of the transition is the boundary of flow state transition.After the fracture height continues to decrease to a certain extent below the boundary of the flow state transition fracture height,the form of gas flow gradually changes to the ideal Darcy flow,thus the transition interval of the gas flow state in the closing process of fracture can be determined.Based on the three-dimensional microconvex body scanning of the fracture surface,the material properties of fracture and properties of fluid in the fracture,a method for calculating the boundary of flow state transition is established.The experimental test and theoretical calculation show that the limit of the fracture height for the transition from pipe flow to Darcy flow is about twice the sum of the maximum height of the microconvex bodies on the upper and lower sides of the fracture.展开更多
We established a model for the chemical reaction kinetics of biomass pyrolysis via the hightemperature thermal cracking of liquid products. We divided the condensable volatiles into two groups, based on the characteri...We established a model for the chemical reaction kinetics of biomass pyrolysis via the hightemperature thermal cracking of liquid products. We divided the condensable volatiles into two groups, based on the characteristics of the liquid prdoducts., tar and biomass oil. The effects of temperature, residence time, particle size, velocity, pressure, and other parameters on biomass pyrolysis and high-temperature tar cracking were investigated numerically, and the results were compared with experimental data. The simulation results showed a large endothermic pyrolysis reaction effect on temperature and the reaction process. The pyrolysis reaction zone had a constant temperature period in several layers near the center of large biomass particles. A purely physical heating process was observed before and after this period, according to the temperature index curve.展开更多
The authors consider a phase field model for Darcy flows with discontinuous data in porous media;specifically,they adopt the Hele-Shaw-Cahn-Hillard equations of[Lee,Lowengrub,Goodman,Physics of Fluids,2002]to model fl...The authors consider a phase field model for Darcy flows with discontinuous data in porous media;specifically,they adopt the Hele-Shaw-Cahn-Hillard equations of[Lee,Lowengrub,Goodman,Physics of Fluids,2002]to model flows in the Hele-Shaw cell through a phase field formulation which incorporates discontinuities of physical data,namely density and viscosity,across interfaces.For the spatial approximation of the problem,the authors use NURBS—based isogeometric analysis in the framework of the Galerkin method,a computational framework which is particularly advantageous for the solution of high order partial differential equations and phase field problems which exhibit sharp but smooth interfaces.In this paper,the authors verify through numerical tests the sharp interface limit of the phase field model which in fact leads to an internal discontinuity interface problem;finally,they show the efficiency of isogeometric analysis for the numerical approximation of the model by solving a benchmark problem,the so-called"rising bubble"problem.展开更多
This paper considers the problem of hydrodynamics and thermal boundary layers of Darcy flow over horizontal surface embedded in a porous medium. The solutions of such problems for the cases of uniform surface temperat...This paper considers the problem of hydrodynamics and thermal boundary layers of Darcy flow over horizontal surface embedded in a porous medium. The solutions of such problems for the cases of uniform surface temperature and variable surface temperature have been studied and analysed in many papers. This paper, however, attempts to find similarity solutions for the Darcy flow problem with a convective boundary condition at the plate surface. It is found that the solution is possible when the heat transfer coefficient is proportional to x<sup>–2/3</sup>. The numerical solutions thus obtained are analyzed for a range of values of the parameter characterizing the hot fluid convection process. Analytical expressions are provided for local surface heat flux and total surface heat transfer rate while the flow variables are discussed graphically.展开更多
In this paper, we develop an a-priori error analysis of a new unified mixed finite element method for the coupling of fluid flow with porous media flow in R<sup><em>N</em></sup>, <em>N<...In this paper, we develop an a-priori error analysis of a new unified mixed finite element method for the coupling of fluid flow with porous media flow in R<sup><em>N</em></sup>, <em>N</em> ∈ {2,3}, on isotropic meshes. Flows are governed by the Stokes and Darcy equations, respectively, and the corresponding transmission conditions are given by mass conservation, balance of normal forces, and the Beavers-Joseph-Saffman law. The approach utilizes a modification of the Darcy problem which allows us to apply a variant nonconforming Crouzeix-Raviart finite element to the whole coupled Stokes-Darcy problem. The well-posedness of the finite element scheme and its convergence analysis are derived. Finally, the numerical experiments are presented, which confirm the excellent stability and accuracy of our method.展开更多
The fabrication of superhydrophobic surfaces and the studies on water flow characteristics therein are of great significance to many industrial areas as well as to science and technology development. Experiments were ...The fabrication of superhydrophobic surfaces and the studies on water flow characteristics therein are of great significance to many industrial areas as well as to science and technology development. Experiments were car- ried out to investigate slip characteristics of water flowing in circular superhydrophobic microtubes within lam- inar flow region. The superhydrophobic microtubes of stainless steel were fabricated with chemical etching- fluorination treatment. An experimental setup was designed to measure the pressure drop as function of water flow rate. For comparison, superhydrophilic tubes were also tested. Poiseuille number Po was found to be smaller for the superhydrophobic microtubes than that for superhydrophilic ones. The pressure drop reduc- tion ranges from 8% to 31%. It decreases with increasing Reynolds number when Re 〈 900, owing to the transition from Cassie state to Wenze] state. However, it is almost unchanged with further increasing Re after Re 〉 900. The slip length in superhydrophobic microtubes also exhibits a Reynolds number dependence similarly to the pressure drop reduction. The relation between slip length and Darcy friction factor is theoretically analyzed with consideration of surface roughness effect, which was testified with the experimental results.展开更多
Permeability is an important index in reservoir evaluation,oil and gas accumulation control,and production effi ciency.At present,permeability can be obtained through several methods.However,these methods are not suit...Permeability is an important index in reservoir evaluation,oil and gas accumulation control,and production effi ciency.At present,permeability can be obtained through several methods.However,these methods are not suitable for tight sandstone in general because the pore type in tight sandstone is mainly secondary pores and has the characteristics of low porosity and permeability,high capillary pressure,and high irreducible water saturation.Mud invasion depth is closely related to permeability during drilling.In general,the greater the permeability,the shallower the mud invasion depth,and the smaller the permeability,the deeper the mud invasion depth.Therefore,this paper builds a model to predict the permeability of tight sandstone using mud invasion depth.The model is based on the improvement of the Darcy flow equation to obtain permeability using mud invasion depth inversion of array induction logging.The influence of various permeability factors on the model is analyzed by numerical simulation.The model is used to predict the permeability of tight sandstone in the south of the Ordos Basin.The predicted permeability is highly consistent with the core analysis permeability,which verifi es the reliability of the method.展开更多
Multiple fractured horizontal well(MFHW) is widely applied in the development of shale gas. To investigate the gas flow characteristics in shale, based on a new dual mechanism triple continuum model, an analytical sol...Multiple fractured horizontal well(MFHW) is widely applied in the development of shale gas. To investigate the gas flow characteristics in shale, based on a new dual mechanism triple continuum model, an analytical solution for MFHW surrounded by stimulated reservoir volume(SRV) was presented. Pressure and pressure derivative curves were used to identify the characteristics of flow regimes in shale. Blasingame type curves were established to evaluate the effects of sensitive parameters on rate decline curves, which indicates that the whole flow regimes could be divided into transient flow, feeding flow, and pseudo steady state flow. In feeding flow regime, the production of gas well is gradually fed by adsorbed gases in sub matrix, and free gases in matrix. The proportion of different gas sources to well production is determined by such parameters as storability ratios of triple continuum, transmissibility coefficients controlled by dual flow mechanism and fracture conductivity.展开更多
A nonconforming rectangular finite element method is proposed to solve a fluid structure interaction problem characterized by the Darcy-Stokes-Brinkman Equation with discontinuous coefficients across the interface of ...A nonconforming rectangular finite element method is proposed to solve a fluid structure interaction problem characterized by the Darcy-Stokes-Brinkman Equation with discontinuous coefficients across the interface of different structures.A uniformly stable mixed finite element together with Nitsche-type matching condi-tions that automatically adapt to the coupling of different sub-problem combinations are utilized in the discrete algorithm.Compared with other finite element methods in the literature,the new method has some distinguished advantages and features.The Boland-Nicolaides trick is used in proving the inf-sup condition for the multi-domain discrete problem.Optimal error estimates are derived for the coupled prob-lem by analyzing the approximation errors and the consistency errors.Numerical examples are also provided to confirm the theoretical results.展开更多
基金Supported by the Science and Technology Cooperation Project of the CNPC-SWPU Innovation Alliance(2020CX010401)。
文摘Low-speed flow experiments in which ultra-fine copper tubes are used to simulate micro-fractures in carbonate strata are conducted to analyze the variations of gas flow state in fractures of different fracture heights,determine flow state transition limit and transition interval,and establish the calculation method of flow state transition limit.The results show that the ideal Hagen-Poiseuille flow is the main form of gas flow in large fractures.Due to the decrease of fracture height,the gas flow in the fracture changes from Hagen-Poiseuille flow with ideal smooth seam surface to non-Hagen-Poiseuille flow,and the critical point of the transition is the boundary of flow state transition.After the fracture height continues to decrease to a certain extent below the boundary of the flow state transition fracture height,the form of gas flow gradually changes to the ideal Darcy flow,thus the transition interval of the gas flow state in the closing process of fracture can be determined.Based on the three-dimensional microconvex body scanning of the fracture surface,the material properties of fracture and properties of fluid in the fracture,a method for calculating the boundary of flow state transition is established.The experimental test and theoretical calculation show that the limit of the fracture height for the transition from pipe flow to Darcy flow is about twice the sum of the maximum height of the microconvex bodies on the upper and lower sides of the fracture.
文摘We established a model for the chemical reaction kinetics of biomass pyrolysis via the hightemperature thermal cracking of liquid products. We divided the condensable volatiles into two groups, based on the characteristics of the liquid prdoducts., tar and biomass oil. The effects of temperature, residence time, particle size, velocity, pressure, and other parameters on biomass pyrolysis and high-temperature tar cracking were investigated numerically, and the results were compared with experimental data. The simulation results showed a large endothermic pyrolysis reaction effect on temperature and the reaction process. The pyrolysis reaction zone had a constant temperature period in several layers near the center of large biomass particles. A purely physical heating process was observed before and after this period, according to the temperature index curve.
基金partially funded by the INdAM-GNCS Project 2017“Modellistica numerica di fenomeni idro/geomeccanici per la simulazione di eventi sismici”
文摘The authors consider a phase field model for Darcy flows with discontinuous data in porous media;specifically,they adopt the Hele-Shaw-Cahn-Hillard equations of[Lee,Lowengrub,Goodman,Physics of Fluids,2002]to model flows in the Hele-Shaw cell through a phase field formulation which incorporates discontinuities of physical data,namely density and viscosity,across interfaces.For the spatial approximation of the problem,the authors use NURBS—based isogeometric analysis in the framework of the Galerkin method,a computational framework which is particularly advantageous for the solution of high order partial differential equations and phase field problems which exhibit sharp but smooth interfaces.In this paper,the authors verify through numerical tests the sharp interface limit of the phase field model which in fact leads to an internal discontinuity interface problem;finally,they show the efficiency of isogeometric analysis for the numerical approximation of the model by solving a benchmark problem,the so-called"rising bubble"problem.
文摘This paper considers the problem of hydrodynamics and thermal boundary layers of Darcy flow over horizontal surface embedded in a porous medium. The solutions of such problems for the cases of uniform surface temperature and variable surface temperature have been studied and analysed in many papers. This paper, however, attempts to find similarity solutions for the Darcy flow problem with a convective boundary condition at the plate surface. It is found that the solution is possible when the heat transfer coefficient is proportional to x<sup>–2/3</sup>. The numerical solutions thus obtained are analyzed for a range of values of the parameter characterizing the hot fluid convection process. Analytical expressions are provided for local surface heat flux and total surface heat transfer rate while the flow variables are discussed graphically.
文摘In this paper, we develop an a-priori error analysis of a new unified mixed finite element method for the coupling of fluid flow with porous media flow in R<sup><em>N</em></sup>, <em>N</em> ∈ {2,3}, on isotropic meshes. Flows are governed by the Stokes and Darcy equations, respectively, and the corresponding transmission conditions are given by mass conservation, balance of normal forces, and the Beavers-Joseph-Saffman law. The approach utilizes a modification of the Darcy problem which allows us to apply a variant nonconforming Crouzeix-Raviart finite element to the whole coupled Stokes-Darcy problem. The well-posedness of the finite element scheme and its convergence analysis are derived. Finally, the numerical experiments are presented, which confirm the excellent stability and accuracy of our method.
基金the National Natural Science Foundation of China(20476014,51376030)
文摘The fabrication of superhydrophobic surfaces and the studies on water flow characteristics therein are of great significance to many industrial areas as well as to science and technology development. Experiments were car- ried out to investigate slip characteristics of water flowing in circular superhydrophobic microtubes within lam- inar flow region. The superhydrophobic microtubes of stainless steel were fabricated with chemical etching- fluorination treatment. An experimental setup was designed to measure the pressure drop as function of water flow rate. For comparison, superhydrophilic tubes were also tested. Poiseuille number Po was found to be smaller for the superhydrophobic microtubes than that for superhydrophilic ones. The pressure drop reduc- tion ranges from 8% to 31%. It decreases with increasing Reynolds number when Re 〈 900, owing to the transition from Cassie state to Wenze] state. However, it is almost unchanged with further increasing Re after Re 〉 900. The slip length in superhydrophobic microtubes also exhibits a Reynolds number dependence similarly to the pressure drop reduction. The relation between slip length and Darcy friction factor is theoretically analyzed with consideration of surface roughness effect, which was testified with the experimental results.
基金supported by the National Natural Science Foundation of China project(No.41504103 and No.41804097).
文摘Permeability is an important index in reservoir evaluation,oil and gas accumulation control,and production effi ciency.At present,permeability can be obtained through several methods.However,these methods are not suitable for tight sandstone in general because the pore type in tight sandstone is mainly secondary pores and has the characteristics of low porosity and permeability,high capillary pressure,and high irreducible water saturation.Mud invasion depth is closely related to permeability during drilling.In general,the greater the permeability,the shallower the mud invasion depth,and the smaller the permeability,the deeper the mud invasion depth.Therefore,this paper builds a model to predict the permeability of tight sandstone using mud invasion depth.The model is based on the improvement of the Darcy flow equation to obtain permeability using mud invasion depth inversion of array induction logging.The influence of various permeability factors on the model is analyzed by numerical simulation.The model is used to predict the permeability of tight sandstone in the south of the Ordos Basin.The predicted permeability is highly consistent with the core analysis permeability,which verifi es the reliability of the method.
基金Project(2011ZX05015)supported by Important National Science and Technology Specific Projects of the "Twelfth Five-years" Plan Period,China
文摘Multiple fractured horizontal well(MFHW) is widely applied in the development of shale gas. To investigate the gas flow characteristics in shale, based on a new dual mechanism triple continuum model, an analytical solution for MFHW surrounded by stimulated reservoir volume(SRV) was presented. Pressure and pressure derivative curves were used to identify the characteristics of flow regimes in shale. Blasingame type curves were established to evaluate the effects of sensitive parameters on rate decline curves, which indicates that the whole flow regimes could be divided into transient flow, feeding flow, and pseudo steady state flow. In feeding flow regime, the production of gas well is gradually fed by adsorbed gases in sub matrix, and free gases in matrix. The proportion of different gas sources to well production is determined by such parameters as storability ratios of triple continuum, transmissibility coefficients controlled by dual flow mechanism and fracture conductivity.
基金supported by the US NSF grant DMS-1522768,CNFS grants.11371199,11471166.
文摘A nonconforming rectangular finite element method is proposed to solve a fluid structure interaction problem characterized by the Darcy-Stokes-Brinkman Equation with discontinuous coefficients across the interface of different structures.A uniformly stable mixed finite element together with Nitsche-type matching condi-tions that automatically adapt to the coupling of different sub-problem combinations are utilized in the discrete algorithm.Compared with other finite element methods in the literature,the new method has some distinguished advantages and features.The Boland-Nicolaides trick is used in proving the inf-sup condition for the multi-domain discrete problem.Optimal error estimates are derived for the coupled prob-lem by analyzing the approximation errors and the consistency errors.Numerical examples are also provided to confirm the theoretical results.