This thesis evaluates two variants of Volumeof fluid method, namely Piecewise Linear Interface Construction (PLIC) and Multidimensional Tangent of Hyperbola for Interface Capturing (MTHINC) used in two-fluid simulations for interface tracking.Simulations of a single rising bubble in different conditions are performed aimed to assess the accuracy of the methods, and to check their convergence with grid size. The results show that both methods demonstrate convergence with mesh refinement, while PLIC captures the transient part of rising process of bubble more accurately than MTHINC in 2D case. / Denna avhandling utvärderar två varianter av metoden Volume of fluid för flerfas strömningar, nämligen Piecewise Linear Interface Construction (PLIC) och Multidimensional Tangent of Hyperbola for Interface Capturing (MTHINC) som används i tvåvätskesimuleringar för gränssnittsspårning.Simuleringar av en enda stigande bubbla under olika förhållanden utförs för att bedöma metodernas noggrannhet och för att kontrollera deras konvergens med gridsstorlek. Resultaten visar att båda metoderna visar konvergens med grid förfining, medan PLIC fångar den övergående delen av stigande bubbelprocess mer exakt än MTHINC i 2D-fall.

External loop airlift reactors are widely used for biochemical applications such as syngas fermentation and wastewater treatment. To further understand the inherent gas-liquid flow physics within the reactors, computational modeling and simulations of hydrodynamics for air-water external loop airlift reactors were investigated. The gas-liquid flow dynamics in a bubble column were simulated using a FORTRAN code developed by Los Alamos National Laboratory, CFDLib, which employs an Eulerian-Eulerian ensemble averaged method. A two-dimensional Cartesian coordinate system was used to conduct an extensive grid resolution study; it was found that grid cells smaller than the bubble diameter produced unstable solutions. Next, closure models for drag force and turbulent viscosity were investigated for a simple bubble column geometry. The effects of using a bubble pressure model and two drag coefficient models, the White model and the Schiller-Naumann model, were investigated. The bubble pressure model performed best for homogeneous (low velocity) flows and the Schiller-Naumann model was best for all flow regimes. Based on the studies for bubble column flows, an external loop airlift reactor was simulated using both two- and three-dimensional coordinates and results for gas holdup and riser velocity agreed better with experimental data for the 3D simulations. It was concluded that when performing 2D and 3D simulations, care must be taken when specifying the effective bubble diameter size, especially at high flow rates. Population balance models (PBM) for bubble break-up and coalescence were implemented into CFDLib, validated with experiments, and simulated for the external loop airlift reactor at high inlet superficial gas velocities. The PBM predictions for multiple bubble sizes were comparable with the single bubble size simulations; however, the PBM simulations better predicted the formation of the gas bubble in the downcomer. The 3D PBM simulations also gave better predictions for the average bubble diameter size in the riser. It was concluded that a two-dimensional domain is adequate for gas-liquid flow simulations of a simple bubble column geometry, whereas three-dimensional simulations are required for the complex airlift reactor geometry. To conclude, a two-fluid Eulerian-Eulerian model coupled with a PBM is needed for quantitative as well as physical predictions of gas-liquid external loop airlift reactor flows at high inlet superficial gas velocities. / Ph. D.

Biomass is considered a biorenewable alternative energy resource that can potentially reduce the use of natural gas and provide low cost power production or process heating needs. Biomass hydrodynamics in a fluidized bed are extremely important to industries that are using biomass material in gasfication processes to yield high quality producer gas. However, biomass particles are typically difficult to fluidize due to their peculiar shape and a second inert material, such as sand, is typically added to the bed. The large differences in size and density between the biomass and inert particles lead to nonuniform distribution of the biomass within the fluidized bed, and particle interactions and mixing become major issues. The main goal of this research was to use CFD as a tool for modeling and analyzing the hydrodynamic behavior of biomassas a single material or as part of a mixture in a fluidized bed. The first part of this research focused on the characterization of biomass particles in a fluidized bed and validation of a numerical model with experimental results obtained from pressure measurements and CT and X-ray radiograph images. For a 2D fluidized bed of glass beads, the pressure drop, void fraction and mean bed height expansion were in quantitative agreement between the experiments and simulations using Syamlal-O'Brien and Gidaspow drag models. It was encouraging that the Gidaspow model predictions were in close agreement because the model does not require knowing the minimum fluidization as an input. Ground walnut shells were used to represent biomass because the material fluidizes uniformly and is classified as a Geldart type B particle. Two-dimensional simulations of ground walnut shells were analyzed to determine parameters that cannot easily be measured experimentally. The parametric study for ground walnut shell indicated that the material can be characterized with a medium sphericity (~0.6) and a relatively large coefficient of restitution (~0.85). In the second part of this work numerical simulations of a ground walnut shell fluidizing bed with side air injection were compared to CT data for the gas-solid distribution to demonstrate the quantitative agreement for bed fluidization. The findings showed that 2D simulations overpredicted the fluidized bed expansion and the results did not demonstrate a uniformly fluidizing bed. The 3D simulations compared well for all cases. This study demonstrates the importance of using a 3D model for a truly 3D flow in order to capture the hydrodynamics of the fluidized bed for a complicated flow and geometry. Finally, CFD modeling of pressure fluctuations was performed on sand and cotton-sand fluidized beds operating at inlet velocities ranging from 1.0-9.0Umf with the objective of predicting characteristic features of bubbling, slugging, and turbulent fluidization regimes. It was determined that the fluidized bed can be modeled using MUSCL discretization and the Ahmadi turbulence model. Three-dimensional sand fluidized beds were simulated for different fluidization regimes. Fluidized beds for all the regimes behaved as second-order dynamic systems. Bubbling fluidized beds showed one broad peak with a maximum at 2.6 Hz while slugging and turbulent showed two distinct peaks. It was observed that the peak at low frequency increased in magnitude as the flow transitioned from a slugging to a turbulent fluidization regime. CFD simulations of fluidized beds with the purpose of studying pressure fluctuations have demonstrated to be a useful tool to obtain hydrodynamic information that will help determine the fluidization regime. Prediction of slugging and turbulent fluidization regimes using CFD have not been reported to date. The work presented here is the first of its kind and can be an important advantage when designing a reactor and evaluating different operation conditions without the need to test them in a pilot plant or a prototype. / Ph. D.

High-speed craft are in high demand in the maritime industry, for example, in maintenance operations for offshore structures, for search and rescue, for patrolling operations, or as leisure craft to deliver speed and excitement. Design and operation of high-speed craft are often governed by the hydrodynamic phenomena of slamming, which occur when the craft impact the wave surface. Slamming loads affect the high-speed craft system; the crew, the structure and various sub-systems and limit the operation. To meet the ever-increasing demands on safety, economy and reduced environmental impact, there is a need to develop more efficient high-speed craft. This progression is however limited by the prevailing semi-empirical design methods for high-speed planing craft structures. These methods provide only a basic description of the involved physics, and their validity has been questioned. This thesis contributes to improving the conditions for designing efficient highspeed craft by focusing on two key topics: evaluation and development of the prevailing design methods for high-speed craft structures, and development towards structural design based on first principles modeling of the slamming process. In particular a methodological framework that enables detailed studies of the slamming phenomena using numerical simulations and experimental measurements is synthesized and evaluated. The methodological framework involves modeling of the wave environment, the craft hydromechanics and structural mechanics, and statistical characterization of the response processes. The framework forms the foundation for an extensive evaluation and development of the prevailing semi-empirical design methods for high-speed planing craft. Through the work presented in this thesis the framework is also shown to be a viable approach in the introduction of simulation-based design methods based on first principles modeling of the involved physics. Summarizing, the presented methods and results provide important steppingstones towards designing more efficient high-speed planing craft. / <p>QC 20160907</p>

La dynamique instable des jets tournants est étudiée, en tenant compte des variations de masse volumique au sein de l'écoulement. Un code de simulation numérique directe permettant de résoudre les équations de Navier-Stokes à masse volumique variable a été développé, en utilisant une méthode originale et efficace pour résoudre le champs de pression. Analytiquement, deux modes instables bidimensionnels ont été mis en évidence, et sont identifiés comme des modes de Couette-Taylor et de Rayleigh-Taylor, ainsi qu'un troisième mode tridimensionnel, du à un couplage de vitesse. La dynamique instable de cet écoulement résulte d'une compétition entre ces trois modes, et les simulations numériques montrent que ces modes perdurent non linéairement. Ensuite, le comportement spatio-temporel de cette instabilité est étudiée par simulation numérique directe, et il a été montré qu'il existe une transition vers des modes absolument instables, sous l'effet du rapport de densité s ainsi que du taux de rotation q. Cette dynamique est également étudiée expérimentalement au travers de plusieurs méthodes de mesures, et la présence de mode globaux auto-entretenus est mise en évidence qui sont en bon accord avec les résultats numériques. Finalement, le phénomène de l'éclatement tourbillonnaire est étudié, et montre le rôle prépondérant de la viscosité réelle. En effet, l'éclatement tourbillonnaire est un mécanisme permettant de soulager le système de l'intensification de la vorticité, au travers de la viscosité, alors qu'il n'apparaît pas en traitant les équations d'Euler tronquées. / The unstable dynamics of a swirling jet flow is studied, including density variations within the flow. A direct numerical simulation method was developed to solve variable density Navier-Stokes equations, using an accurate and efficient pressure solver. Analitically, two unstable bi-dimensionnal modes are highlighted, and are identified as Couette-Taylor and Rayleigh-Taylor modes. A three-dimensionnal mode is also highlighted, wich is created by the shear. Numerical simulations show that those modes are nonlinearly persistant. Then, the spatio-temporal instability behaviour is studied numerically, and show that the instability undergoes to a convective/absolute transition with density ratio s and rotation rate q. This dynamic is also studied experiementally through different methods, and Global selfsustained modes are highlighted wich are in ggod agreement with numerical results. Finally, the vortex breakdown phenomenon is studied, and show the crucial role of real viscosity. Indeed, the vorticity intensification is relaxed through the viscosity effect, while it is not treating the truncated Euler Equations.

La contribution des processus de fusion de galaxies lors de l'assemblage cosmologique de masses est encore mal comprise. Dans ce contexte, l'étude du support dynamique des galaxies est un moyen permettant de contraindre les différents scénarios d'évolution.Dans une première partie, je présente ma contribution à l'analyse cinématique des galaxies de l'échantillon MASSIV. Je définis une technique d'analyse qui vise à quantifier la fraction d'objets susceptibles d'être le produit d'une fusion.J'introduis dans une seconde partie une méthode de construction de conditions initiales de galaxies idéalisées.Dans une troisième partie, je présente la définition de l'échantillon de simulations MIRAGE, reproduisant les propriétés physiques des galaxies jeunes. Les simulations montrent en particulier l'absence de flambées de formation stellaire durant les fusions ce qui suggère un mécanisme de saturation.Une analyse des propriétés des complexes granulaires issus d'instabilités gravitationnelles est réalisée dans la quatrième partie, soulignant la compatibilité du scénario granules âgés avec les observations en terme de taux d'éjection de gaz et d'âge de populations stellaires. Je mets en évidence que le mécanisme de fusion de granules permet de contrôler la croissance de la masse du bulbe, d'aplatir le profil central de densité du halo de matière sombre ainsi que d'éjecter de grandes quantités de gaz.Dans une cinquième partie, je complète ce travail en réalisant une étude comparative des données cinématiques MASSIV à plus de 4000 pseudo-observations construites à partir des simulations de l'échantillon MIRAGE afin de déterminer la capacité à détecter des signatures de fusion. / The contribution of the fusion processes to the cosmological mass assembly is still poorly understood. In this context, the study of the dynamical support of galaxies is a way to constrain different evolution scenarios.In the first part, I introduce my contribution to the kinematical analysis of galaxies MASSIV. I define an analysis that aims at quantifying the fraction of objects that are likely to be the product of a recent merger.I introduce in the second part a method to build initial conditions of idealized galaxies.In the third part, I present the definition of the simulation sample MIRAGE, designed to reproduce the physical properties of young galaxies. The simulations show in particular the absence of a star formation bursts in mergers of fragmented and turbulent disks, suggesting a saturation mechanism.An analysis of the properties of clumps that are generated by gravitational instabilities is performed in the fourth part of this thesis, highlighting the compatibility of long-lived clumps scenarios with observations in terms of outflowing material and stellar population ages. I also highlight that the gas rich clump fusion mechanism is able to control the bulge mass growth, to erode the central profile of the dark matter halo and to drive massive gas outflows in the disk plane.Finally, in the fifth part, I complete this work by carrying out a comparative study of the MASSIV kinematical data to a set of more than 4000 pseudo-observations at z=1.7 built from simulations of the MIRAGE sample to determine the ability to detect galaxy merger signatures.

Lorsqu'un phénomène de détonation survient dans des milieux confinés comme un bâtiment présentant un système de galeries, l'onde de choc générée se propage et peut engendrer des dégâts matériels et humains considérables. En effet, l’onde de choc en se propageant dans ces systèmes de canaux confinés va interagir avec les obstacles qu’elle rencontre (humains et matériels), et va provoquer de très fortes élévations de la pression localement. Ce phénomène peut être très destructeur pour les structures et létale pour les humains. L'objectif de cette étude, expérimentale et numérique, est donc d'étudier la propagation d'une onde de choc dans différents systèmes de canaux afin de comprendre son comportement et de proposer des solutions pour en atténuer les dégâts. L’étude s’articule autour de deux grands axes. Le premier étudie de manière très académique la propagation d’une onde de choc au travers de cinq configurations différentes ; deux sans variation de section (un coude à 45$^{\circ}$ et un coude à 90$^{\circ}$), et trois avec variation de section (un divergent, un élargissement brusque et une bifurcation en « Y »). L’influence de la singularité sur la pression réfléchie en fond de configuration est étudiée ainsi que les mécanismes physiques complexes qui se produisent derrière l’onde de choc le long de la singularité. Le deuxième axe étudie l’atténuation de la pression réfléchie par l’ajout de piège (cavités) le long du système étudié. Une étude paramétrique montre que la position ainsi que la taille et la forme du piège ont un impact sur l’atténuation de l’onde de choc. / In the search for protection from explosions phenomena, a variety of underground shelters were studied to minimize the risks related to the propagation of shock waves in closed areas. Indeed, the blast effect could be really destructive for equipment and humans encountering the shock wave propagation, because of the high elevation of the local pressure it generates. Thus, the propagation of a shock wave through various canal systems is both experimentally and numerically studied. This study is based on two topics. The first part is focused on the study of the propagation of a shock wave through five various configurations ; two whithout aera change (a 45$^{\circ}$ bend duct and a 90$^{\circ}$ bend), and three configurations whith a sectional enlargement (a divergent, an abrupt area change and a « Y » bifurcation). The impact of the geometry on the end wall reflected pressure is studied as well as some complex physical mechanisms that occur behind the incident shock wave. The second part explores the mitigation of the pressure level in the device with the addition of traps (cavities) along the studied configuration. A parametric study based on the shape and size of the cavities, shows that these parameters have a considerable impact on the pressure level in the duct system.

Thermodynamic analysis of solid oxide cells Petr Vágner The thesis deals with continuum thermodynamic modeling and analysis of phe- nomena in solid oxide electrochemical cells. A general description of the evo- lution of charged mixtures using partial mass densities, momentum density, entropy density, electric induction, magnetic field, polarization, and magnetiza- tion based on the GENERIC framework is formulated. The formulation is used to recover the Landau-Lifshitz magnetization relaxation model, the Single Re- laxation Time model for dielectrics, and the generalized Poisson-Nernst-Planck model. The latter model is consequently linked to the second part, where a novel double layer model of an yttria-stabilized zirconia interface is formulated within non-equilibrium thermodynamics. The model is solved for numerically in the time domain, and cyclic voltammetry of the system is analyzed. The last part of the thesis demonstrates the limits of Exergy Analysis on a simple solid oxide hydrogen fuel cell model with non-isothermal boundary. It is demon- strated that the minimization of entropy production does not necessarily lead to the maximization of the electric power for certain optimization scenarios. The thesis consists of a compilation of published and unpublished results of the author.

Cette étude s’intéresse à la modélisation du comportement thermoviscoplastique ainsi qu’à la rupture d’un acier de construction « S355NL » et d’un acier à blindage naval. Une campagne expérimentale a permis d’observer la réponse mécanique de ces matériaux sur une large gamme de condition de sollicitation, notamment leurs sensibilités à la vitesse de déformation ( de 〖10〗^(-3) à ~〖10〗^4 s^(-1)) et à la température (de -100 à 200 °C). Ces résultats ont été exploités lors de l’identification de lois de comportement permettant de représenter précisément le comportement mécanique de ces matériaux. Afin d’améliorer la description de la sensibilité à la vitesse de déformation, un nouvelle approche a également été développée. Une fois ces lois déterminées, elles ont été utilisées dans des simulations numériques modélisant l’impact d’un projectile sur ces aciers. Les résultats obtenus ont été comparés avec les essais expérimentaux correspondant afin d’estimer la fiabilité des données et des hypothèses utilisées lors des simulations. / The present study deals with the modelling of the thermo-viscoplastic behaviour and the failure of a structural steel “S355NL” and a naval armour steel. An experimental campaign have been performed to observe the mechanical response of these material over a wide range of conditions, especially their sensitivities to the strain rate (from 〖10〗^(-3) to ~〖10〗^4 s^(-1)) and to the temperature (from -100 to 200 °C). The obtained results have been used to identify the parameters of several constitutive relations. A new approach have been developed to improve the description of the strain rate sensitivity. These relations have been implemented in numerical simulation to model the impact of a projectile on these steels. The obtained results have been compared with corresponding experimental tests in order to assess the reliability of the data and the hypothesis used in the simulations.

The local universe is dominated by quiescent galaxies with little or no ongoing star formation. Once star formation has been shut down in a possible quasar phase, energy feedback in the form of outflows from a supermassive black hole (SMBH) accreting at a low rate is one of the leading candidates for heating up or removing gas accreted afterwards. In this work we performed hydrodynamic simulations of radiatively inefficient accretion flows around supermassive black holes and investigate the generation of outflows by the accretion flow. We found that hot accretion flows can produce powerful subrelativistic winds that carry considerable amounts of energy away and they can provide feedback inside the host galaxy. / O universo local é dominado por galáxias quiescentes com pouca, ou nenhuma, taxa atual de formação estelar. Uma vez que a formação estelar é suprimida numa possível fase de quasar do núcleo ativo da galáxia, o feedback de energia na forma de outflows do buraco negro supermassivo acretando a pequenas taxas é um dos principais candidatos a aquecer ou remover o gás do sistema. Nesse trabalho, executamos simulações hidrodinâmicas de escoamentos acretivos radiativamente ineficientes ao redor de buracos negros supermassivos e investigamos a geração de outflows pelo escoamento. Encontramos que escoamentos acretivos quentes podem produzir ventos subrelativísticos poderosos que carregam para fora consideráveis quantidades de energia e que podem providenciar feedback dentro da galáxia hospedeira.