Abstract: In this study, the Reynolds-Stress-Navier-Stokes framework is utilized to investigate the flow inside the diesel injector nozzle. The flow is assumed to be multiphase as the formation of vapor by pressure drop is visualized. For pressure and velocity linkage, the coupled algorithm is used. Since the cavitation phenomenon inherently is unsteady, the quasi-steady approach is utilized for saving time and resources in the current study. Schnerr-Sauer cavitation model is used, which was capable of predicting flow behavior both at the initial and final steps of the cavitation process. Two different turbulent models were used in this study to clarify which one is more capable in predicting cavitation inception and super-cavitation. It was found that K-ε was more compatible with the Shnerr-Sauer cavitation model; therefore, the mentioned model is used for the rest of this study.
Abstract: Reducing sloshing is one of the major challenges in industries where transporting of liquid is involved. The present study investigates the sloshing effect for different liquid levels of 50% of the tank capacity. CFD simulation for two different baffle configurations has been carried out using a time-based multiphase Volume of fluid (VOF) scheme. Baffles were introduced to examine the sloshing effect inside the tank. Results were compared against the baseline case to assess the effectiveness of baffles; maximum liquid height over the period of the simulation was considered as the parameter for measuring the sloshing effect inside the tank. It was found that the addition of baffles reduced the sloshing effect inside the tank as compared to the baseline model.
Abstract: The aim of this experimental and numerical study is to analyze the effects of acoustic streaming generated by 40 kHz ultrasonic waves on heat transfer in forced convection, with and without 40 PPI aluminum metal foam. Preliminary dynamic and thermal studies were done with COMSOL Multiphase, to see heat transfer enhancement degree by inserting a 40PPI metal foam (10 × 2 × 3 cm) on a heat sink, after having determined experimentally its permeability and Forchheimer's coefficient. The results obtained numerically are in accordance with those obtained experimentally, with an enhancement factor of 205% for a velocity of 0.4 m/s compared to an empty channel. The influence of 40 kHz ultrasound on heat transfer was also tested with and without metallic foam. Results show a remarkable increase in Nusselt number in an empty channel with an enhancement factor of 37,5%, while no influence of ultrasound on heat transfer in metal foam presence.
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: A Computational Fluid Dynamics (CFD) technique has been applied to simulate the gas-liquid phase in double stirred tank of Rushton impeller. Eulerian-Eulerian model was adopted to simulate the multiphase with standard correlation of Schiller and Naumann for drag co-efficient. The turbulence was modeled by using standard k-ε turbulence model. The present CFD model predicts flow pattern, local gas hold-up, and local specific area. It also predicts local kLa (mass transfer rate) for single impeller. The predicted results were compared with experimental and CFD results of published literature. The predicted results are slightly over predicted with the experimental results; however, it is in reasonable agreement with other simulated results of published literature.
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: The purpose of this study is to present a mathematical phrase for calculating the dispersion rate of spilled oil in water column under non-breaking waves. In this regard, a multiphase numerical model is applied for which waves and oil phase were computed concurrently, and accuracy of its hydraulic calculations have been proven. More than 200 various scenarios of oil spilling in wave waters were simulated using the multiphase numerical model and its outcome were collected in a database. The recorded results were investigated to identify the major parameters affected vertical oil dispersion and finally 6 parameters were identified as main independent factors. Furthermore, some statistical tests were conducted to identify any relationship between the dependent variable (dispersed oil mass in the water column) and independent variables (water wave specifications containing height, length and wave period and spilled oil characteristics including density, viscosity and spilled oil mass). Finally, a mathematical-statistical relationship is proposed to predict dispersed oil in marine waters. To verify the proposed relationship, a laboratory example available in the literature was selected. Oil mass rate penetrated in water body computed by statistical regression was in accordance with experimental data was predicted. On this occasion, it was necessary to verify the proposed mathematical phrase. In a selected laboratory case available in the literature, mass oil rate penetrated in water body computed by suggested regression. Results showed good agreement with experimental data. The validated mathematical-statistical phrase is a useful tool for oil dispersion prediction in oil spill events in marine areas.
Abstract: In the thermochemical water splitting process by Cu-Cl cycle, oxygen gas is produced by an endothermic thermolysis process at a temperature of 530oC. Oxygen production reactor is a three-phase reactor involving cuprous chloride molten salt, copper oxychloride solid reactant and oxygen gas. To perform optimal performance, the oxygen reactor requires accurate control of heat transfer to the molten salt and decomposing solid particles within the thermolysis reactor. In this paper, the scale up analysis of the oxygen reactor that is heated by an internal helical tube is performed from the perspective of heat transfer. A heat balance of the oxygen reactor is investigated to analyze the size of the reactor that provides the required heat input for different rates of hydrogen production. It is found that the helical tube wall and the service side constitute the largest thermal resistances of the oxygen reactor system. In the analysis of this paper, the Cu-Cl cycle is assumed to be heated by two types of nuclear reactor, which are HTGR and CANDU SCWR. It is concluded that using CANDU SCWR requires more heat transfer rate by 3-4 times than that when using HTGR. The effect of the reactor aspect ratio is also studied and it is found that increasing the aspect ratio decreases the number of reactors and the rate of decrease in the number of reactors decreases by increasing the aspect ratio. Comparisons between the results of this study and pervious results of material balances in the oxygen reactor show that the size of the oxygen reactor is dominated by the heat balance rather than the material balance.
Abstract: The thermochemical copper-chlorine (Cu-Cl) cycle is considered as a sustainable and efficient technology for a hydrogen production, when linked with clean-energy systems such as nuclear reactors or solar thermal plants. In the Cu-Cl cycle, water is decomposed thermally into hydrogen and oxygen through a series of intermediate reactions. This paper investigates the thermal scale up analysis of the three phase oxygen production reactor in the Cu-Cl cycle, where the reaction is endothermic and the temperature is about 530 oC. The paper focuses on examining the size and number of oxygen reactors required to provide enough heat input for different rates of hydrogen production. The type of the multiphase reactor used in this paper is the continuous stirred tank reactor (CSTR) that is heated by a half pipe jacket. The thermal resistance of each section in the jacketed reactor system is studied to examine its effect on the heat balance of the reactor. It is found that the dominant contribution to the system thermal resistance is from the reactor wall. In the analysis, the Cu-Cl cycle is assumed to be driven by a nuclear reactor where two types of nuclear reactors are examined as the heat source to the oxygen reactor. These types are the CANDU Super Critical Water Reactor (CANDU-SCWR) and High Temperature Gas Reactor (HTGR). It is concluded that a better heat transfer rate has to be provided for CANDU-SCWR by 3-4 times than HTGR. The effect of the reactor aspect ratio is also examined in this paper and is found that increasing the aspect ratio decreases the number of reactors and the rate of decrease in the number of reactors decreases by increasing the aspect ratio. Finally, a comparison between the results of heat balance and existing results of mass balance is performed and is found that the size of the oxygen reactor is dominated by the heat balance rather than the material balance.
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: This paper presents a speed estimation scheme based
on second-order sliding-mode Super Twisting Algorithm (STA) and
Model Reference Adaptive System (MRAS) estimation theory for
Sensorless control of multiphase induction machine. A stator current
observer is designed based on the STA, which is utilized to take the
place of the reference voltage model of the standard MRAS
algorithm. The observer is insensitive to the variation of rotor
resistance and magnetizing inductance when the states arrive at the
sliding mode. Derivatives of rotor flux are obtained and designed as
the state of MRAS, thus eliminating the integration. Compared with
the first-order sliding-mode speed estimator, the proposed scheme
makes full use of the auxiliary sliding-mode surface, thus alleviating
the chattering behavior without increasing the complexity. Simulation
results show the robustness and effectiveness of the proposed
scheme.
Abstract: Multiphase Induction Machine (IM) is normally
controlled using rotor field oriented vector control. Under phase(s)
loss, the machine currents can be optimally controlled to satisfy
certain optimization criteria. In this paper we discuss the performance
of double manifold sliding mode observer (DM-SMO) in Sensorless
control of multiphase induction machine under unsymmetrical
condition (one phase loss). This observer is developed using the IM
model in the stationary reference frame. DM-SMO is constructed by
adding extra feedback term to conventional single mode sliding mode
observer (SM-SMO) which proposed in many literature. This leads to
a fully convergent observer that also yields an accurate estimate of
the speed and stator currents. It will be shown by the simulation
results that the estimated speed and currents by the method are very
well and error between real and estimated quantities is negligible.
Also parameter sensitivity analysis shows that this method is rather
robust against parameter variation.
Abstract: Liquid-Liquid Equilibrium (LLE) data are measured
for the ternary mixtures of water + 1-butanol + butyl acetate and
quaternary mixtures of water + 1-butanol + butyl acetate + glycerol at
atmospheric pressure at 313.15 K. In addition, isothermal
vapor–liquid–liquid equilibrium (VLLE) data are determined
experimentally at 333.15 K. The region of heterogeneity is found to
increase as the hydrophilic agent (glycerol) is introduced into the
aqueous mixtures. The experimental data are correlated with the
NRTL model. The predicted results from the solution model with the
model parameters determined from the constituent binaries are also
compared with the experimental values.
Abstract: The effects of flame-holder position, the ratio of flame
holder diameter to combustion chamber diameter and injection angle
on fuel propulsive droplets sizing and effective mass fraction have
been studied by a cold flow. We named the mass of fuel vapor inside
the flammability limit as the effective mass fraction. An empty
cylinder as well as a flame-holder which are a simulator for duct
combustion has been considered. The airflow comes into the cylinder
from one side and injection operation will be done by four nozzles
which are located on the entrance of cylinder. To fulfill the
calculations a modified version of KIVA-3V code which is a
transient, three-dimensional, multiphase, multi component code for
the analysis of chemically reacting flows with sprays, is used.