Search results

The purpose of this study is to design numerical simulations of bubbly flow around a cylinder to better understand the characteristics of flow around a rigid obstacle.

The bubbly flow around a circular cylinder was numerically simulated using a semi-Lagrangian–Lagrangian method composed of a vortex-in-cell method for the liquid phase and a Lagrangian description of the gas phase. Additionally, a penalization method was applied to account for the cylinder inside the flow. The slip condition of the bubbles on the cylinder’s surface was enforced, and the outflow conditions were applied to the liquid flow at the far field.

The simulation clarified the characteristics of a bubbly flow around a circular cylinder. The bubbles were shown to move around and separate from both sides of the cylinder, because of entrainment by the liquid shear layers. Once the bubbly flow fully developed, the bubbles distributed into groups and were dispersed downstream of the cylinder. A three-dimensional vortex structure of various scales was also shown to form downstream, whereas a quasi-stable two-dimensional vortex structure was observed upstream. Overall, the proposed method captured the characteristics of a bubbly flow around a cylinder well.

A semi-Lagrangian–Lagrangian approach was applied to simulate a bubbly flow around a circular cylinder. The simulations provided the detail features of these flow phenomena.

The purpose of this paper is to propose a prediction method of gas-liquid two-phase flow patterns and reveal the flow characteristics in the suction chamber of a centrifugal pump.

A transparent model pump was experimentally studied, and the gas-liquid two-phase flow in the pump was numerically simulated based on the Eulerian–Eulerian heterogeneous flow model. The numerical simulation method was verified from three aspects: the flow pattern in the suction chamber, the gas spiral length and the external characteristics of the pump. The two-phase flow in the suction chamber was studied in detail by using the numerical simulation method.

There are up to eight flow patterns in the suction chamber. However, at a certain rotational speed, only six flow patterns are observed at the most. At some rotational speeds, only four flow patterns appear. The gas spiral length has little relationship with the gas flow rate. It decreases with the increase of the liquid flow rate and increases with the increase of the rotational speed. The spiral flow greatly increases the turbulence intensity in the suction chamber.

A method for predicting the flow pattern was proposed. Eight flow patterns in the suction chamber were identified. The mechanism of gas-liquid two-phase flow in the suction chamber was revealed. The research results have reference values for the stable operation of two-phase flow pumps and the optimization of suction chambers.

This study aims to provide discussions of the numerical method and the bubbly flow characteristics of an annular bubble plume.

The bubbles, released from the annulus located at the bottom of the domain, rise owing to buoyant force. These released bubbles have diameters of 0.15–0.25 mm and satisfy the bubble flow rate of 4.1 mm³/s. The evolution of the three-dimensional annular bubble plume is numerically simulated using the semi-Lagrangian–Lagrangian (semi-L–L) approach. The approach is composed of a vortex-in-cell method for the liquid phase and a Lagrangian description of the gas phase.

First, a new phenomenon of fluid dynamics was discovered. The bubbly flow enters a transition state with the meandering motion of the bubble plume after the early stable stage. A vortex structure in the form of vortex rings is formed because of the inhomogeneous bubble distribution and the fluid-surface effects. The vortex structure of the flow deforms as three-dimensionality appears in the flow before the flow fully develops. Second, the superior abilities of the semi-L–L approach to analyze the vortex structure of the flow and supply physical details of bubble dynamics were demonstrated in this investigation.

The semi-L–L approach is applied to the simulation of the gas–liquid two-phase flows.

The purpose of this study is to investigate the monodisperse cavitation of bubbly mixture flow for water and hydrogen mixture flows through a nozzle having a stenosis on the wall.

Two flow regions, namely, quasi-statically stable and quasi-statically unstable increase in the bubble radius, are considered. Different oscillating periods of bubbles in downstream corresponding to various values of Reynolds number are taken into account. The Range–Kutta method is used to tackle nonlinear coupled system of governing equations.

It is observed that for the larger values of Reynolds number, the void fraction at the upstream section, even at small values, yields instabilities at the downstream. Consequently, owing to sudden increase in the velocity, the bubbles strike the wall with high speed that eventually remove the existing stenosis. This process can be considered as an effective cardiac surgery for arteries with semi-blockage.

Original research work and to the best of author’s knowledge, this model is reported for the first time.

The purpose of this paper is to conduct a review of types of tomographic systems that have been widely researched within the past 10 years. Decades of research on non-invasively and non-intrusively visualizing and monitoring gas-liquid multi-phase flow in process plants in making sure that the industrial system has high quality control. Process tomography is a developing measurement technology for industrial flow visualization.

A review of types of tomographic systems that have been widely researched especially in the application of gas-liquid flow within the past 10 years was conducted. The sensor system operating fundamentals and assessment of each tomography technology are discussed and explained in detail.

Potential future research on gas-liquid flow in a conducting vessel using ultrasonic tomography sensor system is addressed.

The authors would like to undertake that the above-mentioned manuscript is original, has not been published elsewhere, accepted for publication elsewhere or under editorial review for publication elsewhere and that my Institute’s Universiti Teknologi Malaysia representative is fully aware of this submission.

Numerical simulations of a multistage multiphase pump at different operating conditions were performed to study the variational characteristics of flow parameters for each impeller. The simulation results were verified against the experimented results. Because of the compressibility of the gas, inlet volume flow rate qi and inlet flow angle ß_i for each impeller decrease gradually from the first to the last stage. The volume flow rate at the entrance of the pump q, rotational speed n and inlet gas volume fraction (IGVF) affect the characteristics of qi and ß_i.

The hydraulic design features of the impellers in the multistage multiphase pump are obtained based on the flow parameter characteristics of the pump. Using the hydraulic setup features, stage-by-stage design of the multistage multiphase pump for a nominal IGVF has been conducted.

The numerical simulation results show that hydraulic loss in impellers of the optimized pump is substantially reduced. Furthermore, the hydraulic efficiency of the optimized pump increases by 3.29 per cent, which verifies the validation of the method of stage-by-stage design.

Under various operating conditions, q_i and ß_i decrease gradually from the first to the fifth stage because of the compressibility of the gas. For this characteristic, the fluid behavior varies at each stage of the pump. As such, it is necessary to design impellers stage by stage in a multistage rotodynamic multiphase pump.

These results will have substantial effect on various practical operations in the industry. For example, in the development of subsea oilfields, the conventional conveying equipment, which contains liquid-phase pumps, compressors and separators, is replaced by multiphase pumps. Multiphase pumps directly transport the mixture of oil, gas and water from subsea oilwells through a single pipeline, which can simplify equipment usage, decrease backpressure of the wellhead and save capital costs.

Characteristics of a multistage multiphase pump under different operating conditions were investigated along with features of the inlet flow parameters for every impeller at each compression stage. Our simulation results have established that the change in the inlet flow parameters of every impeller is mainly because of the compressibility of the gas. The operational parameters q, n and IGVF all affect the characteristics of qi and ß_i. However, the IGVF has the most prominent effect. Lower values of IGVF have an insignificant effect on the gas compressibility. Higher values of IGVF have a significant effect on the gas compressibility. All these characteristics affect the hydraulic design of the impellers for a multistage multiphase pump. In addition, the machining precision should also be considered. Considering all these factors, when IGVF is lower than 10 per cent, all the impellers in the pump can be designed uniformly. When IGVF varies from 10 to 30 per cent, the first two stages should be designed separately, and the latter stages are uniform starting with the second stage. When IGVF varies from 30 to 50 per cent, the first three stages should be designed separately, and the latter stages are going to be similar to the third stage. An additional increase in IGVF results in degeneration of the differential pressure of the pump, which will reduce the compressibility of the gas. As such, it can be deduced that only the first three stages should be designed separately, and the latter stages will be similar to the third stage. In addition, for the pump working under a lower volume flow rate than 25 m³/h, the first three stages should be designed individually while keeping the geometrical structure of the subsequent stages the same as the third stage.

The modeling of interphase forces plays a significant role in the numerical simulation of gas–liquid flow in a rotodynamic multiphase pump, which deserves detailed study.

Numerical analysis is conducted to estimate the influence of interphase forces, including drag force, lift force, virtual mass force, wall lubrication force and turbulent dispersion force.

The results show that the magnitude of the interphase forces can be sorted by: drag force > virtual mass force > lift force > turbulent dispersion force > wall lubrication force. The relations between interphase forces and velocity difference of gas–liquid flow and also the interphase forces and gas volume fraction are revealed. The distribution characteristics of interphase forces in the passages from impeller inlet to diffuser outlet are illustrated and analyzed. According to the results, apart from the drag force, the virtual mass force, lift force and turbulent dispersion force are required, whereas wall lubrication force can be neglected for numerical simulation of gas–liquid flow in a rotodynamic multiphase pump. Compared with the conventional numerical method which considers drag force only, the relative errors of predicted pressure rise and efficiency based on the proposed numerical method in account of four major forces can be reduced by 4.95 per cent and 3.00 per cent, respectively.

The numerical analysis reveals the magnitude and distribution of interphase forces inside multiphase pump, which is meaningful for the simulation and design of multiphase pump.

The purpose of this paper is to compute the pressure drop through sudden expansions and contractions for two‐phase flow of oil/water emulsions.

Two‐phase computational fluid dynamics (CFD) calculations, using Eulerian–Eulerian model, are employed to calculate the velocity profiles and pressure drops across sudden expansions and contractions. The pressure losses are determined by extrapolating the computed pressure profiles upstream and downstream of the expansion/contraction. The oil concentration is varied over a wide range of 0‐97.3 percent by volume. The flow field is assumed to be axisymmetric and solved in two dimensions. The two‐dimensional equations of mass, momentum, volume fraction and turbulent quantities along with the boundary conditions have been integrated over a control volume and the subsequent equations have been discretized over the control volume using a finite volume technique to yield algebraic equations which are solved in an iterative manner for each time step. The realizable per phase k‐ ε turbulent model is considered to update the fluid viscosity with iterations and capture the individual turbulence in both the phases.

The contraction and expansion loss coefficients are obtained from the pressure loss and velocity data for different concentrations of oil–water emulsions. The loss coefficients for the emulsions are found to be independent of the concentration and type of emulsions. The numerical results are validated against experimental data from the literature and are found to be in good agreement.

The present computation could not use the surface tension forces and the energy equation due to huge computing time requirement.

The present computation could compute realistically the two‐phase pressure drop through sudden expansions and contractions by using a two‐phase Eulerian model and hence this model can be effectively used for industrial applications where two‐phase flow comes into picture.

The original contribution of the paper is in the use of the state‐of‐the‐art Eulerian two‐phase flow model to predict the velocity profile and pressure drop through industrial piping systems.

Interphase forces between the gas and liquid phases determine many phenomena in bubbly flow. For the interphase forces in a multiphase rotodynamic pump, the magnitude analysis was carried out within the framework of two-fluid model. The purpose of this paper is to clarify the relative importance of various interphase forces on the mixed transport process, and the findings herein will be a base for the future study on the mechanism of the gas blockage phenomenon, which is the most challenging issue for such pumps.

Four types of interphase forces, i.e. drag force, lift force, virtual mass force and turbulent dispersion force (TDF) were taken into account. By comparing with the experiment in the respect of the head performance, the effectiveness of the numerical model was validated. In conditions of different inlet gas void fractions, bubble diameters and rotational speeds, the magnitude analyses were made for the interphase forces.

The results demonstrate that the TDF can be neglected in the running of the multiphase rotodynamic pump; the drag force is dominant in the impeller region and the outlet extended region. The sensitivity analyses of the bubble diameter and the rotational speed were also performed. It is found that larger bubble size is accompanied by smaller predicted drag but larger predicted lift and virtual mass, while the increase of the rotational speed can raise all the interphase forces mentioned above.

This paper has revealed the magnitude information and the relative importance of the interphase forces in a multiphase rotodynamic pump.

The purpose of this paper is to focus on the application of the Trefftz method to the calculation of the two-dimensional (2D) temperature field in the boiling refrigerant flow through an asymmetrically heated vertical minichannel with a rectangular cross-section. The considerations were limited to determining the temperature of the continuous phase – liquid for bubbly and bubbly-slug flow. The numerical solution found with the Trefftz methods was compared with the simplified solution. For nucleate boiling, heat transfer coefficient at the heating foil – liquid contact was determined.

The Trefftz method was used to determine 2D temperature distributions for the glass pane, the heating foil and the boiling liquid. The temperature fields were approximated by the sum of the particular solution and the linear combination of suitable Trefftz functions. Coefficients of linear combination were computed using experimental data, including heating foil temperature measurements obtained with the liquid-crystal method and experimentally determined void fraction. The computations were based on the Trefftz method supplemented with the adjustment calculus.

The way of solving direct and inverse problems of heat conduction in solid bodies (isolating glass, heating foil) and in liquids (boiling refrigerant flowing through the minichannel) was presented. For the first time, both 2D temperature fields for the heating foil and the boiling liquid were calculated while simultaneously using the Trefftz method. The known temperature values of the foil and liquid allowed the calculation of the heat transfer coefficient and the heat flux at the heating foil-liquid contact. Adjustment calculus implemented into the Trefftz method was used to smooth the measurement data and to reduce their errors.

The approach proposed in the paper can be applied to determining 2D temperature field, heat flux and heat transfer coefficient in direct and inverse problems concerning two-phase flowing miniature compact heat exchangers.

The paper presents a novel implementation of the Trefftz method to simultaneous solving an inverse problem in the heating foil and the contacting flowing liquid.