Understanding how swirl flows can be applied to processes in the oil and gas industry and how problems might hinder them, are the focus of this thesis. Three application areas were identified: wet gas metering, liquid loading in gas wells and erosion at pipe bends due to sand transport. For all three areas, Computational Fluid Dynamics (CFD) simulations were performed. Where available, experimental data were used to validate the CFD results. As a part of this project, a new test loop was conceived for the investigation of sand erosion in pipes. The results obtained from CFD simulations of two-phase (air-water) flow through a pipe with a swirl-inducing device show that generating swirl flow leads to separation of the phases and creates distinct flow patterns within the pipe. This effect can be used in each of the three application areas of interest. For the wet gas metering application, a chart was generated, which suggests the location of maximum liquid deposition downstream of the swirling device used in the ANUMET meter. This will allow taking pressure and phase fraction measurements (from which the liquid flow rate can be determined) where they are most representative of the flow pattern assumed for the ANUMET calculation algorithms. For the liquid loading application, which was taken as an upscaling of the dimensions investigated for the wet gas metering application, the main focus was on the liquid hold-up. This parameter is defined as the ratio of the flowing area occupied by liquid to the total area. Results obtained with CFD simulations showed that as the water rate increases, the liquid hold-up increases, implying a more effective liquid removal. Thus, it was concluded that the introduction of a swirler can help unload liquid from a gas well, although no investigation was carried out on the persistance of the swirl motion downstream of the device. For the third and final application, the erosion at pipe bends due to sand transport, the main focus was to check the erosion rate on the pipe wall with and without the introduction of a swirler. The erosion rate was predicted by CFD simulations. The flow that was investigated consisted of a liquid phase with solid particles suspended in it. The CFD results showed a significant reduction in erosion rate at the pipe walls when the swirler was introduced, which could translate into an extended working life for the pipe. An extensive literature review performed on this topic, complemented by the CFD simulations, showed the need for a dedicated multiphase test loop for the investigation of sand erosion in horizontal pipes and at bends. The design of a facility of this type is included in this thesis. The results obtained with this work are very encouraging and provide a broad perspective of applications of swirl flows and CFD for the oil and gas industry.

Lean premixed combustors are highly susceptible to combustion instabilities, caused by the coupling between heat release fluctuations and combustor acoustics. In order to predict the conditions under which these instabilities occur and their limit cycle amplitudes, understanding of the amplitude dependent response of the flame to acoustic excitation is required. Extensive maps of the flame response were obtained as a function of perturbation amplitude, frequency, and flow velocity. These maps illustrated substantial nonlinearity in the perturbation velocity - heat release relationship, with complex topological dependencies that illustrate folds and kinks when plotted in frequency-amplitude-heat release space. A detailed analysis of phase locked OH PLIF images of acoustically excited swirl flames was used to identify the key controlling physical processes and qualitatively discuss their characteristics. The results illustrate that the flame response is not controlled by any single physical process but rather by several simultaneously occurring processes which are potentially competing, and whose relative significance depends upon forcing frequency, amplitude of excitation, and flame stabilization dynamics. An in-depth study on the effect of acoustic forcing on the turbulent flame properties was conducted in a turbulent Bunsen flame using PIV measurements. The results showed that the flame brush thickness and the local consumption speed were modulated in the presence of acoustic forcing. These results will not only be a useful input to help improve combustion dynamics predictions but will also help serve as validation data for models.

This work addresses the issue of transverse combustion instabilities in annular gas turbine combustor geometries. While modern low-emissions combustion strategies have made great strides in reducing the production of toxic emissions in aircraft engines and power generation gas turbines, combustion instability remains one of the foremost technical challenges in the development of next generation combustor technology. To that end, this work investigates the response of a swirling flow and swirl-stabilized flame to a transverse acoustic field is using a variety of high-speed laser techniques, especially high-speed particle image velocimetry (PIV) for detailed velocity measurements of this highly unsteady flow phenomenon. A description of the velocity-coupled transverse instability mechanism is explained with companion measurements describing each of the velocity disturbance pathways. Dependence on acoustic frequency, amplitude, and field symmetry is discussed. Significant emphasis is placed on the response of a swirling flow field to a transverse acoustic field. Details of the dynamics of the vortex breakdown bubble and the shear layers are explained using a wide variety of measurements for both non-reacting and reacting flow cases. This thesis concludes with an overview of the impact of this work and suggestions for future research in this area.

<p>This thesis work aimed at increasing the knowledge regarding phenomena occurring when gas is injected using a top-blown lance on to a bath. All results are based on physical modeling studies carried out both using low and high gas flow rates and nozzle diameters ranging from 0.8 mm to 3.0 mm. At the low gas flow rates, the penetration depth in the bath was studied. The experiments focused on studying the effect of nozzle diameter, lance height and gas flow rate on the penetration depth. It was found that the penetration depth increases with decreasing nozzle diameter, decreasing lance height and with increasing gas flow rate. The results were also compared with previous work. More specifically, it was studied how the previous published empirical relationships fitted the current experimental data. It was found that the relationships of Banks [1], Davenport [2], Chatterjee [3] and Qian [4] agreed well with the experimental data of this investigation for nozzle diameters of 2.0 mm and 3.0 mm. However, for smaller nozzle diameters there were considerable deviations. Therefore, a new correlation heuristically derived from energy conservation consideration was suggested and showed better agreement for small nozzle diameters.</p><p>The experiments carried out at higher gas flow rates focused on the study of swirl motion. The effects of nozzle diameter, lance height, gas flow rate and aspect ratio on the swirl motion were investigated. The amplitude and period of the swirl as well as the starting time and the damping time of the swirl were determined. The amplitude was found to increase with an increased nozzle diameter and gas flow rate, while the period had a constant value of about 0.5 s for all nozzle diameters, gas flow rates and lance heights. The starting time for the swirl motion was found to decrease with an increased gas flow, while the damping time was found to be independent of gas flow rate, nozzle diameter, lance height and ratio of depth to diameter.</p>

<p>Three-dimensional turbulent flow in a cylindrical hydrocyclone is considered and studied by means of computational fluid dynamics using software packages CFX and Fluent. The aim has been to identify the methods that can be used for accurate simulation of the flow in three-dimensional configurations in hydrocyclones at high swirl numbers.</p><p>As a starting point, swirling pipe flows created by tangential inlets, where detailed experimental data were available in literature, were considered. It was found that the velocity profiles for the flow with a swirl number of 2.67 could be predicted accurately using a Reynolds stress model and an accurate numerical discretization on a fine-enough mesh. At a higher swirl number, 7.84, under-prediction in the tangential velocity profiles was observed; however the prediction of the axial velocity profiles was satisfactory. The validated methods were then used to simulate the flow in a cylindrical hydrocyclone at a swirl number as large as 21. The calculated tangential velocity profiles were compared against experimental data measured with a pitometer. Acceptable agreements were recorded except near the geometric axis of the cyclone. Due to the lack of the aircore in the numerical model, disagreements near the axis of the cyclone could be expected to some extent.</p><p>Numerical experiments performed in the present work indicated that the RNG k-ε model is not likely to be capable to predict highly swirling flows accurately and a Reynolds stress model is required. For three-dimensional models, where the computing capacity and the available memory set strong restrictions on the computational mesh, optimizing the maximum mesh resolution available play an important role on the accuracy and stability of the solution procedure. The most stable results in the present study were found using the Reynolds stress model proposed by Launder et al. on an as regular and structured mesh as possible using a higher order discretization scheme in Fluent. Therefore, the meshing capabilities of the pre-processor, the available turbulence models and the accuracy of the numerical methods must be considered in parallel. Acceptable results were also generated using the Baseline Reynolds stress model implemented in CFX, however, only with a transient procedure which was likely to be more time-consuming.</p><p>Present simulations present a complex flow structure in the cylindrical cyclone with a double axial flow reversal. The effect of such a flow pattern on the fractionation of the fibres with small differences in density needs to be investigated in future studies.</p>

A study was carried out into the computational fluid dynamic simulation of bio-oil combustion. Measurements were taken in an empirical burner to obtain information regarding the flow behaviour. A surrogate fuel was developed to mimic the unique chemical and physical properties of bio-oil combustion. The resulting computational model of the burner domain and surrogate fuel was compared with empirical data. The bio-oil model displayed a good agreement with the data in terms of the combustion behaviour, but was limited by the uncertain flow solution associated with the burner used.

A study was carried out into the computational fluid dynamic simulation of bio-oil combustion. Measurements were taken in an empirical burner to obtain information regarding the flow behaviour. A surrogate fuel was developed to mimic the unique chemical and physical properties of bio-oil combustion. The resulting computational model of the burner domain and surrogate fuel was compared with empirical data. The bio-oil model displayed a good agreement with the data in terms of the combustion behaviour, but was limited by the uncertain flow solution associated with the burner used.

The aim of this thesis is to increase engine efficiency by studying the flow structures created in an engine cylinder during the intake phase and the effect of the subsequent compression. The invention of the combustion engine has enabled three centuries of economic growth fueled by energy stored as hydrocarbons. However, during the latter part of the twentieth century negative consequences on health and environment of the combustion engine were observed. In order to reduce emissions without increasing fuel consumption, improved knowledge of all physical processes occurring in the engine are necessary. The aim of this thesis is to increase the understanding of the flow prior to combustion, which can lead to reduced engine emissions and fuel consumption. Intake flow structures are studied using large eddy simulations and experiments on a steady swirl test rig. Flow acceleration was observed to reduce the swirl coefficient, and higher swirl coefficient was found during valve closing as compared to during valve opening. This implies that the rotation is stronger during the later part of the intake then what has been previously assumed. In addition, the computations show that the volume above the valves has a profound effect on the swirl created during the intake. To take this into account a novel way of calculating the swirl number was suggested. This approach gives a lower swirl number as compared to the commonly used Thien methodology. The effects of compression are studied using simulations of predefined flow structures undergoing compression. The peak turbulence levels were found to be increasing with tumble number and decreasing with swirl. It was noted that compression increased the turbulent fluctuations in the cylinder axis leading to anisotropic turbulence and that a small tilt angle was observed to have a significant effect on swirl homogeneity at top dead center.  In this thesis, a new methodology was proposed and validated for calculation of in-cylinder turbulence for a flat piston. The results of the thesis enhance the understanding of the dynamic effects encountered during intake as well recognizing that a small tumble component has a strong effect on the flow structures prior to combustion. These results can be used to improve the simplified computational methods used to optimize the engine. / <p>QC 20150420</p>

Short annular diffusers were essential components for turbomachines that have been used to expand the air entering the compressor, as interstage ducts between gas generators and power turbines, and on the exhaust gases exiting the turbine. The industrial community was interested and invested in improving diffuser design that was challenging owing to the unfavourable fluid flow effects. Efficient design of fluid flow devices was possible through the complementary use of experimental testing and computational fluid dynamics (CFD). A numerical shape optimization study was undertaken to determine preferential annular diffuser configurations. Experimental data were compared against CFD that simulated the steady-state Reynolds-averaged Navier-Stokes equations with two-equation turbulence models. This investigation reached equivalent conclusions with respect to the influences associated with diffuser geometry and swirl. Vorticity effects caused by square tabs, that were not as well understood, were investigated. The tabs were effective in reducing the central toroidal recirculation zone created by a swirling flow, but at a static pressure penalty for the area ratio, AR<2.73, diffusers tested. Results identified several shortcomings in the CFD that typically over-estimated pressure recovery and outlet velocity uniformity; however, properly qualitatively predicted wall pressure distributions and outlet velocity profiles. The use of CFD on modest grids, with preference given to the realizable k-epsilon turbulence model, for annular diffusers that have length to inlet height ratio of 12 and at least AR=2.73 with up to 20-degrees inlet swirl was encouraged as a design tool. / Thesis (Ph.D, Mechanical and Materials Engineering) -- Queen's University, 2014-05-29 09:03:16.591

Des crises vibratoires ont été constatées dans plusieurs centrales thermiques d’EDF opérant avec du fioul lourd, certaines ayant entraîné l’arrêt du foyer. Ce travail traite des instabilités de combustion pouvant se déclencher dans ce type de système où le combustible liquide est injecté avec de la vapeur d’eau et où l’écoulement d’air est mis en rotation. Ces phénomènes vibratoires résultent d’un couplage résonant entre la dynamique de la combustion et l’acoustique du foyer. La réponse acoustique des flammes diphasiques non-prémélangées swirlées reste largement méconnue et est difficilement analysable sur le foyer réel. L’objectif de ce travail est donc d’étudier la stabilité des chaudières EDF à partir de l’analyse de la réponse d’une flamme diphasique non-prémélangée swirlée issue d’un injecteur générique et soumise à des perturbations de la vitesse acoustique. Cette réponse est déterminée sur un dispositif original (DIFAV) équipé d’un swirler et d’un injecteur bi-fluides fonctionnant à la vapeur d’eau et au dodécane. Ce système est constitué des principaux éléments des brûleurs utilisés sur les centrales thermiques EDF à une échelle 1/7000. Le dispositif est conçu pour facilement modifier la géométrie de la tête d’injection, les conditions d’injection de combustible et de vapeur et ainsi contrôler le spray généré. Des visualisations à la sortie d’une buse d’injection montrent l’influence de la topologie de l’écoulement diphasique dans l’injecteur sur la taille des gouttes mesurées dans le spray. Des mesures de taille et de vitesse des gouttes lorsque le rapport des débits de vapeur et de combustible (GLR) est modifié sont réalisées. Ces données comparées à des modèles ont permis d’estimer l’évolution de la taille des gouttes générées par l’injecteur qui équipe les centrales thermiques lorsque le GLR varie. Une analyse modale du foyer DIFAV et d’un modèle simplifié de la chaudière réelle est ensuite menée. Les fréquences propres et les taux d’amortissement du foyer DIFAV sont déterminés expérimentalement en soumettant le système à une modulation acoustique externe. Un modèle acoustique simplifié composé de trois cavités couplées représentatif du brûleur DIFAV est également développé. Des simulations acoustiques réalisées avec COMSOL Multiphysics sur une coupe transverse d’une chaudière générique représentative de la chaudière industrielle permettent d’identifier trois modes à basses fréquences établis entre les plenums et la chambre de combustion qui sont susceptibles d’être instables. La sensibilité de ces modes à la géométrie du foyer et aux conditions limites est étudiée. La réponse de la flamme générique lorsqu’elle est soumise à des modulations acoustiques de l’écoulement d’air en amont du brûleur est ensuite mesurée sur le banc DIFAV pour différents niveaux d’excitation et deux topologies de flamme lorsque les conditions d’injection sont modifiées. Les mécanismes qui pilotent l’évolution du gain de l’une des fonctions de transfert généralisées (FDF) de la flamme sont étudiés à l’aide de visualisations en moyenne de phase de l’écoulement et de mesures des vitesses axiale et azimutale de l’écoulement d’air au cours d’un cycle de modulation. Une forte sensibilité de la phase de la FDF à l’amplitude des perturbations acoustiques est observée. Un adimensionnement par le nombre de Strouhal basé sur la vitesse débitante et la longueur efficace de la flamme est proposé pour transposer ces FDFs sur le brûleur réel. Une analyse de stabilité du foyer DIFAV est réalisée en intégrant les FDF au modèle acoustique afin de déterminer les cycles limites des oscillations lorsque la longueur de la chambre de combustion varie. Ces calculs sont comparés aux fréquences des instabilités auto-entretenues mesurées aux cycles limites dans le foyer DIFAV. [...] / Vibratory crises have been observed in EDF thermal power plants operating with heavy fuel oil. Such instabilities may lead to shutdown and damage the boiler. This work deals with combustion instabilities that can take place in boilers equipped with steam-assisted atomizers and where the airflow is swirled. These vibratory phenomena result from a resonant coupling between the combustion dynamics and the boiler acoustics. Analyses of combustion dynamics of non-premixed swirling spray flames remain rare and are difficult to realize on the real system. The objective of this work is to analyze the stability of EDF boilers using the response of generic non-premixed swirling spray flames submitted to acoustic velocity disturbances. This response is determined on an original device (DIFAV) equipped with a swirling vane and a twin-fluid atomizer operated with steam and dodecane. This burner is equipped with the main elements of those used in the thermal power plant, but has a reduced scale of 1/7000. The influence of the injector geometry and of the operating conditions on the spray generated by the injector can be studied. Spray visualizations at the outlet of the injector reveal the relationship between the topology of the two-phase flow in the injector and the measured droplet size. Measurements of the droplet diameter and velocity as a function of the gas-to-liquid ratio (GLR) have been performed at the outlet of the injector. These data have been compared to models and were used to estimate the evolution of the droplets diameter as a function of the GLR generated by the industrial injector. A modal analysis of the DIFAV combustor is then carried out and a simplified acoustic model made of three coupled cavities is developed. The natural frequencies and damping rates of the DIFAV combustor are determined experimentally when it is submitted to acoustic modulation. Acoustic simulations are performed with COMSOL Multiphysics on a simplified geometrical model of the industrial boiler. Three low frequency modes established between the plenums and the combustion chamber have been identified and may be unstable. Their sensitivity to modifications of the boiler geometry and boundary conditions are studied. Flame responses subjected to acoustic modulations of the airflow rate are then measured on the DIFAV combustor for several amplitudes and two flames topologies obtained at globally lean condition. Phase-conditioned flame visualizations and measurements of swirl number fluctuations during an acoustic forcing cycle are conducted to explain the mechanisms that control the evolution of gain of the Flame Describing Function (FDF). A high sensitivity of the phase of the FDF to the amplitude of the acoustic disturbance is observed. The Strouhal number based on the airflow velocity and the effective length of the flame is used to transpose these FDF on the industrial burner. FDF are integrated in the acoustic model of the DIFAV setup to carry out a stability analysis and predict the limit cycle oscillations as a function of the combustion chamber length. These calculations are compared to frequencies of self-sustained instability measured at the limit cycles in the DIFAV combustor. A reasonable agreement is obtained showing the validity of the stability analysis for the non-premixed two-phase flames investigated based on the knowledge of their FDF. Finally, a stability analysis of the EDF boiler is conducted with the COMSOL Multiphysics model by including the acoustic flame response of the industrial burner in the simulation. This FDF is deducted from the dimensionless FDF measured on the generic burner. The Rayleigh criterion is used to analyze the stability of the combustor as a function of the flame length for different boundary conditions. Indications are given to improve the stability of the EDF boiler.