O fenômeno de injeção de vapor em tanques de alívio é de relevância no projeto de reatores de água leve, sejam eles do tipo reator de água pressurizada (PWR) ou reator de água fervente (BWR). Este fenômeno permite a rápida absorção do vapor injetado em massa de água, por meio de sua condensação, uma vez que este vapor pode conter contaminantes químicos ou radiológicos que não permitem o seu descarte diretamente no ambiente. Desta forma, facilita-se a coleta do vapor produzido por descarga de vapor da água do resfriamento do reator, radiologicamente contaminada, e evita-se o que projeto de dispositivos e equipamentos necessite considerar a elevada pressão do vapor. A rapidez com que se dá a condensação é fruto de processos físicos que ocorrem na interface de vapor e água e que ainda não possuem modelo analítico e numérico definido. Em 1972 um modelo semi-empírico foi proposto, o qual, desde então, vem evoluindo. Não obstante, até o presente momento, não há modelo definitivo que se proponha a abranger toda extensão das condições experimentais. Estes modelos são fortemente dependentes do fluxo de massa que atravessa a interface de vapor e água, entretanto, até a presente data, não há expressão que determine este fluxo de massa, de tal forma que o valor de 275 Kg/m2/s vem sendo assumido como \"representativo da ordem de grandeza do fenômeno\" até o presente momento. Neste trabalho, é proposto um método de cálculo analítico do fluxo de massa, considerando-se como premissa a isentropia da injeção, e o desenvolvimento da 1ª e 2ª leis da Termodinâmica. Ainda, o fenômeno é analisado experimentalmente, por meio da análise dos dados produzidos no experimento do Circuito Termo Hidráulico de 150 bar (Loop 150), realizado nas dependências do CENTRO TECNOLÓGICO DA MARINHA EM SÃO PAULO. Por fim, um modelo numérico em software comercial foi desenvolvido para complementar a análise. Os resultados obtidos comprovam que a formulação isentrópica do fluxo de massa corrige de maneira satisfatória o fluxo de massa constante utilizado até então nos modelos semi-empíricos. Tal comprovação se deu através de análise numérica e da confrontação com dados experimentais obtidos na literatura. / The phenomenon of vapor injection in relief tanks presents relevance in the design of light water reactors, be they of the type pressurized water reactor (PWR) or boiling water reactor (BWR). This phenomenon allows the rapid absorption of the vapor injected in a mass of water, by condensation. Since this vapor may contain chemical or radiological contaminants that do not allow its discharge directly in the environment, it must be collected. The condensation avoids the design of devices and equipment, which need to consider the high vapor pressure, and allows the vapor to be collected. The rapidity with which the condensation occurs is the result of physical processes that occur at the interface of steam and water. These processes do not yet have a defined numerical and analytical model. In 1972 a semi-empirical model was proposed, which has, since then, evolved. Nevertheless, up to the present moment, there is no definitive model that intends to cover every extension of the experimental conditions. These models are strongly dependent on the mass flow through the steam and water interface, however, up to date, there is no expression that determines this mass flow. For the sake of this, the value of 275 kg / m2 / s has been assumed as \"representative of the order of magnitude of the phenomenon\" up to the present moment. In this work, a method of analytical calculation of mass flow is proposed, considering as premise the isotropy of the injection, and the development of the 1st and 2nd laws of thermodynamics. Still, the phenomenon is analyzed experimentally, by means of the analysis of the data produced in the experiment of Loop 150, realized in dependencies of the CENTRO TECNOLÓGICO DA MARINHA EM SÃO PAULO. Finally, a numerical model in commercial software was developed to complement the analysis. The result is proven with an isentropic mass flow formulation, which satisfactorily corrected the mass used in the former semi-empirical models. Such verification was performed through a series of data and confrontation with experimental data in the literature.

Ce sujet concerne un travail pluridisciplinaire mathématique et astrophysique. Le but de cette thèse est l'étude des modèles d'hydrodynamique radiative dont l'application est bien évidemment très vaste en physique et astrophysique. Les modèles M1-multigroupes sont explorés pour décrire le transfert radiatif sans faire à priori d'hypothèse sur la profondeur optique du milieu. L'intérêt qui découle directement de ce travail est le développement du code d'hydrodynamique radiative HADES 2D permettant le calcul massivement parallèle. Il autorise des simulations dans des configurations astrophysiques réalistes en termes de nombre de Mach et de contraste de densité et de température entre les différents milieux. Nous nous sommes concentrés sur deux applications intéressantes : les jets d'étoiles jeunes et les chocs radiatifs dont les premières simulations seront présentées. / This topic is a multidisciplinary work between mathematics and astrophysics. The aim of this thesis is the study of radiation hydrodynamic models of which application is obviously very broad in physics and astrophysics. M1-multigroup models are explored to describe the radiative transfer without a priori assumption on the optical depth of the medium. The interest ensuing directly from this work is the development of a radiation hydrodynamic code, namely HADES 2D, for massively parallel computing. It allows simulations in realistic astrophysical configurations in terms of Mach number, density and temperature contrasts between different environments. We focused on two interesting applications: the jets from young stars and the radiative shocks of which first simulations will be presented.

Dans cette thèse on s’intéresse à la simulation numérique d’écoulements des matériaux compressibles, voir fluides et solides élastiques. Les matériaux considérés sont décrits avec un modèle monolithique eulérian, fermé avec une loi d’état hyperélastique qui considère les différents comportements des matériaux. On propose un nouveau schéma de relaxation qui résout les écoulements compressibles dans des différents régimes, avec des nombres de Mach très petits jusqu’à l’ordre 1. Le schéma a une formulation générale qui est la même pour tous le matériaux considérés, parce que il ne dépend pas directement de la loi d’état. Il se base sur une discrétisation complétement implicite, facile à implémenter grâce à la linéarité de l’opérateur de transport du système de relaxation. La discrétisation en espace est donnée par la combinaison de flux upwind et centrés, pour retrouver la correcte viscosité numérique dans les différents régimes. L’utilisation de mailles cartésiennes pour les cas 2D s’adapte bien à une parallélisation massive, qui permet de réduire drastiquement le temps de calcul. De plus, le schéma a été adapté pour la résolution sur des mailles quadtree, pour implémenter l’adaptativité de la maille avec des critères entropiques. La dernière partie de la thèse concerne la simulation numérique d’écoulements multi-matériaux. On a proposé une nouvelle méthode d’interface “sharp”, en dérivant les conditions d’équilibre en implicite. L’objectif est la résolution d’interfaces physiques dans des régimes faiblement compressibles et avec un nombre de Mach faible, donc les conditions multi-matériaux sont couplées au schéma implicite de relaxation. / In this thesis we are concerned with the numerical simulation of compressible materials flows, including gases, liquids and elastic solids. These materials are described by a monolithic Eulerian model of conservation laws, closed by an hyperelastic state law that includes the different behaviours of the considered materials. A novel implicit relaxation scheme to solve compressible flows at all speeds is proposed, with Mach numbers ranging from very small to the order of unity. The scheme is general and has the same formulation for all the considered materials, since a direct dependence on the state law is avoided via the relaxation. It is based on a fully implicit time discretization, easily implemented thanks to the linearity of the transport operator in the relaxation system. The spatial discretization is obtained by a combination of upwind and centered schemes in order to recover the correct numerical viscosity in different Mach regimes. The scheme is validated with one and two dimensional simulations of fluid flows and of deformations of compressible solids. We exploit the domain discretization through Cartesian grids, allowing for massively parallel computations (HPC) that drastically reduce the computational times on 2D test cases. Moreover, the scheme is adapted to the resolution on adaptive grids based on quadtrees, implementing adaptive mesh refinement techinques. The last part of the thesis is devoted to the numerical simulation of heterogeneous multi-material flows. A novel sharp interface method is proposed, with the derivation of implicit equilibrium conditions. The aim of the implicit framework is the solution of weakly compressible and low Mach flows, thus the proposed multi-material conditions are coupled with the implicit relaxation scheme that is solved in the bulk of the flow.

This work aims at the numerical simulation of compressible fluid/deformable structure interactions. In particular, we have developed a partitioned coupling algorithm between a Finite Volume method for the compressible fluid and a Discrete Element method capable of taking into account fractures in the solid. A survey of existing fictitious domain methods and partitioned algorithms has led to choose an Embedded Boundary method and an explicit coupling scheme. We first showed that the Discrete Element method used for the solid yielded the correct macroscopic behaviour and that the symplectic time-integration scheme ensured the preservation of energy. We then developed an explicit coupling algorithm between a compressible inviscid fluid and an undeformable solid. Mass, momentum and energy conservation and consistency properties were proved for the coupling scheme. The algorithm was then extended to the coupling with a deformable solid, in the form of a semi-implicit scheme. Finally, we applied this method to unsteady inviscid flows around moving structures: comparisons with existing numerical and experimental results demonstrate the excellent accuracy of our method

Submitted by Emanuel Varela Cardoso (emanuel.varela@ufcg.edu.br) on 2018-06-12T18:29:02Z No. of bitstreams: 1 EDNA GOMES DA SILVA – TESE (PPGEP) 2016.pdf: 5406028 bytes, checksum: f01a907e44b485aba5813ac32c84fd35 (MD5) / Made available in DSpace on 2018-06-12T18:29:02Z (GMT). No. of bitstreams: 1 EDNA GOMES DA SILVA – TESE (PPGEP) 2016.pdf: 5406028 bytes, checksum: f01a907e44b485aba5813ac32c84fd35 (MD5) Previous issue date: 2016-08-19 / Este trabalho tem como objetivo estudar teórico e experimentalmente a secagem combinada (convectiva e micro-ondas) de sólidos esferoidais prolatos. Um modelo matemático bidimensional transiente em coordenadas esferoidais prolatas foi desenvolvido para predizer a transferência de calor e massa e simular a distribuição do teor de umidade e da temperatura no interior do sólido, considerando condição de contorno convectiva na superfície, termo de geração de calor por onda eletromagnética e ao longo do processo, propriedades termo físicas e volume constantes. O método dos volumes finitos foi utilizado para a solução numérica das equações governantes, usando uma malha regular. Para validação da metodologia teórica foram realizados experimentos de secagem por micro-ondas com grãos de arroz em casca in natura (variedade BRSMG Conai) para três níveis de potências absorvidas: 0,192; 0,491 e 0,694 W. Resultados do efeito da razão de aspecto do sólido, números de Biot de transferência de massa e de calor, densidade de potência e fator de atenuação da onda eletromagnética sobre a cinética de secagem e aquecimento, e distribuições de temperatura e teor de umidade no interior do sólido são apresentados e analisados. A comparação entre os dados teóricos e experimentais do teor de umidade e da temperatura dos grãos de arroz em casca resultou num bom ajuste e possibilitou a estimativa dos coeficientes de difusão de massa e de transferência de massa e calor convectivos. Verificou-se que, o aumento na potência absorvida proporcionou menor tempo de secagem e maior temperatura superficial do grão de arroz, e que os coeficientes de transporte de calor e de massa preditos tendem a aumentar com o aumento da potência aplicada. / This work aims to study theoretical and experimentally the combined drying (convective and microwave) of prolate spheroidal solids. A transient two-dimensional mathematical model in prolate spheroidal coordinates it was developed to predict the heat and mass transfer and simulate the moisture content and temperature distributions inside the solid considering convective boundary condition at the surface, microwave heat generation term, thermo physical properties and volume constants. The finite volume method was used for the numerical solution of the governing equations, using a regular grid. To validate the theoretical methodology microwave drying experiments were carried out with rice grains in fresh bark (variety BRSMG Conai) in three levels of absorbed power: 0.192; 0.491 and 0.694 W. Results of the effect of the solid aspect ratio, Biot numbers of mass and heat transfer, power density and attenuation factor of the electromagnetic wave on the drying and heating kinetics, and temperature and moisture content distributions within the solid are presented and analyzed. The comparison between theoretical and experimental data of moisture content and temperature of the rice grain resulted in a good fit and estimative at the mass diffusion coefficient and convective heat and mass transfer coefficients. It was found that the increase in absorbed power provided shorter drying time and higher surface temperature of the rice grain , and that the predicted heat and mass transport coefficients tend to increase with increased power applied.

Les discrétisations adaptatives sont importantes dans les problèmes de fluxcompressible/incompressible puisqu'il est souvent nécessaire de résoudre desdétails sur plusieurs niveaux, en permettant de modéliser de grandes régionsd'espace en utilisant un nombre réduit de degrés de liberté (et en réduisant letemps de calcul).Il existe une grande variété de méthodes de discrétisation adaptative, maisles grilles cartésiennes sont les plus efficaces, grâce à leurs stencilsnumériques simples et précis et à leurs performances parallèles supérieures.Et telles performance et simplicité sont généralement obtenues en appliquant unschéma de différences finies pour la résolution des problèmes, mais cetteapproche de discrétisation ne présente pas, au contraire, un chemin faciled'adaptation.Dans un schéma de volumes finis, en revanche, nous pouvons incorporer différentstypes de maillages, plus appropriées aux raffinements adaptatifs, en augmentantla complexité sur les stencils et en obtenant une plus grande flexibilité.L'opérateur de Laplace est un élément essentiel des équations de Navier-Stokes,un modèle qui gouverne les écoulements de fluides, mais il se produit égalementdans des équations différentielles qui décrivent de nombreux autres phénomènesphysiques, tels que les potentiels électriques et gravitationnels. Il s'agitdonc d'un opérateur différentiel très important, et toutes les études qui ontété effectuées sur celui-ci, prouvent sa pertinence.Dans ce travail seront présentés des approches de différences finies et devolumes finis 2D pour résoudre l'opérateur laplacien, en appliquant des patchsde grilles superposées où un niveau plus fin est nécessaire, en laissant desmaillages plus grossiers dans le reste du domaine de calcul.Ces grilles superposées auront des formes quadrilatérales génériques.Plus précisément, les sujets abordés seront les suivants:1) introduction à la méthode des différences finies, méthode des volumes finis,partitionnement des domaines, approximation de la solution;2) récapitulatif des différents types de maillages pour représenter de façondiscrète la géométrie impliquée dans un problème, avec un focussur la structure de données octree, présentant PABLO et PABLitO. Le premier estune bibliothèque externe utilisée pour gérer la création de chaque grille,l'équilibrage de charge et les communications internes, tandis que la secondeest l'API Python de cette bibliothèque, écrite ad hoc pour le projet en cours;3) la présentation de l'algorithme utilisé pour communiquer les données entreles maillages (en ignorant chacune l'existence de l'autre) en utilisant lesintercommunicateurs MPI et la clarification de l'approche monolithique appliquéeà la construction finale de la matrice pour résoudre le système, en tenantcompte des blocs diagonaux, de restriction et de prolongement;4) la présentation de certains résultats; conclusions, références.Il est important de souligner que tout est fait sous Python comme framework deprogrammation, en utilisant Cython pour l'écriture de PABLitO, MPI4Py pour lescommunications entre grilles, PETSc4py pour les parties assemblage et résolutiondu système d'inconnues, NumPy pour les objets à mémoire continue.Le choix de ce langage de programmation a été fait car Python, facile àapprendre et à comprendre, est aujourd'hui un concurrent significatif pourl'informatique numérique et l'écosystème HPC, grâce à son style épuré, sespackages, ses compilateurs et pourquoi pas ses versions optimisées pour desarchitectures spécifiques. / Adaptive discretizations are important in compressible/incompressible flow problems since it is often necessary to resolve details on multiple levels,allowing large regions of space to be modeled using a reduced number of degrees of freedom (reducing the computational time).There are a wide variety of methods for adaptively discretizing space, but Cartesian grids have often outperformed them even at high resolutions due totheir simple and accurate numerical stencils and their superior parallel performances.Such performance and simplicity are in general obtained applying afinite-difference scheme for the resolution of the problems involved, but this discretization approach does not present, by contrast, an easy adapting path.In a finite-volume scheme, instead, we can incorporate different types of grids,more suitable for adaptive refinements, increasing the complexity on thestencils and getting a greater flexibility.The Laplace operator is an essential building block of the Navier-Stokes equations, a model that governs fluid flows, but it occurs also in differential equations that describe many other physical phenomena, such as electric and gravitational potentials, and quantum mechanics. So, it is a very importantdifferential operator, and all the studies carried out on it, prove itsrelevance.In this work will be presented 2D finite-difference and finite-volume approaches to solve the Laplacian operator, applying patches of overlapping grids where amore fined level is needed, leaving coarser meshes in the rest of the computational domain.These overlapping grids will have generic quadrilateral shapes.Specifically, the topics covered will be:1) introduction to the finite difference method, finite volume method, domainpartitioning, solution approximation;2) overview of different types of meshes to represent in a discrete way thegeometry involved in a problem, with a focuson the octree data structure, presenting PABLO and PABLitO. The first one is anexternal library used to manage each single grid’s creation, load balancing and internal communications, while the second one is the Python API ofthat library written ad hoc for the current project;3) presentation of the algorithm used to communicate data between meshes (beingall of them unaware of each other’s existence) using MPI inter-communicators and clarification of the monolithic approach applied building the finalmatrix for the system to solve, taking into account diagonal, restriction and prolongation blocks;4) presentation of some results; conclusions, references.It is important to underline that everything is done under Python as programmingframework, using Cython for the writing of PABLitO, MPI4Py for the communications between grids, PETSc4py for the assembling and resolution partsof the system of unknowns, NumPy for contiguous memory buffer objects.The choice of this programming language has been made because Python, easy to learn and understand, is today a significant contender for the numerical computing and HPC ecosystem, thanks to its clean style, its packages, its compilers and, why not, its specific architecture optimized versions.

Cette thèse aborde la résolution des problèmes d'interaction fluide-structure par un algorithme consistant en un couplage entre deux solveurs : un pour le fluide et un pour la structure. Pour assurer la cohérence entre les maillages fluide et structure, on considère également une discrétisation de chaque domaine par volumes finis. En raison des difficultés de décomposition du domaine en sous-domaines, nous considérons pour chaque environnement un algorithme parallèle de multi-splitting (ou multi-décomposition) qui correspond à une présentation unifiée des méthodes de sous-domaines avec ou sans recouvrement. Cette méthode combine plusieurs applications de points fixes contractantes et nous montrons que, sous des hypothèses appropriées, chaque application de points fixes est contractante dans des espaces de dimensions finies normés par des normes hilbertiennes et non-hilbertiennes. De plus, nous montrons qu'une telle étude est valable pour les résolutions parallèles synchrones et plus généralement asynchrones de grands systèmes linéaires apparaissant lors de la discrétisation des problèmes d'interaction fluide-structure et peut être étendue au cas où le déplacement de la structure est soumis à des contraintes. Par ailleurs, nous pouvons également considérer l’analyse de la convergence de ces méthodes de multi-splitting parallèles asynchrones par des techniques d’ordre partiel, lié au principe du maximum discret, aussi bien dans le cadre linéaire que dans celui obtenu lorsque les déplacements de la structure sont soumis à des contraintes. Nous réalisons des simulations parallèles pour divers cas test fluide-structure sur différents clusters, en considérant des communications bloquantes et non bloquantes. Dans ce dernier cas nous avons eu à résoudre une difficulté d'implémentation dans la mesure où une erreur irrécupérable survenait lors de l'exécution ; cette difficulté a été levée par introduction d’une méthode assurant la terminaison de toutes les communications non bloquantes avant la mise à jour du maillage. Les performances des simulations parallèles sont présentées et analysées. Enfin, nous appliquons la méthodologie présentée précédemment à divers contextes d'interaction fluide-structure de type industriel sur des maillages non structurés, ce qui constitue une difficulté supplémentaire. / This thesis deals with the solution of fluid-structure interaction problems by an algorithm consisting in the coupling between two solvers: one for the fluid and one for the structure. In order to ensure the consistency between fluid and structure meshes, we also consider a discretization of each domain by finite volumes. Due to the difficulties of decomposing the domain into sub-domains, we consider a parallel multi-splitting algorithm for each environment which represents a unified presentation of sub-domain methods with or without overlapping. This method combines several contracting fixed point mappings and we show that, under appropriate assumptions, each fixed point mapping is contracting in finite dimensional spaces normalized by Hilbertian and non-Hilbertian norms. In addition, we show that such a study is valid for synchronous parallel solutions and more generally asynchronous of large linear systems arising from the discretization of fluidstructure interaction problems and can be extended to cases where the displacement of the structure is subject to constraints. Moreover, we can also consider the analysis of the convergence of these asynchronous parallel multi-splitting methods by partial ordering techniques, linked to the discrete maximum principle, both in the linear frame and in the one obtained when the structure's displacements are subjected to constraints. We carry out parallel simulations for various fluidstructure test cases on different clusters considering blocking and non-blocking communications. In the latter case, we had to solve an implementation problem due to the fact that an unrecoverable error occurred during execution; this issue has been overcome by introducing a method to ensure the termination of all non-blocking communications prior to the mesh update. Performances of parallel simulations are presented ans analyzed. Finally, we apply the methodology presented above to various fluid-structure interaction cases on unstructured meshes, which represents an additional difficulty.

Submitted by Emanuel Varela Cardoso (emanuel.varela@ufcg.edu.br) on 2018-07-25T19:08:37Z No. of bitstreams: 1 HUGO CARVALHO AMORIM – DISSERTAÇÃO (PPGEM) 2016.pdf: 4230356 bytes, checksum: 63473c25aa9b0817eed244e6b1fc8738 (MD5) / Made available in DSpace on 2018-07-25T19:08:37Z (GMT). No. of bitstreams: 1 HUGO CARVALHO AMORIM – DISSERTAÇÃO (PPGEM) 2016.pdf: 4230356 bytes, checksum: 63473c25aa9b0817eed244e6b1fc8738 (MD5) Previous issue date: 2016-07-29 / Frutas são alimentos compostos normalmente por 80% a 90% de água, tornando-as extremamente perecíveis. Dessa forma, métodos de conservação devem ser aplicados para prolongar sua vida útil e conservar características de aceitação pelo consumidor. Assim, este trabalho teve como objetivo estudar a transferência de calor durante os processos de resfriamento e congelamento de banana utilizando geometria esferoidal prolata. Para a modelagem matemática, utilizou-se a equação da conservação de energia escrita em coordenadas esferoidais prolatas. A solução numérica da equação governante foi realizada através do método dos volumes finitos com uma formulação totalmente implícita. Um programa computacional foi escrito para simular o processo, através de linguagem computacional do software Mathematica, afim de se obter o histórico de temperatura no centro da banana e o perfil de temperatura, sendo os dados comparados a valores experimentais de congelamento de banana prata com casca, visando obter o coeficiente de transferência de calor convectivo e condutividade térmica da banana. Avaliou-se a influência das dimensões da fruta na cinética de resfriamento e nas frentes de congelamento. Concluiu-se que para a curva de resfriamento da banana, a modelagem pôde prever com bastante exatidão o período de resfriamento mas não o de pós-congelamento. Verificou-se que quanto menor o tamanho do produto, mais rápido é o processo e quanto maior a razão de aspecto de forma, maiores serão os gradientes de temperatura na ponta do produto. As taxas de resfriamento sofreram maior influência quando variou-se as dimensões mantendo um aspecto de forma constante do que variando as dimensões e mudando também o aspecto de forma. Por fim, verificou-se que quanto maior o aspecto de forma, mais não uniforme se dá o processo de transferência de calor na banana. / Fruits are a type of food that are typically made by 80% to 90% water, making them extremely perishable. Thus, conservation methods should be applied to prolong their shelf life and consumer acceptance characteristics. This work aims to study the heat transfer during the cooling and freezing processes of fruits with a prolate spheroid shape. The mathematical modeling consisted of energy conservation equation written in spheroidal prolates coordinates. The numerical solution of the governing equation was performed using the finite volume method with a fully implicit formulation. The process was simulated by computer-based language within the Mathematica software, to obtain the temperature history at the center of the fruit and the temperature profile, the data being compared with experimental values of freezed banana prata, to obtain the transfer coefficient convective-heat and thermal conductivity of the fruit. We evaluated the influence of fruit dimensions of the cooling kinetics and the freezing fronts. It was found that the lower the product size, the faster the process, and the higher the aspect ratio of shape, higher will be the temperature gradient at the top of the product.

Le comportement des matériaux multiphasiques couvre une multitude de phénomènes qui suscitent un grand intérêt dans le domaine scientifique et professionnel. Les propriétés mécaniques de ces types de matériau trouvent leur origine dans les phases dont ils sont composés, leur distribution et interaction. Un nouveau modèle hydrodynamique couplé est présenté dans ce travail de thèse, à appliquer à l'analyse de l'hydrodynamique des milieux granulaires saturés. Le modèle associe la méthode des éléments discrets (DEM) pour la modélisation de la phase solide, avec une formulation en volumes finis, à l'échelle des pores (PFV), du problème de l'écoulement. Une importance particulière est donné à la description de l'interaction entre les phases, avec la détermination des forces fluides à appliquer sur chacune des particule, tout en assurant un coût de calcul abordable, qui permet la modélisation de plusieurs milliers des particules en trois dimensions. Le milieux est considéré saturé par un fluide incompressible. Les pores et leur connectivité est basée sur une triangulation régulière des assemblages. L'analogie de cette formulation avec la théorie classique de Biot est présenté. Le modèle est validé par la comparaison des résultats numériques obtenus pour un problème de consolidation d'un sol granulaire avec la solution analytique de Terzaghi. Une approche pour analyser l'hydrodynamique d'un sédiment granulaire est finalement présenté. La reproduction du phénomène de liquéfaction d'un sol est également présentée. / The behaviour of multiphase materials covers a wide range of phenomena of interest to both scientists and engineers. The mechanical properties of these materials originate from all component phases, their distribution and interaction. A new coupled hydromechanical model is presented in this work, to be applied to the analysis of the hydrodynamics of saturated granular media. The model associates the discrete element method (DEM) for the solid phase, and a pore-scale finite volume (PFV) formulation of the flow problem. The emphasis of this model is, on one hand, the microscopic description of the interaction between phases, with the determination of the forces applied on solid particles by the fluid; on the other hand, the model involve affordable computational costs, that allow the simulation of thousands of particles in three dimensional models. The medium is assumed to be saturated of an incompressible fluid. Pore bodies and their connections are defined locally through a regular triangulation of the packings. The analogy of the DEM-PFV model and the classic Biot's theory of poroelasticity is discussed. The model is validated through comparison of the numerical result of a soil consolidation problem with the Terzaghi's analytical solution. An approach to analyze the hydrodynamic of a granular sediment is finally presented. The reproduction of the phenomenon of soil liquefaction is analysed and discussed.

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