La modélisation des phénomènes physiologiques induits par un écoulement, tels que l'écoulement sanguin au travers d'une sténose ou l'écoulement d'air lors de la production de parole, repose souvent sur des théories quasi-unidimensionnelles ou bi-dimensionnelles. Cependant, il est établi que le développement des couches limites dépend de la section transversale. Le but de cette thèse est de modéliser, simuler et caractériser l'importance potentielle de la section transversale sur les écoulements laminaires, contrôlés en pression, en l'absence ou en présence d'une constriction. Des coordonnées de translation sont utilisées pour obtenir des solutions pour des écoulement visqueux au travers d'une section de forme arbitraire. Cette paramétrisation est appliquée à la résolution des équations physiques pour des formes à deux et à trois dimensions. Un modèle d'écoulement simplifié quasi-tridimensionnel, qui prend en compte les pertes dissipatives par convection, la viscosité et la forme de la section est présenté et appliqué à la description de l'écoulement le long d'une sténose. Des données expérimentales et issues de simulations numériques sont collectées afin de caractériser l'influence de la forme de la section transversale dans le cas d'une constriction. simulation numérique sont comparées. / Physical models of physiological flow-induced phenomena, such as blood flow through a stenosis or air flow during human speech production, often rely on a quasi-one-dimensional or two-dimensional flow model, so that details of the cross section shape are neglected. Nevertheless, boundary layer development is known to depend on the cross section shape. The aim of this thesis is to model, simulate and characterize the potential impact of the cross section shape for pressure-driven laminar channel flow without and with constriction. Stretched coordinates are introduced to obtain viscous flow solutions for channels with an arbitrary cross section. The proposed cross section shape parametrization is applied to solve physical equations for two-dimensional and three-dimensional shapes. A simplified quasi-three-dimensional flow model, which accounts for kinetic losses, viscosity and the cross section shape, is presented and applied to describe the flow through a stenosis. Finally, flow data are gathered experimentally and numerically in order to characterize the influence of the cross section shape in the case of a constricted channel. Modeled, experimental and numerical data are compared.

Made available in DSpace on 2014-06-11T19:23:38Z (GMT). No. of bitstreams: 0 Previous issue date: 2010-08-30Bitstream added on 2014-06-13T20:30:31Z : No. of bitstreams: 1 rodrigues_tt_me_ilha.pdf: 2332044 bytes, checksum: 6e277868f01ea21a46dfbf5c827109a7 (MD5) / Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) / Fundação de Ensino Pesquisa e Extensão de Ilha Solteira (FEPISA) / A compreensão do escoamento em válvulas de compressores herméticos alternativos é de fundamental importância para introduzir modificações no projeto delas de maneira a aumentar a performance dos compressores, e por fim, dos ciclos de refrigeração. A válvula do compressor é um dispositivo ímpar, umas vez que seu funcionamento se dá pela ação da pressão exercida pelo escoamento, caracterizando um problema de forte interação fluido- estrutura. O uso da modelagem numérica através das ferramentas da mecânica dos fluidos computacional (CFD) tem se destacado como a alternativa mais dinâmica para o estudo do fenômeno. O trabalho desenvolvido foi voltado para o estudo numérico do escoamento através do difusor radial, o qual é um modelo simplificado da válvula, com o emprego do Método da Fronteira Imersa com Modelo Físico Virtual para a modelagem do disco superior do difusor (palheta). O ponto forte desta metodologia é que a representação de regiões sólidas é feita pelo cálculo de um campo de força, o qual é introduzido nas equações das células na vizinhança do sólido. Este procedimento dispensa o uso de malhas que se adaptam ao corpo, possibilitando o uso de malhas cartesianas convencionas para modelar geometrias complexas e móveis. A metodologia foi acoplada com a solução das equações governantes do escoamento em coordenadas cilíndricas através do Método dos Volumes Finitos. Inicialmente, a metodologia foi validada, utilizando como dados de referência resultados provenientes de estudos numéricos e experimentais, e foi avaliada a influência dos parâmetros do procedimento na qualidade final dos resultados. Na segunda etapa foram desenvolvidos estudos preliminares referentes ao movimento do disco superior, com a imposição artificial dos processos de abertura e fechamento da válvula. Os resultados obtidos mostraram que a metodologia adotada... / The fully understanding of the flow through automatic valves of alternative hermetic compressors is essentiall to introduce modifications in its project aiming the improvement of the compressor performance and, also, the refrigeration cycle. The compressor valve is a singular device, once its operation is ruled by the flow pressure, characterizing a case with a strong fluid-structure interaction. The using of numerical tools trough the methods of computational fluid dynamics (CFD) has gained especial attention due to its flexibility to study the phenomenon. The present work was developed to study numerically the flow through the radial diffuser, which is a simplified model of the valve, with the employment of the Immersed Boundary Method with Virtual Physical Model to modeling the superior disk (valve reed). The main advantage of this methodology is that the modeling of solid boundaries is performed with the calculus of a force field, which is introduced in the cells equations nearby the solid. This procedure dispenses the using of body-fitted meshes, enabling the adoption of conventional Cartesian meshes to model complex and moving geometries. The methodology was coupled with the solution of the governing equations in cylindrical coordinates though the Finite Volume Method. Firstly, the methodology was validated, confronting the results obtained with data from numerical and experimental studies, where the influence of the main parameters in the quality of the final results was evaluated. In the second step were developed preliminary studies concerning the movement of the superior disk, whose opening and closing movements were artificially imposed. The results obtained showed that the adopted methodology is quit promising and flexible, and can be employed in more refined studies to the comprehension of the flow through the valve regarding the fluid-structure interaction that rules the problem

Made available in DSpace on 2014-06-11T19:23:39Z (GMT). No. of bitstreams: 0 Previous issue date: 2009-09-08Bitstream added on 2014-06-13T19:09:39Z : No. of bitstreams: 1 lacerda_jf_me_ilha.pdf: 4631781 bytes, checksum: cad33ad7c5cac9de6cec08fd94b419e7 (MD5) / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Nesse trabalho resolve-se numericamente o escoamento em difusores radiais, os quais têm sido utilizados como modelo para o estudo do escoamento em válvulas automáticas de compressores de refrigeração. As equações governantes do escoamento, escritas no sistema de coordenadas cilíndricas, são resolvidas utilizando um código numérico baseado no Método dos Volumes Finitos. O Método da Fronteira Imersa, em conjunto com o Modelo Físico Virtual, foi implementado no código numérico e utilizado para representar a região sólida imersa no escoamento. Inicialmente, o código numérico foi utilizado para resolver o problema do escoamento em torno de um cilindro de base quadrada, como parte do processo de validação do código. O confronto dos resultados numéricos com dados da literatura indicou a validação parcial do código. Posteriormente, realizou-se um estudo preliminar do comportamento da solução do escoamento no difusor radial em relação a diversos parâmetros geométricos e de simulação numérica, com o objetivo de identificar a configuração numérica capaz de fornecer, simultaneamente, resultados satisfatórios com o menor custo computacional. Usando esta configuração, o código numérico foi validado através da comparação dos resultados da distribuição de pressão sobre o disco frontal (palheta) do difusor com dados experimentais da literatura, para duas distâncias entre disco frontal e disco inferior (assento), s=0,02 e 0,025 cm, e números de Reynolds variando entre 1500 e 2500. As comparações entre esses resultados mostraram que a metodologia utilizada é adequada para estudar o problema. Finalmente, a geometria do disco inferior foi modificada através da inserção de um chanfro com três ângulos de inclinação (30, 45 e 60º), com o objetivo de avaliar sua influência sobre o comportamento da distribuição de pressão sobre o disco frontal... / In this work the flow in radial diffusers, which have been used as a model to study the flow in refrigeration compressors valves, is solved numerically. The governing equations, written in cylindrical coordinates, are solved using a numerical code based on the Finite Volume Method. The Immersed Boundary Method, with the Virtual Physical Method, was implemented in the numerical code and used to represent the solid region immersed in the flow. Initially, the numerical code was used to solve the flow around a square base cylinder, as part of the code validation. The comparison of the numerical results with literature data indicated the partial validation of the code. Afterwards it was performed a preliminary study of the behavior of the flow solution in the radial diffuser relating to several geometrical and numerical parameters, with the objective of identifying a configuration capable of providing, simultaneously, satisfactory results with the smaller computational cost. Using this configuration, the numerical code was validated through the comparison with experimental pressure distribution on the frontal disk (reed) for two gaps between the frontal disk and inferior disk (seat), s=0.020 and 0.025 cm, and Reynolds numbers varying between 1500 and 2500. These comparisons have shown that the implemented methodology is suitable to study this problem. Finally, the inferior disk geometry was modified by inserting a chamfer with three inclination angles (30, 45 e 60º), with the objective of evaluating its influence on the pressure distribution on the frontal disk. The decreasing of the total pressure gradient through the flow for increasing inclination angles indicates less amount of energy to drive the flow. This is an important result that can be used to design refrigeration compressors with lower compression power

Orientador: José Luiz Gasche / Banca: Sérgio Said Mansur / Banca: Aristeu da Silveira Neto / Resumo: Nesse trabalho resolve-se numericamente o escoamento em difusores radiais, os quais têm sido utilizados como modelo para o estudo do escoamento em válvulas automáticas de compressores de refrigeração. As equações governantes do escoamento, escritas no sistema de coordenadas cilíndricas, são resolvidas utilizando um código numérico baseado no Método dos Volumes Finitos. O Método da Fronteira Imersa, em conjunto com o Modelo Físico Virtual, foi implementado no código numérico e utilizado para representar a região sólida imersa no escoamento. Inicialmente, o código numérico foi utilizado para resolver o problema do escoamento em torno de um cilindro de base quadrada, como parte do processo de validação do código. O confronto dos resultados numéricos com dados da literatura indicou a validação parcial do código. Posteriormente, realizou-se um estudo preliminar do comportamento da solução do escoamento no difusor radial em relação a diversos parâmetros geométricos e de simulação numérica, com o objetivo de identificar a configuração numérica capaz de fornecer, simultaneamente, resultados satisfatórios com o menor custo computacional. Usando esta configuração, o código numérico foi validado através da comparação dos resultados da distribuição de pressão sobre o disco frontal (palheta) do difusor com dados experimentais da literatura, para duas distâncias entre disco frontal e disco inferior (assento), s=0,02 e 0,025 cm, e números de Reynolds variando entre 1500 e 2500. As comparações entre esses resultados mostraram que a metodologia utilizada é adequada para estudar o problema. Finalmente, a geometria do disco inferior foi modificada através da inserção de um chanfro com três ângulos de inclinação (30, 45 e 60º), com o objetivo de avaliar sua influência sobre o comportamento da distribuição de pressão sobre o disco frontal... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: In this work the flow in radial diffusers, which have been used as a model to study the flow in refrigeration compressors valves, is solved numerically. The governing equations, written in cylindrical coordinates, are solved using a numerical code based on the Finite Volume Method. The Immersed Boundary Method, with the Virtual Physical Method, was implemented in the numerical code and used to represent the solid region immersed in the flow. Initially, the numerical code was used to solve the flow around a square base cylinder, as part of the code validation. The comparison of the numerical results with literature data indicated the partial validation of the code. Afterwards it was performed a preliminary study of the behavior of the flow solution in the radial diffuser relating to several geometrical and numerical parameters, with the objective of identifying a configuration capable of providing, simultaneously, satisfactory results with the smaller computational cost. Using this configuration, the numerical code was validated through the comparison with experimental pressure distribution on the frontal disk (reed) for two gaps between the frontal disk and inferior disk (seat), s=0.020 and 0.025 cm, and Reynolds numbers varying between 1500 and 2500. These comparisons have shown that the implemented methodology is suitable to study this problem. Finally, the inferior disk geometry was modified by inserting a chamfer with three inclination angles (30, 45 e 60º), with the objective of evaluating its influence on the pressure distribution on the frontal disk. The decreasing of the total pressure gradient through the flow for increasing inclination angles indicates less amount of energy to drive the flow. This is an important result that can be used to design refrigeration compressors with lower compression power / Mestre

O escoamento em torno de corpos tem sido objeto de estudo de muitos pesquisadores e é muito explorado experimental e computacionalmente, devido a sua grande aplicabilidade na engenharia. No entanto, simular computacionalmente este tipo de escoamento requer uma atenção especial ao escolher o tipo malha a ser utilizado. Em muitos casos faz-se necessário o uso de uma malha que se adapte ao contorno do obstáculo, o que pode ocasionar um aumento no esforço computacional. Um maneira de contornar este problema é a utilização do Método das Fronteiras Imersas, que possibilita o uso de malha cartesiana na simulação computacional do escoamento em torno de obstáculos. Isso é possível através da adição de um termo forçante nas equações que modelam o escoamento, e assim as forças que agem sobre o contorno do corpo são transferidas diretamente para a malha. O objetivo deste trabalho de mestrado foi implementar o método das Fronteiras Imersas e simular o escoamento em torno de um cilindro circular em repouso, movimentando-se na mesma direção do escoamento, na direção perpendicular ao escoamento, ou rotacionando em torno do próprio eixo. As simulações computacionais possibilitaram a captura do fenômeno de Atrelagem Síncrona, caracterizado pela sincronia entre a frequência de desprendimento natural de vórtices e a frequência de oscilação do mesmo. O Método das Fronteiras Imersas mostrou um ótimo desempenho quando comparado a resultados experimentais e numéricos encontrados na literatura / The flow around bodies have been studied by many researchers. Both experimental and computational approaches have been extensively explored in researches on flow around bodies and have been applied in many engeneering problems. However, to choose an appropriate type of mesh to perform computational simulations of this type of problem requires special attention. In many cases, it is necessary to use a mesh that is able to conform to the boundary if a given obstacle. The need to perform this adaptation may increase the computational effort. The Immersed Boundary Method enables the use of cartesian meshes to perform computational simulations of flows around obstacles. The idea of this method is to add a forcing term in the equations that model the flow. Thus, the forces applied on the body boundaries are directly transfered to the mesh. The aim of this work was to perform a computational implementation of the Immersed Boundary Method to simulate the flow over a oscilating circular cylinder. This oscilation may be inline with the flow, cross-flow, or rotating. The computational simulations enabled the capture of the lock-in phenomena, which consists of the syncronization between the vortex shedding frequency and the cylinder oscilation frequency. The results obtained from the computational simulations using the Immersed Boundary Method were in good agreement with the numerical and experimental results found in the literature

Dans cette thèse, une nouvelle méthode de frontières immergées a été développée pour la simulation d'interaction fluide-structure, appelée la méthode de frontières immergées mobiles (en langage anglo-saxon: MIBM). L'objectif principal de cette nouvelle méthode est de déplacer arbitrairement les solides à géométrie complexe dans un fluide visqueux incompressible, sans remailler le domaine fluide. Cette nouvelle méthode a l'avantage d'imposer la condition de non-glissement à l'interface d'une manière exacte via une force sans introduire des constantes artificielles modélisant la structure rigide. Cet avantage conduit également à la satisfaction de la condition CFL avec un pas de temps plus grand. Pour un calcul précis de la force induite par les frontières mobiles, un système linéaire a été introduit et résolu par la méthode de gradient conjugué. La méthode proposée peut être intégrée facilement dans des solveurs résolvant les équations de Navier-Stokes. Dans ce travail la MIBM a été mise en œuvre en couplage avec un solveur fluide utilisant une méthode de projection adaptée pour obtenir des solutions d'ordre deux en temps et en espace. Le champ de pression a été obtenu par l'équation de Poisson qui a été résolue à l'aide de la méthode du gradient conjugué préconditionné par la méthode multi-grille. La combinaison de ces deux méthodes a permis un gain de temps considérable par rapport aux méthodes classiques de la résolution des systèmes linéaires. De plus le code de calcul développé a été parallélisé sur l'unité graphique GPU équipée de la bibliothèque CUDA pour aboutir à des hautes performances de calcul. Enfin, comme application de nos travaux sur la MIBM, nous avons étudié le couplage "fort" d'interaction fluide-structure (IFS). Pour ce type de couplage, un schéma implicite partitionné a été adopté dans lequel les conditions à l'interface sont satisfaites via un schéma de type "point fixe". Pour réduire le temps de calcul inhérent à cette application, un nouveau schéma de couplage a été proposé pour éviter la résolution de l'équation de Poisson durant les itérations du "point fixe". Cette nouvelle façon de résoudre les problèmes IFS a montré des performances prometteuses pour des systèmes en IFS complexe. / In this thesis a novel non-body conforming mesh formulation is developed, called the moving immersed boundary method (MIBM), for the numerical simulation of fluid-structure interaction (FSI). The primary goal is to enable solids of complex shape to move arbitrarily in an incompressible viscous fluid, without fitting the solid boundary motion with dynamic meshes. This novel method enforces the no-slip boundary condition exactly at the fluid-solid interface with a boundary force, without introducing any artificial constants to the rigid body formulation. As a result, large time step can be used in current method. To determine the boundary force more efficiently in case of moving boundaries, an additional moving force equation is derived and the resulting system is solved by the conjugate gradient method. The proposed method is highly portable and can be integrated into any fluid solver as a plug-in. In the present thesis, the MIBM is implemented in the fluid solver based on the projection method. In order to obtain results of high accuracy, the rotational incremental pressure correction projection method is adopted, which is free of numerical boundary layer and is second order accurate. To accelerate the calculation of the pressure Poisson equation, the multi-grid method is employed as a preconditioner together with the conjugate gradient method as a solver. The code is further parallelized on the graphics processing unit (GPU) with the CUDA library to enjoy high performance computing. At last, the proposed MIBM is applied to the study of two-way FSI problem. For stability and modularity reasons, a partitioned implicit scheme is selected for this strongly coupled problem. The interface matching of fluid and solid variables is realized through a fixed point iteration. To reduce the computational cost, a novel efficient coupling scheme is proposed by removing the time-consuming pressure Poisson equation from this fixed point interaction. The proposed method has shown a promising performance in modeling complex FSI system.

[pt] O escoamento de cápsulas suspensas em uma fase líquida através de canais e capilares micrométricos representa um problema complexo que ocorre em diferentes aplicações, de glóbulos vermelhos em hemodinâmica até escoamento em meios porosos. Em aplicações de meios porosos, a compreensão da dinâmica na microescala é fundamental para avaliar o comportamento macroscópico do escoamento. Canais e capilares com constrição podem ser usados para modelar uma garganta conectando dois poros adjacentes. O escoamento de uma cápsula suspensa através de tais modelos foi analisado para avaliar as características do escoamento considerando os efeitos inerciais (isto é, número de Reynolds finito), incluindo a máxima diferença de pressão necessária para empurrar uma cápsula através da constrição em função do raio da cápsula, a tensão inicial e o material da membrana, geometria do canal e do capilar, assim como as condições de escoamento. De fato, neste estudo, a resposta da pressão é fundamental para avaliar o efeito de bloqueio da cápsula. As fases líquidas internas e externas foram descritas pelas equações de Navier-Stokes, enquanto que a dinâmica da membrana da cápsula foi modelada por uma estrutura flexível 1-D tipo mola. O problema de interação fluido-estrutura foi resolvido usando o método de elementos finitos acoplado ao método de fronteira imersa. Os resultados mostraram a redução da mobilidade da fase contínua devido à presença da cápsula através da constrição. Tais resultados podem ser usados para projetar microcápsulas para bloquear caminhos preferenciais de fluxo da água no processo de deslocamento de óleo em meios porosos. / [en] The flow of capsules suspended in a liquid phase through small channels and capillaries poses a complex problem presented in different applications, from red blood cells on hemodynamics to flow in porous media. In applications of porous media, the understanding of microscale dynamics is fundamental to assess the macroscopic flow behavior. Constricted channels and capillaries can be used to model a pore throat connecting two adjacent pore bodies. The flow of a suspended capsule through such models was analyzed to evaluate the flow characteristics considering inertial effects (i.e. finite Reynolds numbers), including the maximum pressure difference required to push a capsule through the constriction as a function of capsule radius, initial membrane tension, membrane material, channel and capillary geometries, as well as flow conditions. In fact, in this study, the pressure response is fundamental in order to assess the capsule blocking mechanism. Inner and outer liquid phases were described by the Navier-Stokes equations and capsule membrane dynamics was modeled by a 1-D spring-like flexible structure. The fluid-structure interaction problem was solved using the finite element method coupled with the immersed boundary method. Results showed the mobility reduction of the continuous phase due to the presence of a capsule as it flows through the constriction. Such results can be used to design microcapsules to block preferential water flow paths in oil displacement process in porous media.

Em diversos escoamentos sobre superfícies há a presença de protuberâncias, como por exemplo rebites, parafusos e juntas. Estas protuberâncias podem influenciar a camada limite, acelerando a transição do escoamento do estado laminar para o estado turbulento. Em alguns casos isto pode ser indesejável, já que o escoamento turbulento implica necessariamente em uma força de atrito maior do que aquela referente ao escoamento laminar. Existem alguns aspectos neste tipo de escoamento que ainda não estão bem compreendidos. O objetivo deste trabalho é estudar a influência de uma rugosidade isolada no escoamento sobre uma superfície. Este estudo contribui para se entender o que ocorre em casos de maior complexidade. O estudo é de natureza computacional, em que se utiliza simulação numérica direta das equações de Navier-Stokes. A técnica de fronteiras imersas é utilizada para representar a rugosidade no escoamento sobre a superfície. O código numérico é verificado por meio do método de soluções manufaturadas. Comparações entre resultados experimentais, da teoria de estabilidade linear e numéricos também são utilizados para a validação do código. Resultados obtidos com diferentes alturas de rugosidade e variações no gradiente de pressão permitiram analisar a influência de elemento rugoso tridimensional em escoamentos de camada limite. / The presence of protuberances on surfaces, for example, rivets, screws and gaskets, can influence the boundary layer by accelerating the transition from laminar flow to turbulent flow. In some cases this may be undesirable, since the turbulent flow involves frictional forces greater than the ones at the laminar regime. There are some aspects of the flow in the boundary layer perturbed by a single roughness element that are not well understood. The aim of this work is to study the influence of an isolated roughness on the boundary layer. This study is a step towards to the understanding of what can happen in more complex cases. The nature of this study is computational, therefore a Direct Numerical Simulation code is used. The immersed boundary method is used to represent the roughness in the flow on the surface. The numerical code is verified via theMethod ofManufactured Solutions. Comparisons between experimental data, Linear Stability Theory and numerical results are also used for the validation of the code. Results obtained with different roughness heights and variations in the pressure gradient allowed the analysis of the influence of a three-dimensional roughness element in boundary layer flows.

La propreté inclusionnaire reste un enjeu majeur en élaboration des métaux par voie liquide. La flottation, principal procédé retenu en métallurgie secondaire pour éliminer les particules d’inclusions, consiste à injecter des bulles de gaz au sein du réacteur. Lors de leur ascension, les bulles vont capter les plus grosses inclusions et favoriser la collision et l’agrégation des particules. Dans le but de quantifier les phénomènes influents à l’échelle des inclusions sur la dynamique d’agrégation entre deux inclusions à proximité des bulles, un modèle numérique 3D a été développé. L’écoulement local est modélisé par un cisaillement plan permanent et résolu par une méthode de Boltzmann sur réseau. Le couplage entre les particules et le fluide a été assuré par une méthode de frontière immergée permettant de calculer la perturbation hydrodynamique engendrée par la présence des particules et de mettre à jour les interactions entre particules pour leur suivi lagrangien. Les simulations numériques réalisées ont mis en évidence que les effets hydrodynamiques ont une influence non négligeable sur le comportement des inclusions. Des sections efficaces de collision ont pu être extraites, à partir desquelles ont été calculés des noyaux d’agrégation, données macroscopiques rendant compte des effets à petite échelle. Une première application de ce travail a été menée avec le calcul des fréquences d’agrégation d’un train de bulle dans un réacteur canal à partir de résultats de simulations DNS. Les noyaux d’agrégation ont également été exploités en vue de déterminer, à partir de résultats RANS de l’hydrodynamique d’une poche d’acier, l’évolution de la concentration d’inclusions par un bilan de population global. / Inclusion cleanliness remains a major challenge faced in process metallurgy in liquid phase. Flotation, the main process used in secondary metallurgy to remove inclusions, consists in injecting gas bubbles into the reactor. Rising gas bubbles entrap the biggest inclusions at their surface or in their wake. Besides, they promote collision and aggregation among particles. A 3D numerical model has been developed in order to quantify the roles of the prevailing phenomena on aggregation dynamics between inclusions in the vicinity of bubbles. At inclusion (mesoscopic) scale, the turbulent flow is locally modeled by a steady plane shear flow which is solved using a lattice-Boltzmann method. The coupling between both liquid and solid phases is ensured using an immersed boundary method. This method resolves the hydrodynamic perturbation induced by particles, and hence their interactions that are, in turn, used to update their Lagrangian tracking. The conducted numerical simulations bring out the influence of hydrodynamic effects on inclusion behavior. Collision cross sections have been determined from which ensuing aggregation kernels have been calculated. Such cross sections could provide macroscopic models to represent local particle dynamics. A first application of these results is presented to calculate aggregation frequencies in bubble swarms in a channel flow reactor that was simulated using DNS. Additionally, evolution of inclusion populations in molten steel has been determined from RANS simulation of a liquid steel ladle by means of a global population balance implementing the aggregation kernels determined in the present work.

This thesis presents a numerical method for the phase-resolving simulation of rigid particles and deformable bubbles in viscous, magnetohydrodynamic flows. The presented approach features solid robustness and high numerical efficiency. The implementation is three-dimensional and fully parallel suiting the needs of modern high-performance computing. In addition to the steps towards magnetohydrodynamics, the thesis covers method development with respect to the immersed boundary method which can be summarized in simple words by From rigid spherical particles to deformable bubbles. The development comprises the extension of an existing immersed boundary method to non-spherical particles and very low particle-to-fluid density ratios. A detailed study is dedicated to the complex interaction of particle shape, wake and particle dynamics. Furthermore, the representation of deformable bubble shapes, i.e. the coupling of the bubble shape to the fluid loads, is accounted for. The topic of bubble interaction is surveyed including bubble collision and coalescence and a new coalescence model is introduced. The thesis contains applications of the method to simulations of the rise of a single bubble and a bubble chain in liquid metal with and without magnetic field highlighting the major effects of the field on the bubble dynamics and the flow field. The effect of bubble coalescence is quantified for two closely adjacent bubble chains. A framework for large-scale simulations with many bubbles is provided to study complex multiphase phenomena like bubble-turbulence interaction in an efficient manner.