Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

Coalescence of secondary dispersions in fibrous beds

SUMMARY A study has been made of the coalescence of secondary dispersions in a fibrous bed. The literature pertaining to the formation, hydrodynamic behaviour and methods of separation of droplets less than one hundred micrometres in diameter has been reviewed with particular reference to fibrous bed coalescers. The main operating parameters were identified as inlet drop size distribution, phase ratio, superficial velocity and the thickness and voidage of the bed . A recirculatory rig with interchangeable fibrous bed pads was designed and operated with toluene-water dispersions generated by a combination of centrifugal pumps . Inlet drop sizes were analysed using a Coulter Counter and outlet drops were sized photographically. A novel technique, involving conductivity measur ements at different planes in the bed, was developed to measure hold up distribution. Single phase flow and two phase flow pressure drops were correlated by a Blake-Kozeny type equation. Exit drop size was independent of inlet drop size distribution and phase ratio but a function of superficialvelocity and packing thickness. Average bed hold up was independent of inlet drop size distribution and phase ratio, but decreased with increase in superficial velocity. Hold up was not evenly distributed in the bed, the highest value occurred at the inlet followed by a sharp -2 drop at approximately 1.2 x 10 m. Hold up remained constant throughout the rest of the bed until the exit plane, where it increased. From the results, a mechanism is postulated involving: (a) Capture of the inlet drops followed by interdrop coalescence until an equilibrium value is reached. (b) Equilibrium size droplets flowing as rivulets through the intermediate portion of the bed, and (c) Each rivulet forms droplets at the exit face, which detach by a 'drip point' mechanism.

Coalescence of secondary dispersions

A study has been made of the coalescence of secondary dispersions in beds of woven meshes. The variables investigated were superficial velocity, bed depth, mesh geometry and fibre material; the effects of presoaking the bed in the dispersed phase before operation were also considered. Equipment was design~d to generate a 0.1% phase ratio toluene in water dispersion whose mean drop size was determined using a Coulter Counter. The coalesced drops were sized by photography and a novel holographic technique was developed to evaluate the mean diameter of the effluent secondary drops. Previous models describing single phase flow in porous media are reviewed and it was found that the experimental data obtained in this study is best represented by Keller's equation which is based on a physical model similar to the internal structure of the meshes. Statistical analysis of two phase data produced a correlation, for each mesh tested, relating the pressure drop to superficial velocity and bed depth. The flow parameter evaluated from the single phase model is incorporated into a theoretical comparison of drop capture mechanisms which indicated that direct and indirect interception are predominant. The resulting equation for drop capture efficiericy is used to predict the initial, local drop capture rate in a coalescer. A mathematical description of the saturation profiles was formulated and verified by average saturation data. Based 6n the Blake-Kozeny equation, an expression is derived analytically to predict the two phase pressure drop using the parameters which characterise the saturation profiles. By specifying the local saturation at the inlet face for a given velocity, good agreement between experimental pressure drop data and the model predictions was obtained.

Coalescence of secondary dispersions in packed beds

A study has been made of the coalescence of secondary dispersions in beds of monosized glass ballotini. The variables investigated were superficial velocity, bed depth, ballotini size and dispersed phase concentration. Equipment was designed to generate a toluene ln water dispersion with phase ratios from 0.1 - 1.0 v/v % and whose mean drop size was determined using a Coulter Counter. The coalesced drops were sized by photography and the mean diameter of the effluent drops was determined using a Malvern Particle Size Analyser. Previous models describing single phase flow in porous media are reviewed and it was found that the experimental data obtained in this study is best represented by the Carman-Kozeny equations. Relative permeability correlations were used to predict the saturation profiles across the bed from measured two phase pressure drop data. Theoretical comparison of drop capture mechanisms indicated that direct and indirect interception are predominant. The total capture efficiency for the bed can also be evaluated using Spielman and Fitzpatrick's correlation.The resulting equation is used to predict the initial, local drop capture rate in a coalescer. A mathematical description of the saturation profiles is formulated and verified by the saturation profiles obtained by relative permeability. Based on the Carman-Kozeny equation, an expression is derived analytically to .predict the two phase pressure drop using the parameters which characterise the saturation profiles. By specifying the local saturation at the inlet face for a given velocity and phase ratio, good agreement between experimental pressure drop data and the model predictions was obtained. An attempt to predict the exit drop size has been made using an analogy for flow through non cylindrical channels.

The separation of secondary oil-water dispersions in particulate beds

The literature relating to haze formation, methods of separation, coalescence mechanisms, and models by which droplets <100 μm are collected, coalesced and transferred, have been reviewed with particular reference to particulate bed coalescers. The separation of secondary oil-water dispersions was studied experimentally using packed beds of monosized glass ballotini particles. The variables investigated were superficial velocity, bed depth, particle size, and the phase ratio and drop size distribution of inlet secondary dispersion. A modified pump loop was used to generate secondary dispersions of toluene or Clairsol 350 in water with phase ratios between 0.5-6.0 v/v%.Inlet drop size distributions were determined using a Malvern Particle Size Analyser;effluent, coalesced droplets were sized by photography. Single phase flow pressure drop data were correlated by means of a Carman-Kozeny type equation. Correlations were obtained relating single and two phase pressure drops, as (ΔP2/μc)/ΔP1/μd) = kp Ua Lb dcc dpd Cine A flow equation was derived to correlate the two phase pressure drop data as, ΔP2/(ρcU2) = 8.64*107 [dc/D]-0.27 [L/D]0.71 [dp/D]-0.17 [NRe]1.5 [e1]-0.14 [Cin]0.26  In a comparison between functions to characterise the inlet drop size distributions a modification of the Weibull function provided the best fit of experimental data. The general mean drop diameter was correlated by: q_p q_p p_q /β      Γ ((q-3/β) +1) d qp = d fr  .α        Γ ((P-3/β +1 The measured and predicted mean inlet drop diameters agreed within ±15%. Secondary dispersion separation depends largely upon drop capture within a bed. A theoretical analysis of drop capture mechanisms in this work indicated that indirect interception and London-van der Waal's mechanisms predominate. Mathematical models of dispersed phase concentration m the bed were developed by considering drop motion to be analogous to molecular diffusion.The number of possible channels in a bed was predicted from a model in which the pores comprised randomly-interconnected passage-ways between adjacent packing elements and axial flow occured in cylinders on an equilateral triangular pitch. An expression was derived for length of service channels in a queuing system leading to the prediction of filter coefficients. The insight provided into the mechanisms of drop collection and travel, and the correlations of operating parameters, should assist design of industrial particulate bed coalescers.

Estudo experimental sobre o comportamento do escoamento bifásico água - ar para a medição de vazão em placa de orifício

The measurement of flow through the prediction of differential pressure is widely used in industrial day-to-day, this happens mainly due to the fact that it is used for various types of fluids, such as gas flow and liquid with viscosity distinct even flow of fluids with particles in suspension. The suitability of this equipment for measuring mass flow in two-phase flow is of paramount importance for technological development and reliability of results. When it comes to two-phase flow the relationship between the fluids and their interactions are of paramount importance in predicting the flow. In this paper, we propose the use of concentric orifice plate used in small diameter pipes of 25.4 mm order where a two-phase flow flows between water-air. The measurement of single-phase flow was made with the use of data in NBR 5167-1 which was used to Stolz equation for measuring discharge coefficient. In the two-phase flow was used two correlations widely used in the prognosis of mass flow, the pattern of Zhang (1992) and the model of Chisholm (1967), to the homogeneous flow model. It was observed that the behavior found in Zhang model are consistent more realistic way the mass flow of two-phase flow, since the model Chisholm extrapolate the parameters for the downstream pressure P2, the orifice plate, and the rated discharge coefficient. The use of the change in pressure drop P1-P2 and discharge coefficient, led to a better convergence of the values obtained for the two-phase air-water stream.

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle Hows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained.

Non-Linear Scale Interactions in a Forced Turbulent Boundary Layer

This thesis explores the dynamics of scale interactions in a turbulent boundary layer through a forcing-response type experimental study. An emphasis is placed on the analysis of triadic wavenumber interactions since the governing Navier-Stokes equations for the flow necessitate a direct coupling between triadically consist scales. Two sets of experiments were performed in which deterministic disturbances were introduced into the flow using a spatially-impulsive dynamic wall perturbation. Hotwire anemometry was employed to measure the downstream turbulent velocity and study the flow response to the external forcing. In the first set of experiments, which were based on a recent investigation of dynamic forcing effects in a turbulent boundary layer, a 2D (spanwise constant) spatio-temporal normal mode was excited in the flow; the streamwise length and time scales of the synthetic mode roughly correspond to the very-large-scale-motions (VLSM) found naturally in canonical flows. Correlation studies between the large- and small-scale velocity signals reveal an alteration of the natural phase relations between scales by the synthetic mode. In particular, a strong phase-locking or organizing effect is seen on directly coupled small-scales through triadic interactions. Having characterized the bulk influence of a single energetic mode on the flow dynamics, a second set of experiments aimed at isolating specific triadic interactions was performed. Two distinct 2D large-scale normal modes were excited in the flow, and the response at the corresponding sum and difference wavenumbers was isolated from the turbulent signals. Results from this experiment serve as an unique demonstration of direct non-linear interactions in a fully turbulent wall-bounded flow, and allow for examination of phase relationships involving specific interacting scales. A direct connection is also made to the Navier-Stokes resolvent operator framework developed in recent literature. Results and analysis from the present work offer insights into the dynamical structure of wall turbulence, and have interesting implications for design of practical turbulence manipulation or control strategies.

Pressure drops of single and two-phase flows through T-type microchannel mixers

In this paper, preliminary experimental results are presented on pressure drop characteristics of single and two-phase flows through two T-type rectangular microchannel mixers with hydraulic diameters of 528 and 333 mum, respectively. It is shown that both N-2 and water single-phase laminar flows in microchannels, with consideration of experimental uncertainties, are consistent with classic theory, if additional effects, such as entrance effects that will interfere with the interpretation of experimental results, are eliminated by carefully designing the experiments. The obtained pressure drop data of N-2-water two-phase flow in micromixers are analyzed and compared with existing flow pattern-independent models. It is found that the Lockhart-Martinelli method generally underpredicts the frictional pressure drop. Thereafter, a modified correlation of C value in the Chisholm's equation based on linear regression of experimental data is proposed to provide a better prediction of the two-phase frictional pressure drop. Also among the homogeneous flow models investigated, the viscosity correlation of McAdams indicates the best performance in correlating the frictional pressure drop data (mean deviations within +/-20% for two micromixers both). Finally it is suggested that systematic studies are still required to accurately predict two-phase frictional performance in microchannels. (C) 2004 Elsevier B.V. All rights reserved.

Continuous flow hydroformylation using supported ionic liquid phase catalysts with carbon dioxide as a carrier.

A supported ionic liquid phase (SILP) catalyst prepared from [PrMIM][Ph2P(3-C6H4SO3)] (PrMIM = 1-propyl-3-methylimidazolium), [Rh(CO)(2)(acac)] (acacH = 2,4-pentanedione) [OctMIM]NTf2 (OctMIM = 1-n-octyl-3-methylimidazolium, Tf = CF3SO2) and microporous silica has been used for the continuous flow hydroformylation of 1-octene in the presence of compressed CO2. Statistical experimental design was used to show that the reaction rate is neither much affected by the film thickness (IL loading) nor by the syngas: substrate ratio. However, a factor-dependent interaction between the syngas: substrate ratio and film thickness on the reaction rate was revealed. Increasing the substrate flow led to increased reaction rates but lower overall yields. One of the most important parameters proved to be the phase behaviour of the mobile phase, which was studied by varying the reaction pressure. At low CO2 pressures or when N-2 was used instead of CO2 rates were low because of poor gas diffusion to the catalytic sites in the SILP. Furthermore, leaching of IL and Rh was high because the substrate is liquid and the IL had been designed to dissolve in it. As the CO2 pressure was increased, the reaction rate increased and the IL and Rh leaching were reduced, because an expanded liquid phase developed. Due to its lower viscosity the expanded liquid allows better transport of gases to the catalyst and is a poorer solvent for the IL and the catalyst because of its reduced polarity. Above 100 bar (close to the transition to a single phase at 106 bar), the rate of reaction dropped again with increasing pressure because the flowing phase becomes a better and better solvent for the alkene, reducing its partitioning into the IL film. Under optimised conditions, the catalyst was shown to be stable over at least 40 h of continuous catalysis with a steady state turnover frequency (TOF, mol product (mol Rh)(-1)) of 500 h(-1) at low Rh leaching (0.2 ppm). The selectivity of the catalyst was not much affected by the variation of process parameters. The linear: branched (1:b) ratios were ca. 3, similar to that obtained using the very same catalyst in conventional organic solvents.

Flow patterns and pressure drop of ionic liquid-water two-phase flows in microchannels

The two-phase flow of a hydrophobic ionic liquid and water was studied in capillaries made of three different materials (two types of Teflon, FEP and Tefzel, and glass) with sizes between 200µm and 270µm. The ionic liquid was 1-butyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}amide, with density and viscosity of 1420kgm and 0.041kgms, respectively. Flow patterns and pressure drop were measured for two inlet configurations (T- and Y-junction), for total flow rates of 0.065-214.9cmh and ionic liquid volume fractions from 0.05 to 0.8. The continuous phase in the glass capillary depended on the fluid that initially filled the channel. When water was introduced first, it became the continuous phase with the ionic liquid forming plugs or a mixture of plugs and drops within it. In the Teflon microchannels, the order that fluids were introduced did not affect the results and the ionic liquid was always the continuous phase. The main patterns observed were annular, plug, and drop flow. Pressure drop in the Teflon microchannels at a constant ionic liquid flow rate, was found to increase as the ionic liquid volume fraction decreased, and was always higher than the single phase ionic liquid value at the same flow rate as in the two-phase mixture. However, in the glass microchannel during plug flow with water as the continuous phase, pressure drop for a constant ionic liquid flow rate was always lower than the single phase ionic liquid value. A modified plug flow pressure drop model using a correlation for film thickness derived for the current fluids pair showed very good agreement with the experimental data. © 2013 Elsevier Ltd.

Distributed generation modelling for unbalanced three-phase power flow calculations in smart grids

Distributed Generators (DG) are generally modeled as PQ or PV buses in power flow studies. But in order to integrate DG units into the distribution systems and control the reactive power injection it is necessary to know the operation mode and the type of connection to the system. This paper presents a single-phase and a three-phase mathematical model to integrate DG in power flow calculations in distribution systems, especially suited for Smart Grid calculations. If the DG is in PV mode, each step of the power flow algorithm calculates the reactive power injection from the DG to the system to keep the voltage in the bus in a predefined level, if the DG is in PQ mode, the power injection is considered as a negative load. The method is tested on two well known test system, presenting single-phase results on 85 bus system, and three-phase results in the IEEE 34 bus test system. © 2011 IEEE.

Three dimensional finite element simulation of polymer melting and flow in a single-screw extruder : optimization of screw channel geometry

Single-screw extrusion is one of the widely used processing methods in plastics industry, which was the third largest manufacturing industry in the United States in 2007 [5]. In order to optimize the single-screw extrusion process, tremendous efforts have been devoted for development of accurate models in the last fifty years, especially for polymer melting in screw extruders. This has led to a good qualitative understanding of the melting process; however, quantitative predictions of melting from various models often have a large error in comparison to the experimental data. Thus, even nowadays, process parameters and the geometry of the extruder channel for the single-screw extrusion are determined by trial and error. Since new polymers are developed frequently, finding the optimum parameters to extrude these polymers by trial and error is costly and time consuming. In order to reduce the time and experimental work required for optimizing the process parameters and the geometry of the extruder channel for a given polymer, the main goal of this research was to perform a coordinated experimental and numerical investigation of melting in screw extrusion. In this work, a full three-dimensional finite element simulation of the two-phase flow in the melting and metering zones of a single-screw extruder was performed by solving the conservation equations for mass, momentum, and energy. The only attempt for such a three-dimensional simulation of melting in screw extruder was more than twenty years back. However, that work had only a limited success because of the capability of computers and mathematical algorithms available at that time. The dramatic improvement of computational power and mathematical knowledge now make it possible to run full 3-D simulations of two-phase flow in single-screw extruders on a desktop PC. In order to verify the numerical predictions from the full 3-D simulations of two-phase flow in single-screw extruders, a detailed experimental study was performed. This experimental study included Maddock screw-freezing experiments, Screw Simulator experiments and material characterization experiments. Maddock screw-freezing experiments were performed in order to visualize the melting profile along the single-screw extruder channel with different screw geometry configurations. These melting profiles were compared with the simulation results. Screw Simulator experiments were performed to collect the shear stress and melting flux data for various polymers. Cone and plate viscometer experiments were performed to obtain the shear viscosity data which is needed in the simulations. An optimization code was developed to optimize two screw geometry parameters, namely, screw lead (pitch) and depth in the metering section of a single-screw extruder, such that the output rate of the extruder was maximized without exceeding the maximum temperature value specified at the exit of the extruder. This optimization code used a mesh partitioning technique in order to obtain the flow domain. The simulations in this flow domain was performed using the code developed to simulate the two-phase flow in single-screw extruders.

Heat transfer performance investigation of nanofluids flow in pipe

Different types of base fluids, such as water, engine oil, kerosene, ethanol, methanol, ethylene glycol etc. are usually used to increase the heat transfer performance in many engineering applications. But these conventional heat transfer fluids have often several limitations. One of those major limitations is that the thermal conductivity of each of these base fluids is very low and this results a lower heat transfer rate in thermal engineering systems. Such limitation also affects the performance of different equipments used in different heat transfer process industries. To overcome such an important drawback, researchers over the years have considered a new generation heat transfer fluid, simply known as nanofluid with higher thermal conductivity. This new generation heat transfer fluid is a mixture of nanometre-size particles and different base fluids. Different researchers suggest that adding spherical or cylindrical shape of uniform/non-uniform nanoparticles into a base fluid can remarkably increase the thermal conductivity of nanofluid. Such augmentation of thermal conductivity could play a more significant role in enhancing the heat transfer rate than that of the base fluid. Nanoparticles diameters used in nanofluid are usually considered to be less than or equal to 100 nm and the nanoparticles concentration usually varies from 5% to 10%. Different researchers mentioned that the smaller nanoparticles concentration with size diameter of 100 nm could enhance the heat transfer rate more significantly compared to that of base fluids. But it is not obvious what effect it will have on the heat transfer performance when nanofluids contain small size nanoparticles of less than 100 nm with different concentrations. Besides, the effect of static and moving nanoparticles on the heat transfer of nanofluid is not known too. The idea of moving nanoparticles brings the effect of Brownian motion of nanoparticles on the heat transfer. The aim of this work is, therefore, to investigate the heat transfer performance of nanofluid using a combination of smaller size of nanoparticles with different concentrations considering the Brownian motion of nanoparticles. A horizontal pipe has been considered as a physical system within which the above mentioned nanofluid performances are investigated under transition to turbulent flow conditions. Three different types of numerical models, such as single phase model, Eulerian-Eulerian multi-phase mixture model and Eulerian-Lagrangian discrete phase model have been used while investigating the performance of nanofluids. The most commonly used model is single phase model which is based on the assumption that nanofluids behave like a conventional fluid. The other two models are used when the interaction between solid and fluid particles is considered. However, two different phases, such as fluid and solid phases is also considered in the Eulerian-Eulerian multi-phase mixture model. Thus, these phases create a fluid-solid mixture. But, two phases in the Eulerian-Lagrangian discrete phase model are independent. One of them is a solid phase and the other one is a fluid phase. In addition, RANS (Reynolds Average Navier Stokes) based Standard κ-ω and SST κ-ω transitional models have been used for the simulation of transitional flow. While the RANS based Standard κ-ϵ, Realizable κ-ϵ and RNG κ-ϵ turbulent models are used for the simulation of turbulent flow. Hydrodynamic as well as temperature behaviour of transition to turbulent flows of nanofluids through the horizontal pipe is studied under a uniform heat flux boundary condition applied to the wall with temperature dependent thermo-physical properties for both water and nanofluids. Numerical results characterising the performances of velocity and temperature fields are presented in terms of velocity and temperature contours, turbulent kinetic energy contours, surface temperature, local and average Nusselt numbers, Darcy friction factor, thermal performance factor and total entropy generation. New correlations are also proposed for the calculation of average Nusselt number for both the single and multi-phase models. Result reveals that the combination of small size of nanoparticles and higher nanoparticles concentrations with the Brownian motion of nanoparticles shows higher heat transfer enhancement and thermal performance factor than those of water. Literature suggests that the use of nanofluids flow in an inclined pipe at transition to turbulent regimes has been ignored despite its significance in real-life applications. Therefore, a particular investigation has been carried out in this thesis with a view to understand the heat transfer behaviour and performance of an inclined pipe under transition flow condition. It is found that the heat transfer rate decreases with the increase of a pipe inclination angle. Also, a higher heat transfer rate is found for a horizontal pipe under forced convection than that of an inclined pipe under mixed convection.

Observations on transition in plane bubble plumes

Flow visualization studies of plane laminar bubble plumes have been conducted to yield quantitative data on transition height, wavelength and wave velocity of the most unstable disturbance leading to transition. These are believed to be the first results of this kind. Most earlier studies are restricted to turbulent bubble plumes. In the present study, the bubble plumes were generated by electrolysis of water and hence very fine control over bubble size distribution and gas flow rate was possible to enable studies with laminar bubble plumes. Present observations show that (a) the dominant mode of instability in plane bubble plumes is the sinuous mode, (b) transition height and wavelength are related linearly with the proportionality constant being about 4, (c) wave velocity is about 40 % of the mean plume velocity, and (d) normalized transition height data correlate very well with a source Grashof number. Some agreement and some differences in transition characteristics of bubble plumes have been observed compared to those for similar single-phase flows.