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.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

This work presents the coupling between two different formulations applied to flow simulation and analysis of wind rotors, integral and differential formulations. First, for the integral formulation is defined a control volume where the variables problem are defined, as well as the necessaries working hypothesis, then a proposed mathematical modeling is defined. Simulations through NACA airfoils, using Computational Dynamic Fluids, are performed in order to evaluate drag and lift coefficients, to be used in the integral formulation. The Navier-Stokes equations are solved in house and the Smagorinsky turbulence model with Van Driest damping function is retained. The computational code is implemented with structured cartesian mesh, where the airfoil is modeled using the Immersed Boundary Methodology. The results of simulation through a NACA0012 airfoil are shown for several attack angles and Re = 10000. Results of energetic efficiency are presented for a horizontal axis wind turbine using the integral formulation, where the coefficients are given by differential formulations. / Este trabalho apresenta o acoplamento entre as duas formulações diferentes aplicadas à simulação do escoamento e análise de rotores eólicos, formulações integral e diferencial. Primeiramente para a formulação integral é definido um volume de controle onde as variáveis do problema são definidas, bem como as hipóteses simplificadoras necessárias, para então ser proposta uma modelagem matemática. Simulações do escoamento em aerofólios NACA, utilizando Dinâmica dos Fluidos Computacional, são realizadas a fim de determinar os coeficientes de arrasto e sustentação, os quais foram utilizados na formulação integral. As equações de Navier-Stokes são resolvidas em um código computacional e o modelo de turbulência de Smagorinsky com função de amortecimento de Van Driest é utilizado. O código computacional é implementado com uma malha cartesiana estruturada, e o aerofólio é modelado utilizando a Metodologia da Fronteira Imersa. Os resultados da simulação através de um aerofólio NACA0012 são mostrados para vários ângulos de ataque e Re = 10000. Os resultados de eficiência energética são apresentados para uma turbina eólica de eixo horizontal utilizando a formulação integral, onde os coeficientes são dados pelas formulações diferenciais. / Mestre em Engenharia Mecânica

La méthode LS-STAG est une méthode de type frontière immergée/cut-cell pour le calcul d’écoulements visqueux incompressibles qui est basée sur la méthode MAC pour grilles cartésiennes décalées, où la frontière irrégulière est nettement représentée par sa fonction level-set, résultant en un gain significatif en ressources informatiques par rapport aux codes MFN commerciaux utilisant des maillages qui épousent la géométrie. La version 2D est maintenant bien établie et ce manuscrit présente son extension aux géométries 3D avec une symétrie translationnelle dans la direction z (configurations extrudées 3D). Cette étape intermédiaire sera considérée comme la clé de voûte du solveur 3D complet, puisque les problèmes de discrétisation et d’implémentation sur les machines à mémoire distribuée sont abordés à ce stade de développement. La méthode LS-STAG est ensuite appliquée à divers écoulements newtoniens et non-newtoniens dans des géométries extrudées 3D (conduite axisymétrique, cylindre circulaire, conduite cylindrique avec élargissement brusque, etc.) pour lesquels des résultats de références et des données expérimentales sont disponibles. Le but de ces investigations est d’évaluer la précision de la méthode LS-STAG, d’évaluer la polyvalence de la méthode pour les applications d’écoulement dans différents régimes (fluides newtoniens et rhéofluidifiants, écoulement laminaires stationnaires et instationnaires, écoulements granulaires) et de comparer ses performances avec de méthodes numériques bien établies (méthodes non structurées et de frontières immergées) / The LS-STAG method is an immersed boundary/cut-cell method for viscous incompressible flows based on the staggered MAC arrangement for Cartesian grids where the irregular boundary is sharply represented by its level-set function. This approach results in a significant gain in computer resources compared to commercial body-fitted CFD codes. The 2D version of LS-STAG method is now well-established and this manuscript presents its extension to 3D geometries with translational symmetry in the z direction (3D extruded configurations). This intermediate step will be regarded as the milestone for the full 3D solver, since both discretization and implementation issues on distributed memory machines are tackled at this stage of development. The LS-STAG method is then applied to Newtonian and non-Newtonian flows in 3D extruded geometries (axisymmetric pipe, circular cylinder, duct with an abrupt expansion, etc.) for which benchmark results and experimental data are available. The purpose of these investigations is to evaluate the accuracy of LS-STAG method, to assess the versatility of method for flow applications at various regimes (Newtonian and shear-thinning fluids, steady and unsteady laminar to turbulent flows, granular flows) and to compare its performance with well-established numerical methods (body-fitted and immersed boundary methods)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In this work, the combined multi-direct forcing and immersed boundary method (IBM) were presented to simulate uid-structure interaction problems. The multi-direct forcing is used aim at satisfying the no-slip condition in the immersed boundary. For the numerical simulations was used a multi-purpose computer code that is being developed in the MFlab - Fluid Mechanics Laboratory of Federal University of Uberl^andia. Tests are made to validate the numerical schemes and routines were implemented to simulate uid-structures interaction problems. Furthermore, computational tools are developed to construct and manage and optimize the use of a Beowulf cluster where all the parallel simulations presented in this work were done. The Method of Manufactured Solutions has been used for order-of-accuracy verication in the computational uid dynamics code. Two uid-structure interaction problems were studied using this methodology. The rst is a ow over a sphere for some Reynolds numbers. The results were compared to empirical results, obtaining satisfactory approximations. The second one is a immersed simple pendulum. For this problem the results are in agreement with physics. Indeed, these are preliminar results. New tests must be done to make progress in the methodology. Improvements are proposed in the IBM, in the uid-structure model, in the turbulence model, in the method used to discretize the uid domain. It is also proposed to apply the methodology to real problems as risers and valves. / O presente trabalho tem como principal objetivo a aplicação do método multifoçagem (MMF) para solução numérica tridimensional de problemas de interação uidoestrutura, buscando-se garantir a condição de não-escorregamento na região da fronteira imersa. Para as simulações numéricas foi utilizado um código computacional multipropósito em desenvolvimento no MFlab - Laboratório de Mecânica dos Fluidos da Universidade Federal de Uberlândia. Foram feitas modificações nesse código para que se pudesse validá-lo para solução de problemas com fronteira imersa e foi implementada uma rotina para solução de um problema de interação uido-estrutura total. Além disso, foi desenvolvido um pacote de ferramentas computacionais que possibilitou instalar e melhorar o desempenho de um cluster do tipo Beowulf utilizado para o desenvolvimento das simulações num eriças em paralelo do presente trabalho. Utilizando o Método das Soluções Manufaturadas foram obtidas soluções sintetizadas para as equações de Navier-Stokes, o que possibilitou obter a ordem de convergência numérica do código computacional para problemas contínuos e a validação deste código para problemas envolvendo corpos imersos ao combinar a o método das soluções manufaturadas com a metodologia de fronteira imersa. Na sequência foi solucionado o problema de escoamento ao redor de uma esfera parada, cujos resultados foram comparados com referencias empíricas, obtendo-se boa aproximação. Ainda para esse caso foi feita a avalição da norma L2 para as soluções num eriças obtidas nos pontos lagrangianos verificando a garantia da condição de não-escorregamento e feita uma análise da inuência dos número de ciclos utilizados no método multi-forçagem. Foi vericado que a solução numérica obtida depende do número de ciclos o que faz com que seja necessário se estabelecer um critério de convergência para este método. Um segundo problema de interação uido-estrutura total foi estudado. Consiste em um pêndulo simples imerso em um uido que parte de uma dada posição angular inicial e oscila em torno da sua posição de equilíbrio, até parar. Para esse caso foram feitas análises quantitativas. Os resultados são preliminares mas coerentes com a física do problema, indicando que a metodologia é adequada para solução deste tipo de problema. / Mestre em Engenharia Mecânica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Modern engineering increasingly requires the comprehension of phenomena related to combustion, aeroacustics, turbulence transition, among others. For these purposes the Computational Fluids Dynamics (CFD) requires the used high order methods. One of these methods is the Fourier pseudo-spectral method, that provides an excellent numerical accuracy, and with the use of the Fast Fourier Transform algorithm (FFT), it presents a low computational cost in comparison to anothers high-order methods. Another important issue is the projection method of the pression term, which does not require the pressure computation from the Navier-Stokes equations. The procedure to calculate the pression field is usually the most onerous in classical methodologies. Nevertheless, the pseudo-spectral method can be only applied to periodic boundary flows, thus limiting its use. Aiming to solve this restriction, a new methodology is proposed at the present work, which has the objective of simulating nonperiodic flows using the Fourier pseudo-spectral method. For this purpose the immersed boundary method, that represents the boundary conditions through a force field imposed at Navier-Stokes equations is used. As a test to this new methodology, a classic problem of Computational Fluid Dynamics, The Lid Driven Cavity was simulated. The obtained results are promising and demonstrate the possibility to simulating non-periodic flows making use of the Fourier pseudo-spectral method. / Para compreender fenômenos relacionados à combustão, aeroacústica, transição a turbulência entre outros, a Dinâmica de Fluídos Computacional (CFD) utiliza os métodos de alta ordem. Um dos mais conhecidos é o método pseudo-espectral de Fourier, o qual alia: alta ordem de precisão na resolução das equações, com um baixo custo computacional. Este está ligado à utilização da FFT e do método da projeção do termo da pressão, o qual desvincula os cálculos da pressão da resolução das equações de Navier-Stokes. O procedimento de calcular o campo de pressão, normalmente é o mais oneroso nas metodologias convencionais. Apesar destas vantagens, o método pseudo-espectral de Fourier só pode ser utilizado para resolver problemas com condições de contorno periódicas, limitando o seu uso no campo da dinâmica de fluídos. Visando resolver essa restrição uma nova metodologia é proposta no presente trabalho, que tem como objetivo simular escoamentos não-periódicos utilizando o método pseudo-espectral de Fourier. Para isso, é utilizada a metodologia da Fronteira Imersa, a qual representa as condições de contorno de um escoamento através de um campo de força imposto nas equações de Navier-Stokes. Como teste, para essa nova metodologia, foi simulada uma cavidade com tampa deslizante (Lid Driven Cavity), problema clássico da mecânica de fluídos, que objetiva validar novas metodologias e códigos computacionais. Os resultados obtidos são promissores e demostram que é possível simular um escoamento não-periódico com o método pseudo-espectral de Fourier. / Mestre em Engenharia Mecânica

Este trabalho consiste no desenvolvimento de ferramentas computacionais para análise não linear geométrica de interação fluido-casca utilizando o Método dos Elementos Finitos (MEF). O algoritmo para dinâmica dos fluidos é explícito e a integração temporal é baseada em linhas características. O código computacional é capaz de simular as equações de Navier-Stokes para escoamentos compressíveis tanto na descrição Euleriana como na descrição Lagrangeana-Euleriana arbitrária (ALE), na qual é possível prescrever movimentos para a malha do fluido. A estrutura é modelada em descrição Lagrangeana total através de uma formulação de MEF para análise dinâmica não linear geométrica de cascas baseada no teorema da mínima energia potencial total escrito em função das posições nodais e vetores generalizados e não em deslocamentos e rotações. Essa característica evita o uso de aproximações de grandes rotações. Dois modelos de acoplamentos são desenvolvidos. O primeiro modelo, ideal para problemas onde a escala de deslocamentos não é muito grande comparada com as dimensões do domínio do fluido, é baseado na descrição ALE e o acoplamento entre as duas diferentes malhas é feito através do mapeamento das posições locais dos nós do contorno do fluido sobre os elementos de casca e vice-versa, evitando a necessidade de coincidência entre os nós da casca e do fluido. A malha do fluido é adaptada dinamicamente usando um procedimento simples baseado nas posições e velocidades nodais da casca. O segundo modelo de acoplamento, ideal para problemas com grande escala de deslocamentos tais como estruturas infláveis, considera a casca imersa na malha do fluido e consiste em um procedimento robusto baseado em curvas de nível da função distância assinalada do contorno, o qual integra o algoritmo Lagrangeano de casca com o Fluido em descrição Euleriana, sem necessidade de movimentação da malha do fluido, onde a representação computacional do fluido se resume a uma malha não estruturada maior ou igual ao domínio inicial do fluido e a interface fluido-casca dentro da malha do fluido é identificada por meio de curvas de nível da função distância assinalada do contorno. Ambos os modelos são testados através de exemplos numéricos mostrando robustez e eficiência. Finalmente, como uma sugestão para o futuro desenvolvimento desta pesquisa, iniciaram-se estudos relativos a funções B-splines. O uso desse tipo de funções deverá resolver problemas de estabilidade relativos a oscilações espúrias devidas ao uso de polinômios de Lagrange para a representação de descontinuidades. / This work consists of the development of computational tools for nonlinear geometric fluid-shell interaction analysis using the Finite Element Method (FEM). The fluid solver is explicit and its time integration based on characteristics. The computational code is able to simulate the Navier-Stokes equations for compressible flows written in the Eulerian description as well as in the arbitrary Lagrangian-Eulerian (ALE) description, enabling movements prescription for the fluid mesh. The structure is modeled in a total Lagrangian description, using a FEM formulation to deal with geometrical nonlinear dynamics of shells based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors, not displacements and rotations, avoiding the use of large rotation approximations. Two partitioned coupling models are developed. The first model, ideal for simulations where the displacements scale is not very large compared to the fluid domain, is based on the ALE description and the coupling between the two different meshes is done by mapping the fluid boundary nodes local positions over the shell elements and vice-versa, avoiding the need for matching fluid and shell nodes. The fluid mesh is adapted using a simple approach based on shell nodal positions and velocities. The second model, ideal for problems with large scales of displacements such as inflatable structures, is based on immersed boundary and consists of a robust level-set based approach that integrates the Lagrangian shell finite and the Eulerian finite element high speed fluid flow solver, with no need for mesh adaptation, where the fluid representation relies on a fixed unstructured mesh larger or equal to the initial fluid domain and the fluid-shell interface inside the fluid mesh is tracked with level sets of a boundary signed distance function. Both models are tested with numerical examples, showing efficiency and robustness. Finally, as a suggestion for future development of this research, we started studies relatives to B-Spline functions. The use of this kind of functions should solve stability problems related to spurious oscillations due to the use of Lagrange polynomials for representing discontinuities.

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.

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.