Abstract: We present results on the initial formation of ripples from an initially flattened erodible bed. We use direct numerical simulations (DNS) of turbulent open channel flow over a fixed sinusoidal bed coupled with hydrodynamic stability analysis. We use the direct forcing immersed boundary method to account for the presence of the sediment bed. The resolved flow provides the bed shear stress and consequently the sediment transport rate, which is needed in the stability analysis of the Exner equation. The approach is different from traditional linear stability analysis in the sense that the phase lag between the bed topology, and the sediment flux is obtained from the DNS. We ran 11 simulations at a fixed shear Reynolds number of 180, but for different sediment bed wavelengths. The analysis allows us to sweep a large range of physical and modelling parameters to predict their effects on linear growth. The Froude number appears to be the critical controlling parameter in the early linear development of ripples, in contrast with the dominant role of particle Reynolds number during the equilibrium stage.
Abstract: Efficiency of the cooling process for cryogenic
propellant boiling in engine cooling channels on space applications is
relentlessly affected by the phase change occurs during the boiling.
The effectiveness of the cooling process strongly pertains to the
type of the boiling regime such as nucleate and film. Geometric
constraints like a non-transparent cooling channel unable to use
any of visualization methods. The ultrasonic (US) technique as a
non-destructive method (NDT) has therefore been applied almost
in every engineering field for different purposes. Basically, the
discontinuities emerge between mediums like boundaries among
different phases. The sound wave emitted by the US transducer is
both transmitted and reflected through a gas-liquid interface which
makes able to detect different phases. Due to the thermal and
structural concerns, it is impractical to sustain a direct contact
between the US transducer and working fluid. Hence the transducer
should be located outside of the cooling channel which results in
additional interfaces and creates ambiguities on the applicability
of the present method. In this work, an exploratory research is
prompted so as to determine detection ability and applicability of
the US technique on the cryogenic boiling process for a cooling
cycle where the US transducer is taken place outside of the channel.
Boiling of the cryogenics is a complex phenomenon which mainly
brings several hindrances for experimental protocol because of
thermal properties. Thus substitute materials are purposefully selected
based on such parameters to simplify experiments. Aside from
that, nucleate and film boiling regimes emerging during the boiling
process are simply simulated using non-deformable stainless steel
balls, air-bubble injection apparatuses and air clearances instead
of conducting a real-time boiling process. A versatile detection
algorithm is perennially developed concerning exploratory studies
afterward. According to the algorithm developed, the phases can be
distinguished 99% as no-phase, air-bubble, and air-film presences.
The results show the detection ability and applicability of the US
technique for an exploratory purpose.
Abstract: A semi-implicit phase field method for droplet evolution is proposed. Using the phase field Cahn-Hilliard equation, we are able to track the interface in multiphase flow. The idea of a semi-implicit finite difference scheme is reviewed and employed to solve two nonlinear equations, including the Navier-Stokes and the Cahn-Hilliard equations. The use of a semi-implicit method allows us to have larger time steps compared to explicit schemes. The governing equations are coupled and then solved by a GMRES solver (generalized minimal residual method) using modified Gram-Schmidt orthogonalization. To show the validity of the method, we apply the method to the simulation of a rising droplet, a leaky dielectric drop and the coalescence of drops. The numerical solutions to the phase field model match well with existing solutions over a defined range of variables.
Abstract: For the past decades, CO2 flooding has been used as a successful method for enhanced oil recovery (EOR). However, high mobility ratio and fingering effect are considered as important drawbacka of this process. Low temperature injection of CO2 into high temperature reservoirs may improve the oil recovery, but simulating multiphase flow in the non-isothermal medium is difficult, and commercial simulators are very unstable in these conditions. Furthermore, to best of authors’ knowledge, no experimental work was done to verify the results of the simulations and to understand the pore-scale process. In this paper, we present results of investigations on injection of low temperature CO2 into a high-pressure high-temperature micromodel with injection temperature range from 34 to 75 °F. Effect of temperature and saturation changes of different fluids are measured in each case. The results prove the proposed method. The injection of CO2 at low temperatures increased the oil recovery in high temperature reservoirs significantly. Also, CO2 rich phases available in the high temperature system can affect the oil recovery through the better sweep of the oil which is initially caused by penetration of LCO2 inside the system. Furthermore, no unfavorable effect was detected using this method. Low temperature CO2 is proposed to be used as early as secondary recovery.
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: In this paper, the average heat transfer characteristics
for a cross flow cylinder of 16 mm diameter in a vertical pipe has
been studied for single-phase flow (water/oil) and multicomponent
(non-boiling) flow (water-air, water-oil, oil-air and water-oil-air). The
cylinder is uniformly heated by electrical heater placed at the centre
of the element. The results show that the values of average heat
transfer coefficients for water are around four times the values for oil
flow. Introducing air as a second phase with water has very little
effect on heat transfer rate, while the heat transfer increased by 70%
in case of oil. For water–oil flow, the heat transfer coefficient values
are reflecting the percentage of water up to 50%, but increasing the
water more than 50% leads to a sharp increase in the heat transfer
coefficients to become close to the values of pure water. The
enhancement of heat transfer by mixing two phases may be attributed
to the changes in flow structure near to cylinder surface which lead to
thinner boundary layer and higher turbulence. For three-phase flow,
the heat transfer coefficients for all cases fall within the limit of
single-phase flow of water and oil and are very close to pure water
values. The net effect of the turbulence augmentation due to the
introduction of air and the attenuation due to the introduction of oil
leads to a thinner boundary layer of oil over the cylinder surface
covered by a mixture of water and air bubbles.
Abstract: The offshore facilities condition of Pertamina Hulu Energi Offshore North West Java (PHE ONWJ) varies greatly from place to place, depending on the characteristics of the presently installed facilities. In some locations, such as ESA platform, gas trap is mainly caused by the occurrence of flash gas phenomenon which is known as mechanical-physical separation process of multiphase flow. Consequently, the presence of gas trap at main oil line would accumulate on certain areas result in a reduced oil stream throughout the pipeline. Any presence of discrete gaseous along continuous oil flow represents a unique flow condition under certain specific volume fraction and velocity field. From gas lift source, a benefit line is used as a motive flow for ejector which is designed to generate a syphon effect to minimize the gas trap phenomenon. Therefore, the ejector’s exhaust stream will flow to the designated point without interfering other systems.
Abstract: This paper presents a Euler-Lagrange model of the water-particles multiphase flows in a Vorsyl separator where particles with different densities are separated. A series of particles with their densities ranging from 760 kg/m3 to 1380 kg/m3 were fed into the Vorsyl separator with water by means of tangential inlet. The simulation showed that the feed materials acquired centrifugal force which allows most portion of the particles with a density less than water to move to the center of the separator, enter the vortex finder and leave the separator through the bottom outlet. While the particles heavier than water move to the wall, reach the throat area and leave the separator through the side outlet. The particles were thus separated and particles collected at the bottom outlet are pure and clean. The influence of particle density on separation efficiency was investigated which demonstrated a positive correlation of the separation efficiency with increasing density difference between medium liquid and the particle. In addition, the influence of the split ratio on the performance was studied which showed that the separation efficiency of the Vorsyl separator can be improved by the increase of split ratio. The simulation also suggested that the Vorsyl separator may not function when the feeding velocity is smaller than a certain critical feeding in velocity. In addition, an increasing feeding velocity gives rise to increased pressure drop, however does not necessarily increase the separation efficiency.
Abstract: Relative permeabilities are practical factors that are used to correct the single phase Darcy’s law for application to multiphase flow. For effective characterisation of large-scale multiphase flow in hydrocarbon recovery, relative permeability and capillary pressures are used. These parameters are acquired via special core flooding experiments. Special core analysis (SCAL) module of reservoir simulation is applied by engineers for the evaluation of these parameters. But, core flooding experiments in shale core sample are expensive and time consuming before various flow assumptions are achieved for instance Darcy’s law. This makes it imperative for the application of coreflooding simulations in which various analysis of relative permeabilities and capillary pressures of multiphase flow can be carried out efficiently and effectively at a relative pace. This paper presents a Sendra software simulation of core flooding to achieve to relative permeabilities and capillary pressures using different correlations. The approach used in this study was three steps. The first step, the basic petrophysical parameters of Marcellus shale sample such as porosity was determined using laboratory techniques. Secondly, core flooding was simulated for particular scenario of injection using different correlations. And thirdly the best fit correlations for the estimation of relative permeability and capillary pressure was obtained. This research approach saves cost and time and very reliable in the computation of relative permeability and capillary pressures at steady or unsteady state, drainage or imbibition processes in oil and gas industry when compared to other methods.
Abstract: One of the most important parameters in petroleum reservoirs is the pressure distribution along the reservoir, as the pressure varies with the time and location. A popular method to determine the pressure distribution in a reservoir in the unsteady state regime of flow is applying Darcy’s equation and solving this equation numerically. The numerical simulation of reservoirs is based on these numerical solutions of different partial differential equations (PDEs) representing the multiphase flow of fluids. Pressure profile has obtained in a one dimensional system solving Darcy’s equation explicitly. Changes of pressure profile in three situations are investigated in this work. These situations include section length changes, step time changes and time approach to infinity. The effects of these changes in pressure profile are shown and discussed in the paper.
Abstract: The performance of high-resolution schemes is investigated for unsteady, inviscid and compressible multiphase flows. An Eulerian diffuse interface approach has been chosen for the simulation of multicomponent flow problems. The reduced fiveequation and seven equation models are used with HLL and HLLC approximation. The authors demonstrated the advantages and disadvantages of both seven equations and five equations models studying their performance with HLL and HLLC algorithms on simple test case. The seven equation model is based on two pressure, two velocity concept of Baer–Nunziato [10], while five equation model is based on the mixture velocity and pressure. The numerical evaluations of two variants of Riemann solvers have been conducted for the classical one-dimensional air-water shock tube and compared with analytical solution for error analysis.
Abstract: Downward turbulent bubbly flows in pipes were
modeled using computational fluid dynamics tools. The
Hydrodynamics, phase distribution and turbulent structure of twophase
air-water flow in a 57.15 mm diameter and 3.06 m length
vertical pipe was modeled by using the 3-D Eulerian-Eulerian
multiphase flow approach. Void fraction, liquid velocity and
turbulent fluctuations profiles were calculated and compared against
experimental data. CFD results are in good agreement with
experimental data.
Abstract: This paper presents an experimental investigation using Acoustic Emission (AE) technology to monitor sand transportation in multiphase flow. The investigations were undertaken on three-phase (air-water-sand) flow in a horizontal pipe where the superficial gas velocity (VSG) had a range of between 0.2msˉ¹ to 2.0msˉ¹ and superficial liquid velocity (VSL) had a range of between 0.2msˉ¹ to 1.0msˉ¹. The experimental findings clearly show a correlation exists between AE energy levels, sand concentration, superficial gas velocity (VSG), and superficial liquid velocity (VSL).
Abstract: A robust AUSM+ upwind discretisation scheme has been developed to simulate multiphase flow using consistent spatial discretisation schemes and a modified low-Mach number diffusion term. The impact of the selection of an interfacial pressure model has also been investigated. Three representative test cases have been simulated to evaluate the accuracy of the commonly-used stiffenedgas equation of state with respect to the IAPWS-IF97 equation of state for water. The algorithm demonstrates a combination of robustness and accuracy over a range of flow conditions, with the stiffened-gas equation tending to overestimate liquid temperature and density profiles.
Abstract: Most CT reconstruction system x-ray computed
tomography (CT) is a well established visualization technique in
medicine and nondestructive testing. However, since CT scanning
requires sampling of radiographic projections from different viewing
angles, common CT systems with mechanically moving parts are too
slow for dynamic imaging, for instance of multiphase flows or live
animals. A large number of X-ray projections are needed to
reconstruct CT images, so the collection and calculation of the
projection data consume too much time and harmful for patient. For
the purpose of solving the problem, in this study, we proposed a
method for tomographic reconstruction of a sample from a limited
number of x-ray projections by using linear interpolation method. In
simulation, we presented reconstruction from an experimental x-ray
CT scan of a Aluminum phantom that follows to two steps: X-ray
projections will be interpolated using linear interpolation method and
using it for CT reconstruction based upon Ordered Subsets
Expectation Maximization (OSEM) method.
Abstract: In this research the separation efficiency of deoiling hydrocyclone is evaluated using three-dimensional simulation of multiphase flow based on Eulerian-Eulerian finite volume method. The mixture approach of Reynolds Stress Model is also employed to capture the features of turbulent multiphase swirling flow. The obtained separation efficiency of Colman's design is compared with available experimental data and showed that the separation curve of deoiling hydrocyclones can be predicted using numerical simulation.
Abstract: Multiphase flow transport in porous medium is very common and significant in science and engineering applications. For example, in CO2 Storage and Enhanced Oil Recovery processes, CO2 has to be delivered to the pore spaces in reservoirs and aquifers. CO2 storage and enhance oil recovery are actually displacement processes, in which oil or water is displaced by CO2. This displacement is controlled by pore size, chemical and physical properties of pore surfaces and fluids, and also pore wettability. In this study, a technique was developed to measure the pressure profile for driving gas/liquid to displace water in pores. Through this pressure profile, the impact of pore size on the multiphase flow transport and displacement can be analyzed. The other rig developed can be used to measure the static and dynamic pore wettability and investigate the effects of pore size, surface tension, viscosity and chemical structure of liquids on pore wettability.
Abstract: This paper reports the results of an experimental study
conducted to characterise the gas-liquid multiphase flows
experienced within a vertical riser transporting a range of gas-liquid
flow rates. The scale experiments were performed using an
air/silicone oil mixture within a 6 m long riser. The superficial air
velocities studied ranged from 0.047 to 2.836 m/ s, whilst
maintaining a liquid superficial velocity at 0.047 m/ s. Measurements
of the mean cross-sectional and time average radial void fraction
were obtained using a wire mesh sensor (WMS). The data were
recorded at an acquisition frequency of 1000 Hz over an interval of
60 seconds. For the range of flow conditions studied, the average
void fraction was observed to vary between 0.1 and 0.9. An analysis
of the data collected concluded that the observed void fraction was
strongly affected by the superficial gas velocity, whereby the higher
the superficial gas velocity, the higher was the observed average void
fraction. The average void fraction distributions observed were in
good agreement with the results obtained by other researchers. When
the air-silicone oil flows were fully developed reasonably symmetric
profiles were observed, with the shape of the symmetry profile being
strongly dependent on the superficial gas velocity.
Abstract: Considering non-ideal behavior of fluids and its effects on hydrodynamic and mass transfer in multiphase flow is very essential. Simulations were performed that takes into account the effects of mass transfer and mixture non-ideality on hydrodynamics reported by Irani et al. In this paper, by assuming the density of phases to be constant and Raullt-s law instead of using EOS and fugacity coefficient definition, respectively for both the liquid and gas phases, the importance of non-ideality effects on mass transfer and hydrodynamic behavior was studied. The results for a system of octane/propane (T=323 K, P =445 kpa) also indicated that the assumption of constant density in simulation had major role to diverse from experimental data. Furthermore, comparison between obtained results and the previous report indicated significant differences between experimental data and simulation results with more ideal assumptions.