Two-phase Flow Analysis Research Articles

RGB Mapping: A Dynamic Approach for Flow Pattern Identification and Classification

Thermal-hydraulic modeling and simulation are heavily dependent on the knowledge of flow regimes, which sort multiphase flow patterns into classifications with static boundaries that display characteristic features of flow patterns and phase distributions. Flow regime maps heretofore have been primarily developed mainly for vertical or horizontal pipe orientations by capturing images of flows at varying conditions and assigning them to different flow regimes based on visualization studies. As such, this method is inherently subjective. Furthermore, classifying two-phase flows into regimes with static boundaries is not physical, as two-phase flows undergo gradual changes as they transition from one regime to another. Therefore, such a static approach is not only unphysical but also inaccurate, resulting in the erroneous analysis of two-phase flows. Additionally, subjective observations at such flow conditions are prone to disagreement and inconsistencies in labeling. These difficulties are exacerbated in inclined two-phase flows. In this study, a novel two-phase flow pattern map employing the dynamic red, green, and blue (RGB) color model is presented. The model contains the embedded information of three normalized time-averaged statistics from signals obtained by a dual-ring impedance meter. A database of impedance meter signals from a 25.4-mm-inner-diameter adiabatic air-water inclinable test facility is used to compute statistics for each experimental flow condition. Multiple statistical measures are evaluated, and a principal component analysis is performed for feature selection. By mapping each parameter to intensities of red, green, and blue, a dynamic RGB flow pattern map can be produced in which groups of similar colors are representative of various flow regimes and the dynamic transition characteristics between the flow regimes are characterized physically. Therefore, this new method makes viable the observation of how flow patterns gradually change with respect to angle, as well as the fuzzification or quantification of the uncertainty surrounding experimental transition boundaries.

Analysis of flow boiling incipience models in computational fluid dynamics

Recent advances in the computational fluid dynamics and heat transfer (CFD/HT) modeling of flow boiling have enabled a more detailed analysis and understanding of thermal management in relatively complex microsystems. Most advanced modeling approaches allow the transient analysis of two-phase flow and phase-change processes, providing the capability to post-process in detail the temperature, phase and pressure distributions in microscale applications with the use of moderate computational resources. Although these models have significantly evolved in recent years, there are still several idealizations and assumptions that have yet to be improved or better justified, the boiling incipience treatment being a commonly questioned topic. In the present investigation, a boiling incipience model is incorporated into an existing CFD/HT model for the analysis of flow boiling mechanisms, where the superheat temperature is derived from fundamental thermodynamic relations and compared to the original CFD/HT model that does not account for such effects. Results are compared using a non-trivial case of a silicon micro-cooling layer with variable density of pin fins that has been demonstrated to be an efficient thermal control device in high-power applications. The dielectric fluid HFE-7200 is used as the coolant in flow boiling conditions. The two-phase flow regimes and heat transfer results for both models are compared.

Simulation analysis of gas-solid two-phase flow and reaction in a novel catalytic cracking reactor

A gas-solid two-phase flow reaction model was developed to evaluate the performance of a novel catalytic cracking reactor in the residue-to-chemicals (RTC) process. Coupled TFM, EMMS drag, and 12-lump kinetic models were employed to simulate the flow reaction process. The reaction temperature and product yields at the reactor exit aligned well with real industrial RTC reactor test results, indicating the reliability and effectiveness of the coupling model. Catalyst particles formed an internal circulation structure in the reactor, increasing the instantaneous catalyst-to-oil ratio. The reactor exhibited low axial and radial temperature gradients with higher reaction temperatures, accelerating the catalytic cracking rate and enhancing the selectivity for high-value products. The effects of the catalyst-to-oil ratio, catalyst inlet temperature, and feedstock mass flow rate on the flow reaction process were optimized. Results showed that both the catalyst-to-oil ratio and catalyst inlet temperature influenced the overall catalyst velocity distribution. Within the limit range, increasing the catalyst-to-oil ratio increased feedstock conversion, gasoline, LPG, dry gas, and coke yields, but decreased diesel yield. Conversely, increasing the feedstock mass flow rate showed opposite trends. Different trends in product yields were observed with varying inlet catalyst temperatures, with optimal product distribution at 973.15 K.

Uncertainties in measurements of bubbly flows using phase-detection probes

The analysis of bubbly two-phase flows is challenging due to their turbulent nature and the need for intrusive phase-detection probes. However, accurately characterizing these flows is crucial for safely designing critical infrastructure such as dams and their appurtenant structures. The combination of dual-tip intrusive phase-detection probes with advanced signal processing algorithms enables the assessment of pseudo-instantaneous 1-D velocity time series; for which the limitations are not fully fathomed. In this investigation, we theoretically define four major sources of error, which we quantify using synthetically generated turbulent time series, coupled with the simulated response of a phase-detection probe. Based on the analysis of 1010 simulated bubble trajectories, our findings show that typical high-velocity flows in hydraulic structures hold up to 15% error in the mean velocity estimations and up to 35% error in the turbulence intensity estimations for the most critical conditions, typically occurring in the proximity of the wall. Based on thousands of simulations, our study provides a novel data-driven tool for the estimation of these baseline errors (bias and uncertainties) in real-word phase-detection probe measurements of bubbly flows (air concentrations c<40%).

ARIANT Assessments for CANDU Analysis Using RD-14M Experiments

ARIANT (AlgoRIthm for Analysis of Network Thermalhydraulics) is a Canadian Nuclear Laboratories system thermal-hydraulic code for the modeling and analysis of two-phase flow and heat transfer for pressurized heavy water reactors, light water pressurized water reactors, and advanced reactor applications. This paper presents ARIANT models and simulations of RD-14M experiments, including small-break loss-of-coolant accidents, large-break loss-of-coolant accidents, loss-of-flow accidents, station blackout, and natural circulation, that are representative of accident scenarios in a CANDU reactor. ARIANT predictions of pressures, flow rates, temperatures, and void fractions are compared against the steady-state and transient data over the course of the tests. The results show that ARIANT predicted the key parameters with reasonable accuracy, as well as the overall behavior of the five transient events. These assessments support ARIANT’s applicability to the corresponding CANDU design-basis accidents and demonstrate ARIANT as an alternative to existing system thermal-hydraulic codes for CANDU safety analysis.

Non-similar analysis of two-phase hybrid nano-fluid flow with Cattaneo-Christov heat flux model: a computational study

The purpose of this research elaborates the behaviour of the two-phase hybrid nano-fluid flow of viscous fluid through the stretching sheet. Impacts of heat transfer, for MHD two-dimensional steady flow over stretching surface with Cattaneo-Christov (CC) model are discussed. Bio-convective, thermal radiation and convective conditions are considered. Porous surface and velocity slip effects are also considered. Darcy-Forchheimer relation characterizes porous medium. The governing equations are solved through MATLAB. Plots of temperature and velocity fields are discussed via various estimations of emerging variables. Results findings shows the velocity profile decreases in case of both simple and hybrid nano-fluids for higher estimations of magnetic parameter and local inertia coefficient parameter. Fluid temperature reduces for maximum thermal relaxation parameter and possess higher values against radiation parameter. Microorganism profile reduces bio-convection Lewis number but increases for the magnetic parameter. Prominent results from this study include determining the effective parameters to enhance heat transfer efficiency and defining the circumstances in which the nano-fluid performs better than traditional fluids. These discoveries offer significant perspectives on the utilization of Cattaneo-Christov model with hybrid nano-fluid flow in comprehensive thermal management systems and establish a basis for subsequent practical and theoretical investigations within the domain.

Analysis and optimization of LN2 two-phase flow in CRAFT NNBI cryopump

As an important component of Negative ion based Neutral Beam Injector (NNBI), the cryopump mainly provides a suitable vacuum environment for beam generation and transmission. In the paper, Liquid nitrogen (LN2) pipe structure models of two cryopumps were established for simulation and experimental studies. Thermal analysis of the thermal radiation shielding baffle (LN2 cooling) was carried out by Ansys steady-state thermal analysis software, while Fluent was used to simulate the two-phase flow in the LN2 pipeline, then the pressure drops, temperature, velocity, gas-phase volume fraction, and other parameters of the two pipeline structures were analyzed and compared. The simulation results show that the optimized structure B can reduce the pressure drop by 3 bar, the outlet temperature of structure A is 78.8 K, the outlet temperature of structure B is 79.2 K, the temperature difference is 0.4 K. The outlet velocity increases from 2.067 m/s to 4.947 m/s for Structure A, and from 1.096 m/s to 6.614 m/s for Structure B. The experiment results show that the outlet temperature of structure A is 79.4 K, meanwhile the outlet temperature of structure B is 80.1 K, the optimized structure B can reduce the pressure drop about 3 bar.

Detailed 3D URANS analysis of two-phase flow in an airlift pump

An airlift pump is a vertical tube that utilizes the buoyant effects of a gas to lift a liquid. Unlike a standard mechanical pump, the liquid flow rate through the airlift pump is not directly controlled; rather, it depends on the supplied gas flow rate, the tube length and diameter, and the relative height of the liquid supply free surface (submergence ratio). The present study uses the commercial CFD code ANSYS CFX to model the isothermal, 3D, transient flow in an airlift pump using water and air. The model applies pressure boundary conditions at both ends of the tube and specifies the mass flow rate of air through multiple openings in the side of the tube. The bottom of the tube is an inlet of water only and the outlet is a two-phase flow opening. A time-dependent, homogeneous, VOF two-phase RANS CFD modelling approach is used with the air treated as an ideal gas. This work found that a complete 3D domain was necessary for consistent prediction of the airlift performance and physically realistic two-phase flow structures. Statistical analysis of the two-phase flow structures was applied to characterize airlift pump instability and better understand the physics of the airlift pump.

Flow pattern and pressure drop analysis of two-phase segmented flow driven by a thermomagnetic pump

The two-phase segmented flow is a flow pattern in which two immiscible fluids alternate in a flow channel with distinct segments. In this study, a thermomagnetic pump operated by heat and an external magnetic field was used to drive the flow and the temperature-sensitive magnetic fluid (TSMF) was chosen as the continuous-phase fluid to generate a two-phase segmented flow. The dispersed-phase fluid was selected as air or deionized water driven by a syringe pump. The two-phase segmented flow experiments were conducted to investigate the relation between segment length and the dispersed-continuous phase flow rate ratio at the air flow rate from 0.90 ml/h to 4.49 ml/h for air-TSMF experiment and water flow rate from 0.36 ml/h to 0.72 ml/h for water-TSMF experiment. The ratio was proportional to the bubble/droplet length and inversely proportional to the slug length. The air/water length per channel width can be up to 3.28 and 1.79, respectively, at the highest air/water flow rate in the current experiment. Moreover, the pressure drop was measured to increase with the dispersed-continuous phase flow rate ratio. The maximum pressure drop was observed as 428.2 Pa in the water-TSMF experiment with the highest water flow rate and the maximum pressure drop was 314.7 Pa for the air-TSMF experiment. These findings provide valuable insights for controlling the segment length in the two-phase segmented flow driven by a thermomagnetic pump and it could be used as precise control of liquid/gas fuel delivery in a micro-propulsion system in small satellites.

Numerical simulation of sand erosion in Venturi tube under different flow characteristics and sand content

Abstract In the solid-liquid two-phase flow, there is flow erosion of solid relative to the surrounding wall or other components, which has attracted the attention of many scholars. Wherever flows are involved, erosion tends to be a problem. Especially in rivers with large sediment content, hydraulic machinery including pumps, water turbines, reversible water pump turbines and pipelines in a long time of flow erosion, it is easy to cause the material damage of the flow parts, affect the head and efficiency to a certain extent, may also lead to the increase of the gap between the flow parts resulting in noise and other hazards. However, in different flow conditions, the location of erosion is different, which makes it difficult to prevent erosion. As a common measuring tool, Venturi tube has a certain representative geometric shape. The erosion prediction of Venturi tube by computational fluid dynamics can provide a certain reference for other structures. In this study, the solid-liquid two-phase flow analysis method based on CFD was used to analyze the flow erosion under different velocity and sediment content. Through numerical simulation, it is found that the erosion problem in Venturi tube is closely related to its flow characteristics.

Multiscale entropy applications for complexity analysis of two-phase flow

Multiscale entropy applications for complexity analysis of two-phase flow

A fully coupled compositional wellbore/reservoir model for predicting liquid loading in vertical and inclined gas wells

Liquid loading presents a formidable challenge for mature gas wells, often resulting in substantial economic losses. Traditional research has predominantly centered on the analysis of gas-liquid two-phase flow within the wellbore to predict critical gas velocity or rate, aiding in identifying the onset of liquid loading. This study introduces a fully coupled compositional wellbore-reservoir simulator designed to detect liquid loading in both vertical and inclined gas wells. Leveraging the drift-flux model to evaluate flow pattern transitions, this simulator employs pressure or rate constraints at the wellhead as boundary conditions. It comprehensively captures the flow dynamics in both the wellbore and reservoir, unveiling significant changes in gas production rate, water production rate, gas velocity, flow regime, and the reserved position of the liquid film under liquid-loaded conditions. Moreover, the accumulation of liquid at the bottom hole leads to increased reservoir pressure and gas saturation near the wellbore. The simulator predicts a typical unstable production period, emphasizing its crucial role in implementing effective strategies to mitigate liquid loading. This paper investigates the capability of the coupled wellbore-reservoir model to characterize transient liquid loading phenomena from a systematic perspective. The proposed model can function as a real-time tool for predicting the status of liquid loading in gas wells.

Permeability anisotropy analysis of two-phase flow during hydrate dissociation process

Natural gas hydrate (NGH) is a versatile and energy source. But in the actual extraction process, NGH in the reservoir will be out of phase equilibrium conditions and dissociate. This will induce alterations in permeability characteristics of reservoir and will exert a significant influence on the assessment of hydrate reservoir capacity. In the research, an innovative CT triaxial system was used to obtain the micro-structure of sediment specimen during hydrate dissociation. The anisotropy of absolute permeability (k), capillary pressure (Pc) and relative permeability (kr) were also analyzed using the pore network model (PNM). As demonstrated by the results, k increases gradually in all three directions during the dissociation. Pc gradually decreases. The degree of anisotropy of both Pc and gas-water relative permeability (krg and krw) gradually decreases. The impact of hydrate dissociation on the kr anisotropy during primary gas flooding surpasses that observed during secondary water flooding. As hydrate saturation decreases, the krw increases and the krg decreases in the horizontal directions. The krw does not change significantly and the krg increases in the vertical direction. Along the length of the sample, the distribution of hydrate and pore space exhibits heterogeneity across the sections. This inhomogeneous distribution will lead to anisotropy in permeability.

Spectral analysis of dispersed multiphase flows in the presence of fluid interfaces

Spectral analysis of dispersed two-phase flows is highly desirable to reveal the interplay of the various flow scales, much larger or much smaller than the size of the dispersed bodies. This is a challenging task as the matching conditions at the body interfaces generate singularities in the fields describing the two-phase mixture. The nature of these singularities and their consequences on the spectra are theoretically analyzed for bubble or droplet flows. Results of direct numerical simulations are reported and spatial spectra of the mixture velocity, the flow forces and their power are examined. The regular part of the spectral densities of energy production, dissipation and transfers between scales are separated from their singular part. The resulting spectral energy balance, free of the footprint of the singularities, is found in agreement with coarse-grained simulations where the interfaces are filtered out before solving the Navier–Stokes equations. These results pave the way for the spectral analysis of more complex turbulent dispersed flows.

Assessment of MARS-KS prediction capability for natural circulation flow in passive heat removal system

Considering that system analysis codes are used for the evaluation of the performance of Passive Safety Systems (PSSs), it is important to investigate the capability of the system analysis code to reliably predict the heat transfer and natural circulation flow, which are the main phenomena governing the performance of a PSS. Since MARS-KS has been widely validated for heat transfer models, this study focuses on evaluating its capability to predict the single and two-phase pressure drops and natural circulation flow. The straight pipe simulation results indicate that the pressure drop predictions are reliable within ±5 % error margin for the single-phase flow and the errors of pressure drop up to −30 % for the two-phase flow. Through single-phase natural circulation flow analysis, it is concluded that the use of the appropriate K-factor modeling based on the flow regimes is important since the natural circulation flow rate in MARS-KS is mainly affected by the form loss factor modeling. With two-phase natural circulation flow analysis, this study emphasizes the behavior of the system could change significantly depending on the two-phase wall friction and pressure loss modeling. With the analysis results, modeling considerations for the PSS performance evaluation with the system analysis codes are proposed.

CFD Applications to Pressurized Thermal Shock-Related Phenomena

In pressurized water reactor accident scenarios, the injection of water from the emergency core cooling system (ECCS) (ECC injection) might induce a pressurized thermal shock (PTS), affecting the reactor pressure vessel (RPV) integrity. Therefore, PTS is a vital research issue in reactor safety, and its analysis is essential for evaluating the integrity of RPVs, which determines the reactor life. The PTS analysis comprises a coupled analysis between thermal–hydraulic and structural analyses. The thermal–hydraulic approach is particularly crucial, and reliable computational fluid dynamic (CFD) simulations should play a vital role in the future because predicting the temperature gradient of the RPV wall requires data on the transient temperature distribution of the downcomer (DC). Since one-dimensional codes cannot predict the complex three-dimensional flow features during ECC injection, PTS is one reactor safety issue where CFD simulation can benefit from complement evaluations with thermal–hydraulic system analysis codes. This study reviewed from the viewpoint of the turbulence models most affecting PTS analysis based on papers published since 2010 on single- and two-phase flow CFD simulation for the experiment on PTS performed in the Rossendorf coolant mixing model (ROCOM), transient two-phase flow (TOPFLOW), upper plenum test facility (UPTF), and large-scale test facility (LSTF). The results revealed that in single-phase flow CFD simulation, where knowledge and experience are sufficient, various turbulence models have been considered, and many analyses using large eddy simulation (LES) have been reported. For two-phase flow analysis of air–water conditions, interface capturing/tracking methods were used in addition to two-fluid models. The standard k–ε and shear stress transport (SST) k–ω models were still in the validated phase, and various turbulence models have yet to be fully validated. In the two-phase flow analysis of steam–water conditions, many studies have used two-fluid models and Reynolds-averaged Navier-Stoke (RANS), and NEPTUNE_CFD, in particular, has been reported to show excellent prediction performance based on years of accumulated validation.

Comparison of Refractive Index Matching Techniques and PLIF40 Measurements in Annular Flow.

This paper investigates non-invasive techniques for annular two-phase flow analysis, focusing on liquid film characterization to understand the interfacial phenomena that are crucial for heat and mass transfer. Limited methods allow the study of the temporal and spatial evolution of liquid film, such as Planar Laser-Induced Fluorescence (PLIF). However, this method possesses optical challenges, leading to the need for improved techniques to mitigate refraction and reflection, such as Refractive Index Matching (RIM). This study utilizes an experimental annular flow facility to analyze both RIM and non-RIM PLIF over a range of liquid Reynolds numbers from 4200 to 10,400. Three configurations-PLIF RIM90, PLIF RIM40, and PLIF nRIM40-are compared from both qualitative and quantitative perspectives. In the quantitative analysis, key variables of the liquid film are measured, namely mean film thickness, disturbance wave height, and frequency. Variations in the analyzed variables indicate minor deviations, which are not likely to be caused by the technique used. However, all three methodologies exhibited errors that are estimated to be within a maximum of 10%, with a mean value of approximately 8%.

Effect of slug characteristics on the nonlinear dynamic response of a long flexible fluid-conveying cylinder

Hydrocarbon flows in a marine riser may appear in multiple phases with varying flow patterns, among which, slug flows are known to exhibit complex flow characteristics due to fluctuations in multiphase mass, velocities, and pressure changes. Bulk of the literature concerning internal single and two-phase flow induced vibrations have largely focused on the use of a linearized tension beam model. In this study, the fluid-structure interaction phenomena of a submerged long flexible cylinder conveying two-phase slug flows is investigated taking into account the geometric and hydrodynamic nonlinearities. A semi-empirical theoretical model which consists of nonlinear structural equations of coupled out-of-plane, in-plane and axial structural motions is presented for the analysis and prediction of two-phase slug flow induced vibrations (SIV). The model includes equations involving centrifugal and Coriolis forces to capture mass variations in the slug flow regime. A numerical space-time finite difference approach combined with a frequency domain analysis is used for analyzing the highly nonlinear three-dimensional responses. The model assumes constant geometric properties across the span of the cylinder and utilizes an idealized slug unit concept, wherein slug properties are considered to be fully developed and undisturbed by pipe oscillations. Model validations are performed through comparisons with published experimental internal flow-induced vibration results. Parametric study captures the effects of several key slug characteristics on the vibration responses of a fluid-conveying cylinder. The results demonstrate amplitude-modulation response, mean displacements, three-dimensional cylinder displacements due to geometric nonlinear couplings and the influence of internal flow velocities. Higher dominant modes and oscillation frequencies were observed when the cylinder experiences large amplitude motions at specific slug formations. A new dimensionless parameter is introduced to predict scenarios of high amplitude oscillations for long flexible cylinders carrying two-phase slug flows.

Analysis and Calculation of Two-Phase Flow in Industrial Pipelines Using Regime Classification

The analysis shows that using regime classification for individual phases, it is possible to obtain a more reliable design scheme for two phase flows. In contrast to the flow structure, the concept of mode is more accessible for calculating an industrial pipeline. In this case, transient processes can be determined and calculated. The analysis shows that in field conditions it is very difficult to determine transition boundaries for flow structures. However, taking into account the structural approach, it is very difficult to determine the boundaries of the flow, and based on the laws of hydrodynamics of a homogeneous liquid and gas, it is possible to determine the transition of one regime to another using the Reynolds number. This technique is based on the theory of the great Russian scientist Academician A.P. Krylov. Let us note that this theory has not yet lost its specialty in the fishing industry. Based on numerous data, a method for calculating the gas-liquid mixture in horizontal pipes is proposed. Based on the laws of homogeneous liquid and gas, an equation for the distribution of speed, flow and pressure loss in a pipeline is obtained. The analysis shows that, based on the laws of homogeneous liquid and gas, it is possible to propose a method for calculating the hydrodynamic parameters of a two-phase flow, such as liquid and gas in a round pipe. This technique is based on the theory of Academician A.P. Krylov, which characterizes the reliability of this technique.

Numerical analysis of coarse particle two-phase flow in deep-sea mining vertical pipe transport with forced vibration

Exploring the dynamics and flow characteristics of coarse particles within vertical lifting pipes is essential to ensure safe and efficient pipeline transportation. This study combines computational fluid dynamics (CFD) and the discrete element method (DEM) to investigate the behavior of coarse particle solid–liquid two-phase flow in a vertical pipe utilized for deep-sea mining under forced vibration and compares it with that in stationary pipes. This study examined the effects of pipe amplitude, vibration frequency, and initial internal flow velocity separately. The amplitude values tested were 0D (static), 1D, 2D, 3D, and 4D, where D represents pipe diameter. The vibration frequencies of 0 Hz (static), 0.25 Hz, 0.5 Hz, 1 Hz, and 2 Hz were tested, respectively. Finally, initial velocities of 1 m/s, 2 m/s, 3 m/s, and 4 m/s were tested. Compared with stationary pipes, particle-induced turbulence gradually disrupts the flow field, leading to a more dispersed particle volume distribution during pipe vibration. An increase in the amplitude and vibration frequency of the pipe results in an augmentation of the pressure drop within the pipe. By comparing the results of the vibrating pipe with those of the stationary pipe, the effect of the pipe vibration gradually weakened as the initial velocity increased.