Gas Void Research Articles

Co-Attention-Based Cross-Stitch Network for Parameter Prediction of Two-Phase Flow

Gas-liquid two-phase flow is pervasive in industrial processes, and accurate measurement of its parameters can effectively improve industrial efficiency. Void fraction, a key parameter in gas-liquid two-phase flow, plays a crucial role in mixture density determination and flow structure analysis. In this paper, a novel multi-task learning method is proposed to predict the gas void fraction. Firstly, we conduct vertical upward gas-liquid two-phase flow experiments to obtain fluid signals relying on the four-sector distributed conductance sensor. Then, we pair the measurement of gas void fraction and its classification and design a Co-Attention based Cross-Stitch Network (CACSNet) to process these two closely related tasks simultaneously. In CACSNet, two parallel embedding branches are adopted to extract task-specific features. There are three feature extractors in each branch to capture in-depth representations of temporal and spatial. To learn an optimal combination of these task-specific features and generate shared representations, we propose the Cross-stitch Embedding with Attention Modules (CEAM), which is incorporated the Channel Co-Attention Module (CCAM) and the Temporal Co-Attention Module (TCAM). Finally, we compare our CACSNet with single-task baseline and other competitive methods, which proves that our framework achieves superior performance in gas void fraction prediction. We also conduct detailed model ablation and parameter analysis to illustrate the efficiency of our proposed structure.

Experimental investigations of gas–liquid two-phase flow in a horizontal mini pipe: Flow regime, void friction, and frictional pressure drops

Experimental studies of air–water two-phase flows in a mini tube with an inner diameter of 3.12 mm have been conducted in terms of the flow pattern, void fraction, and pressure gradients. The fluid velocities range from 0.065 to 21.78 m/s and from 0.109 to 1.835 m/s, respectively. A right-angle prism is applied to capture images simultaneously from two perpendicular directions. Three-dimensional gas–liquid interfaces are reconstructed with the obtained images from two directions. The gas void fraction is investigated using this method, which is experimentally validated with the quick closing valves method. Three flow patterns are obtained in the present work, which is bubbly, slug, and annular flows. The Probability Density Functions of the cross-sectional void fraction show that the intermittent flow has two peaks and the annular flow has only one peak. Volumetric void fraction is affected by the flow pattern and the flow rate. However, at very high gas velocities, volumetric void fraction is independent of the liquid rate. A new frictional pressure drop correlation is proposed based on the experimental data.

Simulation study on multiphase flow considered seawater salinity in wellbore for deep-water bearing hydrate formation drilling

The gas-liquid-solid multiphase flow caused by decomposition of hydrate cuttings is a new problem facing the drilling of marine hydrate formations. In this paper, a gas-liquid-solid non-isothermal transient flow model is proposed considering the coupling relationships between hydrate decomposition, and heat and mass transfer in multiphase flow. The accuracy of the model in predicting gas void fraction, friction pressure and wellbore temperature is verified by using indoor experimental data and field measured data. Using this model, the influence of seawater salinity on the gas-liquid-solid flow behaviors is investigated. Numerical simulation results show that the decomposition area and decomposition rate of hydrate in the wellbore gradually increase with the increase of seawater salinity. Moreover, the higher the salinity of seawater, the greater the reduction in bottomhole pressure and the higher the decomposed gas void fraction reaching the wellhead. Increasing wellhead backpressure and lowering inlet temperature of seawater can effectively control the amount of hydrate decomposed in the wellbore. To ensure wellhead safety and reduce the rate of hydrate decomposition, the wellhead backpressure should be greater than 3.5 MPa. Compared with the wellhead backpressure of 2 MPa, the gas void fraction at the wellhead decreases 1.75 times when the wellhead backpressure is 5 MPa. When the seawater inlet temperature is in the range of 15 °C–30 °C, the seawater salinity control in the range of 8‰∼19‰ can manage the maximum gas phase volume fraction at the wellhead within 10%. Compared with the inlet temperature of 15 °C, the hydrate decomposition rate at the wellhead is 1.6 times higher when the inlet temperature is 30 °C.

Magnetic field modulated dynamics of partial discharges in defects of high voltage insulating materials

This paper presents original measurement methodology and detection approach to determine the influence of the magnetic field on the partial discharge (PD) dynamics. The application areas refer to insulation systems of both grid and industrial network devices and emerging segments such as high-speed rail, electric vehicles or more electric aircrafts. Conventionally PD measurements are executed only in electric field, however the interaction of magnetic and electric fields influences the dynamics of PDs. The measurement technique allowed to detect quantitatively the effect of magnetic fields on PDs in two representative arrangements: in gaseous void in dielectric material and in corona point-plane setup. Measurements in both configurations have revealed amplification of PDs intensity. The quantitative comparison of PDs evolution in the magnetic field is a novel aspect shown in this paper. Combination of phase resolved images and time-sequence intensity diagrams obtained at magnetic field switching allowed to visualize and determine quantitatively this impact. This effect is attributed to the extension of the charged particle trajectory and the amplification of the electron energy due to acceleration. Thus, the investigated impact of a magnetic field can be perceived as an additional element that influences the PD dynamics.

An experimental study on the liquid phase properties of a bubble column operated in the homogeneous and in the heterogeneous regimes

We present an experimental study on the liquid phase properties in a bubble column. The reactor, that has a diameter of 40 cm, can be operated in both the homogeneous and the heterogeneous regimes, reaching superficial gas velocities and void fractions up to 23.2 cm/s and 30%, respectively.Our measurements, obtained with a Pavlov tube, have a temporal resolution of 20 Hz, and cover different superficial gas velocities and axial and radial positions. The moderate temporal resolution achieved allowed to quantify the liquid velocity via probability density functions, standard deviations, spectra and spatial and temporal integral scales.Spectral analysis was performed to study the validity of bubble-induced turbulence models. While such models were developed for the homogeneous regime, we show that they can be applied, with some limitations, to the heterogeneous one. In particular, a -3 power law is found in the liquid velocity spectra, with integral scales related to the column’s diameter.

Large Eddy Simulations of Turbulent Gas-Liquid Flows in a Diverging Horizontal Channel Using a Hybrid Multiphase Approach

Abstract Centrifugal pumps conveying gas–liquid flows are often designed based on pure liquid flows due to limited understanding of gas formations within a gas–liquid flow. Computational fluid dynamics (CFD) can provide insights into the flow characteristics, yet standard multiphase models are not able to describe a wide range of spatial scales. In this study, a hybrid multiphase approach for turbulent flows based on an Eulerian-Eulerian solver with volume-of-fluid (VOF) enhancements is applied, which can handle multiple flow regimes (continuous gas/liquid, disperse bubbles) at once. Large eddy simulations (LES) of gas–liquid-flows through a diverging horizontal channel (diffuser) are performed to evaluate the possibilities and limitations of this approach. Three representative cases, each containing different flow characteristics, are computed using a fixed bubble size approach. The results show that the hybrid approach is able to capture all flow characteristics and predict the gas void size and position precisely. However, it can be observed that the performance of the approach significantly depends on the specified bubble size.

The Lunar Breccia Dhofar 1442: Noble Gases, Nitrogen, and Carbon Released by Stepwise Combustion and Crushing

AbstractStepwise crushing and combustion methods were applied to study the KREEP-rich lunar breccia Dhofar 1442, the clastic material of which is cemented by porous matrix. The stepwise crushing released significant amount of gases of extraterrestrial origin from gas voids. Argon, nitrogen, and carbon are simultaneously released by stepwise combustion at 1100°С. The simultaneous high-temperature degassing of these gases, as well as the coincidence of С/N ratio and nitrogen and carbon contents in high-temperature combustion steps with those of crushing indicate that the gas carriers are voids in high-temperature phases (in particular, minerals, glasses), which are decomposed/melted at these temperatures. Helium and neon are released from the same positions at lower temperatures. The isotopic composition of neon obtained by stepwise combustion and crushing corresponds to the composition of fractionated solar wind. The fraction of argon in the first crushing steps is higher than that of any other of studied gases. The 40Ar/36Ar in the trapped lunar argon is ~18, which is not consistent with empirical model implying that 40Аr is implanted from lunar atmosphere (McKay et al., 1986; Eugster et al., 2001; Joy et al., 2011). We believe that the entrapment of volatile elements in gas voids of the meteorite Dhofar 1442 was caused by the redistribution of gases from one structural sites into others during impact events that accompanied the cratering, in particular, leading to the formation of the impact melt breccia Dhofar 1442. The trapped gases of the meteorite Dhofar 1442 contain not only typical volatile components (solar, radiogenic, cosmogenic, re-implanted 40Ar) of lunar breccias, but also nitrogen and carbon formed through the oxidation of organic matter of metamorphosed chondrites, which are present in the breccia. With increasing number of strokes and, correspondingly, a degree of crushing, the elemental ratios change. A slight decrease of 4He/20Ne ratio during crushing is likely related to the different diffusion ability and permeability of helium relative to neon under temperature influence and/or to the heterogeneous distribution of these gases in voids of different size. The 4He/36Ar, 20Ne/36Ar, 14N/36Ar, and 12С/36Ar ratios increase by factors of 10–100 during crushing. This can be explained by the combination of dynamically different processes leading to the argon fractionation relative to other gases and uneven redistribution of gases from different positions in voids of different sizes during impact metamorphism.

Numerical simulation of bubbly jets in crossflow using OpenFOAM

This paper conducted a computational fluid dynamics study of bubbly jets (not bubble plumes due to pure gas injection) in crossflow to explore the hydrodynamics that are still unknown. A three-dimensional model was developed, calibrated, and validated by coupling the Euler–Euler two-fluid model with unsteady Reynolds-averaged Navier–Stokes approach in OpenFOAM. The results showed that the modeled gas void fraction, bubble velocity, water jet centerline trajectory, and jet expansion agree well with the experimental data. The vertical distribution of turbulent kinetic energy evolves from mono-peak to dual-peaks as the jet penetrates farther for the bubbly jet due to the interactions between bubbles and ambient water flow. Water velocity distribution was examined at cross sections of both the air- and water-phases of bubbly jets in crossflow, and counter-rotating vortex pairs can be clearly observed for both phases. Generally, the center-plane maximum concentration decreases in the crossflow direction. Compared to pure water jets, bubbly jets are stretched wider in the vertical direction due to the lift of bubbles, and thus, dilution is larger. Interestingly, the vorticity at water jet cross sections of bubbly jets evolves from two vertical “kidney-shapes” to two axisymmetric “thumb-up-shapes.” Moreover, effects of ambient crossflow on bubbly jet behaviors were systematically examined. As the crossflow velocity increases, the locations of maximum concentration, maximum velocity magnitude, maximum vorticity magnitude, as well as water jet centerline, all tend to be lower for bubbly jets.

Experimental study of the critical sand starting velocity of gas-water-sand flow in an inclined pipe

The purpose of this paper is to study the critical sand starting velocity and transformation law of flow pattern based on gas–water–sand three-phase flow in an inclined pipe. Firstly, the indoor simulation experiment system of gas–water–sand three-phase flow was used to test the conversion law of flow pattern based upon the different gas void fraction. Secondly, the influence of slug bubbles on sand migration was investigated according to distinctive hole deviation angles, gas void fraction and sand concentration. Finally, the critical sand starting velocity was tested based on dissimilar hole deviation angles, gas void fraction, sand concentration and sand particle size, and then the influence of the above-mentioned key parameters on the sand starting velocity was debated based on the force analysis of the sand particles. The experimental results illustrated that when the gas void fraction was less than 5%, it was bubbly flow. When it increased from 5% to 30%, the bubbly flow and slug flow coexisted. When it was between 30% and 50%, the slug flow and agitated flow coexisted. When it reached 50%, it was agitated flow. Providing that the hole deviation angle was 90°, the phenomenon of overall migration and wave-like migration on the surface of sand bed was observed. On the contrary, the phenomenon of rolling and jumping migration was recognized. The critical sand starting velocity was positively correlated with the hole deviation angle and sand particle size, but negatively associated with the gas void fraction and sand concentration. This research can provide a certain reference for sand-starting production in the field of petroleum engineering.

Frictional Pressure Drop for Gas–Liquid Two-Phase Flow in Coiled Tubing

Coiled tubing (CT) is widely used in drilling, workover, completion, fracturing and stimulation in the field of oil and gas exploration and development. During CT operation, the tubing will present a gas–liquid two-phase flow state. The prediction of frictional pressure drop for fluid in the tube is an important part of hydraulic design, and its accuracy directly affects the success of the CT technique. In this study, we analyzed the effects of the gas void fraction, curvature ratio and fluid inlet velocity on frictional pressure drop in CT, numerically. Experimental data verified simulated results. Flow friction sensitivity analysis shows the frictional pressure drop reaches its peak at a gas void fraction of 0.8. The frictional pressure gradient increases with the increase in curvature ratio. As the strength of secondary flow increases with the increase in inlet velocity, the increased trend of gas–liquid two-phase flow friction is aggravated. The correlation of friction factor for gas–liquid two-phase flow in coiled tubing is developed by regression analysis of simulation results. The research results can support high quality CT hydraulics design, through which the success of CT operations can be guaranteed.

Numerical Investigation of Flow Characteristics for Gas–Liquid Two–Phase Flow in Coiled Tubing

Coiled tubing (CT) is widely used for horizontal well fracturing, squeeze cementing, and sand and solid washing in the oil and gas industry. During CT operation, a gas–liquid two-phase flow state appears in the tubing. Due to the secondary flow, this state produces a more extensive flow-friction pressure loss, which limits its application. It is crucial to understand the gas–liquid flow behavior in a spiral tube for frictional pressure drop predictions in the CT technique. In this study, we numerically investigated the velocity distribution and phase distribution of a gas–liquid flow in CT. A comparison of experimental data and simulated results show that the maximum average error is 2.14%, verifying the accuracy of the numerical model. The gas and liquid velocities decrease first and then rise along the axial direction due to the effect of gravity. Due to the difference in the gas and liquid viscosity, i.e., the flow resistance of the gas and liquid is different, the gas–liquid slip velocity ratio is always greater than 1. The liquid velocity exhibits a D-shaped step distribution at different cross-sections of spiral tubing. The secondary-flow intensity, caused by radial velocity, increases along the tubing. Due to the secondary-flow effect, the zone of the maximum cross-section velocity is off-center and closer to the outside of the tube. However, under the combined action of centrifugal force and the density difference between gas and liquid, the variation in the gas void fraction along the tubing is relatively stable. These research results are helpful in understanding the complex flow behavior of gas–liquid two-phase flow in CT.

Study on Distribution Law of Gas Phase and Cavitation in the Pressurization Unit of Helical Axial Flow Multiphase Pump

Due to the irregular change of gas void fraction (GVF) in multiphase pumps, the pressure distribution in the pump is often uneven, which leads to the formation of low-pressure area and thus the occurrence of cavitation. In order to study the gas phase and cavitation distribution in the impeller region of a multiphase pump under different cavitation stages and GVF conditions, this study used numerical calculations as the main method and experimental verification as a secondary method to investigate the cavitation phenomenon in the pump under different stages and GVF conditions. The results showed that at different stages, both the volume fraction and the covering area of the gas phase were reduced to a certain extent with the increase in blade height, and the distribution law of the gas phase on the blade changed with the development of the cavitation stage, especially on the blade surface. At different GVFs, cavitation first occurred at the inlet of the blade SS and then extended along the blade streamline from the inlet to outlet, with the volume fraction and distribution of cavitation gradually increasing and then extending to the blade PS. The results showed that the presence of the gas phase inhibited the development of cavitation in the multiphase pump to some extent, and the cavitation performance of the multiphase pump was better in the presence of the gas phase than in pure water conditions. The results of this study provide a theoretical basis for improving the cavitation performance of multiphase pumps.

Investigation on multi-physical field simulations of blade ECM using vertical flow

Electrochemical machining (ECM), an advanced manufacturing technology, is widely used in aero-engine blades machining. In traditional ECM, the electrolyte flows through the inter-electrode gap (IEG), generating hydrogen bubbles and heat, which affect the conductance and thus influence the machining quality. This paper focuses on the effect of bubble movement on the flow field and the machining quality of ECM. A novel vertical flow mode of electrolyte is proposed according to the bubbles dynamics analysis. Multi-physical fields simulations of blade ECM using vertical and horizontal flows were carried out. With an initial gas void fraction, flow rate, and temperature at the inlet of 0, 19.7 m/s, and 302.65 K, respectively, in both flow modes, the vertical flow reduces the gas void fraction, flow rate, and temperature at the outlet by 2.4%, 0.4 m/s, and 0.6 K, and increases the conductance by 0.47 S/m. Thus, the vertical flow of the electrolyte is beneficial in reducing the gas void fraction and controlling the temperature rise, while enhancing the conductance. Then, the corresponding experiments using a vertical flow were carried out. The maximum machining deviation ranges from 3.4 to 75.6 μm and surface roughness Ra < 0.35 μm. The machining quality is high and the variation observed in the experiments is consistent with the simulation results, the validity and correctness of the simulations are verified. Thus, the vertical flow mode proposed in this paper is appropriate, can be used for other complex structures in ECM.

Numerical Investigation of Gas-Liquid Flow in a Multiphase Pump with Special Emphasis on the Effect of Tip Leakage Vortex on the Gas Flow Pattern

In this paper, the gas-liquid flow is comprehensively analyzed under different inlet gas void fractions, and the effect of tip leakage vortex (TLV) on the gas flow pattern in multiphase pumps is revealed. The results show that the gas flow pattern in an impeller is closely related to the centrifugal force, low-pressure region, and vortex motion. Most gas is present near the hub and suction surface of the blade as well as in the TLV. The two- and three-dimensional spatiotemporal evolution of the gas is presented, and the gas motion during the inception, development, and dissipation of TLV is revealed. It is reflected that the gas volume fraction is the highest at the TLV core and gradually weakens along the radial direction with the vortex core at the center. Additionally, the TLV energy dissipation is closely related to the gas and pressure difference, and strong energy dissipation occurs in the jet-wake flow.

Creation of hollow silica-fiberglass soft ceramics for thermal insulation

Hollow-structured materials show promise in thermal insulation because the shells encapsulating gaseous voids can interrupt heat transport pathways. Here, we present two low-cost routes to fabricate hollow silica nanoshells, via gas-phase and liquid-phase methods. The gas-phase synthesis method generates hollow shells by a droplet surface precipitation mechanism in a flame aerosol reactor. The liquid-phase synthesis route forms hollow shells by removal of a carbon template, which is produced by hydrothermal reaction of glucose. Both approaches (gas- and liquid-phase) provide hollow silica with amorphous structure, low thermal conductivity (0.023 and 0.026 W m−1 K−1), small particle size (442 and 383 nm), thin shell (35 and 36 nm), and low density (0.015 and 0.033 g cm−3). We employed high shear mechanical mixing to fabricate hollow silica-fiberglass composite ceramics. The resulting three-dimensional network provides the ceramics with robust mechanical elasticity and fire-retardancy while maintaining low thermal conductivity, dramatically outperforming an analogous material using commercial silica gel in place of the hollow nanoshells. Our findings provide two practical routes to synthesize hollow silica, either of which can be used to manufacture a class of hollow shell-fiber nanocomposite soft ceramics for energy-saving applications.

Gas–Liquid Interaction Characteristics in a Multiphase Pump under Different Working Conditions

In this study, we analyze gas–liquid interaction characteristics using a heterogeneous two-fluid model to investigate the influence of interphase force on multiphase pump performance. Two-phase transport platforms are used in oil and gas development to eliminate the need for separation equipment and reduce costs. Full-channel numerical calculations were conducted for an axial-flow multiphase pump based on different inlet gas void fractions (IGVFs) and flow rates. The results indicate that the interaction force of each phase is relatively large in the rotor–stator interference region, and the drag, lift, virtual mass, and turbulent dispersion forces increase with an increase in IGVF or when deviating from the design condition (Q = 50 m3/h). The interphase forces (resistance, lift, virtual mass force, and turbulent dispersion) increase considerably in the impeller passage and minimally in the guide blade passage. Under the conditions of small and high flows, the force of each phase changes considerably in the impeller and diffuser passages, respectively. Furthermore, the turbulent kinetic energy in the flow passage corresponds to the change trend of the interphase force, indicating that the interphase force causes energy loss inside a multiphase pump. These results provide essential information for the optimization of the hydraulic design of multiphase pumps.

Influence factors and prediction model of enstrophy dissipation from the tip leakage vortex in a multiphase pump.

In a multiphase pump, tip clearance is the required distance between the blade tip and the pump body wall of the impeller, forming tip leakage vortex (TLV), causing unstable flow and energy dissipation. In the present work, the enstrophy dissipation theory is innovatively applied to quantitatively study the energy dissipation of the TLV. The flow rate, tip clearance, and inlet gas void fraction (IGVF) play a crucial role in affecting the enstrophy dissipation of the TLV. The results show that increasing flow rate, tip clearance, and IGVF significantly exacerbate the TLV pattern and raise the TLV scale, which gradually raises volume enstrophy dissipation and decreases wall enstrophy dissipation. As the flow rate increases, the separation angle between the primary TLV trajectory and the blade gradually decreases, and widely dispersing the enstrophy dissipation near the shroud. However, as the tip clearance increases, the tip separated vortex scale increases and extends to the suction surface, raising the velocity gradient. Besides, as the IGVF increases, the secondary TLV develops from a continuous sheet vortex to a scattered strip vortex, increasing the significantly increasing the enstrophy dissipation. Considering the flow rate, tip clearance, and IGVF as independent variables, simple and multiple nonlinear regression models have the ability to predict the enstrophy dissipation of the TLV accurately.

Characterization of the gassy sediment layer in shallow water using an acoustical method: Lake Kinneret as a case study

AbstractRemote characterization and parameterization of gassy sediments have significant environmental importance for quantifying the global methane budget and assessing its impact on climate change. Acoustic techniques that have been developed hold advantages over direct sediment sampling (e.g., using pressurized and frozen cores), as they permit comparatively quick and cost‐effective assessments over the large bottom areas. This paper proposes a non‐invasive acoustic method that allows simultaneous assessments of the free gas content () and the thickness () of a gassy layer in the aquatic surface sediments. The method is based on amplitude measurements and frequency analysis of the bottom reflection coefficient in the wide frequency band (300–3500 Hz). The spatial variability of and in freshwater Lake Kinneret (Israel) is studied, where the upper sedimentary layer is characterized by a high organic matter content, high methane production rates, and a large . The assessed values of and varied from 0.1% to 0.6%, and from 20 to 40 cm, respectively, depending on the location of measurements. These results are in reasonable agreement with gas void fractions measured directly in frozen sediment cores, where the depth‐averaged varied from to . The suggested methodology should have considerable practical implementation for remote spatiotemporal monitoring of shallow gassy sediments in aquatic ecosystems.

Investigation on the performance of a helico-axial multiphase pump under slug flow

The helico-axial multiphase pump is often used for gas-liquid mixture transportation in offshore oilfields, and slug flow is the main reason for the unstable operation of the pump. Aimed for slug flow condition, a self-designed three-stage multiphase pump is set to the object to perform unsteady simulations and fluid-structure interaction calculations, and the inlet gas void fraction (IGVF) is set from 20% to 80%. The results show that affected by the flow from the slug, the gas-liquid two-phase flow pattern in the multiphase pump changes sharply, resulting in severe fluctuations in the differential pressure, spindle torque and deformation of the multiphase pump. The gas-phase enters the first-impeller along the suction blade surface when affected by Taylor bubbles, while the second and third-stage impellers gas-phases are in the form of small air masses flow into the impeller along the pressure blade surface. The deformation trend of impeller torque, differential pressure and the main pump spindle is similar to that of trigonometric function, while the fluctuation of torque is more intense, and the shape variable of spindle increases with the inflow of liquid plug, and the maximum deformation amount increases by 10.9% at high GVF relative to IGVF.

Numerical study on gas–liquid two phase flow characteristic of multistage electrical submersible pump by using a novel multiple-size group (MUSIG) model

Electrical submersible pumps (ESPs) face enormous challenges in the petroleum industry while handling gas–liquid two-phase flow. The major difficulty is caused by the accumulation of gas bubbles inside ESP-impellers, which results in mild to severe degradation in pump performance. Therefore, to analyze the influence of gas entrainment and bubble size, a combination of experimental and numerical analysis is performed on a five-stage mixed-flow ESP in the present study. The experiments are first conducted to analyze the performance of ESP under pure water conditions at different rotating speeds, followed by the gas–liquid two-phase flow experiments that are performed at constant rotating speed (1475 r/min) and for a wide range of inlet gas void fractions (IGVFs). For numerical calculations, a novel multiple-size group (MUSIG) model is applied in ANSYS CFX to analyze the performance and different flow patterns in ESP in different IGVFs and understand the coalescence and breakup phenomena of gas bubbles in the impeller flow passage. The simulation results from the MUSIG model are compared with the Euler–Euler two-fluid model and test results. The MUSIG model can more accurately predict the changes in the performance and internal flow-field of ESP under two-phase flow conditions. Moreover, when the MUSIG model is used to calculate the two-phase flow of the ESP, the first-stage impeller has a higher head than other stages because the flow inside the second and other stages is affected by the disoriented flow coming from the first-stage diffuser and other return channels. Furthermore, this study gives an insight into the comprehensive application of the novel MUSIG model for complex turbo-machine designs such as ESP.