Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 me...Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.展开更多
Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing addit...Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the ...The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the two-phase flow pumps, i.e., low overall efficiency and severe abrasion. In this study, the three-dimensional, steady, incompressible, and turbulent solid-liquid two-phase flows in a low-specific-speed centrifugal pump are numerically simulated and analyzed by using a computational fluid dynamics （CFD） code based on the mixture model of the two-phase flow and the RNG k-~ two-equation turbulence model, in which the influences of rotation and curvature are fully taken into account. The coupling between impeller and volute is implemented by means of the frozen rotor method. The simulation results predicted indicate that the solid phase properties in two-phase flow, especially the concentration, the particle diameter and the density, have strong effects on the hydraulic performance of the pump. Both the pump head and the efficiency are reduced with increasing particle diameter or concentration. However, the effect of particle density on the performance is relatively minor. An obvious jet-wake flow structure is presented near the volute tongue and becomes more remarkable with increasing solid phase concentration. The suction side of the blade is subject to much more severe abrasion than the pressure side. The obtained results preliminarily reveal the characteristics of solid-liquid two-phase flow in the centrifugal pump, and are helpful for improvement and empirical correction in the hydraulic design of centrifugal pumps.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to T...A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to TECPLOT as the post-processor, the simulation results isvisualized. The results show the main flow characteristics: There exists backflow and aberrantvelocities at inlet area and a relative velocity slip between two phases; A jet-wake flow pattern isdiscerned around the shroud-suction side area; The relative velocity vector of solid phase iscloser to the pressure surface than that of liquid phase and the trend is more obvious with theincrease of diameter; The kinetic energy of turbulence k and the dissipation rate e reach theirpeaks at the corner of pressure and suction surface. The simulation results show a good agreementwith the experimental flow features in the impeller channel, which prove the turbulent model used isvalid and provide a theoretical design basis to non-clogging pumps.展开更多
Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid pa...Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.展开更多
Taking an oxygen enriched side-blown furnace as the prototype,a hydraulic model was established according to the similarity principle.The influence of three factors on the gas-liquid two-phase flow was analyzed,i.e.th...Taking an oxygen enriched side-blown furnace as the prototype,a hydraulic model was established according to the similarity principle.The influence of three factors on the gas-liquid two-phase flow was analyzed,i.e.the airflow speed,the submerged depth and the downward angle of the nozzle.A numerical simulation of the hydraulic model was carried out trying to find the suitable turbulence model which can describe the side-blown two-phase flow correctly by comparing the simulation results with the experimental data.The experiment shows that the airflow speed has a great influence on the flow of the water.The submerged depth of the nozzle has a relatively smaller influence on the penetration depth and the surface fluctuation height in the liquid phase.When the nozzle is at a downward angle of 15°,the penetration depth and the surface fluctuation height are reduced.It is concluded that the numerical results with the realizable k-εturbulence model are the closest to the experiment for the penetration depth,the surface fluctuation height and the bubble scale.展开更多
Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model a...Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.展开更多
Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consum...Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.展开更多
According to the Lagrange discrete phase model of multiphase flow, mathematical model of the abrasive water-jet （AWJ） nozzle based on the multi-phase movement was introduced, then the boundary conditions was determi...According to the Lagrange discrete phase model of multiphase flow, mathematical model of the abrasive water-jet （AWJ） nozzle based on the multi-phase movement was introduced, then the boundary conditions was determined and the liquid-solid turbulence which is isothermal, can not be compressed and steadystate in the cone-cylinder nozzles of the export of the pre-mixed AWJ was simulated applying the software FLUENT. The results showed that： the axial velocity and dynamic pressure of the continuous phase in the nozzle were axial symmetry notable, and at the axis had a extreme point; abrasive accelerated at two points, in front contractive segment, the rate increased rapidly, and in the back straight one, the speed accelerated slowly. The length of the cylinder is 100 mm, the diameter of the nozzle is 8 mm, and the angle of the cone is 15°. There is a extreme point of the rate at the point 10 mm in the established model. The results of simulation laid the foundation for optimizing the nozzle structure, improving efficiency and developing the nozzle.展开更多
The flow focusing nozzle is a new type of nozzle that performs effective atomization of the discrete phase by means of high-speed motion of the continuous phase.The flow pattern and its morphological changes have a si...The flow focusing nozzle is a new type of nozzle that performs effective atomization of the discrete phase by means of high-speed motion of the continuous phase.The flow pattern and its morphological changes have a significant effect on the atomization, but the influence of different parameters on the morphological change of the flow pattern remains unclear.The flow focusing pattern and morphological changes in the two-phase flow inside the nozzle were simulated numerically, based on the volume of fluid method.The results demonstrate that the ratio of the nozzle-to-capillary distance and capillary diameter, the gas–liquid velocity ratio, and capillary diameter have significant effects on the flow pattern.When the ratio of the nozzle-to-capillary distance H and capillary diameter D increases, or the capillary diameter D increases, the flow pattern tends to transform into a laminar form; however, when the gas–liquid velocity ratio V increases, the flow pattern tends to transform into a turbulence form.Furthermore, we define the cone-shaped expansion rate, cone-shaped focusing rate,and cone angle in order to study the morphological changes in the cone shape inside the nozzle.The results indicate that the morphological change of the cone shape and flow pattern transformation is interrelated.When the cone shape tends to be unstable, the flow pattern changes towards flow blurring, whereas, a stable cone indicates that the flow tends to exhibit a droplet pattern.展开更多
Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid...Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid velocities.The flow rates in the two branches were measured accurately to determine how the two considered phases distribute in the two outlets.The experimental results have shown that when the two outlet pressures are asymmetric,the two-phase flow always tends to flow into the outlet which has a lower pressure.As the inlet liquid velocity increases,however,the two-phase flow gradually tends to split evenly.Compared with the experiment results,the pressure difference between the two outlets can be determined more accurately by means of numerical simulation.The trends of experimental results and simulations are in very good agreement.展开更多
Using a gas-solid two-phase model(a discrete phase model),the authors investigated the flow field inside the multi-channel nozzle for surface nanocrystallization(SNC)induced by the ultrasonic particulate peening(USPP)...Using a gas-solid two-phase model(a discrete phase model),the authors investigated the flow field inside the multi-channel nozzle for surface nanocrystallization(SNC)induced by the ultrasonic particulate peening(USPP).By computation,the velocity fields of both the gas and the solid phases were simulated and the track of the solid phase was analyzed in detail.It can be found that the velocities of the two phases are able to reach an ultrasonic level;meanwhile,the dispersion width of the solid phase at the nozzle exit is less than that of the gas phase.When particle diameters are less than 5 μm,there is a decreasing trend in the dispersion width of the solid phase with an increase in particle diameters.The trend becomes stable as the particle diameters are greater than 5 μm;in the meantime,the distribution of solid particles is near the axis of the jet flow.The optimal standoff distance between the nozzle and the substrate in the process of USPP is about 120 mm.Simulation results can help improve the design of mass-production-oriented multi-channel nozzles for SNC induced by USPP.展开更多
Numerical simulation on the flow,heat transfer and cracking reactions in commercial fluid catalyticcracking(FCC)riser reactors were carried out employing the developed turbulent gas-solid two-phase flow-reac-tion mode...Numerical simulation on the flow,heat transfer and cracking reactions in commercial fluid catalyticcracking(FCC)riser reactors were carried out employing the developed turbulent gas-solid two-phase flow-reac-tion model for FCC riser reactors given in Part Ⅰ of the present paper.Detailed information about the turbulentflow fields in the riser reactor obtained revealed the basic characteristics of the gas-solid two-phase turbulentflows when heat transfer and catalytic cracking reactions were co-existing in the riser.Results showed that thedistributions of the flow,the turbulence kinetic energy and the catalyst particle concentration are not uniform inthe axial,radial and tangential directions.The most complicated part of the riser reactor is the feed injectingzone.The complicated configuration of the turbulent gas-solid two-phase flows would exert a great influence onthe results of interphase heat transfer and cracking reactions.展开更多
The gas-solid two-phase flous of the precalciner were simulated by different multiphase models,such as mixture model,the Enderium model,including mixture and dispersed,and discrete phase model（DPM）,The results of th...The gas-solid two-phase flous of the precalciner were simulated by different multiphase models,such as mixture model,the Enderium model,including mixture and dispersed,and discrete phase model（DPM）,The results of the different multiphase models were analyzed and compared.showing the rationality of the diffusion and mixture of the cenment raic meals and coal poroder some extent Moreover,the results also shose the rationality of the given inlets parameters of actual process of the precalciner.展开更多
Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical sol...Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.展开更多
Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale...Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.展开更多
In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering comp...In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.展开更多
文摘Shale gas production involves complex gas-water two-phase flow,with flow patterns in proppant-filled fractures playing a critical role in determining production efficiency.In this study,3D geometric models of 40/70 mesh ceramic particles and quartz sand proppant clusters were elaborated using computed tomography(CT)scanning.These models were used to develop a numerical simulation framework based on the lattice Boltzmann method(LBM),enabling the investigation of gas-water flow behavior within proppant-filled fractures under varying driving forces and surface tensions.Simulation results at a closure pressure of 15 MPa have revealed that ceramic particles exhibit a simpler and more porous internal structure than quartz sand of the same size.Under identical flow conditions,ceramic proppants demonstrate higher fluid replacement efficiency.Replacement efficiency increases with higher porosity,greater driving force,and lower surface tension.Furthermore,fluid displacement is strongly influenced by pore geometry:flow is faster in straighter and wider channels,with preferential movement through larger pores forming dominant flow paths.The replacement velocity exhibits a characteristic time evolution,initially rapid,then gradually decreasing,correlating positively with the development of these dominant channels.
基金National Key Research and Development Program of China(2022YFB4600902)Shandong Provincial Science Foundation for Outstanding Young Scholars(ZR2024YQ020)。
文摘Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金National Natural Science Foundation (Project number 40604013).
文摘Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
基金supported by National Natural Science Foundation of China (Grant No. 21076198)Zhejiang Provincial Natural Science Foundation of China (Granted No. R1100530)National Basic Research Program of China (973 Program,Grant No. 2009CB724303)
文摘The flow with solid-liquid two-phase media inside centrifugal pumps is very complicated and the relevant method for the hydraulic design is still immature so far. There exist two main problems in the operation of the two-phase flow pumps, i.e., low overall efficiency and severe abrasion. In this study, the three-dimensional, steady, incompressible, and turbulent solid-liquid two-phase flows in a low-specific-speed centrifugal pump are numerically simulated and analyzed by using a computational fluid dynamics （CFD） code based on the mixture model of the two-phase flow and the RNG k-~ two-equation turbulence model, in which the influences of rotation and curvature are fully taken into account. The coupling between impeller and volute is implemented by means of the frozen rotor method. The simulation results predicted indicate that the solid phase properties in two-phase flow, especially the concentration, the particle diameter and the density, have strong effects on the hydraulic performance of the pump. Both the pump head and the efficiency are reduced with increasing particle diameter or concentration. However, the effect of particle density on the performance is relatively minor. An obvious jet-wake flow structure is presented near the volute tongue and becomes more remarkable with increasing solid phase concentration. The suction side of the blade is subject to much more severe abrasion than the pressure side. The obtained results preliminarily reveal the characteristics of solid-liquid two-phase flow in the centrifugal pump, and are helpful for improvement and empirical correction in the hydraulic design of centrifugal pumps.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
文摘A mathematical model is set to evaluate the 3-D dense solid-liquid two-phaseturbulent flow in a non-clogging mud pump, the flow feature in the impeller channel is simulatedwith the tool of IPSA. Meanwhile, resort to TECPLOT as the post-processor, the simulation results isvisualized. The results show the main flow characteristics: There exists backflow and aberrantvelocities at inlet area and a relative velocity slip between two phases; A jet-wake flow pattern isdiscerned around the shroud-suction side area; The relative velocity vector of solid phase iscloser to the pressure surface than that of liquid phase and the trend is more obvious with theincrease of diameter; The kinetic energy of turbulence k and the dissipation rate e reach theirpeaks at the corner of pressure and suction surface. The simulation results show a good agreementwith the experimental flow features in the impeller channel, which prove the turbulent model used isvalid and provide a theoretical design basis to non-clogging pumps.
基金supported by the Fund of Innovation Research Group of National Natural Science Foundation of China (Grant NO.5052160450323001)Major Program of National Natural Science Foundation of China (Grant No.50536020)
文摘Chokes are one of the most important components of downhole flow-control equipment. The particle erosion mathematical model, which considers particle-particle interaction, was established and used to simulate solid particle movement as well as particle erosion characteristics of the solid-liquid two-phase flow in a choke. The corresponding erosion reduction approach by setting ribs on the inner wall of the choke was advanced. This mathematical model includes three parts： the flow field simulation of the continuous carrier fluid by an Eulerian approach, the particle interaction simulation using the discrete particle hard sphere model by a Lagrangian approach and calculation of erosion rate using semiempirical correlations. The results show that particles accumulated in a narrow region from inlet to outlet of the choke and the dominating factor affecting particle motion is the fluid drag force. As a result, the optimization of rib geometrical parameters indicates that good anti-erosion performance can be achieved by four ribs, each of them with a height （H） of 3 mm and a width （B） of 5 mm equaling the interval between ribs （L）.
基金Project(2018YFC1901606)supported by the National Key R&D Program of China.
文摘Taking an oxygen enriched side-blown furnace as the prototype,a hydraulic model was established according to the similarity principle.The influence of three factors on the gas-liquid two-phase flow was analyzed,i.e.the airflow speed,the submerged depth and the downward angle of the nozzle.A numerical simulation of the hydraulic model was carried out trying to find the suitable turbulence model which can describe the side-blown two-phase flow correctly by comparing the simulation results with the experimental data.The experiment shows that the airflow speed has a great influence on the flow of the water.The submerged depth of the nozzle has a relatively smaller influence on the penetration depth and the surface fluctuation height in the liquid phase.When the nozzle is at a downward angle of 15°,the penetration depth and the surface fluctuation height are reduced.It is concluded that the numerical results with the realizable k-εturbulence model are the closest to the experiment for the penetration depth,the surface fluctuation height and the bubble scale.
基金Project supported by the Scientific Research Foundation of Hebei University of Technology of China
文摘Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.
文摘Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.
基金Sponsored by the Key Project of Chinese Ministry of Education(108135)Program for New Century Excellent Talentsin University(NCET-07-0799)+2 种基金National Eleventh Five-Year Key Science & Technology Project(2006BAK03B0303)the Fok Ying-Tong Education Foundation for Young Teachers in the Higher Education Institutions of China(111053)Beijing Science & Technology New StarPlan(2006A81)
文摘According to the Lagrange discrete phase model of multiphase flow, mathematical model of the abrasive water-jet （AWJ） nozzle based on the multi-phase movement was introduced, then the boundary conditions was determined and the liquid-solid turbulence which is isothermal, can not be compressed and steadystate in the cone-cylinder nozzles of the export of the pre-mixed AWJ was simulated applying the software FLUENT. The results showed that： the axial velocity and dynamic pressure of the continuous phase in the nozzle were axial symmetry notable, and at the axis had a extreme point; abrasive accelerated at two points, in front contractive segment, the rate increased rapidly, and in the back straight one, the speed accelerated slowly. The length of the cylinder is 100 mm, the diameter of the nozzle is 8 mm, and the angle of the cone is 15°. There is a extreme point of the rate at the point 10 mm in the established model. The results of simulation laid the foundation for optimizing the nozzle structure, improving efficiency and developing the nozzle.
基金Supported by the National Natural Science Foundation of China(51776016,51606006)Beijing Natural Science Foundation(3172025,3182030)+2 种基金National Key Research and Development Program(2017YFB0103401)National Engineering Laboratory for Mobile Source Emission Control Technology(NELMS2017A10)the Talents Foundation of Beijing Jiaotong University(2018RC017)
文摘The flow focusing nozzle is a new type of nozzle that performs effective atomization of the discrete phase by means of high-speed motion of the continuous phase.The flow pattern and its morphological changes have a significant effect on the atomization, but the influence of different parameters on the morphological change of the flow pattern remains unclear.The flow focusing pattern and morphological changes in the two-phase flow inside the nozzle were simulated numerically, based on the volume of fluid method.The results demonstrate that the ratio of the nozzle-to-capillary distance and capillary diameter, the gas–liquid velocity ratio, and capillary diameter have significant effects on the flow pattern.When the ratio of the nozzle-to-capillary distance H and capillary diameter D increases, or the capillary diameter D increases, the flow pattern tends to transform into a laminar form; however, when the gas–liquid velocity ratio V increases, the flow pattern tends to transform into a turbulence form.Furthermore, we define the cone-shaped expansion rate, cone-shaped focusing rate,and cone angle in order to study the morphological changes in the cone shape inside the nozzle.The results indicate that the morphological change of the cone shape and flow pattern transformation is interrelated.When the cone shape tends to be unstable, the flow pattern changes towards flow blurring, whereas, a stable cone indicates that the flow tends to exhibit a droplet pattern.
基金the National Science and Technology Major Project of China(No.2016ZX05028-004-003).
文摘Dedicated experiments and numerical simulations have been conducted to investigate the splitting characteristics of a gas-liquid two phase flow at a T junction.The experiments were carried out for different gas-liquid velocities.The flow rates in the two branches were measured accurately to determine how the two considered phases distribute in the two outlets.The experimental results have shown that when the two outlet pressures are asymmetric,the two-phase flow always tends to flow into the outlet which has a lower pressure.As the inlet liquid velocity increases,however,the two-phase flow gradually tends to split evenly.Compared with the experiment results,the pressure difference between the two outlets can be determined more accurately by means of numerical simulation.The trends of experimental results and simulations are in very good agreement.
基金supported by the National High-Tech.R&D Program of China(the National 863 plans projects,Grant No.2007AA03Z352)
文摘Using a gas-solid two-phase model(a discrete phase model),the authors investigated the flow field inside the multi-channel nozzle for surface nanocrystallization(SNC)induced by the ultrasonic particulate peening(USPP).By computation,the velocity fields of both the gas and the solid phases were simulated and the track of the solid phase was analyzed in detail.It can be found that the velocities of the two phases are able to reach an ultrasonic level;meanwhile,the dispersion width of the solid phase at the nozzle exit is less than that of the gas phase.When particle diameters are less than 5 μm,there is a decreasing trend in the dispersion width of the solid phase with an increase in particle diameters.The trend becomes stable as the particle diameters are greater than 5 μm;in the meantime,the distribution of solid particles is near the axis of the jet flow.The optimal standoff distance between the nozzle and the substrate in the process of USPP is about 120 mm.Simulation results can help improve the design of mass-production-oriented multi-channel nozzles for SNC induced by USPP.
文摘Numerical simulation on the flow,heat transfer and cracking reactions in commercial fluid catalyticcracking(FCC)riser reactors were carried out employing the developed turbulent gas-solid two-phase flow-reac-tion model for FCC riser reactors given in Part Ⅰ of the present paper.Detailed information about the turbulentflow fields in the riser reactor obtained revealed the basic characteristics of the gas-solid two-phase turbulentflows when heat transfer and catalytic cracking reactions were co-existing in the riser.Results showed that thedistributions of the flow,the turbulence kinetic energy and the catalyst particle concentration are not uniform inthe axial,radial and tangential directions.The most complicated part of the riser reactor is the feed injectingzone.The complicated configuration of the turbulent gas-solid two-phase flows would exert a great influence onthe results of interphase heat transfer and cracking reactions.
文摘The gas-solid two-phase flous of the precalciner were simulated by different multiphase models,such as mixture model,the Enderium model,including mixture and dispersed,and discrete phase model（DPM）,The results of the different multiphase models were analyzed and compared.showing the rationality of the diffusion and mixture of the cenment raic meals and coal poroder some extent Moreover,the results also shose the rationality of the given inlets parameters of actual process of the precalciner.
文摘Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.
基金Supported by the National Natural Science Foundation of China(52374043)Key Program of the National Natural Science Foundation of China(52234003).
文摘Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.
文摘In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.
