The research in this thesis focuses on the investigation of tidal turbines using a small scale horizontal axis tidal turbine and a 2D hydrofoil testing rig, combining experiments with simulations to provide comprehensive results and to better understand some of the variables that affect their performance. The experimental campaigns were carried out at the University of Victoria fluids research lab and the Sustainable Systems Design Lab (SSDL). The experimental testing rigs were re-designed by the author and are now fully automated, including a friendly graphical user interface for easy implementation. Particle image velocimetry (PIV) technique was used as the quantitative flow visualization method to obtain the time-averaged flow fields. This thesis presents three investigations. The first study aims to quantify the impacts of channel blockage, free surface effects and foundations on hydrokinetic turbine performance, using porous discs and an axial flow rotor. The results were used to cross-validate computational fluid dynamics (CFD) simulations. It was found that as wall blockage increases, thrust and power are incremented with and without the inclusion of free surface deformation. Discrepancies between simulations and experimental results on free surface effects compared to a slip wall were obtained and hence further research is recommended and the author gives some advice on how to proceed in this investigation. The second study determines the performance of four hydrofoil candidates over a range of low Reynolds number (Re), delivering useful information that can be applied to low Re energy conversion systems and, specifically in this case, to improve the performance of the small scale tidal turbine at the SSDL lab. The study combines the 2D hydrofoil test rig along with PIV measurements in order to experimentally obtain lift and drag coefficients. The experiments were carried out in the recirculating flume tank over the range of low Re expected for the small scale rotor rig, in order to provide more accurate results to improve rotor blade design. In addition, numerical simulations using XFOIL, a viscid-inviscid coupled method, were introduced to the study. These results were analysed against experiments to find the most suitable parameters for reliable performance prediction. The final results suggested that adding a numerical trip at a certain chordwise distance produced more reliable results. Finally, an experimental study on turbine rotor performance and tip vortex behavior was performed using again the rotor rig and PIV. Blade design and rotor performance were assessed, showing good agreement with Blade Element Momentum (BEM) simulations, particularly at predicting the tip speed ratio corresponding to the maximum power coefficient point. Regarding the wake structure, tip vortex locations (shed from the blade tips) were captured using PIV in the near wake region, showing evidence of wake expansion. The velocity and vorticity fields are also provided to contribute to the development and validation of CFD and potential flow codes. / Graduate / 0548 / 0547 / 0538 / iafranch@uvic.ca

Le présent manuscrit de thèse présente l'étude de différents phénomènes dans un convertisseur d?acier, grâce au développement de deux modèles distincts. Le premier modèle décrit la cavité produite à la surface libre du bain de métal par l'impact du jet d'oxygène supersonique. Il est basé sur le découpage du domaine de calcul en deux régions. Les effets de compressibilité du gaz sont pris en compte uniquement dans la région du jet où la vitesse est élevée, alors que partout ailleurs, le gaz est considéré comme incompressible. La méthode Volume Of Fluid (VOF) est utilisée pour suivre la déformation de la surface libre du bain. Les résultats de simulations sont présentés pour des systèmes bi- et triphasés et comparés à des données expérimentales obtenues dans diverses maquettes froides. L'influence sur la taille et la forme de la cavité de différents paramètres (parmi lesquels les conditions aux limites en sortie de la lance d'injection, le schéma d'advection de la méthode VOF et le modèle de turbulence) est étudiée. Le modèle est ensuite utilisé pour simuler l'interaction entre un jet supersonique d'oxygène et la surface libre d'un bain d'acier liquide dans un convertisseur de taille pilote. Le second modèle se focalise sur l'écoulement du gaz, le transfert de chaleur et la réaction de postcombustion dans la phase gazeuse au-dessus du bain de métal. Il utilise l'algorithme Simple Chemical Reaction Scheme pour décrire le transport des espèces chimiques, et prend en compte l'absorption d'oxygène dans le bain et les transferts thermiques par rayonnement. Les prédictions numériques sont en assez bon accord avec les mesures recueillies dans une expérience de laboratoire et dans un four à l'échelle pilote / The present thesis treats different phenomena taking place in a steelmaking converter through the development of two separate models. The first model describes the cavity produced at the free surface of the metal bath by the high speed impinging oxygen jet. It is based on a zonal approach, where gas compressibility effects are taken into account only in the high velocity jet region while elsewhere the gas is treated as incompressible. The Volume Of Fluid (VOF) method is employed to follow the deformation of the bath free surface. Calculations are presented for two- and three-phase systems and compared against experimental data obtained in various cold model experiments. The influence on the size and shape of the cavity of various parameters (including the jet inlet boundary conditions, the VOF advection scheme and the turbulence modelling) is studied. Next, the model is used to simulate the interaction of a supersonic oxygen jet with the surface of a liquid steel bath in a pilot-scale converter. The second model concentrates on fluid flow, heat transfer and the post-combustion reaction in the gas phase above the metal bath. It uses the Simple Chemical Reaction Scheme approach to describe the transport of the chemical species and takes into account the consumption of oxygen by the bath and thermal radiative transfer. The numerical predictions are in reasonable agreement with measurements collected in a laboratory experiment and in a pilot-scale furnace

Ce travail est consacré au développement de schémas numériques pour approcher les solutions de modèles d'écoulement type “depth averaged”. Dans un premier temps, nous détaillons la construction d'approches Volumes Finis pour le système Shallow Water avec termes sources sur maillages non structurés. En se basant sur une reformulation appropriée des équations, nous mettons en place un schéma équilibré et préservant la positivité de la hauteur d'eau, et suggérons des extensions MUSCL adaptées. La méthode est capable de gérer des topographies irrégulières et exhibe de fortes propriétés de stabilité. L'inclusion des termes de friction fait l'objet d'une analyse poussée, aboutissant à l'établissement d'une propriété type “Asymptotic Preserving” à travers l'amélioration d'un autre récent schéma Volumes Finis. La seconde composante de cette étude concerne les méthodes Elements Finis type Galerkin discontinu. Certaines des idées avancées dans le contexte Volumes Finis sont employées pour aborder le système Shallow Water surmaillages triangulaires. Des résultats numériques sont exposés et la méthode se révèle bien adaptée à la description d'une large variété d'écoulements. Partant de ces observations nous proposons finalement d'exploiter ces caractéristiques pour étendre l'approche à une nouvelle famille d'équations type Green-Nadghi. Des validations numériques sont également proposées pour valider le modèle numérique. / This work is devoted to the development of numerical schemes to approximatesolutions of depth averaged flow models. We first detail the construction of Finite Volume approaches for the Shallow Water system with source terms on unstructured meshes. Based on a suitable reformulation of the equations, we implement a well-balanced and positive preserving approach, and suggest adapted MUSCL extensions. The method is shown to handle irregular topography variations and demonstrates strong stabilities properties. The inclusion of friction terms is subject to a thorough analysis, leading to the establishment of some Asymptotic Preserving property through the enhancement of another recent Finite Volume scheme.The second aspect of this study concerns discontinuous Galerkin Finite Elementmethods. Some of the ideas advanced in the Finite Volume context areemployed to broach the Shallow Water system on triangular meshes. Numericalresults are exposed and the method turns out to be well suited to describe a large variety of flows. On these observations we finally propose to exploit its features to extend the approach to a new family of Green-Nadghi equations. Numerical experiments are also proposed to validate this numerical model.

La rencontre entre une couche limite laminaire se développant sur une plaque plane sous une surface libre et un obstacle rectangulaire émergent mène à l’apparition d’une zone de recirculation complexe appelée tourbillon en fer à cheval à l’amont de l’obstacle. Cette zone de recirculation est composée d’un certain nombre de vortex s’enroulant autour de l’obstacle et pouvant présenter des dynamiques variées. Le but de ce travail est de réaliser une étude expérimentale paramétrique du tourbillon en fer à cheval laminaire et faiblement turbulent. Les mesures sont réalisées par PIV et par trajectographies et des méthodes de détection basés sur la décomposition modale (POD) et la détection des points critiques sont utilisées afin de faciliter l’étude de la dynamique tourbillonnaire. Ce manuscrit se propose de : (i) Déterminer l’évolution des caractéristiques du tourbillon en fer à cheval en fonction des paramètres adimensionnels de l’écoulement. (ii) Mettre en place une typologie des différentes dynamiques exhibées. (iii) Déterminer les mécanismes à l’origine de l’apparition de ces différentes dynamiques à l’aide d’un modèle conceptuel. (iv) Et finalement, étudier l’influence de la longueur de l’obstacle sur le tourbillon en fer à cheval. / A laminar boundary layer developing under a free-surface and interacting with an emerging obstacle lead to the creation of a complex recirculation zone called horseshoe vortex upstream from the obstacle. This recirculation zone is composed of a given number of vortices wrapping around the obstacle with various kind of dynamics. The aim of this work is to perform an experimental parametric study of the laminar, and slightly turbulent horseshoe vortex. Measurements are obtained by PIV and trajectographies and detection method based on modal decomposition (POD) and critical points are used to ease the study of the vortex dynamics. This manuscript intends to : (i) Explain the horseshoe vortex main characteristics evolution with the non-dimensional parameters of the flow. (ii) Define a typology of the different observed vortex dynamics. (iii) Identify the physical mechanisms at the origin of those dynamics through the use of a conceptual model. (iv) And finally, show the influence of the obstacle length on the horseshoe vortex.

Este trabalho apresenta uma técnica de diferenças finitas para resolver a equação constitutiva Fluido de Segunda Ordem para escoamentos tridimensionais com superfície livre. As equações governantes são resolvidas pelo método de diferenças finitas em uma malha deslocada 3D. A superfície livre é modelada por células marcadoras (Marker-and-Cell) e as condições de contorno a superfície livre são empregadas. O método numérico apresentado neste trabalho foi validado pela comparação entre as soluções numéricas obtidas para o escoamento em um tubo com a solução analítica correspondente para Fluidos de Segunda Ordem. Ao fazer refinamento de malha, a convergência do método numérico foi verificada. Resultados numéricos da simulação do problema do inchamento do extrudado para números de Deborah De \'< OU =\' 0:3 são apresentados / This work presents a finite difference method to simulate three-dimensional viscoelastic flow with free surfaces governed by the constitutive equation Second Order Fluid. The governing equations are solved by the finite difference method in a three-dimensional shifted mesh. The free surface of fluid is modeled by the Marker-and-Cell method which allows for the visualization and the location of the free surface of fluid. The full free surface stress conditions are employed. The numerical method developed in this work is validated by comparing the numerical and analytic solutions for the steady state flow of a Second Order Fluid in a pipe. By using mesh refinement convergence results are given. Numerical results of the simulation of the transient extrudate swell of a Second Order Fluid of the Deborah number De \'< OR =\' 0:3 are presented

O método numérico para fluidos incompressíveis desenvolvido no presente estudo é o Moving Particle Semi-Implicit Method (MPS) que enxerga o domínio discretizado em partículas, é baseado em representação lagrangeana e não tem a necessidade de utilização de malhas. O método MPS tem como equações governantes uma forma particular da equação de Navier-Stokes e a equação da continuidade para fluidos incompressíveis e não viscosos. Os métodos de simulação de fluidos mais comumente utilizados são baseados em representação euleriana e utilizam malhas para descrever a geometria do domínio a ser simulado. Devido a essas diferenças, uma das grandes virtudes do método de partículas é a facilidade de investigação de fenômenos altamente não-lineares como o de superfície livre com quebra de ondas, de líquidos no interior de uma embarcação em movimento, de ondas batendo na parte externa do casco de um navio, etc. Em artigos já publicados, resultados de experimentos físicos mostraram boa aderência aos resultados numéricos de simulações realizadas com o método MPS. No presente trabalho, resultados das forças de excitação das simulações com ondas regulares foram comparados com os resultados do programa Wave Analysis MIT (WAMIT) que é um programa consagrado no meio científico. Houve uma boa concordância de resultados entre os dois programas. A otimização do cálculo de vizinhança forneceu uma grande economia de tempo computacional. A maior contribuição deste estudo foi a otimização da função que resolve o sistema linear implementando no código desenvolvido um código paralelizado de uso público existente chamado Portable, Extensible Toolkit for Scientific Computation (PETSc) que proporcionou um bom ganho de desempenho. / A numerical method called Moving Particle Semi-implicit (MPS) method was developed in this study to analyze incompressible fluids. It is a particle method using a lagrangean representation without any grid. The governing equations are the Navier-Stokes equation and continuity equation for incompressible and non-viscous flow. Most of the computational fluid dynamics (CFD) methods are based on eulerian representation and use grids to describe the geometry of the simulated domain. These differences make the MPS method easier to analyze highly nonlinear phenomena as free surface with wave breaking, sloshing, slamming, etc. In previously published articles, results of physical experiments had shown good agreement with the numerical results obtained with MPS method. In the present work, results of exciting forces were compared with the results obtained with a validated program called Wave Analysis MIT (WAMIT). It had a good agreement of results between these two programs. The optimization of the neighborhood calculation function got a good economy of computational time. The greatest contribution of this study was the optimization of the linear system solver. It was made implementing in the developed code a parallelized public code called Portable, Extensible Toolkit for Scientific Computation (PETSc) that provided a good performance profit.

No contexto do método MAC e baseado em esquemas de diferenças finitas, este trabalho apresenta três estudos: i) uma análise de estabilidade, ii) o desenvolvimento de técnicas implícitas e, iii) a construção de métodos de projeção para escoamentos com superfície livre. Na análise de estabilidade, o principal resultado mostra que o método de Crank-Nicolson torna-se condicionalmente estável quando aplicado para uma malha deslocada com a discretiza ção explícita das condições de contorno do tipo Dirichlet. Entretanto, o mesmo método com condições de contorno implícitas é incondicionalmente estável. Para obter métodos mais estáveis, formulações implícitas são desenvolvidas para a equação da pressão na superfície livre, derivada da condição de tensão normal. Esta estratégia resulta no acoplamento dos campos de velocidade e pressão, o que exige a introdução de novos métodos de projeção. Os métodos de projeção assim desenvolvidos resultam em novas metodologias para escoamentos com superfície livre que são apropriados para o tratamento de problemas com baixo número de Reynolds. Além disso, mostra-se que os métodos propostos podem ser aplicados para fluidos viscoelásticos. Novas estratégias são derivadas para obter métodos de projeção de segunda ordem de precisão para escoamentos com superfícies livres. Além dos resultados teóricos sobre a estabilidade de esquemas numéricos, técnicas implícitas e métodos de projeção, testes computacionais são realizados e comparados para consolidação da teoria apresentada. Os resultados numéricos são obtidos no sistema FREEFLOW. A eficiência e robustez das técnicas desenvolvidas neste trabalho são demonstradas na solução de problemas tridimensionais complexos com superfície livre e baixo número de Reynolds, incluindo os problemas do jato oscilante e do inchamento do extrudado / In the context of the MAC method and based on finite difference schemes, this work presents three studies: i) a stability analysis, ii) the development of implicit techniques, and iii) the construction of projection methods for free surface flows. In the stability analysis, the main result shows a precise stability restriction on the Crank-Nicolson method when one uses a staggered grid with Dirichlet explicit boundary conditions. However, the same method with implicit boundary conditions becomes unconditionally stable. In order to obtain more stable methods, implicit formulations are applied for the pressure equation at the free surface, which is derived from the normal stress condition. This approach results in a coupling of the velocity and pressure fields; hence new projection methods for free surface flows need to be developed. The developed projection methods result in new methodologies for low Reynolds number free surface flows. It is also shown that the proposed methods can be applied for viscoelastic fluids. New strategies are derived for obtaining second-order accurate projection methods for free surface flows. In addition to the theoretical results on the stability of numerical schemes, implicit techniques and projection methods, computational tests are carried out and the results compared to consolidate the theory. The numerical results are obtained by the FREEFLOW system. The eficiency and robustness of the techniques in this work are demonstrated by solving complex tridimensional problems involving free surface and low Reynolds numbers, including the jet buckling and the extrudate swell problems

The marine current turbine (MCT) is an exciting proposition for the extraction of renewable tidal and marine current power. However, the numerical prediction of the performance of the MCT is difficult due to its complex geometry, the surrounding turbulent flow and the free surface. The main purpose of this research is to develop a computational tool for the simulation of a MCT in turbulent flow and in this thesis, the author has modified a 3D Large Eddy Simulation (LES) numerical code to simulate a three blade MCT under a variety of operating conditions based on the Immersed Boundary Method (IBM) and the Conservative Level Set Method (CLS). The interaction between the solid structure and surrounding fluid is modelled by the immersed boundary method, which the author modified to handle the complex geometrical conditions. The conservative free surface (CLS) scheme was implemented in the original Cgles code to capture the free surface effect. A series of simulations of turbulent flow in an open channel with different slope conditions were conducted using the modified free surface code. Supercritical flow with Froude number up to 1.94 was simulated and a decrease of the integral constant in the law of the wall has been noticed which matches well with the experimental data. Further simulations of the marine current turbine in turbulent flow have been carried out for different operating conditions and good match with experimental data was observed for all flow conditions. The effect of waves on the performance of the turbine was also investigated and it has been noticed that this existence will increase the power performance of the turbine due to the increase of free stream velocity.

"A Java pre- and post-processing graphical user-oriented interface has been developed by the authors to aid a mill engineer with little or no finite element experience throughout the analysis process of the finishing rolling stands. A case study is presented that uses the commercial finite element code ABAQUS/Explicit to predict roundness and tolerance customer requirements. Other parameters that are determined include spread, crosssectional area, percentage reduction in area, incremental plastic strain, total plastic strain and roll force. All parameters are compared to theoretical models and some are compared to full-scale mill testing."

O presente trabalho tem por objetivo desenvolver uma formulação numérica baseada em Análise Isogeométrica para o estudo de problemas de interação fluido-estrutura (IFE) em aplicações envolvendo corpos rígidos submersos, onde escoamentos incompressíveis de fluidos Newtonianos com superfície livre são considerados. Propõe-se o emprego da Análise Isogeométrica por permitir a unificação entre os procedimentos de pré-processamento e análise, melhorando assim as condições de continuidade das funções de base empregadas tanto na discretização espacial do problema como na aproximação das variáveis do sistema de equações. O sistema de equações fundamentais do escoamento é formado pelas equações de Navier-Stokes e pela equação da conservação de massa, descrita segundo a hipótese de pseudo-compressibilidade, em uma formulação cinemática ALE (Arbitrary Lagrangean- Eulerian). A consideração da superfície livre no escoamento se dá tratando o fluido como um meio bifásico, através do método Level Set. O corpo rígido apresenta não linearidade na rotação e restrições representadas por vínculos elásticos e amortecedores viscosos, sendo a equação de equilíbrio dinâmico resolvida através do método de Newmark. O esquema de acoplamento sólido-fluido adotado é o particionado convencional, que impõe condições de compatibilidade cinemáticas e de equilíbrio sobre a interface sólido-fluido, analisando ambos os meios de maneira sequencial. A discretização das equações governantes é realizada através do esquema explícito de dois passos de Taylor-Galerkin, aplicado no contexto da Análise Isogeométrica. Por fim, são analisados alguns problemas da Dinâmica de Fluidos Computacional, de onde se concluiu que os resultados obtidos são bastante consistentes com os fenômenos envolvidos, com as ferramentas exclusivas da Análise Isogeométrica, como o refinamento k, melhorando a convergência dos resultados. Para escoamentos bifásicos, verificou-se que o método Level Set obteve resultados bastante promissores apresentando, entretanto, uma dissipação numérica excessiva. Propõe-se, para estudos futuros, a elaboração de esquemas numéricos que conservem melhor o volume da fase líquida do escoamento. / The present work aims to development of a numerical formulation based on Isogeometric Analysis for the study of Fluid-Structure Interaction problems in applications involving rigid bodies submerged, considering incompressible Newtonian flows with free surface. The use of the Isogeometric Analysis allows unification between the preprocessing and analysis steps, improving then the continuity of the base functions employed, both in the spatial discretization and approximation of the variables in the system of equations. The fundamental flow equations are formed by the Navier-Stokes and the mass conservation, described by de pseudo-compressibility hypothesis, in an ALE (Arbitrary Lagrangean-Eulerian) kinematic formulation. The free surface consideration of the flow is handled treating the fluid like a two- phase medium, using the Level Set method. The rigid body considers nonlinearity in rotation, and restrictions represented by elastic springs and viscous dampers, with the dynamic equilibrium equation being resolved using the Newmark’s method. The solid-fluid coupling scheme is the conventional partitioned, which imposes kinematics and equilibrium compatibility conditions on the solid-fluid interface, analyzing both mediums in a sequential manner. The governing equations are discretized using the explicit two step Taylor-Galerkin method, applied in an Isogeometric Analisys context. Finally, some Computational Fluid Dinamics problems are analysed, from which it was concluded that the results obtained are quite consistent with phenomena involved, with the unique tools of Isogeometric Analysis, such as k-refinement, improving the convergence of the results. For biphasic flows, it was verified that the Level Set method obtained very promising results, presenting, however, an excessive numerical dissipation. For future studies, it is proposed the elaboration of numerical schemes that better preserve the volume of the liquid phase of the flow.