Ce travail de thèse est consacré à la modélisation et simulation numérique d'écoulements diphasiques liquide-gaz mettant en jeu des transferts de chaleur. La simulation de configurations où la prise en compte des effets de compressibilité de la phase gazeuse est indispensable (micropompes, microactionneurs, etc...) a nécessité l'utilisation d'un modèle original, considérant le liquide incompressible et le gaz compressible sous l'hypothèse faible Mach. Lors de cette thèse, ce modèle a été implémenté dans un code diphasique prenant en compte l'interface à l'aide d'une méthode de front-tracking. Des cas tests ont été développés spécifiquement afin de vérifier la conservation de l'énergie pour des configurations de complexité croissante. Les résultats des cas tests ont permis de mettre en évidence la difficulté à assurer la conservation de l'énergie lorsque l'interface n'est pas discontinue mais lissée, comme c'est le cas dans la méthode de front-tracking standard. Une méthode de traitement d'interface hybride a été proposée, rétablissant le caractère discontinu de l'interface avec la reconstruction d'une fonction indicatrice de phase échelon, tandis que le déplacement de l'interface est assuré d'un pas de temps à l'autre à l'aide du front-tracking. Les résultats obtenus avec cette nouvelle méthode hybride sont très satisfaisants, la méthode hybride permettant d'assurer la conservation de l'énergie et de la masse avec précision dans les simulations.

Les travaux réalisés au cours de cette thèse visent à apporter une aide à la compréhension des mécanismes régissant l'évaporation des films liquides. La physique des problèmes à changement de phase avec ligne de contact est gouvernée par des mécanismes agissant au niveau de ces dernières ainsi qu'au niveau des interfaces gaz-liquide. La description géométrique des interfaces est donc un point clef dans cette étude tant d'un point de vue numérique qu'expérimental. Un outil numérique tridimensionnel de suivi d'interface infiniment fine a été créé afin d'anticiper la simulation numérique de ces problèmes. Cet algorithme repose sur une description quadratique de l'interface qui est une surface du maillage. L’ordre de convergence de cette méthode relatif à différentes caractéristiques (courbure, normale, position) de l'interface ainsi que sa qualité conservative ont été étudiés.Parallèlement à cette approche numérique, une technique expérimentale d'inversion optique a été développée. Celle-ci permet la mesure dynamique du profil d'une interface au voisinage de la ligne de contact dans le cadre de fluide très mouillant. Cette méthode, appliquée dans le cadre d'un problème de film en évaporation quasi axisymétrique a permis, dans le cadre d'une expérience simple, une première quantification des flux de chaleur et de masse échangés au niveau de la ligne de contact. / The tasks fulfilled during this study aim at providing tools for understanding the mechanisms of liquid film evaporation. Evaporation in the presence of a triple lines is ruled by physical mechanisms acting on these singularities as well as on the gas liquid interfaces. Thus, the interface description, under it numerical as experimental aspects, is a key point of this study. A numerical tridimentionnal sharp interface tracking tool has been developed to be used in further numerical simulation of the evaporation problem. In this algorithm, the interface is described as a quadratic grid surface. Its convergence orders relatively to geometrical parameters of the interface (curvature, normal, position) has been studied.In the mean time, an experimental optical inversion technique has been developed and validated. It allows a dynamic measurement of the interface shape in the vicinity of the triple line of highly wetting fluid. This method has been applied to the case of a quasi axisymetric evaporating film and provides a first quantification of the heat and mass transfers generated by the triple line region.

O objetivo desta dissertação de mestrado é estudar técnicas numéricas para simular escoamentos incompressíveis multifásicos e implementar uma ferramenta computacional utilizando malhas hierárquicas e discretizações por diferenças finitas. São apresentados a formulação matemática e o desenvolvimento do método numérico, levando em consideração o caráter multifásico do escoamento. Foi adotado o modelo de força superficial contínua e a representação da interface foi feita pelo método de acompanhamento de fronteira. São expostos todos os testes realizados durante o desenvolvimento da ferramenta para checar cada etapa do método. Finalmente, testes visando verificar o código foram feitos e os resultados obtidos foram considerados satisfatórios para a verificação da ferramenta aqui desenvolvida. / The objective of this masters degree essay is to study numerical techniques to simulate incompressible multiphase flows and to implement a computational tool using hierachical meshes and discretizations by finite diferences. We introduce the mathematical formulation and the development of the numerical method, for the multiphase flow problem. A continuum surface force model is employed with the interface representation by the front tracking method. We show all tests performed to verify each stage of the methods development. Finally, results obtained in classical benchmark flow tests show good agreement with previous published results, corroborating the validity of this newly developed numerical tool.

Métodos de representação de superfícies a partir de pontos não-organizados se mantêm como uma das principais vertentes científicas que aquecem o estado-da-arte em Computação Gráfica e, significativamente, estão sendo reconhecidos como uma ferramenta interessante para definição de interfaces móveis no contexto de simulações numéricas de escoamento de fluidos. Não é difícil encontrar motivos para tais fatos: pelo lado da computação gráfica, por exemplo, a manipulação de conjuntos de pontos massivos com geometrias complexas e sujeitos a informações ruidosas ainda abre margem para novas metodologias. Já no âmbito da mecânica dos fluidos, onde os dados não são originados de \\emph tridimensionais, mas sim de interfaces entre fluidos imiscíveis, mecanismos de representação de superfícies a partir de pontos não-organizados podem apresentar características computacionais e propriedades geométricas que os tornem atrativos para aplicações em simulação de fenômenos físicos. O objetivo principal dessa tese de doutorado foi, portanto, o desenvolvimento de técnicas de representação de superfícies a partir de pontos não-organizados, que sejam capazes de suprir restrições de importantes trabalhos prévios. Nesse sentido, primeiramente focalizamos a elaboração de técnicas baseadas em formulações de mínimos-quadrados-móveis e de uma técnica robusta de partição da unidade implícita adaptativa em duas vias. Além de mecanismos de representação de superfícies a partir de pontos não-organizados, também propusemos um método promissor para representação de interfaces em simulação numérica de escoamento de fluidos multifásicos. Para isso, embasamo-nos numa abordagem Lagrangeana (livre-de-malhas), fundamentada no método dos mínimos-quadrados-móveis algébricos e apresentamos diversos resultados numéricos, estudos de convergências e comparações que evidenciam o potencial dessa metodologia para simulações numéricas de fenômenos físicos. Apesar de a contribuição principal deste trabalho ser o desenvolvimento de métodos para representação de superfícies a partir de pontos não-organizados, a experiência que adquirimos no desenvolvimento dessas técnicas nos conduziu à elaboração de mecanismos para representação de dados volumétricos não-organizados. Por conta disso, apresentamos dois mecanismos de representação a partir de dados volumétricos não-organizados com o intuito de serem aplicáveis a informações oriundas de malhas contendo células arbitrárias, isto é, propusemos a definição de um método de rendering unificado / Surface reconstruction from unorganized points has been one of the most promising scientific research areas in Computer Graphics. In addition, it has been used successfully for the definition of fluid interface in numerical simulation of fluid flow. There are several reasons to that fact: for instance, considering Computer Graphics, we have the handling of out-of-core data from complicated geometries and subject to noisy information that brings out opportunities for the development of new techniques. Further, considering Numerical Fluid Mechanics, where the input data does not come from tridimensional scanners, but from fluid interfaces, schemes that define the surface from unorganized points can offer geometrical and computational properties useful to numerical fluid flow simulation. The main goal of this project was the development of novel techniques for reconstructing surfaces from unorganized points with the capability to overcome the main drawbacks of important previous work. To that end, first we focused on the development of techniques based on moving-least-squares and on a robust twofold partition of unity Implicits. Added to the development of surface reconstruction from unorganized points, we proposed a novel scheme for defining fluid flow interfaces. We approach a meshless Lagrangian based on algebraic moving-least-squares surfaces. In addition, we presented several numerical results, convergence tests and comparisons, which state the power of the method to numerical simulation of physical phenomena. Although our main contributions were focused on surface reconstruction from points, we proposed methods to function reconstruction from unorganized volumetric data. Thus, we present two schemes to represent volumetric data from arbitrary meshes, i.e., a unified rendering scheme

O estudo do comportamento de fluidos é um antigo domínio das ciências da natureza. Ultimamente, fenômenos de engenharia que eram estudados empiricamente passaram a ser estudados com auxílio computacional. A Dinâmica de Fluidos Computacional (DFC) é a área da ciência da computação que estuda métodos computacionais para simulação de escoamento de fluidos, e muitas vezes é a forma mais prática, ou a única, de se observar fenômenos de interesse no escoamento. Este projeto de Mestrado procurou investigar, no âmbito da simulação de um escoamento bifásico, métodos computacionais para representar a interface entre dois fluidos imiscíveis. A separação dos fluidos por meio de uma interface é necessária para assegurar que, propriedades como viscosidade e densidade, específicas de cada fluido, sejam utilizadas corretamente para o cálculo do movimento de seus respectivos fluidos. Desenvolvemos um método lagrangeano sem a utilização de malhas com o objetivo de suprir algumas restrições de trabalhos prévios. Para representar a interface entre os dois fluidos, este método utiliza uma técnica de reconstrução de superfícies baseada em aproximações de superfícies algébricas de alta ordem. Os resultados numéricos reportados neste documento evidenciam o potencial da nossa abordagem / The study of the behaviour of fluids is an ancient field in natural sciences. Recently, engineering phenomena that were empirically studied started to be done with computacional aid. The Computational Fluid Dynamics (CFD) is the area of science that studies computational methods for computer simulation of fluid flow, and often is the most practical way, or the only, to observe phenomena of interest in flow. This Masters degree project sought to investigate, in the context of the simulation of biphasic flows, computational methods to represent the interface between two immiscible fluids. The separation of fluids by the means of an interface is required to ensure that, during the simulation, the physical properties of a fluid, like density and viscosity (specific of each fluid) are properly used in the calculus of the respective fluid motion. We developed a lagrangean method without the use of mesh with the goal of alleviating some of the previous works restrictions. To represent the interface between the two fluids, this method uses a surface reconstruction technique based on approximations of high order algebraic surfaces. The numerical results reported herein show the potential of our approach

<p>U doktorskoj disertaciji posmatrani su Rimanovi problemi kod strogo i slabo hiperboličnih&nbsp;nelinearnih sistema PDJ. U uvodu je predstavljena jednačina zakona održanja u jednoj prostornoj&nbsp;dimenziji i definisani su Ko&scaron;ijevi i Rimanovi problemi. U drugoj glavi, date su osnovne osobine&nbsp;nelinearnih hiperboličnih zakona održanja, uvedeni supojmovi stroge hiperboličnosti i slabog re&scaron;enja&nbsp;zakona održanja. Definisani su Rankin-Igono i entropijski uslovi kao i op&scaron;te re&scaron;enje Rimanovog problema&nbsp;(za dovoljno male početne uslove). U trećoj glavi detaljno je obja&scaron;njena Glimova diferencna &nbsp;&scaron;ema. Metod&nbsp;praćenja talasa predstavljen je u četvrtoj glavi. Pokazano je da se ovom metodom, za dovoljno male&nbsp;početne uslove, dobija stabilno i jedinstveno re&scaron;enje koje u svakom vremenu ima ograničenu totalnu&nbsp;varijaciju. U petoj glavi, posmatrana je jednačina protoka izentropnog gasa u Lagranžovim koordinatama.&nbsp;Uz pretpostavku da je početni uslov ograničen i da ima ograničenu totalnu varijaciju, pokazano je da&nbsp;Ko&scaron;ijev problem ima jedinstveno slabo re&scaron;enje ako je totalna varijacija početnog uslova pomnožena sa &nbsp;0&lt;&epsilon;&lt;&lt; 1 dovoljno mala. Slabo re&scaron;enjedobijeno je metodom praćenja talasa. U glavi &scaron;est ispitana je&nbsp;interakcija dva delta talasa koji su posmatrani kao specijalna vrsta shadowtalasa. U glavi sedam,&nbsp;pokazano je da za proizvoljno velike početne uslove, re&scaron;enje Rimanovog problema jednodimenzionalnog&nbsp;Ojlerovog zakona održanja gasne dinamikepostoji, daje jedinstveno i entropijski dopustivo, uz drugačiju<br />ocenu snaga elementarnih talasa. Data je numerička verifikacija interakcije dva delta talasa kori&scaron;ćenjem&nbsp;metode praćenja talasa.</p> / <p>In this doctoral thesis, Riemann problems for strictly and weakly nonlinear hyperbolic PDE&nbsp;systems were observed. In the introduction, conservation laws in one spatial dimension were presented&nbsp;and the Cauchy and Riemann problems were defined. In the second chapter, the basic properties of&nbsp;nonlinear hyperbolic conservation laws were intorduced, as well as the terms such as strictly hyperbolic&nbsp;system and weak solution of conservation law. Also, Rankine -Hugoniot and entropy conditions were<br />introduced and the general solution to the Riemann problem (for sufficiently small initial conditions) were&nbsp;defined. Glimm&rsquo;s difference scheme was explained in the third chapter. The wave front tracking method&nbsp;was introduced in the fourth chapter. It was shown that, using this method, for sufficiently small initial&nbsp;conditions, it could be obtained a unique solution with bounded total variation for t &ge;0. In the fifth&nbsp;chapter, the Euler equations for isentropic fluid inLagrangian coordinates were observed. Under the&nbsp;assumption that the initial condition was bounded and had bounded total variation, it was shown that the&nbsp;Cauchy problem had a weak unique solution, provided that the total variation of initial condition&nbsp;multiplied by 0&lt;&epsilon;&lt;&lt;1 was sufficiently &nbsp;small. Weak solution was obtained by applying the wave front&nbsp;tracking method. In the sixth chapter, the interaction of two delta shock waves were examined. Delta&nbsp;shock waves were regarded as special kind of shadowwaves. In the chapter seven, it was shown that for&nbsp;arbitrarily large initial conditions, solution to the Riemann problem of one-dimensional Euler&nbsp;conservation laws of gas dynamics existed, it was unique and admissible. New bounds on the strength of&nbsp;elementary waves in the wave front tracking algorithm were given. The numerical verification of two&nbsp;delta shock waves interaction via wave front tracking method was given at the end of the thesis.</p>

No description available.

La simulation de cas réalistes d'écoulements ou de transferts thermiques implique souvent l'utilisation d'obstacles ou d'interfaces de forme complexe. De part leur manque de flexibilité, les maillages structurés ne sont pas initialement adaptés au traitement d'interfaces irrégulières, ces dernières coïncidant rarement avec les lignes du maillage. Afin de permettre à l'approche structurée de traiter des interfaces complexes avec précision, des méthodes dites de domaines fictifs sont nécessaires. La première contribution de cette thèse est une nouvelle méthode de travail sur maillage curviligne structuré qui permet de réutiliser de nombreuses méthodes fonctionnant initialement sur des maillages cartésiens sur maillages curvilignes. Nous avons ensuite mis au point deux nouvelles méthodes de domaines fictifs : la méthode de pénalisation de sous-maille (PSM) pour la gestion des frontières immergées pour les équations elliptiques et de Navier-Stokes et la méthode d'interface immergée algébrique (IIA) pour les problèmes d'interfaces immergées pour les équations elliptiques. L'un des intérêts de ces deux méthodes à l'ordre deux en espace est leur simplicité. Ces différents développements ont finalement été appliqués à des cas de couplage fluide-structure académiques et réalistes (sédimentation d'un cylindre, hydroplanage d'un pneu, écoulements dans une tête de forage et convection naturelle dans la grotte de Lascaux).

Dans cette thèse, on considère plusieurs problèmes issus de la modélisation macroscopique des mouvements de foule. Le premier modèle consiste en une loi de conservation avec un flux discontinu, le second est un système mixte hyperbolique-elliptique et le dernier est une équation non-locale. D'abord, on utilise le modèle de Hughes une dimension pour décrire l'évacuation d'un couloir avec deux sorties. Ce modèle couple une loi de conservation avec un flux discontinu à une équation eikonale. On implémente la méthode de suivi de fronts, qui traite explicitement le comportement de la solution non-classique au point de rebroussement, afin d'obtenir des solutions de référence. Elles serviront à tester numériquement la convergence de schémas aux volumes finis classiques. Ensuite, on modélise le croisement de deux groupes marchant dans des directions opposées avec un système de lois de conservation mixte hyperbolique-elliptique dont le flux dépend des deux densités. Le système perd son hyperbolicité pour certainement valeurs de densité. On assiste à l'apparition d'oscillations persistantes mais bornées, ce qui conduit à la reformulation du problème associé dans le cadre des mesures de probabilités. Finalement, on étudie un modèle non-local de trafic piétonnier en deux dimensions. Le modèle consiste en une loi de conservation dont le flux dépend d'une convolution de la densité. Avec ce modèle, on résout un problème d'optimisation pour une évacuation d'une salle avec une méthode de descente, évaluant l'impact du calcul explicite du gradient de la fonction coût avec la méthode de l'état adjoint plutôt que son approximation par différences finies. / In this thesis, we consider nonclassical problems brought out by the macroscopic modeling of pedestrian flow. The first model consists of a conservation law with a discontinuous flux, the second is a mixed hyperbolic-elliptic system of conservation laws and the last one is a nonlocal equation. In the first chapter, we use the Hughes model in one space-dimension to represent the evacuation of a corridor with two exits. The model couples a conservation law with discontinuous flux to an eikonal equation. We implement the wave front tracking scheme, treating explicitly the solution nonclassical behavior at the turning point, to provide a reference solution, which is used to numerically test the convergence of classical finite volume schemes. In the second chapter, we model the crossing of two groups of pedestrians walking in opposite directions with a system of conservation laws whose flux depends on the two densities. This system loses its hyperbolicity for certain density values. We assist to the rising of persistent but bounded oscillations, that lead us to the recast of the problem in the framework of measure-valued solutions. Finally we study a nonlocal model of pedestrian flow in two space-dimensions. The model consists of a conservation law whose flux depends on a convolution of the density. With this model, we solve an optimization problem for a room evacuation with a descent method, evaluating the impact of the explicit computation of the cost function gradient with the adjoint state method rather than approximating it with finite differences.