Liquid Flow Research Articles

Study on Liquid Hydrogen Leakage and Diffusion Behavior in a Hydrogen Production Station

Liquid hydrogen storage is an important way of hydrogen storage and transportation, which greatly improves the storage and transportation efficiency due to the high energy density but at the same time brings new safety hazards. In this study, the liquid hydrogen leakage in the storage area of a hydrogen production station is numerically simulated. The effects of ambient wind direction, wind speed, leakage mass flow rate, and the mass fraction of gas phase at the leakage port on the diffusion behavior of the liquid hydrogen leakage were investigated. The results show that the ambient wind direction directly determines the direction of liquid hydrogen leakage diffusion. The wind speed significantly affects the diffusion distance. When the wind speed is 6 m/s, the diffusion distance of the flammable hydrogen cloud reaches 40.08 m, which is 2.63 times that under windless conditions. The liquid hydrogen leakage mass flow rate and the mass fraction of the gas phase have a greater effect on the volume of the flammable hydrogen cloud. As the leakage mass flow rate increased from 5.15 kg/s to 10 kg/s, the flammable hydrogen cloud volume increased from 5734.31 m3 to 10,305.5 m3. The installation of a barrier wall in front of the leakage port can limit the horizontal diffusion of the flammable hydrogen cloud, elevate the diffusion height, and effectively reduce the volume of the flammable hydrogen cloud. This study can provide theoretical support for the construction and operation of hydrogen production stations.

Heat transfer in 3D radiative Oldroyd-B fluid flow with irregular heat source and activation energy

Heat transfer in 3D radiative Oldroyd-B fluid flow with irregular heat source and activation energy

Coolability of a Corium Pool in a Debris Bed: Impact of Debris Size, Steam, and Liquid Flow Rate; Tilting Angle; and Pressure on Critical Heat Flux—First Numerical Evaluations

In case of a severe accident in a light water reactor, a debris bed may form in the core and possibly melt, as happened in Three Mile Island Unit 2. Knowledge about the coolability of such a molten pool surrounded by debris is crucial to investigate the possibility of stabilizing a part of the fuel inside the vessel. In particular, it is of primary interest to determine the maximum size of a molten pool surrounded by debris that may be stabilized under water. A key parameter is the maximum heat flux [critical heat flux (CHF)] that may be extracted from the pool boundary by water flowing within the debris bed. A facility was built at IRSN to determine the CHF under various conditions. A heated copper surface simulates the boundary of the pool and is placed in contact with a debris bed (monodisperse steel balls), under water. In this article, we first complete previous experimental work with steam and liquid flow rates, ball diameter, tilting angle of the heated surface, and pressure up to 2.5 bar absolute. A general CHF correlation depending on the tilting angle and the steam flow rate is derived, and an example of the use of this correlation to evaluate the maximum mass of the corium pool that can be stabilized under water is given, for some typical reactor conditions. In the last part of the article, the accuracy of the ASTEC code is evaluated based on first calculations intending to reproduce the experimental impact of liquid and steam flow rates.

Developing a Three-Dimensional Bubble Dynamics Model for Centrifugal Nuclear Thermal Propulsion

Nuclear thermal propulsion (NTP) is a highly efficient type of propulsion system that could yield a specific impulse up to 800 s. High-performance NTP allows for possibilities of various mission types from flyby to interplanetary flight, manned or unmanned. Centrifugal nuclear thermal propulsion (CNTP) is a type of NTP that uses liquid fissile material instead of a solid core in a standard reactor configuration. CNTP requires an interaction between liquid fissile fuel and gaseous propellant to perform heat transfer and generate thrust in a spinning enclosure. Because of the liquid state of the fissile fuel, it is essential to monitor bubble behavior and trajectory. For example, the bubble dwell time within the liquid and hydrogen flow rates can influence the liquid volume fraction and volume, which in turn can affect reactivity. A model is developed and presented that, as it matures, will predict the bubbles’ behavior based on the various types of liquid and gas and their corresponding properties as well as assist in cross-validation with the concurring experiment. The model primarily utilizes the smooth particle hydrodynamics (SPH) approach to simulate the liquid and gas interaction within the design space. The physics include buoyancy and drag on the bubbles. The pressure field within the liquid is modeled using hydrostatics, in which centrifugal motion introduces a radial gradient. Boundary conditions are devised to confine the liquid and gas within the enclosure. In this paper, we present the progress to date in two centrifugal fuel element–relevant subscale devices, the so-called “Antfarm” and “Blender II.” Antfarm is a static stage in a rectangular enclosure where buoyancy is purely gravity driven. Blender II refers to a rectangular enclosure that spins at up to 7000 rpm. Both Antfarm and Blender II are experiments that provide important data for validation against our SPH model. Two liquid and gas combinations are modeled in the present study using the Antfarm setup: water-air, and Galinstan-Nitrogen. The bubbles’ behavior was comparable to the experiment, and the velocity was at the same order of magnitude. The liquid simulation performed for the Blender II model shows that the pressure gradient within the liquid in the radial direction matches that predicted analytically from the centrifugal acceleration. The Blender II liquid model awaits experimental data for validation and verification.

Dynamics and Chaos of Convective Fluid Flow

In this paper, a mathematical model of fluid flows in a convective thermal system is developed, and a five-dimensional dynamical system is developed for the investigation of the convective fluid dynamics. The analytical solutions of periodic motions to chaos of the convective fluid flows are developed for steady-state vortex flows, and the corresponding stability and bifurcations of periodic motions in the five-dimensional dynamical system are studied. The harmonic frequency-amplitude characteristics for periodic flows are obtained, which provide energy distribution in the parameter space. Analytical homoclinic orbits for the convective fluid flow systems are developed for the asymptotic convection through the infinite-many homoclinic orbits in the five-dimensional dynamical system. The dynamics of fluid flows in the convective thermal systems are revealed, and one can use such methodology to predict atmospheric and oceanic phenomena through thermal convections.

Optimization of Dropwise Condensation of Steam over Hybrid Hydrophobic–Hydrophilic Surfaces via Enhanced Statistically Based Heat Transfer Modelization

Steam condensation over a hybrid hydrophobic–hydrophilic surface is modeled via simplified heat transfer modelization. Filmwise condensation is assumed over the hydrophilic region. The standard film model is improved, accounting for the liquid flow rate crossing the hydrophobic–hydrophilic boundaries. A threshold for flooding occurrence is also presented. Dropwise condensation is assumed over the hydrophobic region. Compared to the heat transfer models in the literature, based on the statistical drop size distribution, a novel correlation is used for the size distribution of small droplets. The correlations of both the liquid flow rate crossing the hydrophobic–hydrophilic boundary and the size distribution of small drops are derived via Lagrangian simulations, using an in-house code previously developed and validated by the authors. The heat transfer model is validated with experimental data in the literature involving a hybrid surface, composed by alternate vertical hydrophobic–hydrophilic stripes. Then, the optimization of the hybrid surface geometry is performed in terms of hydrophobic width and hydrophilic width, with the aim of enhancing the heat flux.

Nanoenzyme Hydrogel Film-Based Portable Point-of-Care Testing Platform for Double-Signal Visual Detection of PSA.

The development of the Point-of-Care Testing (POCT) platform that combines convenience and cost-effectiveness is crucial for enabling the visual detection of disease biomarkers. In this work, a POCT platform for the sensitive in situ detection of prostate specific antigen (PSA) with dual-signal output was constructed by functionalizing the Eppendorf (EP) tube. This was achieved through the modification of aptamer hairpin probes (AHPs) on the lid of the EP tube and the assembly of a nanoenzyme hydrogel film on its inner wall. The target could trigger the release of Ag+ by AHP and subsequently activate Ag+-dependent DNAzyme (Ag-DNAzyme). This would initiate the cleavage of the DNA-Au/Pt NP hydrogel network, leading to the release of Au/Pt NPs. The released Au/Pt NPs exhibit both peroxidase (POD)-like and catalase (CAT)-like activity to produce a colorimetric response and induce liquid flow under pressure. Therefore, the target can be measured visually and quantitatively through colorimetric analysis and the measurement of total dissolved solids (TDS) using a pressure-triggered liquid flow device integrated into the platform. The designed platform is distinguished by its simplicity, specificity, cost-effectiveness, and remarkable sensitivity. It allows for the visual detection of PSA within concentration ranges of 0.5-100 ng/L (colorimetric) and 3-100 ng/L (TDS reading), boasting detection limits as low as 0.15 ng/L (colorimetric) and 0.57 ng/L (TDS reading). The strategy of target-triggered nanoenzyme release significantly enhances sensitivity and provides a guiding approach for visual biomarker detection.

Continuous casting process based on an optimized submerged entry nozzle

Continuous casting process is an important part in the steel manufacturing process. The slab production quality is commonly influenced by mould and flow pattern of submerged entry nozzles (SEN). Hence, the present work aimed to attain a more liquid flow rate to increase the slab production quantity within a certain period. Moreover, a novel African Buffalo Algorithm (ABA) is presented in this work to optimise the SEN parameters and the SEN was designed in the MATLAB environment. The developed ABA model optimises the design parameters of SEN such as the size of the nozzle, the shape of the port, length, submergence depth, etc. The fitness of buffalo effectively optimises the SEN parameters and by utilising the optimised design parameters, the SEN was resigned in the MATLAB to identify the proposed framework's effectiveness. The results demonstrated that the liquid flow rate is increased and surface tension gets minimised. When the flow rate of liquid grows, the wave amplitude also upturns. The effect of liquid and air flow rate, port shape and size, maximum wave amplitude, surface tension, etc., were calculated with 0° parallel port. In addition, the presented SEN's results were compared with other existing SEN to determine the optimised SEN's effectiveness. The comparison result demonstrated that the designed optimised SEN has a high liquid flow rate, casting speed, density and viscosity.

Erosion wear characteristics of metal seal floating ball valve in gas–solid–liquid three-phase flow

During high-pressure acid leaching, sulfuric acid slurry acts as the transportation medium in the withdrawal drum for acid pressure leaching. Over time, the metal-sealed float valve at the pipeline outlet may experience internal leakage. To address this issue, the Euler–Lagrange method for fluid dynamics is used to evaluate the corrosion wear characteristics of the valve in a three-phase flow. Specifically, we investigate the erosion and wear properties of the ball valve core as it opened and closed. To further understand how different ball valve openings impact erosion characteristics, we establish a correlation between the corrosion rate and the drag coefficient. Our research findings indicate that as the valve’s relative opening decreases from 80% to 20%, the corrosion wear rate gradually increases. In addition, the number of particles within the valve cavity also follows a similar trend. These erosion results highlight the significant influence that particle quantity has on erosion. Moreover, with an increase in particle size, both the erosion gradient range and maximum erosion rate decrease.

Ultrasensitive detection of interleukin-6 enabled by three-dimensional spatial amplification

Accurately measuring interleukin-6 (IL-6) levels is crucial in both clinical medicine and research due to its role in various physiological and pathological processes. We present proof-of-concept of a novel three-dimensional (3D) spatial amplification method for developing highly sensitive and miniaturized IL-6 biosensors. We designed an IL-6 immunosensor platform utilizing a unique 3D microfluidic structure fabricated with femtosecond-pulsed laser processing. This design features a porous transducer element with over 4500 microcavities, thereby increasing the reaction area by 13-fold compared to traditional 2D designs while maintaining transparency. By combining this 3D structure with a competitive antigen-antibody reaction, the sensor achieved an exceptional limit of detection of 0.27 pg/mL for IL-6 in human blood samples. Additionally, we demonstrated effective control of liquid flow within the porous element using a reduction of pressure and centrifugation speeds. This 3D spatial amplification method offers a promising approach for developing highly sensitive and compact IL-6 biosensors. Such miniaturized sensors hold potential for point-of-care testing devices, enabling convenient and timely IL-6 analysis.

Причины повреждений корпусных конструкций в носовой оконечности судов

The operation of marine vessels, especially large-tonnage ones, is accompanied by the occurrence of serious damage to hull structures within bow in stormy conditions. A new mechanism to explain the catastrophic destruction of the side grillages of large-tonnage ships in the area of the bow is proposed. Simulation of the process of flow around the fore end when it is digged in a wave using SPH technology showed that when the flow speeds around the fore end occur when the vessel moves in developed waves, there is a sharp drop in pressure in certain areas of the bow end. This circumstance can lead to the occurrence of cavitation that is confirmed by modeling the process of flow around the fore end using ANSYS FLUENT software. During stall cavitation, a cavitation cavern of considerable size may arise, the collapse of which is accompanied by the occurrence of hydrodynamic loads that can destroy hull structures. A calculation method is proposed that allows one to estimate the magnitude of hydrodynamic loads acting on the ship’s hull during the collapse of a cavitation cavern. To prevent damage to the bow in storm conditions, a number of constructive measures have been proposed. One of the options for solving the problem is to install a deck fairing in the bow that eliminates the negative influence of structural cavitators in the traditional design when the bow is digged in the wave. It is also possible to install a compressed air supply system to those areas of the bow where vacuum occurs when a liquid flows around when digged in a wave.

Investigations on transport behaviours of rotating packed bed for CO2 capture by chemical absorption using CFD simulation

A rotating packed bed (RPB) has been proposed to overcome the limitations of high operating and capital costs involved in conventional columns which are used for chemical absorption-based carbon dioxide (CO2) capture. The design of large-scale RPBs depends on their transport characteristics. Owing to the difficulties in predicting the transport characteristics of RPBs experimentally, Computational Fluid Dynamics (CFD) simulation has been implemented. In this study, a two-dimensional (2D) volume of fluid (VOF) simulation is conducted to predict transport characteristics of RPB including, hydrodynamics, and physical and reactive mass transfers. To predict hydrodynamics of RPBs accurately, the liquid holdup by CFD simulation is compared against experimental prediction, showing a strong agreement between simulation and experimentation. Various liquid phase flow patterns within RPB are observed and these are cross-referenced with findings in the existing literature. Additionally, the impacts of rotational motion and flow rate are further assessed in the study. For the mass transfer performance, the prediction of the volumetric mass transfer coefficient, kLae of the air-water system from CFD simulation is studied. Two different mass transfer theories such as surface renewal theory and penetration theory are studied to predict kLae. It is observed that the kLae from the penetration theory closely follows the trend of the experimental data. Finally, a 2D simulation of the reactive absorption of CO2 by diethylenetriamine (DETA) in RPB is developed and compared with the existing available data in the literature. The present model underpredicts CO2 capture efficiency in RPB, which needs further improvement in reaction kinetics. In general, the present CFD model has the capability of predicting the transport characteristics of RPB.

New insights into the parameterization of the dry surface layer and its hydrogeochemical mechanism: An experimental study

Knowledge of the parameterization of the dry surface layer (DSL) is essential for evaluating near-surface water flow and water balance in arid and semi-arid areas. Existing studies have parameterized DSL thickness and vapor flow as functions of the soil moisture content (SMC) in the surface layer to predict soil evaporation. However, hydrochemical processes related to DSL development have been ignored, including changes in hydrochemistry, the underlying hydrochemical mechanism, and the role of dissolved substances in the DSL development. Herein, we performed a series of soil evaporation experiments for 260 days and explored the factors influencing DSL development (e.g., soil texture, atmospheric temperature, SMC, solutes). Evaporation experiments were performed using silty loess, sandy loess, and fine sand with a 60-cm water table. Results showed that the cumulative evaporation of silty loess, sandy loess, and fine sand over the experimental period were 1,391.52, 460.10, and 185.53 mm, respectively, which determined by the maximum height of liquid flow continuity. The content of total dissolved solids (TDS) and major ions at the surface soil were significantly higher than the values at deep depths of 5‒55 cm, which largely depend on evaporative water loss. Evolutionary trends of chemical facies in sand media along the liquid water migration were from HCO3-Ca type to SO4·Cl-Na type. This was attributed to mineral dissolution at a depth of 5–55 cm and their transport with liquid water, resulting in the precipitation of salt crystals at the surface soil. Furthermore, a consolidated DSL with a thickness of 3.0–3.5 cm in the sandy loess and a loose DSL with a thickness of 1.5–2.0 cm in fine sand were observed at the end of the experiments. The accumulation of solutes at the surface leads to a reduction in effective porosity and the aggregation of soil particles during continuous drying, which facilitates the consolidation of DSL in sandy loess. This overestimated the DSL thickness, resulting in a difference between the experimental and predicted evaporation rates by Fick's law. Overall, these results highlight the limitations of considering DSL thickness as a function of SMC only, providing new insights into hydrochemical processes and dissolved solutes involving DSL parameterization during continuous soil drying.

Statistical Analysis of Bubble Parameters from a Model Bubble Column with and without Counter-Current Flow

Bubble columns are widely used in numerous industrial processes because of their advantages in operation, design, and maintenance compared to other multiphase reactor types. In contrast to their simple design, the generated flow conditions inside a bubble column reactor are quite complex, especially in continuous mode with counter-current liquid flow. For the design and optimization of such reactors, precise numerical simulations and modelling are needed. These simulations and models have to be validated with experimental data. For this reason, experiments were carried out in a laboratory-scale bubble column using shadow imaging and particle image velocimetry (PIV) techniques with and without counter-current liquid flow. In the experiments, two types of gases—relatively poorly soluble air and well-soluble CO2—were used and the bubbles were generated with three different capillary diameters. With changing gas and liquid flow rates, overall, 108 different flow conditions were investigated. In addition to the liquid flow fields captured by PIV, shadow imaging data were also statistically evaluated in the measurement volume and bubble parameters such as bubble diameter, velocity, aspect ratio, bubble motion direction, and inclination. The bubble slip velocity was calculated from the measured liquid and bubble velocities. The analysis of these parameters shows that the counter-current liquid flow has a noticeable influence on the bubble parameters, especially on the bubble velocity and motion direction. In the case of CO2 bubbles, remarkable bubble shrinkage was observed with counter-current liquid flow due to the enhanced mass transfer. The results obtained for bubble aspect ratio are compared to known correlations from the literature. The comprehensive and extensive bubble data obtained in this study will now be used as a source for the development of correlations needed in the validation of numerical simulations and models. The data are available from the authors on request.

Dynamics of the ohmic heating and chemically reactive time-dependent flow of magnetized Casson fluid inside a rectangular pipe

Chemical reaction and ohmic heating can significantly influence the flow behavior and heat transfer characteristics in industrial processes involving non-Newtonian fluids. These processes include polymer processing, food processing, and manufacturing processes where non-Newtonian fluids are commonly encountered. Therefore, the impacts of chemical reaction and ohmic heating on the unsteady flow of an incompressible magnetized non-Newtonian Casson liquid through a rectangular pipe with a variable pressure gradient is scrutinized in this article. The impacts of heat generation/absorption and viscous dissipation are also incorporated through the energy equation. The numerical simulations of partial differential equations (PDEs) are acquired utilizing the finite difference method with the incorporating relaxation (SOR) algorithm. The impacts of dissimilar imperative parameters on various flow fields are scrutinized through graphs. Outcomes indicate that liquid velocity decreases with enhancing values of the Casson parameter. Furthermore, the liquid temperature has an enhancing behavior with the impact of the heat source/sink parameter. Comparisons of numerical results with previously published results show an excellent agreement. This kind of research may find specific industrial and medical utilizations such as glass manufacturing, crude oil purification, lubrication, paper production, blood transport study in cardiovascular design, etc.

Flow Velocity Computation in Solid–Liquid Two-Phase Flow by a Hybrid Network CNN–RKSVM

As an advanced detection technique, electrical resistive tomography (ERT) has been applied to detect the solid–liquid two-phase flow velocity based on available ERT measurements. The flow velocity computation by ERT must depend on the relative algorithms, including both the cross-correlation (CC) principle and convolutional neural networks (CNNs). However, these two types of algorithms have poor accuracy and generalization under complex measuring conditions and various flow patterns. To address this issue, in this paper, a hybrid network is proposed that combines a CNN with a reproducing kernel-based support vector machine (RKSVM) technique. The features hidden in ERT measurements are extracted using the CNN, and then the flow velocity is computed by the RKSVM in a high-dimensional feature space. According to the ERT measurements in an actual experimental platform, the results show that the hybrid network has higher accuracy and generalization ability for flow velocity computation compared with the existing CC, RKSVM, and CNN methods.

EMHD micropolar fluid flowing through a micro‐structural slipped surface with heat source/sink and chemical reaction

AbstractThe electromagnetohydrodynamic (EMHD) has a vital role due to its importance in aerospace and plasma physics, including energy transformation systems in which interaction between the liquid and magnetic field is crucial. Further, micropolar fluids with microstructural slip is used in the lubrication and liquid crystals. From this observation, the current study is conducted to express an EMHD flow of liquid on a microstructural slipped surface. The effects of a uniform heat source/sink (HS/S) with homogeneous and heterogeneous chemical reactions have been incorporated in energy and mass profiles. The governing equations are converted to ordinary differential equations using proper similarity variables. The converted equations are computed by the implementation of Runge–Kutta–Fehlberg (RKF) 4th 5th order with shooting technique. A list of the essential dimensionless constraints and graphs showing how they are affected are provided. The outcome of the problem shows that the electric parameter will improve the velocity but decline the microrotation profile. Further, both heterogeneous and homogeneous reactions have a decreasing concentration. While the heat distribution rate increases with greater magnetic and electric fields, the surface drag force decreases.

DROPLET MOTION ON HYDROPHOBIC AND HYDROPHILIC SURFACES USING LATTICE BOLTZMANN METHOD WITH BLOOD TEST MICROPAD APPLICATION

Laboratory sciences and disease diagnostics need accelerated medical tests, ease of testing for patients, and lower diagnostic costs. In this respect, one of the most innovative methods is paper-based microfluidic analysis tools, called micro-PAD. Since the chemical-based micro-PADs need sufficient blood to show an accurate blood test, the shape of the micro-PADs, the blood quantity, and the droplet flow to reach the end of the microchannels are important. Therefore, this research showed the liquid droplet motion on a heterogeneous (hydrophobic and hydrophilic) surface using the lattice Boltzmann method, and two star-shaped micro-PADs were modeled to investigate the effect of the droplet size and the droplet deflection from the center of the micro-PAD on liquid flow in microchannels. Based on the obtained results, fluid behavior was similar in two patterns (ten-pointed pattern and four-pointed pattern). This study compared the effects of parallel, convergent, and divergent channels in a four-pointed pattern on fluid velocity. The off-center error of the droplet was considered at the initial moment. The obtained results confirmed the direct relationship between the increased deviation of the center and the lack of uniformity in the fluid penetration. The present simulation results were validated against the laboratory data.

ELECTRO-OSMOTIC EFFECT ON PERISTALTIC FLOW OF PHAN–THIEN–TANNER FLUID IN A PLANAR CHANNEL

Abstract There are several factors that can cause the excessive accumulation of biofluid in human tissue, such as pregnancy, local traumas, allergic responses or the use of certain therapeutic medications. This study aims to further investigate the shear-dependent peristaltic flow of Phan–Thien–Tanner (PTT) fluid within a planar channel by incorporating the phenomenon of electro-osmosis. This research is driven by the potential biomedical applications of this knowledge. The non-Newtonian fluid features of the PTT fluid model are considered as physiological fluid in a symmetric planar channel. This study is significant, as it demonstrates that the chyme in the small intestine can be modelled as a PTT fluid. The governing equations for the flow of the ionic liquid, thermal radiation and heat transfer, along with the Poisson–Boltzmann equation within the electrical double layer, are discussed. The long-wavelength ( $\delta \ll 1$ ) and low-Reynolds-number approximations ( $Re \to 0$ ) are used to simplify the simultaneous equations. The solutions analyse the Debye electronic length parameter, Helmholtz–Smoluchowski velocity, Prandtl number and thermal radiation. Additionally, streamlines are used to examine the phenomenon of entrapment. Graphs are used to explain the influence of different parameters on the flow and temperature. The findings of the current model have practical implications in the design of microfluidic devices for different particle transport phenomena at the micro level. Additionally, the noteworthy results highlight the advantages of electro-osmosis in controlling both flow and heat transfer. Ultimately, our objective is to use these findings as a guide for the advancement of lab-on-a-chip systems.

Modeling and simulation of air-water upward annular flow characteristics in a vertical tube using CFD

Annular flow refers to a special type of two-phase flow pattern in which liquid flows as a thin film at the periphery of a pipe, tube, or conduit, and gas with relatively high velocity flows at the center of the flow section. This gas also includes dispersed liquid droplets. The liquid film flow rate continuously changes inside the tube due to two processes-entrainment and deposition. To determine the liquid holdup, pressure drop, the onset of dryout, and heat transfer characteristics in annular flow, it is important to have proper knowledge of flow characteristics. Especially a better understanding of entrainment fraction is important for the heat transfer and safe operation of two-phase flow systems operating in an annular two-phase flow regime. Therefore, the objective of this work is to develop a computational model for the simulation of the annular two-phase flow regime and assess the various existing models for the entrainment rate. In this work, Computational Fluid Dynamics (CFD) in ANSYS FLUENT has been applied to determine annular flow characteristics such as liquid film thickness, film velocity, entrainment rate, deposition rate, and entrainment fraction for various gas-liquid flow conditions in a vertical upward tube. The gas core with droplets was simulated using the Discrete Phase Model (DPM) which is based on the Eulerian-Lagrangian approach. The Eulerian Wall Film (EWF) model was utilized to simulate liquid film on the tube wall. Three different models of Entrainment rate were implemented and assessed through user-defined functions (UDF) in ANSYS. Finally, entrainment for fully developed flow was determined and compared with the experimental data available in the literature. From the simulations, it was obtained that the Bertodano correlation performed best in predicting entrainment fraction and the results were within the ±30 % limit when compared to experimental data.