A probabilistic approach to assess hydrate formation and design preventive measures

Formation of hydrates is one of the major flow assurance problems faced by the oil and gas industry. Hydrates tend to form in natural gas pipelines with the presence of water and favorable temperature and pressure conditions, generally low temperatures and corresponding high pressures. Agglomeration of hydrates can result in blockage of flowlines and equipment, which can be time consuming to remove in subsea equipment and cause safety issues. Natural gas pipelines are more susceptible to burst and explosion owing to hydrate plugging. Therefore, a rigorous risk-assessment related to hydrate formation is required, which assists in preventing hydrate blockage and ensuring equipment integrity. This thesis presents a novel methodology to assess the probability of hydrate formation and presents a risk-based approach to determine the parameters of winterization schemes to avoid hydrate formation in natural gas pipelines operating in Arctic conditions. It also presents a lab-scale multiphase flow loop to study the effects of geometric and hydrodynamic parameters on hydrate formation and discusses the effects of geometric and hydrodynamic parameters on multiphase development length of a pipeline. Therefore, this study substantially contributes to the assessment of probability of hydrate formation and the decision making process of winterization strategies to prevent hydrate formation in Arctic conditions.

Numerical Investigation on the Heat Transfer Enhancement Using Micro/Nano Phase-Change Particulate Flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re=600-900 and phase change material particle concentrations ¡Ü0.25 , as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (0.5-0.7) slurry flow. By using the two newly-defined parameters, named effectiveness factor and performance index, respectively, it is found that there exists an optimal relation between the channel design parameters, particle volume fraction, Reynolds number, and the wall heat flux. The influence of the particle volume fraction, particle size, and the particle viscosity, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed.

Modelling of water removal during a paper vacuum dewatering process using a Level-Set method

Water removal in paper manufacturing is an energy-intensive process. The dewatering process generally consists of four stages of which the first three stages include mechanical water removal through gravity filtration, vacuum dewatering and wet pressing. In the fourth stage, water is removed thermally, which is the most expensive stage in terms of energy use. In order to analyse water removal during a vacuum dewatering process, a numerical model was created by using a Level-Set method. Various different 2D structures of the paper model were created in MATLAB code with randomly positioned circular fibres with identical orientation. The model considers the influence of the forming fabric which supports the paper sheet during the dewatering process, by using volume forces to represent flow resistance in the momentum equation. The models were used to estimate the dry content of the porous structure for various dwell times. The relation between dry content and dwell time was compared to laboratory data for paper sheets with basis weights of 20 and 50 g/m2 exposed to vacuum levels between 20 kPa and 60 kPa. The comparison showed reasonable results for dewatering and air flow rates. The random positioning of the fibres influences the dewatering rate slightly. In order to achieve more accurate comparisons, the random orientation of the fibres needs to be considered, as well as the deformation and displacement of the fibres during the dewatering process.

Uma nova abordagem do método Level Set baseada em conhecimento a e priori da forma

The analysis of fluid behavior in multiphase flow is very relevant to guarantee system safety. The use of equipment to describe such behavior is subjected to factors such as the high level of investments and of specialized labor. The application of image processing techniques to flow analysis can be a good alternative, however, very little research has been developed. In this subject, this study aims at developing a new approach to image segmentation based on Level Set method that connects the active contours and prior knowledge. In order to do that, a model shape of the targeted object is trained and defined through a model of point distribution and later this model is inserted as one of the extension velocity functions for the curve evolution at zero level of level set method. The proposed approach creates a framework that consists in three terms of energy and an extension velocity function λLg(θ)+vAg(θ)+muP(0)+θf. The first three terms of the equation are the same ones introduced in (LI CHENYANG XU; FOX, 2005) and the last part of the equation θf is based on the representation of object shape proposed in this work. Two method variations are used: one restricted (Restrict Level Set - RLS) and the other with no restriction (Free Level Set - FLS). The first one is used in image segmentation that contains targets with little variation in shape and pose. The second will be used to correctly identify the shape of the bubbles in the liquid gas two phase flows. The efficiency and robustness of the approach RLS and FLS are presented in the images of the liquid gas two phase flows and in the image dataset HTZ (FERRARI et al., 2009). The results confirm the good performance of the proposed algorithm (RLS and FLS) and indicate that the approach may be used as an efficient method to validate and/or calibrate the various existing equipment used as meters for two phase flow properties, as well as in other image segmentation problems.

Sloshing and slamming oscillations in collapsible channel flow

We consider laminar high-Reynolds-number flow through a finite-length planar channel, where a portion of one wall is replaced by a thin massless elastic membrane that is held under longitudinal tension T and subject to an external pressure distribution. The flow is driven by a fixed pressure drop along the full length of the channel. We investigate the global stability of two-dimensional Poiseuille flow using a method of matched local eigenfunction expansions, which is compared to direct numerical simulations. We trace the neutral stability curve of the primary oscillatory instability of the system, illustrating a transition from high-frequency ‘sloshing’ oscillations at high T to vigorous ‘slamming’ motion at low T . Small-amplitude sloshing at high T can be captured using a low-order eigenmode truncation involving four surface-based modes in the compliant segment of the channel coupled to Womersley flow in the rigid segments. At lower tensions, we show that hydrodynamic modes contribute increasingly to the global instability and we demonstrate a change in the mechanism of energy transfer from the mean flow, with viscous effects being destabilising. Simulations of finite-amplitude oscillations at low T reveal a generic slamming motion, in which the the flexible membrane is drawn close to the opposite rigid wall before rapidly recovering. A simple model is used to demonstrate how fluid inertia in the downstream rigid channel segment, coupled to membrane curvature downstream of the moving constriction, together control slamming dynamics.

Simulador de escoamento multifásico em poços de petróleo (SEMPP)

The multiphase flow occurrence in the oil and gas industry is common throughout fluid path, production, transportation and refining. The multiphase flow is defined as flow simultaneously composed of two or more phases with different properties and immiscible. An important computational tool for the design, planning and optimization production systems is multiphase flow simulation in pipelines and porous media, usually made by multiphase flow commercial simulators. The main purpose of the multiphase flow simulators is predicting pressure and temperature at any point at the production system. This work proposes the development of a multiphase flow simulator able to predict the dynamic pressure and temperature gradient in vertical, directional and horizontal wells. The prediction of pressure and temperature profiles was made by numerical integration using marching algorithm with empirical correlations and mechanistic model to predict pressure gradient. The development of this tool involved set of routines implemented through software programming Embarcadero C++ Builder® 2010 version, which allowed the creation of executable file compatible with Microsoft Windows® operating systems. The simulator validation was conduct by computational experiments and comparison the results with the PIPESIM®. In general, the developed simulator achieved excellent results compared with those obtained by PIPESIM and can be used as a tool to assist production systems development

On the Stability of Supersonic Boundary Layers with Injection

The problem of supersonic flow over a 5 degree half-angle cone with injection of gas through a porous section on the body into the boundary layer is studied experimentally. Three injected gases are used: helium, nitrogen, and RC318 (octafluorocyclobutane). Experiments are performed in a Mach 4 Ludwieg tube with nitrogen as the free stream gas. Shaping of the injector section relative to the rest of the body is found to admit a "tuned" injection rate which minimizes the strength of shock waves formed by injection. A high-speed schlieren imaging system with a framing rate of 290 kHz is used to study the instability in the region of flow downstream of injection, referred to as the injection layer. This work provides the first experimental data on the wavelength, convective speed, and frequency of the instability in such a flow. The stability characteristics of the injection layer are found to be very similar to those of a free shear layer. The findings of this work present a new paradigm for future stability analyses of supersonic flow with injection.

Tutte's 5-flow conjecture

The aim of this thesis is to show and put together the results, obtained so far, useful to tackle a conjecture of graph theory proposed in 1954 by William Thomas Tutte. The conjecture in question is Tutte's 5-flow conjecture, which states that every bridgeless graph admits a nowhere-zero 5-flow, namely a flow with non-zero integer values between -4 and 4. We will start by giving some basics on graph theory, useful for the followings, and proving some results about flows on oriented graphs and in particular about the flow polynomial. Next we will treat two cases: graphs embeddable in the plane $\mathbb{R}^2$ and graphs embeddable in the projective plane $\mathbb{P}^2$. In the first case we will see the correlation between flows and colorings and prove a theorem even stronger than Tutte's conjecture, using the 4-color theorem. In the second case we will see how in 1984 Richard Steinberg used Fleischner's Splitting Lemma to show that there can be no minimal counterexample of the conjecture in the case of graphs in the projective plane. In the fourth chapter we will look at the theorems of François Jaeger (1976) and Paul D. Seymour (1981). The former proved that every bridgeless graph admits a nowhere-zero 8-flow, the latter managed to go even further showing that every bridgeless graph admits a nowhere-zero 6-flow. In the fifth and final chapter there will be a short introduction to the Tutte polynomial and it will be shown how it is related to the flow polynomial via the Recipe Theorem. Finally we will see some applications of flows through the study of networks and their properties.

Avaliação dos fluxos das valvas aórtica e pulmonar com ecocardiografia Doppler pulsátil em cães clinicamente sadios

Avaliaram-se com Doppler pulsátil os fluxos das valvas aórtica (AO) e pulmonar (Pul) por meio de análise qualitativa (presença de regurgitações valvares e características do espectro avaliado) e quantitativa, com obtenção de parâmetros ecocardiográficos como velocidades máxima (V. Max.) e média (V. Me.), integral de velocidade (VTI), tempo de aceleração (TA) e ejeção (TE), de 30 cães considerados clinicamente sadios por meio de exames físico, laboratoriais, eletrocardiográfico, ecocardiográfico (modos uni e bidimensional), radiográfico de tórax e mensuração da pressão arterial sistêmica. Obtiveram-se os seguintes resultados para os referidos parâmetros: V.max. AO= 1,22± 19,38m/s; V. Me. AO= 0,72± 0,08m/s; VTI AO= 0,14± 0,02m; TA AO= 38,80± 11,29ms; TE AO= 197,90± 24,77ms; V. Max. Pul= 0,95± 0,18m/s; V. Me. Pul= 0,63± 0,10m/s; VTI Pul= 0,13± 0,02m; TA Pul= 70,97± 18,87ms; TE Pul= 203,70± 28,98ms. Em apenas três animais observou-se regurgitação pulmonar. Alguns parâmetros apresentaram correlação negativa com a variável freqüência cardíaca (VTI AO, TE AO, VTI Pul, TA Pul, TE Pul); outros correlação positiva com a variável peso (VTI AO, TA AO, TE AO, VTI Pul, TE Pul,) e não se observou influência da variável sexo na maioria dos parâmetros avaliados. Na comparação entre os dois fluxos, observaram-se V. Max. AO e V. Me. AO maiores que V. Max. Pul. e V. Me. Pul., respectivamente, VTI AO maior que VTI Pul, e TA AO menor que TA Pul.

Adaptations in autonomic function during exercise training in heart failure

Although neurohumoral excitation is the hallmark of heart failure (HF), the mechanisms underlying this alteration are not entirely known. Abnormalities in several systems contribute to neurohumoral excitation in HF, including arterial and cardiopulmonary baroreceptors, central and peripheral chemoreceptors, cardiac chemoreceptors, and central nervous system abnormalities. Exercise intolerance is characteristic of chronic HF, and growing evidence strongly suggests that exercise limitation in patients with chronic HF is not due to elevated filling pressures or inadequate cardiac output during exercise, but instead due to skeletal myopathy. Several lines of evidence suggest that sympathetic excitation contributes to the skeletal myopathy of HF, since sympathetic activity mediates vasoconstriction at rest and during exercise likely restrains muscle blood flow, arteriolar dilatation, and capillary recruitment, leading to underperfused areas of working muscle, and areas of muscle ischemia, release of reactive oxygen species (ROS), and inflammation. Although controversial, either unmyelinated, metabolite-sensitive afferent fibers, and/or myelinated, mechanosensitive afferent fibers in skeletal muscle underlie the exaggerated sympathetic activity in HF. Exercise training has emerged as a unique non-pharmacological strategy for the treatment of HF. Regular exercise improves functional capacity and quality of life, and perhaps prognosis in chronic HF patients. Recent studies have provided convincing evidence that these benefits in chronic HF patients are mediated by significant reduction in central sympathetic outflow as a consequence of improvement in arterial and chemoreflex controls, and correction of central nervous system abnormalities, and increase in peripheral blood flow with reduction in cytokines and increase in mass muscle.

Turbulence in Hydraulic Jumps: Experimental Measurements

A hydraulic jump is the transition from a supercritical open channel flow to a subcritical regime. It is characterised by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air-water flow measurements were performed in hydraulic jump flows for a range of inflow Froude numbers. The experiments were conducted in a large-size facility using two types of phase-detection intrusive probes: i.e., single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The present study was focused on the turbulence characteristics of hydraulic jumps with partially-developed inflow conditions. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The present results highlighted some influence of the inflow Froude number onto the air bubble entrainment process. At the largest Froude numbers, the advected air bubbles were more thoroughly dispersed vertically, and larger amount of air bubbles were detected in the turbulent shear layer. In the air-water mixing layer, the maximum void fraction and bubble count rate data showed some longitudinal decay function in the flow direction. Such trends were previously reported in the literature. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The present data suggested some longitudinal decay of the turbulence intensity. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of CHANSON (2006,2007). The results provided the integral turbulent time and length scales of the eddy structures advecting the air bubbles in the developing shear layer. The experimental data showed that the auto-correlation time scale Txx was larger than the transverse cross-correlation time scale Txz. The integral turbulence length scale Lxz was a function of the inflow conditions, of the streamwise position (x-x1)/d1 and vertical elevation y/d1. Herein the dimensionless integral turbulent length scale Lxz/d1 was closely related to the inflow depth: i.e., Lxz/d1 = 0.2 to 0.8, with Lxz increasing towards the free-surface. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. A linear relationship was found between the normalized maximum free-surface fluctuation and the inflow Froude number.

Residence time distributions of gas flowing through rotating drum bioreactors

Residence time distribution studies of gas through a rotating drum bioreactor for solid-state fermentation were performed using carbon monoxide as a tracer gas. The exit concentration as a function of time differed considerably from profiles expected for plug flow, plug flow with axial dispersion, and continuous stirred tank reactor (CSTR) models. The data were then fitted by least-squares analysis to mathematical models describing a central plug flow region surrounded by either one dead region (a three-parameter model) or two dead regions (a five-parameter model). Model parameters were the dispersion coefficient in the central plug flow region, the volumes of the dead regions, and the exchange rates between the different regions. The superficial velocity of the gas through the reactor has a large effect on parameter values. Increased superficial velocity tends to decrease dead region volumes, interregion transfer rates, and axial dispersion. The significant deviation from CSTR, plug flow, and plug flow with axial dispersion of the residence time distribution of gas within small-scale reactors can lead to underestimation of the calculation of mass and heat transfer coefficients and hence has implications for reactor design and scaleup. (C) 2001 John Wiley & Sons, Inc.

On the effect of tracer density on axial dispersion in a batch oscillatory baffled column

In this paper, we report our modelling evaluation on the effect of tracer density on axial dispersion in a batch oscillatory baffled column (OBC). Tracer solution of potassium nitrite, its specific density ranged from 1.0 to 1.5, was used in the study, and was injected to the vertical column from either the top or bottom. Local concentration profiles are measured using conductivity probes at two locations along the height of the column. Using the experimental measured concentration profiles together with both 'Tank-in-Series' and 'Plug Flow with Axial Dispersion' models, axial dispersion coefficients were determined and used to describe the effect of specific tracer density on mixing in the OBC. The results showed that the axial dispersion coefficients evaluated by the two models are very similar in both magnitudes and trends, and the range of variations in such coefficients is generally larger for the bottom injection than for the top one. Empirical correlations linking the mechanical energy for mixing, the specific density of tracer and axial dispersion coefficient were established. Using these correlations, we identified the enhancements of up to 269% on axial dispersion for various specific tracer densities. (C) 2002 Elsevier Science B.V. All rights reserved.

Climatology of the Iberia coastal low-level wind jet: weather research forecasting model high-resolution results

Coastal low-level jets (CLLJ) are a low-tropospheric wind feature driven by the pressure gradient produced by a sharp contrast between high temperatures over land and lower temperatures over the sea. This contrast between the cold ocean and the warm land in the summer is intensified by the impact of the coastal parallel winds on the ocean generating upwelling currents, sharpening the temperature gradient close to the coast and giving rise to strong baroclinic structures at the coast. During summertime, the Iberian Peninsula is often under the effect of the Azores High and of a thermal low pressure system inland, leading to a seasonal wind, in the west coast, called the Nortada (northerly wind). This study presents a regional climatology of the CLLJ off the west coast of the Iberian Peninsula, based on a 9km resolution downscaling dataset, produced using the Weather Research and Forecasting (WRF) mesoscale model, forced by 19 years of ERA-Interim reanalysis (1989-2007). The simulation results show that the jet hourly frequency of occurrence in the summer is above 30% and decreases to about 10% during spring and autumn. The monthly frequencies of occurrence can reach higher values, around 40% in summer months, and reveal large inter-annual variability in all three seasons. In the summer, at a daily base, the CLLJ is present in almost 70% of the days. The CLLJ wind direction is mostly from north-northeasterly and occurs more persistently in three areas where the interaction of the jet flow with local capes and headlands is more pronounced. The coastal jets in this area occur at heights between 300 and 400 m, and its speed has a mean around 15 m/s, reaching maximum speeds of 25 m/s.