Two-phase pipe flow occurs frequently in oil&gas industry,nuclear power plants,and CCUS.Reliable calculations of gas void fraction(or liquid holdup)play a central role in two-phase pipe flow models.In this paper w...Two-phase pipe flow occurs frequently in oil&gas industry,nuclear power plants,and CCUS.Reliable calculations of gas void fraction(or liquid holdup)play a central role in two-phase pipe flow models.In this paper we apply the fractional flow theory to multiphase flow in pipes and present a unified modeling framework for predicting the fluid phase volume fractions over a broad range of pipe flow conditions.Compared to existing methods and correlations,this new framework provides a simple,approximate,and efficient way to estimate the phase volume fraction in two-phase pipe flow without invoking flow patterns.Notably,existing correlations for estimating phase volume fraction can be transformed and expressed under this modeling framework.Different fractional flow models are applicable to different flow conditions,and they demonstrate good agreement against experimental data within 5%errors when compared with an experimental database comprising of 2754 data groups from 14literature sources,covering various pipe geometries,flow patterns,fluid properties and flow inclinations.The gas void fraction predicted by the framework developed in this work can be used as inputs to reliably model the hydraulic and thermal behaviors of two-phase pipe flows.展开更多
There is a strong relationship between analytical and numerical heat transfers due to thermodynamically anticipated findings,making thermo-dynamical modeling an effective tool for estimating the ideal melting point of...There is a strong relationship between analytical and numerical heat transfers due to thermodynamically anticipated findings,making thermo-dynamical modeling an effective tool for estimating the ideal melting point of heat transfer.Under certain assumptions,the present study builds a mathematical model of melting heat transport nanofluid flow of chemical reactions and joule heating.Nanofluid flow is described by higher-order partial non-linear differential equations.Incorporating suitable similarity transformations and dimensionless parameters converts these controlling partial differential equations into the non-linear ordinary differential equations and resulting system of nonlinear equations is established.Plotted graphic visualizations in MATLAB allow for an indepth analysis of the effects of distinguishing factors on fluid flow.Innovative applications of the findings include electronic cooling,heat transfer,reaction processes,nuclear reactors,micro heat pipes,and other related fields.If the exponential index increases,however,the thermal profile becomes worse.By comparing the current findings to those already published in the literature for this particular example,we find that they are highly congruent,therefore validating the present work.Every one of the numerical findings exhibits asymptotic behavior by meeting the specified boundary conditions.展开更多
Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on t...Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.展开更多
Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to ...Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.展开更多
Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typica...Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typically governed by gravity, viscous and capillary forces, these factors lead invasive fluids to occupy pore space irregularly and incompletely. Previous studies have demonstrated capillary numbers,describing the viscous and capillary forces, to quantificationally induce evolution of invasion patterns.While the evolution mechanisms of invasive patterns have not been deeply elucidated under the constant capillary number and three variable parameters including velocity, viscosity, and interfacial tension.Our research employs two horizontal visualization systems and a two-phase laminar flow simulation to investigate the tendency of invasive pattern transition by various parameters at the pore scale. We showed that increasing invasive viscosity or reducing interfacial tension in a homogeneous pore space significantly enhanced sweep efficiency, under constant capillary number. Additionally, in the fingering crossover pattern, the region near the inlet was prone to capillary fingering with multi-directional invasion, while the viscous fingering with unidirectional invasion was more susceptible occurred in the region near the outlet. Furthermore, increasing invasive viscosity or decreasing invasive velocity and interfacial tension promoted the extension of viscous fingering from the outlet to the inlet, presenting that the subsequent invasive fluid flows toward the outlet. In the case of invasive trunk along a unidirectional path, the invasive flow increased exponentially closer to the outlet, resulting in a significant decrease in the width of the invasive interface. Our work holds promising applications for optimizing invasive patterns in heterogeneous porous media.展开更多
Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinni...Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.展开更多
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.展开更多
Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple the...Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.展开更多
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.展开更多
Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement m...Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.展开更多
In response to escalating challenges in energy conservation and emission reduction,this study delves into the complexities of heat transfer in two-phase flows and adjustments to combustion processes within coal-fired ...In response to escalating challenges in energy conservation and emission reduction,this study delves into the complexities of heat transfer in two-phase flows and adjustments to combustion processes within coal-fired boilers.Utilizing a fusion of hybrid modeling and automation technologies,we develop soft measurement models for key combustion parameters,such as the net calorific value of coal,flue gas oxygen content,and fly ash carbon content,within theDistributedControl System(DCS).Validated with performance test data,thesemodels exhibit controlled root mean square error(RMSE)and maximum absolute error(MAXE)values,both within the range of 0.203.Integrated into their respective automatic control systems,thesemodels optimize two-phase flow heat transfer,finetune combustion conditions,and mitigate incomplete combustion.Furthermore,this paper conducts an in-depth exploration of the generationmechanismof nitrogen oxides(NO_(x))and low oxygen emission reduction technology in coal-fired boilers,demonstrating a substantial reduction in furnace exit NO_(x) generation by 30%to 40%and the power supply coal consumption decreased by 1.62 g/(kW h).The research outcomes highlight the model’s rapid responsiveness,enabling prompt reflection of transient variations in various economic indicator parameters.This provides a more effective means for real-time monitoring of crucial variables in coal-fired boilers and facilitates timely combustion adjustments,underscoring notable achievements in boiler combustion.The research not only provides valuable and practical insights into the intricacies of two-phase flow heat transfer and heat exchange but also establishes a pioneering methodology for tackling industry challenges.展开更多
Due to the complex high-temperature characteristics of hydrocarbon fuel,the research on the long-term working process of parallel channel structure under variable working conditions,especially under high heat-mass rat...Due to the complex high-temperature characteristics of hydrocarbon fuel,the research on the long-term working process of parallel channel structure under variable working conditions,especially under high heat-mass ratio,has not been systematically carried out.In this paper,the heat transfer and flow characteristics of related high temperature fuels are studied by using typical engine parallel channel structure.Through numeri⁃cal simulation and systematic experimental verification,the flow and heat transfer characteristics of parallel chan⁃nels under typical working conditions are obtained,and the effectiveness of high-precision calculation method is preliminarily established.It is known that the stable time required for hot start of regenerative cooling engine is about 50 s,and the flow resistance of parallel channel structure first increases and then decreases with the in⁃crease of equivalence ratio(The following equivalence ratio is expressed byΦ),and there is a flow resistance peak in the range ofΦ=0.5~0.8.This is mainly caused by the coupling effect of high temperature physical proper⁃ties,flow rate and pressure of fuel in parallel channels.At the same time,the cooling and heat transfer character⁃istics of parallel channels under some conditions of high heat-mass ratio are obtained,and the main factors affect⁃ing the heat transfer of parallel channels such as improving surface roughness and strengthening heat transfer are mastered.In the experiment,whenΦis less than 0.9,the phenomenon of local heat transfer enhancement and deterioration can be obviously observed,and the temperature rise of local structures exceeds 200℃,which is the risk of structural damage.Therefore,the reliability of long-term parallel channel structure under the condition of high heat-mass ratio should be fully considered in structural design.展开更多
This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with ...This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with a filter-based density correction turbulence model and a modified Zwart cavitation model.The study investigates the dynamic cavitation features of the thermal fluid around the hydrofoil at various incoming flow velocities.It systematically elucidates the evolution of cavitation and vortex dynamics corresponding to each velocity condition.The results indicate that with increasing incoming flow velocity,distinct cavitation processes take place in the flow field.展开更多
Gas-liquid flow(GLF),especially slug and annular flows in oil and gas gathering and transportation pipelines,become particularly complex inside elbows and can easily exacerbate pipeline corrosion and damage.In thisstu...Gas-liquid flow(GLF),especially slug and annular flows in oil and gas gathering and transportation pipelines,become particularly complex inside elbows and can easily exacerbate pipeline corrosion and damage.In thisstudy,FLUENT was used to conduct 3D simulations of slug and annular flow in elbows for different velocitiesto assess the ensuing changes in terms of pressure.In particular,the multifluid VOF(Volume of Fraction)modelwas chosen.The results indicate that under both slug and annular flow conditions,the pressure inside the elbow islower than the outside.As the superficial velocity of liquid and gas increase,the pressure and liquid flow velocityat different positions of the elbow also increase,while the secondary flow weakens.Under annular flow conditions,the liquid film on the outer side of the elbow is thicker than that on the inner side,and the liquid velocityin the main liquid film zone is the lowest.展开更多
The third-order flow Gerdjikov–Ivanov(TOFGI)equation is studied,and the Darboux transformation(DT)is used to obtain the determinant expression of the solution of this equation.On this basis,the soliton solution,ratio...The third-order flow Gerdjikov–Ivanov(TOFGI)equation is studied,and the Darboux transformation(DT)is used to obtain the determinant expression of the solution of this equation.On this basis,the soliton solution,rational solution,positon solution,and breather solution of the TOFGI equation are obtained by taking zero seed solution and non-zero seed solution.The exact solutions and dynamic properties of the Gerdjikov–Ivanov(GI)equation and the TOFGI equation are compared in detail under the same conditions,and it is found that there are some differences in the velocities and trajectories of the solutions of the two equations.展开更多
As the proportion of newenergy increases,the traditional cumulant method(CM)produces significant errorswhen performing probabilistic load flow(PLF)calculations with large-scale wind power integrated.Considering the wi...As the proportion of newenergy increases,the traditional cumulant method(CM)produces significant errorswhen performing probabilistic load flow(PLF)calculations with large-scale wind power integrated.Considering the wind speed correlation,a multi-scenario PLF calculation method that combines random sampling and segmented discrete wind farm power was proposed.Firstly,based on constructing discrete scenes of wind farms,the Nataf transform is used to handle the correlation between wind speeds.Then,the random sampling method determines the output probability of discrete wind power scenarios when wind speed exhibits correlation.Finally,the PLF calculation results of each scenario areweighted and superimposed following the total probability formula to obtain the final power flow calculation result.Verified in the IEEE standard node system,the absolute percent error(APE)for the mean and standard deviation(SD)of the node voltages and branch active power are all within 1%,and the average root mean square(AMSR)values of the probability curves are all less than 1%.展开更多
Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphas...Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphaseflow within the furnace.Understanding the flow structure and temperature distribution in this setup is crucial foroptimizing production.In this study,gas-liquid interactions,and temperature profiles under varying air-injectionconditions are examined by means of numerical simulation for a 3.2 m×20 m furnace.The results indicate that thehigh-velocity regions are essentially distributed near the lance within the reaction region and the flue gas outlet,while low-velocity regions are located close to the furnace walls on both side of the reaction region.Dead regionsappear in the sedimentation region,with gas velocities surpassing those of the molten phase.As the injection rateincreases from 0.50 to 0.80 Nm3/s,the stabilization time of the average liquid surface velocity decreases from 2.6 sto 1.9 s,exhibiting a similar trend to the gas holdup.During stabilization,the average liquid surface velocity risesfrom 0.505 to 0.702 m/s.The average turbulent kinetic energy(TKE)of the fluid in the molten bath increases from0.095 to 0.162 m^(2)/s^(2).The proportion of the area distribution with TKE greater than 0.10 m^(2)/s^(2) and the gas holdupat steady state both rise with an increase in the injection quantity.The maximum splashing height of the melt growsfrom approximately 0.756 to 1.154 m,with the affected area expanding from 14.239 to 20.498 m^(2).Under differentworking conditions with varying injection quantities,the average temperature changes in melt zone and flue gaszone of the furnace are small.The temperature in the melt and in the flue-gas zone spans the interval 1200℃–1257℃,and 1073℃–1121℃,respectively.The temperature distribution of the melt and flue gas reveals a patterncharacterized by elevated temperatures in the reaction zone,gradually transitioning to lower temperatures in thesedimentation region.展开更多
Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical s...Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical simulation study on the dynamics of blood flow in a stenosed artery,focusing on the effects of copper and alumina nanoparticles,is conducted.The study employs a 2-dimensional Newtonian blood flow model infused with copper and alumina nanoparticles,considering the influence of a magnetic field,thermal radiation,and various flow parameters.The governing differential equations are first non-dimensionalized to facilitate analysis and subsequently solved using the 4th order collocation method,bvp4c module in MATLAB.This approach obtains velocity and temperature profiles,revealing the impact of relevant parameters crucial in the biomedical field.The findings of this study underscore the significance of understanding blood flow dynamics in stenosed arteries and the potential benefits of utilizing copper and alumina nanoparticles in treatment strategies.The incorporation of nanoparticles introduces novel avenues for enhancing therapeutic interventions,particularly in mitigating the effects of stenosis.The elucidation of velocity and temperature profiles provides valuable insights into the behavior of blood flow under different conditions,thereby informing the development of targeted biomedical applications.The arterial curvature flow parameter influences temperature profiles,with increased parameters promoting more efficient heat dissipation.The elevated values of Prandtl number and thermal radiation parameter showcase the diminished temperature profiles,indicating stronger dominance of momentum diffusion over thermal diffusion and radiative heat transfer mechanism.Sensitivity analysis of the pertinent physical parameters reveals that the Prandtl number has the most significant impact on blood flow dynamics.A statistical analysis of the present results and existing literature has also been included in the study.Overall,this research contributes to advancing our understanding of vascular health and lays the groundwork for innovative approaches in stenosis treatment and related biomedical fields.展开更多
With the widespread implementation of distributed generation(DG)and the integration of soft open point(SOP)into the distribution network(DN),the latter is steadily transitioning into a flexible distribution network(FD...With the widespread implementation of distributed generation(DG)and the integration of soft open point(SOP)into the distribution network(DN),the latter is steadily transitioning into a flexible distribution network(FDN),the calculation of carbon flow distribution in FDN is more difficult.To this end,this study constructs a model for low-carbon optimal operations within the FDN on the basis of enhanced carbon emission flow(CEF).First,the carbon emission characteristics of FDNs are scrutinized and an improved method for calculating carbon flow within these networks is proposed.Subsequently,a model for optimizing low-carbon operations within FDNs is formulated based on the refined CEF,which merges the specificities of DG and intelligent SOP.Finally,this model is scrutinized using an upgraded IEEE 33-node distribution system,a comparative analysis of the cases reveals that when DG and SOP are operated in a coordinated manner in the FDN,with the cost of electricity generation was reduced by 40.63 percent and the cost of carbon emissions by 10.18 percent.The findings indicate that the judicious optimization of areas exhibiting higher carbon flow rates can effectively enhance the economic efficiency of DN operations and curtail the carbon emissions of the overall network.展开更多
基金financial support from the Energize Program between the University of Texas at Austin and Southwest Research InstituteHydraulic Fracturing and Sand Control Industrial Affiliates Program at the University of Texas at Austin for financially supporting this research。
文摘Two-phase pipe flow occurs frequently in oil&gas industry,nuclear power plants,and CCUS.Reliable calculations of gas void fraction(or liquid holdup)play a central role in two-phase pipe flow models.In this paper we apply the fractional flow theory to multiphase flow in pipes and present a unified modeling framework for predicting the fluid phase volume fractions over a broad range of pipe flow conditions.Compared to existing methods and correlations,this new framework provides a simple,approximate,and efficient way to estimate the phase volume fraction in two-phase pipe flow without invoking flow patterns.Notably,existing correlations for estimating phase volume fraction can be transformed and expressed under this modeling framework.Different fractional flow models are applicable to different flow conditions,and they demonstrate good agreement against experimental data within 5%errors when compared with an experimental database comprising of 2754 data groups from 14literature sources,covering various pipe geometries,flow patterns,fluid properties and flow inclinations.The gas void fraction predicted by the framework developed in this work can be used as inputs to reliably model the hydraulic and thermal behaviors of two-phase pipe flows.
文摘There is a strong relationship between analytical and numerical heat transfers due to thermodynamically anticipated findings,making thermo-dynamical modeling an effective tool for estimating the ideal melting point of heat transfer.Under certain assumptions,the present study builds a mathematical model of melting heat transport nanofluid flow of chemical reactions and joule heating.Nanofluid flow is described by higher-order partial non-linear differential equations.Incorporating suitable similarity transformations and dimensionless parameters converts these controlling partial differential equations into the non-linear ordinary differential equations and resulting system of nonlinear equations is established.Plotted graphic visualizations in MATLAB allow for an indepth analysis of the effects of distinguishing factors on fluid flow.Innovative applications of the findings include electronic cooling,heat transfer,reaction processes,nuclear reactors,micro heat pipes,and other related fields.If the exponential index increases,however,the thermal profile becomes worse.By comparing the current findings to those already published in the literature for this particular example,we find that they are highly congruent,therefore validating the present work.Every one of the numerical findings exhibits asymptotic behavior by meeting the specified boundary conditions.
文摘Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.
基金the support of the Opening Fund of State Key Laboratory of Multiphase Flow in Power Engineering(SKLMF-KF-2102)。
文摘Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.
基金supported by the National Natural Science Foundation of China Joint Fund Project (Grant/Award Number: U20B6003)National Natural Science Foundation of China (Grant/Award Number: 52304054)。
文摘Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typically governed by gravity, viscous and capillary forces, these factors lead invasive fluids to occupy pore space irregularly and incompletely. Previous studies have demonstrated capillary numbers,describing the viscous and capillary forces, to quantificationally induce evolution of invasion patterns.While the evolution mechanisms of invasive patterns have not been deeply elucidated under the constant capillary number and three variable parameters including velocity, viscosity, and interfacial tension.Our research employs two horizontal visualization systems and a two-phase laminar flow simulation to investigate the tendency of invasive pattern transition by various parameters at the pore scale. We showed that increasing invasive viscosity or reducing interfacial tension in a homogeneous pore space significantly enhanced sweep efficiency, under constant capillary number. Additionally, in the fingering crossover pattern, the region near the inlet was prone to capillary fingering with multi-directional invasion, while the viscous fingering with unidirectional invasion was more susceptible occurred in the region near the outlet. Furthermore, increasing invasive viscosity or decreasing invasive velocity and interfacial tension promoted the extension of viscous fingering from the outlet to the inlet, presenting that the subsequent invasive fluid flows toward the outlet. In the case of invasive trunk along a unidirectional path, the invasive flow increased exponentially closer to the outlet, resulting in a significant decrease in the width of the invasive interface. Our work holds promising applications for optimizing invasive patterns in heterogeneous porous media.
基金This work is supported by the National Natural Science Foundation of China(No.52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(No.BYESS2023262)Science Foundation of China University of Petroleum,Beijing(No.2462022BJRC004).
文摘Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.
基金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.
基金supported by the NSFC Grant no.12271492the Natural Science Foundation of Henan Province of China Grant no.222300420550+1 种基金supported by the NSFC Grant no.12271498the National Key R&D Program of China Grant no.2022YFA1005202/2022YFA1005200.
文摘Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant Nos.22275092,52102107 and 52372084)the Fundamental Research Funds for the Central Universities(Grant No.30923010920)。
文摘Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.
文摘In response to escalating challenges in energy conservation and emission reduction,this study delves into the complexities of heat transfer in two-phase flows and adjustments to combustion processes within coal-fired boilers.Utilizing a fusion of hybrid modeling and automation technologies,we develop soft measurement models for key combustion parameters,such as the net calorific value of coal,flue gas oxygen content,and fly ash carbon content,within theDistributedControl System(DCS).Validated with performance test data,thesemodels exhibit controlled root mean square error(RMSE)and maximum absolute error(MAXE)values,both within the range of 0.203.Integrated into their respective automatic control systems,thesemodels optimize two-phase flow heat transfer,finetune combustion conditions,and mitigate incomplete combustion.Furthermore,this paper conducts an in-depth exploration of the generationmechanismof nitrogen oxides(NO_(x))and low oxygen emission reduction technology in coal-fired boilers,demonstrating a substantial reduction in furnace exit NO_(x) generation by 30%to 40%and the power supply coal consumption decreased by 1.62 g/(kW h).The research outcomes highlight the model’s rapid responsiveness,enabling prompt reflection of transient variations in various economic indicator parameters.This provides a more effective means for real-time monitoring of crucial variables in coal-fired boilers and facilitates timely combustion adjustments,underscoring notable achievements in boiler combustion.The research not only provides valuable and practical insights into the intricacies of two-phase flow heat transfer and heat exchange but also establishes a pioneering methodology for tackling industry challenges.
文摘Due to the complex high-temperature characteristics of hydrocarbon fuel,the research on the long-term working process of parallel channel structure under variable working conditions,especially under high heat-mass ratio,has not been systematically carried out.In this paper,the heat transfer and flow characteristics of related high temperature fuels are studied by using typical engine parallel channel structure.Through numeri⁃cal simulation and systematic experimental verification,the flow and heat transfer characteristics of parallel chan⁃nels under typical working conditions are obtained,and the effectiveness of high-precision calculation method is preliminarily established.It is known that the stable time required for hot start of regenerative cooling engine is about 50 s,and the flow resistance of parallel channel structure first increases and then decreases with the in⁃crease of equivalence ratio(The following equivalence ratio is expressed byΦ),and there is a flow resistance peak in the range ofΦ=0.5~0.8.This is mainly caused by the coupling effect of high temperature physical proper⁃ties,flow rate and pressure of fuel in parallel channels.At the same time,the cooling and heat transfer character⁃istics of parallel channels under some conditions of high heat-mass ratio are obtained,and the main factors affect⁃ing the heat transfer of parallel channels such as improving surface roughness and strengthening heat transfer are mastered.In the experiment,whenΦis less than 0.9,the phenomenon of local heat transfer enhancement and deterioration can be obviously observed,and the temperature rise of local structures exceeds 200℃,which is the risk of structural damage.Therefore,the reliability of long-term parallel channel structure under the condition of high heat-mass ratio should be fully considered in structural design.
文摘This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with a filter-based density correction turbulence model and a modified Zwart cavitation model.The study investigates the dynamic cavitation features of the thermal fluid around the hydrofoil at various incoming flow velocities.It systematically elucidates the evolution of cavitation and vortex dynamics corresponding to each velocity condition.The results indicate that with increasing incoming flow velocity,distinct cavitation processes take place in the flow field.
基金supported by the Ministry of Industry and Information Technology High Tech Ship Special Project(Grant No.CBG3N21-2-6).
文摘Gas-liquid flow(GLF),especially slug and annular flows in oil and gas gathering and transportation pipelines,become particularly complex inside elbows and can easily exacerbate pipeline corrosion and damage.In thisstudy,FLUENT was used to conduct 3D simulations of slug and annular flow in elbows for different velocitiesto assess the ensuing changes in terms of pressure.In particular,the multifluid VOF(Volume of Fraction)modelwas chosen.The results indicate that under both slug and annular flow conditions,the pressure inside the elbow islower than the outside.As the superficial velocity of liquid and gas increase,the pressure and liquid flow velocityat different positions of the elbow also increase,while the secondary flow weakens.Under annular flow conditions,the liquid film on the outer side of the elbow is thicker than that on the inner side,and the liquid velocityin the main liquid film zone is the lowest.
基金Project supported by the National Natural Science Foundation of China(Grant No.12201329)the Zhejiang Provincial Natural Science Foundation of China(Grant No.LY24A010002)the Natural Science Foundation of Ningbo(Grant No.2023J126)。
文摘The third-order flow Gerdjikov–Ivanov(TOFGI)equation is studied,and the Darboux transformation(DT)is used to obtain the determinant expression of the solution of this equation.On this basis,the soliton solution,rational solution,positon solution,and breather solution of the TOFGI equation are obtained by taking zero seed solution and non-zero seed solution.The exact solutions and dynamic properties of the Gerdjikov–Ivanov(GI)equation and the TOFGI equation are compared in detail under the same conditions,and it is found that there are some differences in the velocities and trajectories of the solutions of the two equations.
基金supported by Basic Science Research Program through the National Natural Science Foundation of China(Grant No.61867003).
文摘As the proportion of newenergy increases,the traditional cumulant method(CM)produces significant errorswhen performing probabilistic load flow(PLF)calculations with large-scale wind power integrated.Considering the wind speed correlation,a multi-scenario PLF calculation method that combines random sampling and segmented discrete wind farm power was proposed.Firstly,based on constructing discrete scenes of wind farms,the Nataf transform is used to handle the correlation between wind speeds.Then,the random sampling method determines the output probability of discrete wind power scenarios when wind speed exhibits correlation.Finally,the PLF calculation results of each scenario areweighted and superimposed following the total probability formula to obtain the final power flow calculation result.Verified in the IEEE standard node system,the absolute percent error(APE)for the mean and standard deviation(SD)of the node voltages and branch active power are all within 1%,and the average root mean square(AMSR)values of the probability curves are all less than 1%.
基金Supported by Yunnan Fundamental Research Projects(Nos.202301AT070469,202301AT070275)Supported by Yunnan Major Scientific and Technological Projects(No.202202AG050002).
文摘Smelting with oxygen bottom blowing is one of the main methods used in the frame of copper pyrometallurgy.With this approach,feed materials and oxygen-enriched air are introduced in reversed order to enhance multiphaseflow within the furnace.Understanding the flow structure and temperature distribution in this setup is crucial foroptimizing production.In this study,gas-liquid interactions,and temperature profiles under varying air-injectionconditions are examined by means of numerical simulation for a 3.2 m×20 m furnace.The results indicate that thehigh-velocity regions are essentially distributed near the lance within the reaction region and the flue gas outlet,while low-velocity regions are located close to the furnace walls on both side of the reaction region.Dead regionsappear in the sedimentation region,with gas velocities surpassing those of the molten phase.As the injection rateincreases from 0.50 to 0.80 Nm3/s,the stabilization time of the average liquid surface velocity decreases from 2.6 sto 1.9 s,exhibiting a similar trend to the gas holdup.During stabilization,the average liquid surface velocity risesfrom 0.505 to 0.702 m/s.The average turbulent kinetic energy(TKE)of the fluid in the molten bath increases from0.095 to 0.162 m^(2)/s^(2).The proportion of the area distribution with TKE greater than 0.10 m^(2)/s^(2) and the gas holdupat steady state both rise with an increase in the injection quantity.The maximum splashing height of the melt growsfrom approximately 0.756 to 1.154 m,with the affected area expanding from 14.239 to 20.498 m^(2).Under differentworking conditions with varying injection quantities,the average temperature changes in melt zone and flue gaszone of the furnace are small.The temperature in the melt and in the flue-gas zone spans the interval 1200℃–1257℃,and 1073℃–1121℃,respectively.The temperature distribution of the melt and flue gas reveals a patterncharacterized by elevated temperatures in the reaction zone,gradually transitioning to lower temperatures in thesedimentation region.
基金funded by Universiti Teknikal Malaysia Melaka and Ministry of Higher Education(MoHE)Malaysia,grant number FRGS/1/2024/FTKM/F00586.
文摘Nanotechnology holds immense importance in the biomedical field due to its ability to revolutionize healthcare on a molecular scale.Motivated by the imperative of enhancing patient outcomes,a comprehensive numerical simulation study on the dynamics of blood flow in a stenosed artery,focusing on the effects of copper and alumina nanoparticles,is conducted.The study employs a 2-dimensional Newtonian blood flow model infused with copper and alumina nanoparticles,considering the influence of a magnetic field,thermal radiation,and various flow parameters.The governing differential equations are first non-dimensionalized to facilitate analysis and subsequently solved using the 4th order collocation method,bvp4c module in MATLAB.This approach obtains velocity and temperature profiles,revealing the impact of relevant parameters crucial in the biomedical field.The findings of this study underscore the significance of understanding blood flow dynamics in stenosed arteries and the potential benefits of utilizing copper and alumina nanoparticles in treatment strategies.The incorporation of nanoparticles introduces novel avenues for enhancing therapeutic interventions,particularly in mitigating the effects of stenosis.The elucidation of velocity and temperature profiles provides valuable insights into the behavior of blood flow under different conditions,thereby informing the development of targeted biomedical applications.The arterial curvature flow parameter influences temperature profiles,with increased parameters promoting more efficient heat dissipation.The elevated values of Prandtl number and thermal radiation parameter showcase the diminished temperature profiles,indicating stronger dominance of momentum diffusion over thermal diffusion and radiative heat transfer mechanism.Sensitivity analysis of the pertinent physical parameters reveals that the Prandtl number has the most significant impact on blood flow dynamics.A statistical analysis of the present results and existing literature has also been included in the study.Overall,this research contributes to advancing our understanding of vascular health and lays the groundwork for innovative approaches in stenosis treatment and related biomedical fields.
基金supported in part by National Natural Science Foundation of China under Grant 52007026.
文摘With the widespread implementation of distributed generation(DG)and the integration of soft open point(SOP)into the distribution network(DN),the latter is steadily transitioning into a flexible distribution network(FDN),the calculation of carbon flow distribution in FDN is more difficult.To this end,this study constructs a model for low-carbon optimal operations within the FDN on the basis of enhanced carbon emission flow(CEF).First,the carbon emission characteristics of FDNs are scrutinized and an improved method for calculating carbon flow within these networks is proposed.Subsequently,a model for optimizing low-carbon operations within FDNs is formulated based on the refined CEF,which merges the specificities of DG and intelligent SOP.Finally,this model is scrutinized using an upgraded IEEE 33-node distribution system,a comparative analysis of the cases reveals that when DG and SOP are operated in a coordinated manner in the FDN,with the cost of electricity generation was reduced by 40.63 percent and the cost of carbon emissions by 10.18 percent.The findings indicate that the judicious optimization of areas exhibiting higher carbon flow rates can effectively enhance the economic efficiency of DN operations and curtail the carbon emissions of the overall network.