La simulation de scenarios d’incendie nécessite de modéliser de nombreux processus complexe, particulièrement la combustion gazeuse d’hydrocarbure incluant la production de suie et les transferts radiatifs dans un écoulement turbulent. La nature turbulente de l’écoulement fait apparaitre des interactions qui doivent être prises en compte entre ces processus. L’objectif de cette thèse est d’implémenter une méthode de transport de la fonction de densité de probabilité afin de modéliser ces interactions de manière précise. En conjonction avec un modèle de flammelettes, le modèle de Lindstedt et un modèle à large-bande k-corrélé, l’équation de transport de la PDF jointe de composition est résolue avec la méthode des Champs Eulérien Stochastiques. Le modèle est validé en simulant 12 flammes turbulentes recouvrant une large gamme de nombre de Reynolds et de propension à former de la suie par les combustibles. Dans un second temps, les effets des interactions rayonnement-turbulence (TRI) sur l’émission de la suie sont étudiés en détails, montrant que la TRI tend à augmenter l’émission radiative de la suie à cause des fluctuations de température, mais que cette augmentation est plus faible pour des nombres de Reynolds élevés ou des quantités de suie plus élevées. Ceci est dû à la corrélation négative entre le coefficient d’absorption des suies et la fonction de Planck. Finalement, l’influence de la corrélation entre la fraction de mélange et le paramètre de non-adiabaticité est étudiée sur une flamme d’éthylène, montrant qu’elle a peu d’effet sur la structure moyenne de flamme mais tend à limiter les fluctuations de température et les pertes radiatives. / The simulation of fire scenarios requires the numerical modeling of various complex process, particularly the gaseous combustion of hydrocarbons including soot production and radiative transfers in a turbulent. The turbulent nature of the flow induces interactions between these processes that need to be taken accurately into account. The purpose of this thesis is to implement a transported Probability Density function method to model these interactions precisely. In conjunction with the flamelet model, the Lindstedt model, and a wide-band correlated-k model, the composition joint-PDF transport equation is solved using the Stochastic Eulerian Fields method. The model is validated by simulating 12 turbulent jet flames covering a large range of Reynolds numbers and fuel sooting propensity. Model prediction are found to be in reasonable agreement with experimental data. Second, the effects of turbulence-radiation interactions (TRI) on soot emission are studied in details, showing that TRI tends to increase soot radiative emission due to temperature fluctuations, but that this increase is smaller for higher Reynolds numbers and higher soot loads. This is due to the negative correlation between soot absorption coefficient and the Planck function. Finally, the effects of taking into account the correlation between mixture fraction and enthalpy defect on flame structure and radiative characteristics are also studied on an ethylene flame, showing that it has weak effect on the mean flame structure but tends to inhibit both temperature fluctuations and radiative loss.

De nos jours, la combustion représente encore un 90% de la production totale d'énergie au monde. La plupart des brûleurs de type industriel utilisent comme carburant des hydrocarbures en forme liquide. Cependant, un grand nombre d'études ont été dédiés aux flammes gazeuses et l'impact du spray liquide est encore loin d'être totalement compris. Le but de cet étude est l'amélioration de la modélisation des deux phénomènes principaux qui ont lieu entre l'atomisation du spray et la combustion, i.e. la dispersion des gouttes par la turbulence gazeuse et le procès d'évaporation dans le contexte de la Simulation Aux Grandes Echelles (SGE) des configurations complexes. Premièrement, l'approche Euler-Euler mésoscopique (Février et al. (2005)), basée sur une moyenne d'ensemble conditionnée et implémentée dans AVBP est amélioré. Le modèle de fermeture (Simonin et al. (2001); Kaufmann (2004)) pour les moments de deuxième ordre qui apparait dans les équations de transport résolues échoue quand appliqué à des configurations cisaillées (Riber (2007)). Plusieurs modèles proposés récemment par Masi (2010) et qui ont été valides a priori dans une configuration de nappe chargée de particules sont validés a posteriori dans la même configuration. Un analyse quantitative sur plusieurs cas avec diffèrent nombres de Stokes, nombres de Reynolds de la phase gazeuse et résolutions du maillage ont permit de retenir un modèle non-linéaire nommé 2EASM3, qui utilise le tenseur de déformations de la phase dispersée comme échelle de temps caractéristique. La deuxième partie a pour but l'amélioration du modèle d'évaporation implémenté dans AVBP. Ce modèle suppose une conduction infinie dans la phase liquide et symétrie sphérique dans la phase gazeuse ainsi que des lois simplifiées pour les propriétés thermodynamiques et de transport. Un nouveau modèle prenant en compte la dépendance de la viscosité du mélange gazeux avec la composition locale, et des nombres de Prandtl et Schmidt fixés par les valeurs à l'équilibre obtenus par moyen d'une simulation prenant en compte des lois complèxes pour les propriétés thermodynamiques et de transport est proposé. Cette nouvelle méthode produit des résultats en bon accord avec les mesures expérimentales pour l'évaporation d'une goutte isolé en une atmosphère d'azote au calme sans pourtant augmenter le cout du calcul. Finalement, l'impacte des nouveaux modèles est analysé dans une SGE de la configuration semi-industrielle MERCATO (García-Rosa (2008)). Bien que les données expérimentales ne soient pas suffisantes pour confirmer les résultats, les distributions de gouttes et de carburant gazeux sont significativement affectés par les modèles, ce qui pourrait avoir un impact directe sur le procès d'allumage. / Combustion still represents about 90% of the energy production in the world. Most industrial burners are fuelled with liquid hydrocarbons. However, most studies have been dedicated to gaseous ßames and the impact of liquid spray is still misunderstood. The purpose of this study is to improve the modelisation of two main phenomena occurring between atomization and combustion, i.e. the droplet dispersion in the turbulent gaseous flow and the evaporation process, in the context of Large Eddy Simulation (LES) of complex configurations. First, the mesoscopic Euler-Euler approach (Février et al. (2005)) based on a conditioned ensemble averaging and implemented in AVBP is improved. The closure model (Simonin et al. (2001), Kaufmann (2004)) for the second-order moments appearing in the transport equations solved fails in mean-sheared configurations (Riber (2007)). Several new models proposed by Masi (2010) and a priori tested in a particle-laden slab are tested a posteriori in the same configuration. A quantitative analysis based on several calculations varying the Stokes number, the gaseous Reynolds number and the grid resolution allows to retain a non-linear model using the particle rate-of-strain tensor as timescale and called 2EASM3. The second part consists in improving the evaporation model implemented in AVBP which assumes infinite conduction in the liquid and spherical symmetry in the gas phase along with simplified thermodynamics and transport properties calculation. A new model is proposed, where the dependence of gaseous mixture viscosity on local composition is accounted for, and the Prandtl and Schmidt numbers are fixed by a reference equilibrium calculation using complex thermodynamics and transport properties. This method shows good agreement with experimental measurements in the configuration of an isolated droplet evaporating in quiescent N2 without further increasing the computational cost. Finally, the impact of the new models is analysed in the LES of the MERCATO semi-industrial configuration (García-Rosa (2008)). Although the experimental data are not sufficient to confirm the results, both the droplet distribution and the fuel mass fraction are significantly affected, which would eventually affect the ignition process.

Une demande mondiale grandissante en gaz naturel pousse à exploiter des ressources de plus en plus acides (concentration des gaz acides pouvant atteindre 20% en volume). Le procédé de traitement de gaz par des solutions d'amines, existant depuis plus de 50 ans, est le plus répandu pour éliminer les composants acides (CO2, H2S, mercaptans…). Ce procédé comporte deux unités principales : une colonne d'absorption gaz liquide à contre-courant pour la séparation des composants acides du gaz et une colonne de régénération du solvant chargé. On recherche une optimisation du design, en particulier de l'absorbeur, en vue de réduire les couts et d'augmenter l'efficacité. Dans la colonne d'absorption, le transfert de matière s'avère limitant par rapport à la thermodynamique. Ainsi l'optimisation du design de cette unité passe par une maitrise des paramètres hydrodynamiques et de transfert des contacteurs utilisés dans les colonnes. Plusieurs études existantes sur unités pilotes proposent des corrélations majoritairement empiriques pour les paramètres critiques de dimensionnement. Par ailleurs leur extrapolation à l'échelle industrielle montre une divergence importante entre les différentes corrélations. Une meilleure compréhension des phénomènes physiques ainsi qu'une identification des paramètres importants pour l'extrapolation est donc requise. Dans ce contexte, la présente thèse a comme objectif principal la compréhension de l'effet de changement d'échelle sur les paramètres hydrodynamiques et l'aire interfaciale d'échange dans le cas des plateaux à clapets. La méthodologie employée dans cette étude se base sur une complémentarité entre une étude expérimentale et le potentiel offert par les outils de simulation numérique. L'étude expérimentale s'est effectuée sur deux colonnes pilotes rectangulaires transparentes ayant deux longueurs de passe différentes. Des mesures de pertes de charge, de hauteur de l'émulsion, de rétention liquide et d'aire interfaciale d'échange ont été réalisées. Des mesures innovantes de profils de l'émulsion gaz-liquide sur un plateau sont également présentées. Les différents résultats expérimentaux ont permis la proposition d'un diagramme hydrodynamique ainsi qu'une compréhension et une analyse phénoménologique cohérente de l'écoulement sur une large gamme de vitesses liquide et gaz. La comparaison entre les deux colonnes a permis, en premier lieu, l'identification des vitesses liquide et gaz pertinentes pour l'extrapolation. Des similitudes de comportement ont été trouvées pour certains paramètres (rétention liquide moyenne, perte de charge clapets, aire interfaciale) offrant ainsi la possibilité de proposer des corrélations basées sur une description phénoménologique dépendant essentiellement de deux nombres adimensionnels que sont le nombre de Froude (comparant l’inertie gaz au poids liquide sur le plateau) et le paramètre de l’écoulement (comparant les deux inerties liquide et gaz). En revanche une influence notable de la longueur de passe est relevée. En particulier des profils de l'émulsion nettement différents entre la petite et la grande colonne ont été observés. Des risques d'extrapolation sont par conséquent pointés dans cette thèse notamment pour des paramètres tels que la hauteur de liquide clair ou la hauteur moyenne de l’émulsion. Dans une deuxième partie, l'intérêt a été porté sur la simulation numérique des écoulements sur les plateaux. L'importance et la complexité de la modélisation du terme d'interaction entre les deux phases sont soulignées. Une approche proposée dans la littérature a été testée et montre la possibilité de l'emploi des simulations CFD comme outil pour une meilleure compréhension des comportements locaux. En outre une approche de modélisation nouvelle est proposée dans une optique de valorisation des outils numériques pour l'extrapolation. / The increasing demand of natural gas has encouraged the exploration of sour ressources (concentration of acid gases (CO2, H2S, mercaptans…) reaching 20% of volume). The most popular gas sweetening process which uses amine solutions has been in use for over 50 years. This process consists of two main units: a counter current gas-liquid absorption column in which acid compounds are removed from the gas and a stripper column for loaded solvent regeneration. An optimisation is needed to reduce the over-sizing costs and enhance efficiency especially for the absorber column. For the absorption column the mass transfer is the limiting phenomenon in comparison to thermodynamics. Thus design optimisation of such units needs full knowledge of both hydrodynamic and mass transfer parameters of the contactors used in the absorbers. Several literature studies carried out on pilot units propose empirical correlations for these critical design parameters. However the extrapolation to industrial scale shows important divergences between existing correlations. It seems that a better understanding of physical phenomena as well as an identification of the key parameters for extrapolation is still needed. Under this context, the present work aims at understanding the scale effect on hydrodynamics and interfacial area on valve trays contactors. The methodology employed relies on the complementarity between experimental studies and the potential offered by numerical simulation tools. The experimental study has been carried out on two rectangular and transparent pilot columns having two different path lengths. Measurements of pressure drops, liquid retention, emulsion height and interfacial area along with innovating measurements of emulsion profiles have been made. Experimental results allowed the establishment of a hydrodynamic diagram as well as the understanding and the phenomenological analysis of the two phase flow over a large scope of gas and liquid velocities. Comparison between the two columns led, at first, to identifying pertinent liquid and gas velocities for extrapolation. Similarities between both columns were noticed for some parameters (liquid retention, valve pressure drop and interfacial area) allowing the proposition of correlations based on phenomenological description of the flow. These correlations involve mainly two non-dimensional numbers which are Froude number (comparing the gas inertia to liquid weight on the tray) and flow parameter (comparing liquid inertia to gas inertia). However an effect of path length has been observed in particular for emulsion profiles. Risks for extrapolation related to this effect have been consequently pointed out for some parameters such as clear liquid height or mean emulsion height. In a second part of the study, the interest was focused on numerical simulation of tray hydrodynamics. Importance and complexity of gas liquid interaction term modelling have been highlighted. Using a proposed approach in literature, CFD simulations with the Eulerian model under the software Ansys Fluent have been made. They proved the possibility of using numerical simulation as a mean for a better understanding of the two phase flow local behaviour. Moreover a new approach is proposed which could be used in order to make numerical tools more valuable for extrapolation.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos. / A two-phase flow computational code taking a hybrid Arbitrary Lagrangian-Eulerian desciption of movement along with the Finite Element Method is used to simulate the dynamics of an incompressible drop jet in crossflow in the primary breakup zone. Dispersed entities are described by means of a front-tracking method which produces zero-thickness interfaces through contiguous element meshing and adaptive refinement techniques. Periodic boundary conditions are implemented in a variationally consistent way for all the scalar fields involved in the presented simulations and amodified version of the pressure field is added to the one-fluid formulation employed in the momentum equation. Three-dimensional direct numerical simulations for different flow configurations of immiscible liquids pertinent to experimental results found in literature. Analyses of the hydrodynamics of the drop jet in crossflow in these configurations considering trajectories, drop shape variations, spectrum of small disturbances, besides additional aspects relating to mesh quality are presented and discussed.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos. / A two-phase flow computational code taking a hybrid Arbitrary Lagrangian-Eulerian desciption of movement along with the Finite Element Method is used to simulate the dynamics of an incompressible drop jet in crossflow in the primary breakup zone. Dispersed entities are described by means of a front-tracking method which produces zero-thickness interfaces through contiguous element meshing and adaptive refinement techniques. Periodic boundary conditions are implemented in a variationally consistent way for all the scalar fields involved in the presented simulations and amodified version of the pressure field is added to the one-fluid formulation employed in the momentum equation. Three-dimensional direct numerical simulations for different flow configurations of immiscible liquids pertinent to experimental results found in literature. Analyses of the hydrodynamics of the drop jet in crossflow in these configurations considering trajectories, drop shape variations, spectrum of small disturbances, besides additional aspects relating to mesh quality are presented and discussed.

Numerical solution of differential equations having multitude of scales in the solution field is one of the most challenging research areas, but highly demanded in scientific and industrial applications. One of the natural approaches for handling such problems is to separate the scales and approximate the solution of the segregated scales with appropriate numerical method. Variational multiscale method (VMS) is a predominant method in the paradigm of scale separation schemes. In our work we have used the VMS technique to develop a numerical scheme for computations of turbulent flows in time-dependent domains. VMS allows separation of the entire range of scales in the flow field into two or three groups, thereby enabling a different numerical treatment for the different groups. In the context of computational fluid dynamics(CFD), VMS is a significant new improvement over the classical large eddy simulation (LES). VMS does away with the commutation errors arising due to filtering in LES. Further, in a three-scale VMS approach the model for the subgrid scale can be contained to only a part of the resolved scales instead of effecting the entire range of resolved scales. The projection based VMS scheme that we have developed gives a robust and efficient method for solving problems of turbulent fluid flows in deforming domains, governed by incompressible Navier {Stokes equations. In addition to the existing challenges due to turbulence, the computational complexity of the problem increases further when the considered domain is time-dependent. In this work, we have used an arbitrary Lagrangian-Eulerian (ALE) based VMS scheme to account for the domain deformation. In the proposed scheme, the large scales are represented by an additional tensor valued space. The resolved large and small scales are computed in a single unified equation, and the effect of unresolved scales is confined only to the resolved small scales, by using a projection operator. The popular Smagorinsky eddy viscosity model is used to approximate the effects of unresolved scales. The used ALE approach consists of an elastic mesh update technique. Moreover, a computationally efficient scheme is obtained by the choice of orthogonal finite element basis function for the resolved large scales, which allows to reformulate the ALE-VMS system matrix into the standard form of the NSE system matrix. Thus, any existing Navier{Stokes solver can be utilized for this scheme, with modifications. Further, the stability and error estimates of the scheme using a linear model of the NSE are also derived. Finally, the proposed scheme has been validated by a number of numerical examples over a wide range of problems.

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.

A interação entre fluidos e estruturas caracteriza um problema multi-físico não linear e está presente numa grande variedade de áreas da engenharia. Este trabalho apresenta o desenvolvi mento de ferramentas computacionais com base no Método dos Elementos Finitos (MEF) para a análise de interação fluido-estrutura (IFE) considerando escoamentos com baixas velocidades. Dada a interdisciplinaridade do tema, se faz necessário o estudo em três diferentes assuntos: a dinâmica das estruturas computacional, a dinâmica dos fluidos computacional, e o problema de acoplamento. No caso da dinâmica das estruturas empregar-se um elemento finito que seja adequado para a simulação de problemas de IFE, que claramente demandam uma análise não linear geométrica, optando-se pelo emprego de uma formulação descrita em posições, a qual evita problemas relativos à aproximação de rotações finitas. Quanto à dinâmica dos fluidos computacional, é empregado um método estável e ao mesmo tempo sensível à movimentação da estrutura, utilizando a descrição Lagrangeana-Euleriana Arbitrária (ALE). Os casos considerados neste trabalho, assim como muitos dos problemas de engenharia, ocorrem com escoamentos em baixas velocidades, implicando na incompressibilidade do fluido, o que demanda, para um método estável, a utilização de elementos que atendam à condição de Ladyzhenskaya-Babuska-Brezzi (LBB). Além disso, é necessário também o emprego de métodos que consigam neutralizar as variações espúrias decorrentes da não-linearidade de possíveis escoamentos com convecção dominante e que surgem com a aplicação do processo clássico de Galerkin. Para superar esse problema, é aplicado o método Streamline-Upwind/Petrov-Galerkin (SUPG), que adiciona difusividade artificial na direção do escoamento, controlando a amplitude dos termos convectivos. No que se refere ao acoplamento fluido-casca, buscam-se modularidade e versatilidade adotando-se o modelo particionado. O modelo de acoplamento implementado garante ainda a utilização de malhas do fluido e da estrutura sem a necessidade de coincidência de nós. / Interaction between fluids and structures characterizes a nonlinear multi-physics problem presente in a wide range of engineering fields. This works presets the development of computational tools based on finite element method (FEM) for fluid-structure interaction (FSI) analysis considering low speed flows (incompressible), as a great part of the engineering problems. Given the topic multidisciplinary nature, it is necessary to study three different subjects: the computational structural dynamics, the computational fluid mechanics and the coupling problem. Regarding structural mechanics, we seek to employ a finite element adequate to FSI simulation, what clearly demands a geometric nonlinear analysis. We chose to employ shell elements with formulation in terms of positions, which avoids problems related to finite rotations approximations. Concerning computational fluid dynamics, we employ a stable method, at same time sensible o structural movements, which is written in the arbitrary Lagrangian-Eulerian (ALE) description. The flow incompressibility demands, for a stable method, the use of elements according to the Ladyzhenskaya-Bbuska-Brezzi (LBB) condition. It is also necessary to employ methods able to neutralize the spurious variations that appears from convection dominated flows when applying the standard Galerking method. In order to overcome this problem, we apply the Streamline-Upwind/Petrov-Galerkin (SUPG) method, which adds artificial diffusivity to the streamline direction, controlling spurious variations. Considering the fluid-shell coupling, we seek modularity and versatility, adopting the partitioned model. The developed coupling model ensure the use of fluid and structure meshes with no need for matching nodes.

Um método lagrangeano-euleriano arbitrário para a resolução de escoamentos dominados por tensão superficial é apresentado neste trabalho. Tais escoamentos são importantes em muitas aplicações, especialmente em canais capilares que frequentemente aparecem em escoamentos em microescala. A resolução deste tipo de escoamento apresenta vários desafios que são abordados neste trabalho. O escoamento é resolvido somente para a fase líquida, com condições de contorno apropriadas para a superfície livre que delimita o líquido e o gás, que é representada por arestas e vértices da malha computacional. Esta se move e se deforma, sendo que sua qualidade é mantida sob controle para não degradar a solução numérica. As equações de Navier-Stokes são discretizadas pelo método de elementos finitos em um referencial arbitrário. O método de incorporação dos efeitos de tensão superficial e linha de contato é explicado em detalhes. Validações comprovam a precisão do método proposto, com comparações através de soluções pseudo-analíticas para casos simples. Finalmente alguns resultados sobre escoamentos em capilares são apresentados / An arbitrary lagrangian-eulerian finite element method to solve surface tension dominated flows is presented. Such flows are important in many applications, particularly in capillary channels, that appear in microscale flows. The resolution of such flows presents several challenges that are addressed in this work. The flow is solved only in the liquid phase, and proper boundary conditions are applied on the free-surface, bounding the liquid and gas, which is explicitly represented by vertices and edges of the computational mesh. The mesh is moved and deformed, but its quality is kept under control in order to control errors in the numerical solution. The Navier-Stokes equations are discretized by standard Galerkin finite element method in an arbitrary reference. Details of the computation of surface tension and contact line effects are presented. The methodology is validated for a number of simple test cases against known pseudo-analytical solutions, and numerical results are presented, showing the robustness and accuracy of the methodology. Finally, some results about surface-tension-driven flows in capillaries are presented

Submitted by Geandra Rodrigues (geandrar@gmail.com) on 2018-04-18T13:51:58Z No. of bitstreams: 1 marcusviniciusdesouzaferraz.pdf: 4103901 bytes, checksum: e4adcd64380c6ba8941b29bcc9d0abfd (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2018-04-19T17:48:41Z (GMT) No. of bitstreams: 1 marcusviniciusdesouzaferraz.pdf: 4103901 bytes, checksum: e4adcd64380c6ba8941b29bcc9d0abfd (MD5) / Made available in DSpace on 2018-04-19T17:48:41Z (GMT). No. of bitstreams: 1 marcusviniciusdesouzaferraz.pdf: 4103901 bytes, checksum: e4adcd64380c6ba8941b29bcc9d0abfd (MD5) Previous issue date: 2018-02-27 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O conhecimento da topografia das superfícies e uma compreensão da interação entre elas é essencial para qualquer estudo que envolva os fenômenos de atrito, desgaste e lubrificação. O estudo da relação entre o atrito e os parâmetros de rugosidade é um problema difícil e de interesse tanto industrial como acadêmico e trabalhos experimentais e teóricos têm mostrado que uma película de fluido entre duas superfícies rugosas em movimento relativo impede o contato sólido - sólido e pode proporcionar atrito muito baixo e desgaste desprezível. A modelagem matemática utilizada neste trabalho é baseada em modelos clássicos, tais como a equação de Reynolds para a descrição dos fenômenos hidrodinâmicos e as formulações de Hertz (1896) e Greenwood e Williamson (1966) para a modelagem do contato das asperezas entre as superfícies rugosas. Para tratar a complexidade das interações entre o fluido e os pares sólidos contactados, a descrição Lagrangiana-Euleriana Arbitrária é apresentada nesta pesquisa. Através do Método dos Elementos Finitos um modelo tridimensional é gerado no Abaqus ®, a fim de identificar as pressões de contato, as tensões tangenciais e normais resultantes e os coeficientes de atrito decorrrentes do deslizamento entre uma superfície texturizada e lubrificada e um plano rígido (em analogia aos modelos de contato clássicos), cujos perfis de rugosidade são construídos a partir de informações da rugosidade média quadrática de superfícies dentárias. São avaliados também a sensibilidade de alguns parâmetros do lubrificante na determinação do coeficiente de atrito e são propostos modelos com condições de contorno distintas. Entretanto, para a verificação destes últimos busca-se reproduzir qualitativamente o resultado encontrado por Lorentz (2013) na investigação numérica de sistemas tribológicos no regime misto de lubrificação. A metodologia aqui proposta emerge como uma alternativa eficaz no campo da Tribologia, na predição do coeficiente de atrito e outras variáveis pertinentes a um fenômeno ainda pouco compreendido. Realiza-se uma análise de sensibilidade dos parâmetros de modelagem, a fim de identificar como os mesmos afetam consideravelmente o comportamento mecânico na interface de contato. / The knowledge of the topography of surfaces and an understanding of the interaction between them is essential for any study involving the phenomena of friction, wear and lubrication. The study of the relationship between friction and roughness parameters is a difficult problem of both industrial and academic interest and experimental and theoretical works have shown that a fluid film between two rough surfaces in relative motion prevents solid - solid contact and can provide very low friction and negligible wear. The mathematical modeling used in this paper is based on classical models, such as the Reynolds equation for the description of the hydrodynamic phenomena and the formulations of Hertz (1896) and Greenwood e Williamson (1966) of the contact between the asperities of rough surfaces. To address the complexity of the interactions between the fluid and the contacted solid pairs, the Lagrangian-Eulerian Arbitrary description is presented in this research. Through the Finite Element Method, a three-dimensional model is generated in Abaqus ®R to identify contact pressures, resulting tangential and normal stresses, and friction coefficients resulting from sliding between a textured and lubricated surface and a rigid plane (in analogy to classic contact models), whose roughness profiles are constructed from information on the quadratic roughness of dental surfaces. The sensitivity of some lubricant parameters in the determination of the coefficient of friction is also evaluated and models with different boundary conditions are proposed. However, for the vefrification of the latter, it is sought to qualitatively reproduce the result found by Lorentz (2013) in the numerical investigation of tribological systems without mixed lubrication regime. A methodology proposed here emerges as an effective alternative in the field of Tribology, in the prediction of the coefficient of friction and other relevant variables to a phenomenon still little understood. A sensitivity analysis of the modeling parameters is performed, in order to identify how they considerably affect the mechanical behavior at the contact interface.