Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

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.

MATHEMATICAL MODEL OF TWO-PHASE FLUID NONLINEAR FLOW IN LOW-PERMEABILITY POROUS MEDIA WITH APPLICATIONS

A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.

BUBBLE COLLAPSE IN SOLID-LIQUID TWO-PHASE FLUID

Equations of motion for bubble collapse in solid-liquid two-phase fluid have been derived, in which the resistance coupling effects between the liquid and solid particles have been considered. The motion of particles during the bubble collapse and the effects of particles on bubble collapse have been calculated and discussed. Qualitative relations between the concentration and the size of the particles and the rate of bubble collapse have been obtaind.

Introduction to investigations of the negative corona and EHD flow in gaseous two-phase fluids

Research interests have recently been directed towards electrical discharges in multi-phase environments.Natural electrical discharges,such as lightning and coronas,occur in the Earth＇s atmosphere,which is actually a mixture of gaseous phase（air） and suspended solid and liquid particulate matters（PMs）.An example of an anthropogenic gaseous multi-phase environment is the flow of flue gas through electrostatic precipitators（ESPs）,which are generally regarded as a mixture of a post-combustion gas with solid PM and microdroplets suspended in it.Electrical discharges in multi-phase environments,the knowledge of which is scarce,are becoming an attractive research subject,offering a wide variety of possible discharges and multi-phase environments to be studied.This paper is an introduction to electrical discharges in multi-phase environments.It is focused on DC negative coronas and accompanying electrohydrodynamic（EHD） flows in a gaseous two-phase fluid formed by air（a gaseous phase） and solid PM（a solid phase）,run under laboratory conditions.The introduction is based on a review of the relevant literature.Two cases will be considered：the first case is of a gaseous two-phase fluid,initially motionless in a closed chamber before being subjected to a negative corona（with the needle-toplate electrode arrangement）,which afterwards induces an EHD flow in the chamber,and the second,of a gaseous two-phase fluid flowing transversely with respect to the needle-to-plate electrode axis along a chamber with a corona discharge running between the electrodes.This review-based introductory paper should be of interest to theoretical researchers and modellers in the field of negative corona discharges in single-or two-phase fluids,and for engineers who work on designing EHD devices（such as ESPs,EHD pumps,and smoke detectors）.

Development of a 100 MPa water–gas two-phase fluid pressurization device

This paper introduces a 100 MPa water gas twophase fluid pressurization device.The device can provide 100 MPa gas pressure and 200 MPa liquid pressure for small volume(<20 mL)high-pressure experimental devices.This device can make the pressure control independent of the temperature control without changing the material components of the system.The resolution of this device in adjusting the pressure is±0.2 MPa in the process of boosting and depressurizing.This pressure boosting device generates very little vibration during work and it can be used in experiments with strict requirements on vibration.As a thermodynamic parameter,pressure has a great influence on matter.In the field of experimental geochemistry,pressure is not only an experimental method and an extreme condition but an important physical parameter independent of temperature and chemical composition.

Melting Flow Analyzation of Radiative Riga Plate Two-Phase Nano-Fluid Across Non-Flatness Plane with Chemical Reaction

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.

Characteristics of lubrication at nanoscale in two-phase fluid system

Thin film lubrication (TFL) is a condition in which the lubricating features between two surfaces in relative motion are determined by the combination of the properties of the surfaces and the lubricant and viscosity of the lubricant. The effects imposed by couple stress on lubrication characteristics cannot be disregarded in this regime where the ordered molecules dominate the fluid field. There are different tensor measures and constitutive equations in this case other than Newtonian case. The lubrication of two-phase (solid phase and liquid phase) fluid is investigated in this paper. The existence of couple stress will enhance the lubricant viscosity and hence increase the film thickness and improve the load-carrying capability. Size-dependent effects can be seen in the lubrication with couple stress, and the thinner the lubricating film is, the more obvious the effect will be.

Simulation and Analysis on the Two-Phase Flow Fields in a Rotating-Stream-Tray Absorber by Using Computational Fluid Dynamics

The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, and a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall; the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.

Three-dimensional Computational Fluid Dynamics Modeling of Two-phase Flow in a Structured Packing Column

Characterizing the complex two-phase hydrodynamics in structured packed columns requires a power- ful modeling tool. The traditional two-dimensional model exhibits limitations when one attempts to model the de- tailed two-phase flow inside the columns. The present paper presents a three-dimensional computational fluid dy- namics （CFD） model to simulate the two-phase flow in a representative unit of the column. The unit consists of an CFD calculations on column packed with Flexipak 1Y were implemented within the volume of fluid （VOF） mathe- matical framework. The CFD model was validated by comparing the calculated thickness of liquid film with the available experimental data. Special attention was given to quantitative analysis of the effects of gravity on the hy- drodynamics. Fluctuations in the liquid mass flow rate and the calculated pressure drop loss were found to be quali- tatively in agreement with the experimental observations.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

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.

Preliminary Study on Corrosion Behaviour of Carbon Steel in Oil–Water Two-Phase Fluids

The corrosion behavior of type 1045（AISI） carbon steel（CS） in oil-water fluids was investigated by considering the surface wetting state as the breakthrough point.The surface water wetting percentage（SWWP）and corrosion weight loss of CS in different oil-water fluids were measured.The morphology and composition of the corrosion films were detected via scanning electron microscopy,electron-probe facescanning technology and Fourier-transform infrared spectroscopy.The results indicated that the corrosion weight loss of CS in the oil-water fluids corresponded with SWWP and that the corrosion process was influenced by the participation of oil.

Characteristics of high-sulfate wastewater treatment by two-phase anaerobic digestion process with Jet-loop anaerobic fluidized bed

A new anaerobic reactor, Jet-loop anaerobic fluidized bed （JLAFB）, was designed for treating high-sulfate wastewater. The treatment characteristics, including the effect of influent COD/SO42 ratio and alkalinity and sulfide inhibition in reactors, were discussed for a JLAFB and a general anaerobic fiuidized bed （AFB） reactor used as sulfate-reducing phase and methane-producing phase, respectively, in two-phase anaerobic digestion process. The formation of granules in the two reactors was also examined. The results indicated that COD and sulfate removal had different demand of influent COD/SO4^2- ratios. When total COD removal was up to 85%, the ratio was only required up to 1.2, whereas, total sulfate removal up to 95% required it exceeding 3.0. The alkalinity in the two reactors increased linearly with the growth of influent alkalinity. Moreover, the change of influent alkalinity had no significant effect on pH and volatile fatty acids （VFA） in the two reactors. Influent alkalinity kept at 400-500 mg/L could meet the requirement of the treating process. The JLAFB reactor had great advantage in avoiding sulfide and free-H2S accumulation and toxicity inhibition on microorganisms. When sulfate loading rate was up to 8. 1 kg/（m^3.d）, the sulfide and free-H2S concentrations in JLAFB reactor were 58.6 and 49.7 mg/L, respectively. Furthermore, the granules, with offwhite color, ellipse shape and diameters of 1.0-3.0 mm, could be developed in JLAFB reactor. In granules, different groups of bacteria were distributed in different layers, and some inorganic metal compounds such as Fe, Ca, Mg etc. were found.

Fluid−Structure Interaction of Two-Phase Flow Passing Through 90° Pipe Bend Under Slug Pattern Conditions

Numerical simulations of evolution characteristics of slug flow across a 90°pipe bend have been carried out to study the fluid−structure interaction response induced by internal slug flow.The two-phase flow patterns and turbulence were modelled by using the volume of fluid(VOF)model and the Realizable k−εturbulence model respectively.Firstly,validation of the CFD model was carried out and the desirable results were obtained.The different flow patterns and the time-average mean void fraction was coincident with the reported experimental data.Simulations of different cases of slug flow have been carried out to show the effects of superficial gas and liquid velocity on the evolution characteristics of slug flow.Then,a one-way coupled fluid-structure interaction framework was established to investigate the slug flow interaction with a 90°pipe bend under various superficial liquid and gas velocities.It was found that the maximum total deformation and equivalent stress increased with the increasing superficial gas velocity,while decreased with the increasing superficial liquid velocity.In addition,the total deformation and equivalent stress has obvious periodic fluctuation.Furthermore,the distribution position of maximum deformation and stress was related to the evolution of slug flow.With the increasing superficial gas velocity,the maximum total deformation was mainly located at the 90°pipe bend.But as the superficial liquid velocity increases,the maximum total deformation was mainly located in the horizontal pipe section.Consequently,the slug flow with higher superficial gas velocity will induce more serious cyclical impact on the 90°pipe bend.

A two-phase type-curve method with multiscale fluid transport mechanisms in hydraulically fractured shale reservoirs

The quantitative understanding of hydraulic fracture(HF)properties guides accurate production forecasts and reserve estimation.Type curve is a powerful technique to characterize HF and reservoir properties from flowback and long-term production data.However,two-phase flow of water and hydrocarbon after an HF stimulation together with the complex transport mechanisms in shale nanopores exacerbate the nonlinearity of the transport equation,causing errors in type-curve analysis.Accordingly,we propose a new two-phase type-curve method to estimate HF properties,such as HF volume and permeability of fracture,through the analysis of flowback data of multi-fractured shale wells.The proposed type curve is based on a semianalytical solution that couples the two-phase flow from the matrix with the flow in HF by incorporating matrix influx,slippage effect,stress dependence,and the spatial variation of fluid properties in inorganic and organic pores.For the first time,multiple fluid transport mechanisms are considered into two-phase type-curve analysis for shale reservoirs.We analyze the flowback data from a multi-fractured horizontal well in a shale gas reservoir to verify the field application of the proposed method.The results show that the fracture properties calculated by the type-curve method are in good agreement with the long-time production data.

Mathematical model of micropolar fluid in two-phase immiscible fluid flow through porous channel

This paper is concerned with the flow of two immiscible fluids through a porous horizontal channel. The fluid in the upper region is the micropolar fluid/the Eringen fluid, and the fluid in the lower region is the Newtonian viscous fluid. The flow is driven by a constant pressure gradient. The presence of micropolar fluids introduces additional rotational parameters. Also, the porous material considered in both regions has two different permeabilities. A direct method is used to obtain the analytical solu- tion of the concerned problem. In the present problem, the effects of the couple stress, the micropolarity parameter, the viscosity ratio, and the permeability on the velocity profile and the microrotational velocity are discussed. It is found that all the physical parameters play an important role in controlling the translational velocity profile and the microrotational velocity. In addition, numerical values of the different flow parameters are computed. The effects of the different flow parameters on the flow rate and the wall shear stress are also discussed graphically.

Evaluation of frictional pressure drop correlations for air-water and air-oil two-phase flow in pipeline-riser system

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.

Numerical Simulations of Wave Impact Forces on the Open-Type Sea Access Road Using A Two-Phase SPH Model

A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analyses presented by Chen et al.(2022),in which the sea access roads are semi-immersed.In this new configuration,the sea access roads are placed above the still water level,therefore the presence of the air phase becomes a relevant issue in the determination of the wave forces acting on the structures.Indeed,the comparison of wave forces on the open-type sea access roads obtained from the single and two-phase SPH models with the experimental results shows that the latter are in much better agreement.So in the numerical simulations,a two-phaseδ-SPH model is adopted to investigate the dynamical problems.Based on the numerical results,the maximum horizontal and uplifting wave forces acting on the sea access roads are analyzed by considering different wave conditions and geometries of the structures.In particular,the presence of the girder is analyzed and the differences in the wave forces due to the air cushion effects which are created below the structure are highlighted.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

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.

Simulation of Two-Phase Flowback Phenomena in Shale Gas Wells

The gas-water two-phaseflow occurring as a result of fracturingfluidflowback phenomena is known to impact significantly the productivity of shale gas well.In this work,this two-phaseflow has been simulated in the framework of a hybrid approach partially relying on the embedded discrete fracture model(EDFM).This model assumes the region outside the stimulated reservoir volume(SRV)as a single-medium while the SRV region itself is described using a double-medium strategy which can account for thefluid exchange between the matrix and the micro-fractures.The shale gas adsorption,desorption,diffusion,gas slippage effect,fracture stress sensitivity,and capillary imbibition have been considered.The shale gas production,pore pressure distribution and water saturation distribution in the reservoir have been simulated.The influences of hydraulic fracture geometry and nonorthogonal hydraulic fractures on gas production have been determined and discussed accordingly.The simulation results show that the daily gas production has an upward and downward trend due to the presence of a large amount of fracturingfluid in the reservoir around the hydraulic fracture.The smaller the angle between the hydraulic fracture and the wellbore,the faster the daily production of shale gas wells decreases,and the lower the cumulative production.Nonplanar fractures can increase the control volume of hydraulic fractures and improve the production of shale gas wells.

Comparative numerical study of single and two-phase models of nanofluid heat transfer in wavy channel

The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） prediction is used for heat transfer and flow prediction of the single phase and three different two-phase models （mixture, volume of fluid （VOF）, and Eulerian）. The heat transfer coefficient, temperature, and velocity distributions are investigated. The results show that the differences between the temperature fie].d in the single phase and two-phase models are greater than those in the hydrodynamic tleld. Also, it is found that the heat transfer coefficient predicted by the single phase model is enhanced by increasing the volume fraction of nanoparticles for all Reynolds numbers; while for the two-phase models, when the Reynolds number is low, increasing the volume fraction of nanoparticles will enhance the heat transfer coefficient in the front and the middle of the wavy channel, but gradually decrease along the wavy channel.