Cette thèse porte sur l’étude des interactions entre les globules rouges (GRs)et l’endothélium, la couche des cellules qui délimite les vaisseaux sanguins.Il a été démontré que l’endothélium et les GRs jouent des rôles actifs dans divers processus du système vasculaire, et leurs interactions produisent un signal bio chimique grâce à des moyens à la fois chimiques (molécules de signalisation) et mécaniques (taux de cisaillement sur les parois). D’abord,nous étudions le rôle des GRs, y compris dans des conditions pathologiques, dans la création de contraintes de cisaillement spatialement et temporellement dynamiques sur l’endothélium. Il a été montré que les contraintes de cisaillement constituaient un élément critique dans le déclenchement d’un signal bio mécanique depuis l’endothélium. Par ailleurs, étant donné qu’il a été montré que les parois des vaisseaux sanguins ondulent en raison des cellules endothéliales individuelles qui le composent, nous avons intégré à notre modélisation cette géométrie. On trouve que cette ondulation affecte la dynamique des GRs dans l’écoulement ainsi que le taux de cisaillement sur les parois. Nous étudions rapidement dans quelle mesure la déformabilité d’un GR affecte sa trajectoire dans un vaisseau ondulé. Pour cela, nous nous inspirons du processus de fonctionnement un appareil de déplacement latéral déterministe (DLD) qui utilise les variations de trajectoires des particules en fonction de leur taille pour les séparer dans l’écoulement. Nous étudions par ailleurs l’effet des suspensions de GRs sur les caractéristiques rhéologiques et les contraintes de cisaillement sur la paroi du vaisseau.Finalement, nous nous adressons à les interaction chimiques en développons un modèle numérique avec la méthode de Boltzmann sur réseaux-limite immergée (LB-IBM) pour résoudre la diffusion et l’advectiond’un soluté libéré par un particule en mouvement et déformable. L’oxygène et l’adénosine triphosphate (ATP) sont toutes les deux libérées par les GRs,se diffusent dans l'écoulement, et sont absorbées par l’endothélium. Ils représentent des facteurs de signalisation critiques pour les processus de l’inflammation et vasodilatation. Nous montrons que la morphologie des GRs affectera le temps de résidence et la dilution des espèces chimiques lorsqu’elles rentreront en contact avec la paroi du vaisseau. Ensemble, ces éléments nous conduisent vers la développement d’un modèle capable de simuler des processus vitaux du système vasculaire qui résultent d’événements locaux de composants individuels. / This thesis is devoted to the study of the interactions between red blood cells (RBCs) and the endothelium, the monolayer of cells lining blood vessels. The endothelium and RBCs have been shown to be active participants in various processes in the vascular system, and their interactions trigger biochemical signalling by mechanical (wall shear stress) and chemical (signalling molecules) means. We first investigate the role of RBCs, including pathological conditions, in creating time- and space-varying shear stress on the endothelium. Shear stress has been shown to be a critical element in biochemical signalling from the endothelium. In addition, as it has been shown that the endothelium is undulating due to the individual endothelial cells comprising it, we take this into account in our model of the geometry of the vessel wall. We find that this undulation affects the dynamics of the RBCs in the flow and the wall shear stress. We briefly explore how the deformability of a single RBC affects its trajectory in undulating channels, inspired by the idea behind deterministic lateral displacement devices (DLDs) which exploit the differing trajectories of particles based upon their sizes to separate them in flow. We also investigate the effect of suspensions of RBCs in undulating channels on rheological properties and wall shear stress. Finally, we address the chemical interactions by building a numerical model with the lattice Boltzmann-immersed boundary method (LB-IBM) to solve advection-diffusion of solute released from moving, deformable particles. Oxygen and adenosine triphosphate (ATP) are both released by RBCs and are advected and diffused in the flow and uptaken by the endothelium and serve as critical signalling factors in inflammation and vasodilation. We find that the morphology of RBCs will affect the residence time and dilution of the chemical species upon contact with the wall. Together, these elements lead us towards the development of a model capable of simulating vital processes in the vascular system which result from local interactions of individual components.

L’objectif de ce travail est d’étudier les capacités de la méthode de Boltzmann sur réseau (LBM) dans un cadre numériquement contraignant : celui de la simulation aéroacoustique en maillage non-uniforme, à très haut nombre de Reynolds et à nombre de Mach non négligeable (Ma > 0.1), appliquée à l’automobile. La problématique industrielle est celle du calcul du bruit intérieur d’origine aérodynamique, dont le calcul du champ de pression pariétal instationnaire sur le vitrage conducteur est la première étape décisive. Il a été constaté qu’un manque de précision sur la faible part acoustique du champ de pression total sur le vitrage, provenant très probablement d’erreurs au niveau des transitions de résolution du maillage, était la cause d’une surestimation du bruit intérieur. Nous présentons d’abord une construction cohérente et unifiée de la méthode de Boltzmann sur réseau à partir de l’équation de Boltzmann, dans un cadre athermal faiblement compressible. Nous étudions ensuite en détail les propriétés aéroacoustiques de la LBM, en parcourant toutes les grandes familles d’opérateurs de collision de la littérature. Une variante de modèle à temps de relaxation multiples, utilisable pour l’aéroacoustique, est présentée et testée. Un modèle alternatif simplifié de filtrage sélectif, rapide et compact, est développé et validé. La problématique des maillages non-uniformes est abordée. Un recensement exhaustif des études LBM menées dans ce cadre dans la littérature montre qu’aucune ne correspond à nos contraintes. Des algorithmes alternatifs aux transitions sont développés. Enfin, des applications industrielles sont réalisées à l’aide des modèles développés dans le mémoire. / The main goal of this work is to study the capacities of the Lattice Boltzmann Method in a constrained numerical framework : that of numerical simulation in automotive aeroacoustics with non-uniform meshes, at high Reynolds number and non egligible Mach number (Ma > 0.1). The industrial problem is the computation of the interior aerodynamic noise, which includes as its first decisive step the computation of the unsteady wall pressure field on the car windows. It was observed that a lack of precision on the weak acoustic part of the total pressure field on the driver-side window, which is most probably due to errors at mesh refinement interfaces, caused an overestimation of the interior noise. We first present a coherent and unified construction of the Lattice BoltzmannMethod from the Boltzmann equation, in an athermal weakly compressible framework. Then, we study in details the aeroacoustic properties of the LBM by reviewingall the main families of collisional operators that exist in the literature. A variant of multiple relaxation time operator that can be used for aeroacoustics is presented and tested. A simplified alternative selective filter, fast and compact, is developped and numerically validated. The problem of non-uniform meshes is discussed. An exhaustive review of the LBM studies that have been carried out within that framework shows that none of them corresponds to our constraints. Alternative transition nodes algorithms are developed. Finally, all the developed models of this work are applied to industrial cases.

No description available.

Cette étude s'intéresse aux risques associés à la formation d'une ATmosphère EXplosive (ATEX) née d'une fuite d'hydrogène et de sa dispersion dans l'air ambiant. La fuite a été modélisée par un jet turbulent à densité variable libre, impactant sur une sphère de diamètre 20mm ou sur une plaque plane. Dans un premier temps, les champs de vitesses et de concentration ont été obtenus expérimentalement en proche sortie grâce à des mesures de Vélocimétrie par Images de Particules (PIV) et de Fluorescence Induite par Plan Laser sur l'acétone (PLIF). La turbulence et le mélange ont été caractérisés pour le cas d'un jet libre ou en présence d'un obstacle. A partir de ces mesures, la structure générale de l'écoulement a été étudiée à partir des champs moyens et fluctuants par comparaison avec les données de la bibliographie. Puis, les données issues des fluctuations ont été analysées statistiquement par l'étude des fonctions de densité de probabilité du scalaire. Ces travaux se sont poursuivis avec la mise en relation des résultats expérimentaux avec ceux obtenues par des simulations numériques DNS (Direct Numerical Simulation) utilisant la méthode Boltzmann sur Réseau (LBM) d'un scalaire passif dans un jet d'air. Cette étude a permis de recueillir et d'analyser des données supplémentaires sur le mélange d'un jet à masse volumique variable libre ou impactant. Ces données sont directement applicables à la maitrise des risques liés aux fuites d'hydrogène. / This study examines the risks associated with the formation of an explosive atmosphere from a hydrogen leak and its dispersion into the air. We considered the leak as a turbulent jet with density variable, free and impinging a $20,mm$ diameter sphere or a flat plate. Firstly, velocity and scalar fields have been measured experimentally in the near field through Particle Image Velocimetry (PIV) and acetone Planar Laser Induced Fluorescence (LIF). Turbulence and mixing have been defined in the case of free jet and impinging jet. From this measurements, the flow structure has been presented from the mean and fluctuating flow measurements by comparison with literature data. Next, the fluctuation scalar fields are studied with the probability density function method. Finally, a comparison has been conducted between the experiments and direct numerical simulation (DNS) of turbulence based on the lattice Boltzmann method (LBM) for passive scalar in air jet. This study is gathering and analyzing data on the mixing of jet with density variable, free and impinging jet. Such data is directly useful to identify and control risks incurred due to hydrogen leak.

L'objet de cette thèse est la modélisation computationnelle de l'Imagerie par Résonance Magnétique (IRM), appliquée à l'imagerie des réseaux vasculaires. Les images sont influencées par la géométrie des vaisseaux mais aussi par le flux sanguin. Par ailleurs, outre la qualité des modèles développés, il est important que les calculs soient performants. C'est pourquoi, le calcul parallèle est utilisé pour gérer ce type de problèmes complexes. Dans cette thèse, trois solutions sont proposées. La première concerne les algorithmes parallèles pour la modélisation des réseaux vasculaires. Des algorithmes dédiés à différentes architectures sont proposés. Le premier est basé sur le modèle de « passage de messages » pour les machines à mémoires distribuées. La parallélisation concerne l'irrigation de nouvelles zones de tissu par les vaisseaux existants. Le deuxième algorithme est dédié aux machines à mémoire partagée. Il parallélise également le processus de perfusion mais des processeurs différents se chargent de gérer les différents arbres vasculaires. Le troisième algorithme est une combinaison des approches précédentes offrant une solution pour les architectures parallèles hybrides. Les algorithmes proposés permettent d'accélérer considérablement la croissance des réseaux vasculaires complexes, ce qui rend possible la simulation de structures vasculaires plus précises, en un temps raisonnable et aide à améliorer le modèle vasculaire et à tester plus facilement différents jeux de paramètres. Une nouvelle approche de modélisation computationnelle des flux en IRM est également proposée. Elle combine le calcul de flux par la méthode de Lattice Boltzmann, la simulation IRM par le suivi temporel de magnétisations locales, ainsi qu'un nouvel algorithme de transport des magnétisations. Les résultats montrent qu'une telle approche intègre naturellement l'influence du flux dans la modélisation IRM. Contrairement aux travaux de la littérature, aucun mécanisme additionnel n'est nécessaire pour considérer les artéfacts de flux, ce qui offre une grande facilité d'extension du modèle. Les principaux avantages de cette méthode est sa faible complexité computationnelle, son implémentation efficace, qui facilitent le lancement des simulations en utilisant différents paramètres physiologiques ou paramètres d'acquisition des images. La troisième partie du travail de thèse a consisté à appliquer le modèle d'imagerie de flux à des réseaux vasculaires complexes en combinant les modèles de vaisseaux, de flux et d'acquisition IRM. Les algorithmes sont optimisés à tous les niveaux afin d'être performants sur des architectures parallèles. Les possibilités du modèle sont illustrées sur différents cas. Cette démarche de modélisation peut aider à mieux interpréter les images IRM grâce à l'intégration, dans les modèles, de connaissances variées allant de la vascularisation des organes jusqu'à la formation de l'image en passant par les propriétés des flux sanguins. / This PhD thesis concerns computer modeling of magnetic resonance imaging (MRI). The main attention is centered on imaging of vascular structures. Such imaging is influenced not only by vascular geometries but also by blood flow which has to been taken into account in modeling. Next to the question about the quality of developed models, the challenge lies also in the demand for high performance computing. Thus, in order to manage computationally complex problems, parallel computing is in use. In the thesis three solutions are proposed. The first one concerns parallel algorithms of vascular network modeling. Algorithms for different architectures are proposed. The first algorithm is based on the message passing model and thus, it is suited for distributed memory architectures. It parallelizes the process of connecting new parts of tissue to existing vascular structures. The second algorithm is designed for shared memory machines. It also parallelizes the perfusion process, but individual processors perform calculations concerning different vascular trees. The third algorithm combines message passing and shared memory approaches providing solutions for hybrid parallel architectures. Developed algorithms are able to substantially speed up the time-demanded simulations of growth of complex vascular networks. As a result, more elaborate and precise vascular structures can be simulated in a reasonable period of time. It can also help to extend the vascular model and to test multiple sets of parameters. Secondly, a new approach in computational modeling of magnetic resonance (MR) flow imaging is proposed. The approach combines the flow computation by lattice Boltzmann method, MRI simulation by following discrete local magnetizations in time and a new magnetization transport algorithm together. Results demonstrate that such an approach is able to naturally incorporate the flow influence in MRI modeling. As a result, in the proposed model, no additional mechanism (unlike in prior works) is needed to consider flow artifacts, what implies its easy extensibility. In combination with its low computational complexity and efficient implementation, the solution is a user-friendly and manageable at different levels tool which facilitates running series of simulations with different physiological and imaging parameters. The goal of the third solution is to apply the proposed MR flow imaging model on complex vascular networks. To this aim, models of vascular networks, flow behavior and MRI are combined together. In all the model components, computations are adapted to be performed at various parallel architectures. The model potential and possibilities of simulations of flow and MRI in complex vascular structures are shown. The model aims at explaining and exploring MR image formation and appearance by the combined knowledge from many processes and systems, starting from vascular geometry, through flow patterns and ending on imaging technology.

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.

Ce travail de thèse de doctorat, effectué au sein du Département de Chimie Appliquée (au Laboratoire de Génie Chimique, LGC), vise à réaliser létude de lécoulement dun fluide dans des géométries complexes rencontrées en Génie Chimique (les empilages aléatoires et structurés, souvent utilisés dans des colonnes dadsorption et de distillation) par la méthode de Boltzmann sur réseau (Lattice Boltzmann Method, LBM) et à développer un outil numérique de simulation sous la forme dun code de calcul tridimensionnel. Dans le but daméliorer les performances de ces éléments complexes, qui conduisent à une efficacité et une sélectivité accrue des procédés se produisant dans les colonnes de distillation et dadsorption, il est donc intéressant dexaminer avec plus de précision la description des phénomènes de transport se produisant à léchelle locale dans ces milieux. Dans un premier temps, les concepts et les équations macroscopiques qui dérivent de lécoulement dun fluide ont été succinctement passés en revue, ainsi que le modèle de turbulence qui a été utilisé. Nous présentons ensuite, les éléments théoriques de la méthode de Boltzmann sur réseau, après quelques détails sur le modèle Gaz sur réseau, dont elle est dérivée. La méthode de Boltzmann sur réseau se démarque des schémas traditionnels en résolvant numériquement une équation basée sur la physique statistique. Les codes des modèles SRTLBE (équation de Boltzmann sur réseau à un seul temps de relaxation) et MRTLBE (équation de Boltzmann sur réseau à plusieurs temps de relaxation) à trois dimensions, de même que des conditions aux limites convenables ont été écrits dans un langage de programmation, puis validées sur des cas types de la littérature : un écoulement laminaire dans un conduit de section carrée et un écoulement de jet turbulent rectangulaire à 3D. Les codes implémentés dans une grille de calcul ont permis de simuler des écoulements de fluides dans des structures très complexes qui ont été obtenues par une technique expérimentale avancée : la microtomographie à rayons X. Létude dune mousse métallique, dun filtre à charbons actifs et dun filtre de capture daérosols ont été présentés. Les résultats numériques ont été comparés avec succès aux mesures expérimentales et à des corrélations de la littérature. Une autre étude a concerné létude de lécoulement à travers deux feuilles de Mellapak 250 Y. Les champs et contours de vitesse ainsi obtenus ont été soigneusement analysés et commentés à la fois en régime laminaire et en régime turbulent. Les résultats hydrodynamiques ont été comparés avec des mesures expérimentales ainsi quavec des valeurs calculées par des corrélations disponibles dans la littérature. Les modèles existant ont été analysés. Enfin, ces simulations (LBM) ont été comparées à celles effectuées dans les mêmes conditions opératoires, avec le code de calcul classique FLUENT.

No description available.

Cette thèse porte sur l'étude numérique des fluctuations extrêmes de la force de traînée exercée par un écoulement turbulent sur un corps immergé.Ce type d'évènement, très rare, est difficile à caractériser par le biais d'un échantillonnage direct, puisqu'il est alors nécessaire de simuler l'écoulement sur des durées extrêmement longues. Cette thèse propose une approche différente, basée sur l'application d'algorithmes d'échantillonnage d'événements rares. L'objectif de ces algorithmes, issus de la physique statistique, est de modifier la statistique d'échantillonnage des trajectoires d'un système dynamique, de manière à favoriser l'occurrence d'événements rares. Si ces techniques ont été appliquées avec succès dans le cas de dynamiques relativement simples, l'intérêt de ces algorithmes n'est à ce jour pas clair pour des dynamiques déterministes extrêmement complexes, comme c'est le cas pour les écoulement turbulents.Cette thèse présente tout d'abord une étude de la dynamique et de la statistique associée aux fluctuations extrêmes de la force de traînée sur un obstacle carré fixe immergé dans un écoulement turbulent à deux dimensions. Ce cadre simplifié permet de simuler la dynamique sur des durées très longues, permettant d'échantillonner un grand nombre de fluctuations dont l'amplitude est assez élevée pour être qualifiée d'extrême.Dans un second temps, l'application de deux algorithmes d’échantillonnage est présentée et discutée.Dans un premier cas, il est illustré qu'une réduction significative du temps de calcul d'extrêmes peut être obtenue. En outre, des difficultés liées à la dynamique de l'écoulement sont mises en lumière, ouvrant la voie au développement de nouveaux algorithmes spécifiques aux écoulements turbulents. / This thesis discusses the numerical simulation of extreme fluctuations of the drag force acting on an object immersed in a turbulent medium.Because such fluctuations are rare events, they are particularly difficult to investigate by means of direct sampling. Indeed, such approach requires to simulate the dynamics over extremely long durations.In this work an alternative route is introduced, based on rare events algorithms.The underlying idea of such algorithms is to modify the sampling statistics so as to favour rare trajectories of the dynamical system of interest.These techniques recently led to impressive results for relatively simple dynamics. However, it is not clear yet if such algorithms are useful for complex deterministic dynamics, such as turbulent flows.This thesis focuses on the study of both the dynamics and statistics of extreme fluctuations of the drag experienced by a square cylinder mounted in a two-dimensional channel flow.This simple framework allows for very long simulations of the dynamics, thus leading to the sampling of a large number of events with an amplitude large enough so as they can be considered extreme.Subsequently, the application of two different rare events algorithms is presented and discussed.In the first case, a drastic reduction of the computational cost required to sample configurations resulting in extreme fluctuations is achieved.Furthermore, several difficulties related to the flow dynamics are highlighted, paving the way to novel approaches specifically designed to turbulent flows.

Submitted by Geandra Rodrigues (geandrar@gmail.com) on 2018-01-11T14:54:52Z No. of bitstreams: 1 andresricardovaldez.pdf: 6592036 bytes, checksum: 23a86a3d84f13bffa421f219e7e4501d (MD5) / Rejected by Fabíola Rubim (fabiola.rubim@ufjf.edu.br), reason: on 2018-01-12T11:05:10Z (GMT) / Submitted by Geandra Rodrigues (geandrar@gmail.com) on 2018-01-12T11:46:32Z No. of bitstreams: 1 andresricardovaldez.pdf: 6592036 bytes, checksum: 23a86a3d84f13bffa421f219e7e4501d (MD5) / Rejected by Adriana Oliveira (adriana.oliveira@ufjf.edu.br), reason: Favor corrigir: Membro da banca: Filho, José Karam on 2018-01-23T14:01:35Z (GMT) / Submitted by Geandra Rodrigues (geandrar@gmail.com) on 2018-01-23T14:06:58Z No. of bitstreams: 1 andresricardovaldez.pdf: 6592036 bytes, checksum: 23a86a3d84f13bffa421f219e7e4501d (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2018-01-23T14:22:22Z (GMT) No. of bitstreams: 1 andresricardovaldez.pdf: 6592036 bytes, checksum: 23a86a3d84f13bffa421f219e7e4501d (MD5) / Made available in DSpace on 2018-01-23T14:22:22Z (GMT). No. of bitstreams: 1 andresricardovaldez.pdf: 6592036 bytes, checksum: 23a86a3d84f13bffa421f219e7e4501d (MD5) Previous issue date: 2017-09-18 / Neste trabalho são apresentados aspectos de modelagem computacional para o estudo de Interação Fluido-Estrutura (FSI). Numericamente, o Método de Lattice Boltzmann (LBM) é usado para resolver a mecânica dos fluidos, em particular as equações de Navier-Stokes incompressíveis. Neste contexto, são abordados problemas de escoamentos complexos, caracterizado pela presença de obstáculos. A imposição das restrições na interface fluido-sólido é feita utilizando princípios variacionais, empregando o Princípio de Balanço de Potências Virtuais (PVPB) para obter as equações de Euler-Lagrange. Esta metodologia permite determinar as dependências entre carregamentos cinematicamente compatíveis e o estado mecânico adotado. Neste sentido, as condições de interface fluido-sólido são abordadas pelo Método de Fronteira Imersa (IBM) visando técnicas computacionais de baixo custo. A metodologia IBM trata o equilíbrio das equações na interface fluido-sólido através da interpolação entre os nós Lagrangianos (sólidos) e os nós Eulerianos (fluidos). Neste contexto, uma modificação desta estratégia que fornece soluções mais precisas é estudada. Para mostrar as capacidades do acoplamento LBM-IBM são apresentados vários experimentos computacionais que demonstram grande fidelidade entre as soluções obtidas e as soluções disponíveis na literatura. / This work presents computational modeling aspects for studying Fluid-Structure Interaction (FSI). The Lattice Boltzmann Method (LBM) is employed to solve the fluid mechanics considering the incompressible Navier-Stokes equations. The flows studied are complex due to the presence of arbitrary shaped obstacles. The obstacles alters the bulk flow adding complexity to the analysis. In this work the Euler-Lagrange equations are obtained employing the Principle of Virtual Power Balance (PVPB). Consequently, the functional dependencies between the mechanical state and every kinematic compatible loadings are established employing variational arguments. This modeling technique allows to study the fluid-solid boundary constraint. In this context the fluid-solid interface is handled employing the Immersed Boundary Method (IBM). The IBM deals with the fluid-solid interface equilibrium equations performing an interpolation of forces between Lagrangian nodes (solid domain) and Eulerian Lattice grid (fluid domain). In this work a different version of this methodology is studied that allows to obtain more accurate solutions. To show the capabilities of the implemented LBM-IBM solver several experiments are done showing the agreement with the benchmarks results available in literature.