Abstract: This study investigated the void fraction characteristics under low superficial gas velocity (Jg) and low superficial fluid velocity (Jf) conditions in a 3x3 rod bundle geometry. Three arrangements of conductivity probes were set to measure the void fraction at various cross-sectional regions, including rod-gap, sub-channel and rod-wall regions. The experimental tests were performed under the flow conditions of Jg = 0-0.236 m/s and Jf = 0-0.142 m/s, and the time-averaged void fractions were recorded at each flow condition. It was observed that while the superficial gas velocity increases, the small bubbles started to cluster together and become big bubbles. As the superficial fluid velocity increases, the local void fractions of the three test regions will get closer and the bubble distribution will be more uniform across the cross section.
Abstract: Use of ultrasound waves is one of the techniques for increasing the mixing and mass transfer in the microdevices. Ultrasound propagation into liquid medium leads to stimulation of the fluid, creates turbulence and so increases the mixing performance. In this study, CFD modeling of two-phase flow in a pitted micromixer equipped with a piezoelectric with frequency of 1.7 MHz has been studied. CFD modeling of micromixer at different velocity of fluid flow in the absence of ultrasound waves and with ultrasound application has been performed. The hydrodynamic of fluid flow and mixing efficiency for using ultrasound has been compared with the layout of no ultrasound application. The result of CFD modeling shows well agreements with the experimental results. The results showed that the flow pattern inside the micromixer in the absence of ultrasound waves is parallel, while when ultrasound has been applied, it is not parallel. In fact, propagation of ultrasound energy into the fluid flow in the studied micromixer changed the hydrodynamic and the forms of the flow pattern and caused to mixing enhancement. In general, from the CFD modeling results, it can be concluded that the applying ultrasound energy into the liquid medium causes an increase in the turbulences and mixing and consequently, improves the mass transfer rate within the micromixer.
Abstract: This study analyzed the transient signals of central region void fraction of air-water two-phase flow in a 3x3 rod bundle. Experimental tests were carried out utilizing a vertical rod bundle test section along with a set of air-water supply/flow control system, and the transient signals of the central region void fraction were collected through the electrical conductivity sensors as well as visualized via high speed photography. By converting the electric signals, transient void fraction can be obtained through the voltage ratios. With a fixed superficial water velocity (Jf=0.094 m/s), two different superficial air velocities (Jg=0.094 m/s and 0.236 m/s) were tested and presented, which were corresponding to the flow conditions of bubbly flows and cap/slug flows, respectively. The time averaged central region void fraction was obtained as 0.109-0.122 with 0.028 standard deviation for the selected bubbly flow and 0.188-0.221with 0.101 standard deviation for the selected cap/slug flow, respectively. Through Fast Fourier Transform (FFT) analysis, no clear frequency peak was found in bubbly flow, while two dominant frequencies were identified around 1.6 Hz and 2.5 Hz in the present cap/slug flow.
Abstract: We present and analyze reliable numerical techniques
for simulating complex flow and transport phenomena related to
natural gas transportation in pipelines. Such kind of problems
are of high interest in the field of petroleum and environmental
engineering. Modeling and understanding natural gas flow and
transformation processes during transportation is important for the
sake of physical realism and the design and operation of pipeline
systems. In our approach a two fluid flow model based on a system
of coupled hyperbolic conservation laws is considered for describing
natural gas flow undergoing hydratization. The accurate numerical
approximation of two-phase gas flow remains subject of strong
interest in the scientific community. Such hyperbolic problems are
characterized by solutions with steep gradients or discontinuities, and
their approximation by standard finite element techniques typically
gives rise to spurious oscillations and numerical artefacts. Recently,
stabilized and discontinuous Galerkin finite element techniques
have attracted researchers’ interest. They are highly adapted to the
hyperbolic nature of our two-phase flow model. In the presentation
a streamline upwind Petrov-Galerkin approach and a discontinuous
Galerkin finite element method for the numerical approximation of
our flow model of two coupled systems of Euler equations are
presented. Then the efficiency and reliability of stabilized continuous
and discontinous finite element methods for the approximation is
carefully analyzed and the potential of the either classes of numerical
schemes is investigated. In particular, standard benchmark problems
of two-phase flow like the shock tube problem are used for the
comparative numerical study.
Abstract: The diffusion-absorption refrigeration cycle consists of a generator bubble pump, an absorber, an evaporator and a condenser, and usually operates with ammonia/water/ hydrogen or helium as the working fluid. The aim of this paper is to study the stability problem a bubble pump. In fact instability can caused a reduction of bubble pump efficiency. To achieve this goal, we have simulated the behaviour of two-phase flow in a bubble pump by using a drift flow model. Equations of a drift flow model are formulated in the transitional regime, non-adiabatic condition and thermodynamic equilibrium between the liquid and vapour phases. Equations resolution allowed to define void fraction, and liquid and vapour velocities, as well as pressure and mixing enthalpy. Ammonia-water mixing is used as working fluid, where ammonia mass fraction in the inlet is 0.6. Present simulation is conducted out for a heating flux of 2 kW/m² to 5 kW/m² and bubble pump tube length of 1 m and 2.5 mm of inner diameter. Simulation results reveal oscillations of vapour and liquid velocities along time. Oscillations decrease with time and with heat flux. For sufficient time the steady state is established, it is characterised by constant liquid velocity and void fraction values. However, vapour velocity does not have the same behaviour, it increases for steady state too. On the other hand, pressure drop oscillations are studied.
Abstract: Two-phase and multi-phase flows are common flow types in fluid mechanics engineering. Among the basic and applied problems of these flow types, two-phase parallel flow is the one that two immiscible fluids flow in the vicinity of each other. In this type of flow, fluid properties (e.g. density, viscosity, and temperature) are different at the two sides of the interface of the two fluids. The most challenging part of the numerical simulation of two-phase flow is to determine the location of interface accurately. In the present work, a coupled interface tracking algorithm is developed based on Arbitrary Lagrangian-Eulerian (ALE) approach using a cell-centered, pressure-based, coupled solver. To validate this algorithm, an analytical solution for fully developed two-phase flow in presence of gravity is derived, and then, the results of the numerical simulation of this flow are compared with analytical solution at various flow conditions. The results of the simulations show good accuracy of the algorithm despite using a nearly coarse and uniform grid. Temporal variations of interface profile toward the steady-state solution show that a greater difference between fluids properties (especially dynamic viscosity) will result in larger traveling waves. Gravity effect studies also show that favorable gravity will result in a reduction of heavier fluid thickness and adverse gravity leads to increasing it with respect to the zero gravity condition. However, the magnitude of variation in favorable gravity is much more than adverse gravity.
Abstract: Recently increasing oil production from petroleum reservoirs is one of the most important issues in the global energy sector. So, in this paper, the recovery of oil by the waterflooding technique from petroleum reservoir are considered. To investigate the aforementioned phenomena, the relative permeability of two immiscible fluids in sand is measured in the laboratory based on the steady-state method. Two sorts of oils, kerosene and heavy oil, and water are pumped simultaneously into a vertical sand column with different pumping ratio. From the change in fractional discharge measured at the outlet, a method for determining the relative permeability is developed focusing on the displacement mechanism in sand. Then, displacement mechanism of two immiscible fluids in the sand is investigated under the Buckley-Leveret frontal displacement theory and laboratory experiment. Two sorts of experiments, one is the displacement of pore water by oil, the other is the displacement of pore oil by water, are carried out. It is revealed that the relative permeability curves display tolerably different shape owing to the properties of oils, and produce different amount of residual oils and irreducible water saturation.
Abstract: Flow instability during gas lift operation is caused by three major phenomena – the density wave oscillation, the casing heading pressure and the flow perturbation within the two-phase flow region. This paper focuses on the causes and the effect of flow instability during gas lift operation and suggests ways to control it in order to maximise productivity during gas lift operations. A laboratory-scale two-phase flow system to study the effects of flow perturbation was designed and built. The apparatus is comprised of a 2 m long by 66 mm ID transparent PVC pipe with air injection point situated at 0.1 m above the base of the pipe. This is the point where stabilised bubbles were visibly clear after injection. Air is injected into the water filled transparent pipe at different flow rates and pressures. The behavior of the different sizes of the bubbles generated within the two-phase region was captured using a digital camera and the images were analysed using the advanced image processing package. It was observed that the average maximum bubbles sizes increased with the increase in the length of the vertical pipe column from 29.72 to 47 mm. The increase in air injection pressure from 0.5 to 3 bars increased the bubble sizes from 29.72 mm to 44.17 mm and then decreasing when the pressure reaches 4 bars. It was observed that at higher bubble velocity of 6.7 m/s, larger diameter bubbles coalesce and burst due to high agitation and collision with each other. This collapse of the bubbles causes pressure drop and reverse flow within two phase flow and is the main cause of the flow instability phenomena.
Abstract: The present paper discusses the prediction of gas-liquid two-phase frictional pressure drop in a 2.12 mm horizontal circular minichannel using Artificial Neural Network (ANN). The experimental results are obtained with air as gas phase and water as liquid phase. The superficial gas velocity is kept in the range of 0.0236 m/s to 0.4722 m/s while the values of 0.0944 m/s, 0.1416 m/s and 0.1889 m/s are considered for superficial liquid velocity. The experimental results are predicted using different Artificial Neural Network (ANN) models. Networks used for prediction are radial basis, generalised regression, linear layer, cascade forward back propagation, feed forward back propagation, feed forward distributed time delay, layer recurrent, and Elman back propagation. Transfer functions used for networks are Linear (PURELIN), Logistic sigmoid (LOGSIG), tangent sigmoid (TANSIG) and Gaussian RBF. Combination of networks and transfer functions give different possible neural network models. These models are compared for Mean Absolute Relative Deviation (MARD) and Mean Relative Deviation (MRD) to identify the best predictive model of ANN.
Abstract: In many practical situations, bubbles are dispersed in a
liquid phase. Understanding these complex bubbly flows is therefore
a key issue for applications such as shell and tube heat exchangers,
mineral flotation and oxidation in water treatment. Although a large
body of work exists for bubbles rising in an unbounded medium,
that of bubbles rising in constricted geometries has received less
attention. The particular case of a bubble sliding underneath an
inclined surface is common to two-phase flow systems. The current
study intends to expand this knowledge by performing experiments
to quantify the streamwise flow structures associated with a single
sliding air bubble under an inclined surface in quiescent water. This
is achieved by means of two-dimensional, two-component particle
image velocimetry (PIV), performed with a continuous wave laser
and high-speed camera. PIV vorticity fields obtained in a plane
perpendicular to the sliding surface show that there is significant bulk
fluid motion away from the surface. The associated momentum of the
bubble means that this wake motion persists for a significant time
before viscous dissipation. The magnitude and direction of the flow
structures in the streamwise measurement plane are found to depend
on the point on its path through which the bubble enters the plane.
This entry point, represented by a phase angle, affects the nature and
strength of the vortical structures. This study reconstructs the vorticity
field in the wake of the bubble, converting the field at different
instances in time to slices of a large-scale wake structure. This is, in
essence, Taylor’s ”frozen turbulence” hypothesis. Applying this to the
vorticity fields provides a pseudo three-dimensional representation
from 2-D data, allowing for a more intuitive understanding of the
bubble wake. This study provides insights into the complex dynamics
of a situation common to many engineering applications, particularly
shell and tube heat exchangers in the nucleate boiling regime.
Abstract: A Multi-dimensional computational fluid dynamics
(CFD) two-phase model was developed with the aim to simulate
the in-core coolant circuit of a pressurized heavy water reactor
(PHWR) of a commercial nuclear power plant (NPP). Due to the
fact that this PHWR is a Reactor Pressure Vessel type (RPV),
three-dimensional (3D) detailed modelling of the large reservoirs of
the RPV (the upper and lower plenums and the downcomer) were
coupled with an in-house finite volume one-dimensional (1D) code
in order to model the 451 coolant channels housing the nuclear fuel.
Regarding the 1D code, suitable empirical correlations for taking into
account the in-channel distributed (friction losses) and concentrated
(spacer grids, inlet and outlet throttles) pressure losses were used.
A local power distribution at each one of the coolant channels
was also taken into account. The heat transfer between the coolant
and the surrounding moderator was accurately calculated using a
two-dimensional theoretical model. The implementation of subcooled
boiling and condensation models in the 1D code along with the use
of functions for representing the thermal and dynamic properties of
the coolant and moderator (heavy water) allow to have estimations
of the in-core steam generation under nominal flow conditions for a
generic fission power distribution. The in-core mass flow distribution
results for steady state nominal conditions are in agreement with the
expected from design, thus getting a first assessment of the coupled
1/3D model. Results for nominal condition were compared with
those obtained with a previous 1/3D single-phase model getting more
realistic temperature patterns, also allowing visualize low values of
void fraction inside the upper plenum. It must be mentioned that the
current results were obtained by imposing prescribed fission power
functions from literature. Therefore, results are showed with the aim
of point out the potentiality of the developed model.
Abstract: In this study, the three-dimensional cavitating
turbulent flow in a complete Francis turbine is simulated using
mixture model for cavity/liquid two-phase flows. Numerical analysis
is carried out using ANSYS CFX software release 12, and standard k-ε
turbulence model is adopted for this analysis. The computational
fluid domain consist of spiral casing, stay vanes, guide vanes, runner
and draft tube. The computational domain is discretized with a threedimensional
mesh system of unstructured tetrahedron mesh. The
finite volume method (FVM) is used to solve the governing equations
of the mixture model. Results of cavitation on the runner’s blades
under three different boundary conditions are presented and
discussed. From the numerical results it has been found that the
numerical method was successfully applied to simulate the cavitating
two-phase turbulent flow through a Francis turbine, and also
cavitation is clearly predicted in the form of water vapor formation
inside the turbine. By comparison the numerical prediction results
with a real runner; it’s shown that the region of higher volume
fraction obtained by simulation is consistent with the region of runner
cavitation damage.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: In the present work, detailed analysis on flow characteristics of a pair of immiscible liquids through horizontal pipeline is simulated by using ANSYS FLUENT 6.2. Moderately viscous oil and water (viscosity ratio = 107, density ratio = 0.89 and interfacial tension = 0.024 N/m) have been taken as system fluids for the study. Volume of Fluid (VOF) method has been employed by assuming unsteady flow, immiscible liquid pair, constant liquid properties, and co-axial flow. Meshing has been done using GAMBIT. Quadrilateral mesh type has been chosen to account for the surface tension effect more accurately. From the grid independent study, we have selected 47037 number of mesh elements for the entire geometry. Simulation successfully predicts slug, stratified wavy, stratified mixed and annular flow, except dispersion of oil in water, and dispersion of water in oil. Simulation results are validated with horizontal literature data and good conformity is observed. Subsequently, we have simulated the hydrodynamics (viz., velocity profile, area average pressure across a cross section and volume fraction profile along the radius) of stratified wavy and annular flow at different phase velocities. The simulation results show that in the annular flow, total pressure of the mixture decreases with increase in oil velocity due to the fact that pipe cross section is completely wetted with water. Simulated oil volume fraction shows maximum at the centre in core annular flow, whereas, in stratified flow, maximum value appears at upper side of the pipeline. These results are in accord with the actual flow configuration. Our findings could be useful in designing pipeline for transportation of crude oil.
Abstract: An appropriate model to predict the size of the droplets
resulting from the break-up with the structures will help in a better
understanding and modeling of the two-phase flow calculations in the
simulation of a reactor core loss-of-coolant accident (LOCA). A
droplet behavior impacting on a hot surface above the Leidenfrost
temperature was investigated. Droplets of known size and velocity
were impacted to an inclined plate of hot temperature, and the
behavior of the droplets was observed by a high-speed camera. It was
found that for droplets of Weber number higher than a certain value,
the higher the Weber number of the droplet the smaller the secondary
droplets. The COBRA-TF model over-predicted the measured
secondary droplet sizes obtained by the present experiment. A simple
model for the secondary droplet size was proposed using the mass
conservation equation. The maximum spreading diameter of the
droplets was also compared to previous correlations and a fairly good
agreement was found. A better prediction of the heat transfer in the
case of LOCA can be obtained with the presented model.
Abstract: The atomization effect is an important factor of the heat transfer of liquid nitrogen spray. In this paper, two kinds of internal-mixing twin-fluid atomizers were design. According to the fracture theory and fluid mechanics, the model is established to simulate atomization effect. The results showed that: Internal-mixing atomizers, with the liquid nitrogen atomization size from 20um to 40um, have superior performance. Y-jet atomizer spray speed is greater than Multi-jet atomizer, and it can improve the efficiency of heat transfer between the liquid nitrogen and its spray object. Multi-jet atomizer atomization cone angle is about 30°, Y-jet atomizer atomization cone angle is about 20°. During atomizer selection, the size of the heat transfer area should be considered.
Abstract: In the present study, computational fluid dynamics
(CFD) simulation has been executed to investigate the transition
boundaries of different flow patterns for moderately viscous oil-water
(viscosity ratio 107, density ratio 0.89 and interfacial tension of 0.032
N/m.) two-phase flow through a horizontal pipeline with internal
diameter and length of 0.025 m and 7.16 m respectively. Volume of
Fluid (VOF) approach including effect of surface tension has been
employed to predict the flow pattern. Geometry and meshing of the
present problem has been drawn using GAMBIT and ANSYS
FLUENT has been used for simulation. A total of 47037 quadrilateral
elements are chosen for the geometry of horizontal pipeline. The
computation has been performed by assuming unsteady flow,
immiscible liquid pair, constant liquid properties, co-axial flow and a
T-junction as entry section. The simulation correctly predicts the
transition boundaries of wavy stratified to stratified mixed flow.
Other transition boundaries are yet to be simulated. Simulated data
has been validated with our own experimental results.
Abstract: In the study the influence of the physical-chemical properties of a liquid, the width of a channel gap and the superficial liquid and gas velocities on the patterns formed during two phase flows in vertical, narrow mini-channels was investigated. The research was performed in the channels of rectangular cross-section and of dimensions: 15 x 0.65 mm and 7.5 x 0.73 mm. The experimental data were compared with the published criteria of the transitions between the patterns of two-phase flows.
Abstract: The paper investigates parallel channel instabilities of
natural circulation boiling water reactor. A thermal-hydraulic model
is developed to simulate two-phase flow behavior in the natural circulation boiling water reactor (NCBWR) with the incorporation of
ex-core components and recirculation loop such as steam separator, down-comer, lower-horizontal section and upper-horizontal section
and then, numerical analysis is carried out for parallel channel
instabilities of the reactor undergoing both in-phase and out-of-phase
modes of oscillations. To analyze the relative effect on stability of the reactor due to inclusion of various ex-core components and
recirculation loop, marginal stable point is obtained at a particular inlet enthalpy of the reactor core without the inclusion of ex-core
components and recirculation loop and then with the inclusion of the
same. Numerical simulations are also conducted to determine the
relative dominance between two modes of oscillations i.e. in-phase and out-of-phase. Simulations are also carried out when the channels
are subjected to asymmetric power distribution keeping the inlet enthalpy same.