O presente trabalho tem por objetivo desenvolver uma formulação numérica baseada em Análise Isogeométrica para o estudo de problemas de interação fluido-estrutura (IFE) em aplicações envolvendo corpos rígidos submersos, onde escoamentos incompressíveis de fluidos Newtonianos com superfície livre são considerados. Propõe-se o emprego da Análise Isogeométrica por permitir a unificação entre os procedimentos de pré-processamento e análise, melhorando assim as condições de continuidade das funções de base empregadas tanto na discretização espacial do problema como na aproximação das variáveis do sistema de equações. O sistema de equações fundamentais do escoamento é formado pelas equações de Navier-Stokes e pela equação da conservação de massa, descrita segundo a hipótese de pseudo-compressibilidade, em uma formulação cinemática ALE (Arbitrary Lagrangean- Eulerian). A consideração da superfície livre no escoamento se dá tratando o fluido como um meio bifásico, através do método Level Set. O corpo rígido apresenta não linearidade na rotação e restrições representadas por vínculos elásticos e amortecedores viscosos, sendo a equação de equilíbrio dinâmico resolvida através do método de Newmark. O esquema de acoplamento sólido-fluido adotado é o particionado convencional, que impõe condições de compatibilidade cinemáticas e de equilíbrio sobre a interface sólido-fluido, analisando ambos os meios de maneira sequencial. A discretização das equações governantes é realizada através do esquema explícito de dois passos de Taylor-Galerkin, aplicado no contexto da Análise Isogeométrica. Por fim, são analisados alguns problemas da Dinâmica de Fluidos Computacional, de onde se concluiu que os resultados obtidos são bastante consistentes com os fenômenos envolvidos, com as ferramentas exclusivas da Análise Isogeométrica, como o refinamento k, melhorando a convergência dos resultados. Para escoamentos bifásicos, verificou-se que o método Level Set obteve resultados bastante promissores apresentando, entretanto, uma dissipação numérica excessiva. Propõe-se, para estudos futuros, a elaboração de esquemas numéricos que conservem melhor o volume da fase líquida do escoamento. / The present work aims to development of a numerical formulation based on Isogeometric Analysis for the study of Fluid-Structure Interaction problems in applications involving rigid bodies submerged, considering incompressible Newtonian flows with free surface. The use of the Isogeometric Analysis allows unification between the preprocessing and analysis steps, improving then the continuity of the base functions employed, both in the spatial discretization and approximation of the variables in the system of equations. The fundamental flow equations are formed by the Navier-Stokes and the mass conservation, described by de pseudo-compressibility hypothesis, in an ALE (Arbitrary Lagrangean-Eulerian) kinematic formulation. The free surface consideration of the flow is handled treating the fluid like a two- phase medium, using the Level Set method. The rigid body considers nonlinearity in rotation, and restrictions represented by elastic springs and viscous dampers, with the dynamic equilibrium equation being resolved using the Newmark’s method. The solid-fluid coupling scheme is the conventional partitioned, which imposes kinematics and equilibrium compatibility conditions on the solid-fluid interface, analyzing both mediums in a sequential manner. The governing equations are discretized using the explicit two step Taylor-Galerkin method, applied in an Isogeometric Analisys context. Finally, some Computational Fluid Dinamics problems are analysed, from which it was concluded that the results obtained are quite consistent with phenomena involved, with the unique tools of Isogeometric Analysis, such as k-refinement, improving the convergence of the results. For biphasic flows, it was verified that the Level Set method obtained very promising results, presenting, however, an excessive numerical dissipation. For future studies, it is proposed the elaboration of numerical schemes that better preserve the volume of the liquid phase of the flow.

O método numérico para fluidos incompressíveis desenvolvido no presente estudo é o Moving Particle Semi-Implicit Method (MPS) que enxerga o domínio discretizado em partículas, é baseado em representação lagrangeana e não tem a necessidade de utilização de malhas. O método MPS tem como equações governantes uma forma particular da equação de Navier-Stokes e a equação da continuidade para fluidos incompressíveis e não viscosos. Os métodos de simulação de fluidos mais comumente utilizados são baseados em representação euleriana e utilizam malhas para descrever a geometria do domínio a ser simulado. Devido a essas diferenças, uma das grandes virtudes do método de partículas é a facilidade de investigação de fenômenos altamente não-lineares como o de superfície livre com quebra de ondas, de líquidos no interior de uma embarcação em movimento, de ondas batendo na parte externa do casco de um navio, etc. Em artigos já publicados, resultados de experimentos físicos mostraram boa aderência aos resultados numéricos de simulações realizadas com o método MPS. No presente trabalho, resultados das forças de excitação das simulações com ondas regulares foram comparados com os resultados do programa Wave Analysis MIT (WAMIT) que é um programa consagrado no meio científico. Houve uma boa concordância de resultados entre os dois programas. A otimização do cálculo de vizinhança forneceu uma grande economia de tempo computacional. A maior contribuição deste estudo foi a otimização da função que resolve o sistema linear implementando no código desenvolvido um código paralelizado de uso público existente chamado Portable, Extensible Toolkit for Scientific Computation (PETSc) que proporcionou um bom ganho de desempenho. / A numerical method called Moving Particle Semi-implicit (MPS) method was developed in this study to analyze incompressible fluids. It is a particle method using a lagrangean representation without any grid. The governing equations are the Navier-Stokes equation and continuity equation for incompressible and non-viscous flow. Most of the computational fluid dynamics (CFD) methods are based on eulerian representation and use grids to describe the geometry of the simulated domain. These differences make the MPS method easier to analyze highly nonlinear phenomena as free surface with wave breaking, sloshing, slamming, etc. In previously published articles, results of physical experiments had shown good agreement with the numerical results obtained with MPS method. In the present work, results of exciting forces were compared with the results obtained with a validated program called Wave Analysis MIT (WAMIT). It had a good agreement of results between these two programs. The optimization of the neighborhood calculation function got a good economy of computational time. The greatest contribution of this study was the optimization of the linear system solver. It was made implementing in the developed code a parallelized public code called Portable, Extensible Toolkit for Scientific Computation (PETSc) that provided a good performance profit.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-03-07T11:34:58Z No. of bitstreams: 1 maiconwilliamniebusodone.pdf: 9396542 bytes, checksum: 0a9d05fa621e639f71fa74726b219707 (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-03-07T15:03:29Z (GMT) No. of bitstreams: 1 maiconwilliamniebusodone.pdf: 9396542 bytes, checksum: 0a9d05fa621e639f71fa74726b219707 (MD5) / Made available in DSpace on 2017-03-07T15:03:29Z (GMT). No. of bitstreams: 1 maiconwilliamniebusodone.pdf: 9396542 bytes, checksum: 0a9d05fa621e639f71fa74726b219707 (MD5) Previous issue date: 2014-08-28 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Problemas de mudança de fase líquido-sólido e/ou sólido-líquido são amplamente encontrados em diversas aplicações da indústria tais como na moldagem de diferentes tipos de peças, na aplicação de sprays térmicos para tratamento de superfícies e na fundição de materiais energéticos tais como o trinitrotolueno (TNT). Neste último caso, contrações de material, trincas e formação de vazios são frequentemente observados, devendo estes comportamentos serem previstos e evitados. Sendo assim, a simulação numérica de tais processos com algoritmos capazes de acompanhar a evolução da superfície livre, e a distribuição da temperatura no corpo durante o processo de fundição é importante. Pretende-se estudar diferentes estratégias na modelagem e solução de problemas desta natureza através das equações de transporte de massa, transferência de calor e de mudança de fase por meio do método dos volumes finitos. Na metodologia usada, os processos de solidificação e fusão são tratados por meio da aproximação do método da entalpia-porosidade baseado em malhas fixas. / Phase change problems as liquid-solid and/or solid-liquid are widely found in several industrial applications such as in molding of different parts, the application of thermal spray treatment and the casting of energy materials such as trinitrotoluene (TNT). For the latter case, material contractions, cracks and empities formation are often observed and these behaviors should be anticipated and avoided. Thus, the numerical simulation of such processes with algorithms capable to follow free surface evolution and body temperature distribution during casting process is important. In this work we study strategies in modeling and solving problems of this nature through the mass transport equations, heat transfer and phase change using the finite volume method. In the methodology used, the melting and solidification process are treated by means of the enthalpy-porosity approximation method based on fixed grids.

Micro-assembly processes suffer from some breaches due to the continuing trend towards an increase in the production capabilities as well as in the size reduction of the components manipulated. Usual manipulating schemes have reached their limit and capillary forces constitute a valuable alternative strategy.<p><p>The goal of this work is to describe the dynamics of liquid bridges in the application of micro-assembly processes. The description is obtained using the Kelvin-Voigt model, with a spring, a damper, and a mass connected in parallel, supported by numerical simulations, analytical approximations and experiments.<p><p>The works is divided into three parts. First we present important aspects of microfluidics, as well as the constitutive equations and an overview of numerical approaches used to describe fluid flow problems with moving interfaces.<p><p>The second part is devoted to the capillary rise case, intended to validate and to compare the numerical approaches to analytical laws and experimental results. The implementation of the slipping and the dynamic contact angles is discussed.<p><p>The last part focuses on the dynamics of the liquid bridge. The liquid bridge is confined between two circular and parallel plates and presents an axial symmetry. The description reveals that the stiffness depends on the surface tension and on the shape of the air/liquid interface, the damping coefficient depends on the viscosity and the volume of liquid and the equivalent mass depends on the density and the volume.<p> / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

No contexto do método MAC e baseado em esquemas de diferenças finitas, este trabalho apresenta três estudos: i) uma análise de estabilidade, ii) o desenvolvimento de técnicas implícitas e, iii) a construção de métodos de projeção para escoamentos com superfície livre. Na análise de estabilidade, o principal resultado mostra que o método de Crank-Nicolson torna-se condicionalmente estável quando aplicado para uma malha deslocada com a discretiza ção explícita das condições de contorno do tipo Dirichlet. Entretanto, o mesmo método com condições de contorno implícitas é incondicionalmente estável. Para obter métodos mais estáveis, formulações implícitas são desenvolvidas para a equação da pressão na superfície livre, derivada da condição de tensão normal. Esta estratégia resulta no acoplamento dos campos de velocidade e pressão, o que exige a introdução de novos métodos de projeção. Os métodos de projeção assim desenvolvidos resultam em novas metodologias para escoamentos com superfície livre que são apropriados para o tratamento de problemas com baixo número de Reynolds. Além disso, mostra-se que os métodos propostos podem ser aplicados para fluidos viscoelásticos. Novas estratégias são derivadas para obter métodos de projeção de segunda ordem de precisão para escoamentos com superfícies livres. Além dos resultados teóricos sobre a estabilidade de esquemas numéricos, técnicas implícitas e métodos de projeção, testes computacionais são realizados e comparados para consolidação da teoria apresentada. Os resultados numéricos são obtidos no sistema FREEFLOW. A eficiência e robustez das técnicas desenvolvidas neste trabalho são demonstradas na solução de problemas tridimensionais complexos com superfície livre e baixo número de Reynolds, incluindo os problemas do jato oscilante e do inchamento do extrudado / In the context of the MAC method and based on finite difference schemes, this work presents three studies: i) a stability analysis, ii) the development of implicit techniques, and iii) the construction of projection methods for free surface flows. In the stability analysis, the main result shows a precise stability restriction on the Crank-Nicolson method when one uses a staggered grid with Dirichlet explicit boundary conditions. However, the same method with implicit boundary conditions becomes unconditionally stable. In order to obtain more stable methods, implicit formulations are applied for the pressure equation at the free surface, which is derived from the normal stress condition. This approach results in a coupling of the velocity and pressure fields; hence new projection methods for free surface flows need to be developed. The developed projection methods result in new methodologies for low Reynolds number free surface flows. It is also shown that the proposed methods can be applied for viscoelastic fluids. New strategies are derived for obtaining second-order accurate projection methods for free surface flows. In addition to the theoretical results on the stability of numerical schemes, implicit techniques and projection methods, computational tests are carried out and the results compared to consolidate the theory. The numerical results are obtained by the FREEFLOW system. The eficiency and robustness of the techniques in this work are demonstrated by solving complex tridimensional problems involving free surface and low Reynolds numbers, including the jet buckling and the extrudate swell problems

Escoamentos viscoelásticos não estacionários com superfícies livres são comuns em muitos processos industriais e diversas técnicas numéricas têm sido empregadas para reproduzir computacionalmente estes processos. A maioria dos modelos empregados utiliza equações diferenciais na definição do tensor de tensões. Porém, para alguns grupos de fluidos complexos, por exemplo, fluidos de Boger, os modelos integrais mostram-se mais capacitados em fornecer uma boa aproximação para os comportamentos não lineares desses fluidos. Este trabalho trata da solução numérica do modelo constitutivo integral KBKZ-PSM para escoamentos transientes bidimensionais com superfícies livres. O método numérico proposto é uma técnica numérica que utiliza diferenças finitas para simular escoamentos com superfícies livres na presença de paredes sólidas. As principais características do método numérico proposto são: solução das equações de conservação de quantidade de movimento e massa utilizando um método semi-implícito; a condição de contorno na superfície livre é acoplada à equação de Poisson, o que garante conservação de massa; a discretização do tempo t é realizada por uma nova técnica numérica; o tensor de Finger é calculado pelo método dos campos de deformação e avançado no tempo pelo método de Euler modificado. Essa nova técnica é verificada em escoamentos cisalhantes e elongacionais. Adicionalmente, uma solução analítica desenvolvida para escoamentos em canais bidimensionais é empregada para verificar e analisar a convergência do método proposto. Com relação a escoamentos com superfícies livres, a convergência é verificada por meio de refinamento de malha nas simulações de um jato incidindo sobre placa rígida e no problema do inchamento do extrudado. Finalmente, o método é aplicado para investigar os problemas jet buckling e inchamento do extrudado de fluidos KBKZ-PSM. / Unsteady viscoelastic free surface flows are common in many industrial processes and a variety of numerical techniques have been employed to simulate these flows. The majority of constitutive models employed are based on differential equations to define the extra stress tensor. However, for some complex fluids, for instance, Boger fluids, integral models are more adequate to approximate the nonlinear behaviour of these fluids. This work deals with the numerical solution of the integral constitutive model KBKZ-PSM for two-dimensional unsteady free surface flows. The proposed numerical method is a numerical technique that employs finite differences to simulate moving free surface flows that interact with solid walls. The main features of the method are: numerical solution of the momentum and mass equations by an implicit method; the pressure condition on the free surface is implicitly coupled with the Poisson equation for obtaining the pressure field from mass conservation; a novel scheme for defining the past times t is employed; the Finger tensor is calculated by the deformation fields method and is advanced in time by the modified Euler method. This new technique is verified by solving shear and uniaxial elongational flows. Moreover, an analytic solution for channel flows is obtained that is used in the verification and convergence analysis of the proposed methodology. For free surface flows, the assessment of convergence lies on the mesh refinement on the simulation of a jet impinging on a flat surface and the extrudade swell problem. Finally, the new method is applied to investigate the jet buckling phenomenon and extrudate swell of KBKZ-PSM fluids.

Laser welding is commonly used in the automotive-, steel- and aerospace industry. It is a highly non-linear and coupled process where the weld geometry is strongly affected by the flow pattern in the melt pool. Experimental observations are challenging since the melt pool and melt flow below the surface are not yet accessible during welding. Improved process control would allow maintaining, or improving, product quality with less material and contribute further to sustainability by reducing production errors. Numerical modelling with Computational Fluid Dynamics, CFD, provides complementary understanding with access to process properties that are not yet reachable with experimental observation. However, the existing numerical models lack predictability when considering the weld shape. The work presented here is the development of a model for conduction laser welding. The solver upon which the model is based is first described in detail. Then different validation cases are applied in order to test specific parts of the physics implemented. Two cases focus on thermocapillary convection in two-phase and three-phase flows with surface deformation. Finally, a third case considers the melt pool flow during conduction mode welding.It is concluded that the convection of fusion enthalpy, which was neglected in former studies, should be included in the model. The implementation of the thermo capillary force is recommended to be consistent with the other surface forces to avoid unphysical solution. Free surface oscillations, known from experimental observations, are also computed numerically. However, further investigation is needed to check that these oscillations are not disturbed b ynumerical oscillations.

Dans cette thèse, on s’intéresse à l’approximation numérique d’un modèle d’écoulement dispersif en eau peu profonde. Les applications visées par ce type de modélisation sont nombreuses (écoulement dans les océans, les rivières, etc) et cette thèse est motivée en particulier par les risques naturels et la production d’énergie renouvelable. Le modèle étudié a été dérivé par moyenne selon la verticale des équations d’Euler et prend en compte la pression non-hydrostatique. Il est alors nécessaire de résoudre un système de type incompressible; ce qui nous amène à résoudre une équation elliptique en pression. Nous proposons une méthode numérique pour résoudre le système dispersif avec topographie pour les modèles 1D et 2D. L’approche développée est basée sur un schéma de type prediction-correction, initialement introduit par Chorin-Temam pour les équations de Navier-Stokes. Nous définissons un cadre générique qui permet de concevoir un schéma valable en 1D et 2D et aussi de pouvoir augmenter l’ordre de précision. Ainsi, nous proposons une formulation variationnelle qui nous permet d’appliquer la méthode des éléments finis avec des choix d’espaces compatibles. Le travail effectué étant destiné à simuler des processus géophysiques réels, la méthode a été conçue pour pouvoir traiter les transitions de sol sec/mouillé et cette propriété a été confirmée par plusieurs tests numériques. Afin de valider la méthode, nous présentons la comparaison entre certaines solutions analytiques et leurs simulations numériques. / In this PhD thesis we are interested in the numerical approximation of a dispersive shallow water system, aimed at modeling the free surface flows (e.g. ocean and rivers) and motivated by applications for natural hazards and sustainable energy resources. This model is a depth-averaged Euler system and takes into account a non-hydrostatic pressure which brings crucial information for understanding the behavior of the flow, particularly when dispersion occur. We develop a numerical method for the one- and the two-dimensional dispersive shallow water system with a topography. The approach is based on a prediction-correction method initially introduced by Chorin-Temam, and we establish a global framework in order to easily increase the order of accuracy of the method. The prediction part leads to solving a shallow water system for which we use finite volume methods, while the correction part leads to solving a mixed problem in velocity and pressure. We propose a variational formulation of the mixed problem which allows us to apply a finite element method with compatible spaces. In this framework we establish compatible boundary conditions between the prediction part and the correction part. The method is performed for the one-dimensional model and for the two-dimensional problem on unstructured grids. In order to make the method practical for real geophysical cases, we have derived a scheme able to treat wet/dry interfaces and to this end we give many examples to test its performance. Moreover, we provide a comparison of simulated solutions with data from laboratory experiments.

Cette thèse porte sur le problème du ballottement d'un liquide dans un réservoir, rencontré dans le transport et le stockage du GNL par des structures flottantes. La prédiction des chargements réels, dus au ballotement sur les parois du réservoir, est souvent basée sur des études expérimentales à petite échelle. La modélisation expérimentale à petite échelle respecte la similitude de Froude et le rapport de densité entre le gaz et le liquide. Cependant, d’autres similitudes sont biaisées comme la similitude par rapport au nombre de Weber et la similitude par rapport au nombre de Reynolds. De plus, les pressions enregistrées montrent une grande variabilité quand le même essai est répété. Dans une première partie, différentes sources physiques responsables de cette variabilité sont discutées, à savoir les instabilités de surface libre, la retombée des gouttes et des jets liquides sur la surface libre, et la production et l'entraînement des bulles dans le liquide. En fait, ces phénomènes sont à l'origine des perturbations de l'écoulement, de la variabilité de la géométrie de la vague et de cette façon des pressions engendrées par cette dernière sur la paroi. D'autres mécanismes de dissipation d’énergie sont identifiés. Ils sont liés aux frottements aux parois et aux déferlements de vagues. Nous montrons que cette dissipation induit un effet mémoire à courte durée pour l’écoulement, permettant de reproduire pour chaque impact la distribution statistique des pics de pression avec une courte durée des excitations. Ces sources de variabilité et ces mécanismes de dissipation dépendent de la tension de surface et de la viscosité du liquide. Ainsi nous étudions dans une deuxième partie, les effets de ces paramètres physiques. Nous montrons que la forme locale de la vague dépend de la tension de surface. Par contre, les effets sur la forme globale de la vague sont négligeables. Plus la tension de surface diminue, plus les pics pression sont faibles. Ce qui est dû aux différents phénomènes liés au développement des ligaments, la fragmentation en gouttes et la génération de la mousse sur la crête de la vague, et à l’entraînement des bulles dans le liquide. Quant à la viscosité du liquide, elle affecte à la fois la forme globale et la forme locale de la vague, là encore les pressions sont modifiées. Cette étude paramétrique permet, dans une troisième partie, d'étudier et comprendre les effets du nombre de Weber et du nombre de Reynolds, en comparant les résultats pour deux échelles différentes 1:40 et 1:20, quand les mêmes fluides sont considérés. De plus, en se basant sur différents cas de comparaison avec la similitude de Reynolds et/ou la similitude de Weber, nous montrons que la double similitude est indispensable pour obtenir une forme de vague avant l'impact indépendante de l'échelle. Cependant, la distribution statistique des pics de pression dépend aussi d’autres nombres adimensionnels à savoir le nombre de Mach du liquide et le nombre de Mach du gaz. / This work focuses on sloshing problem, encountered in the transport and storage of LNG by floating structures. The prediction of real sloshing loads is often based on small-scale experimental studies, respecting the Froude similarity and the density ratio between the gas and the liquid. However, other similarities are biased such as the Weber similarity and the Reynolds similarity. In addition, the recorded pressures show great variability when the same test is repeated. In a first part, different physical sources responsible for this variability are discussed, which are the free surface instabilities, the falling droplets and liquid jets impinging on the free surface, and the liquid entrainment by bubbles. In fact, these phenomena are at the origin of the flow disturbances, the variability of the wave shape, and hence its pressures on the wall. Other dissipation mechanisms are identified. They are related to wall frictions and breaking waves. Thanks to this energy dissipation, we show that the flow is characterized by a short-term memory, making it possible to reproduce for each impact its statistical distribution of pressure peaks with a short duration of excitations. These sources of variability and dissipation mechanisms depend on the surface tension and the viscosity of the liquid. Thus, we study, in a second part, these physical parameters. We show that the local wave shape depends on the surface tension. However, its effects on the global wave shape are negligible. Besides, when the surface tension is reduced, the statistical pressures are reduced. This is due to various phenomena related to the development of liquid ligaments, their fragmentation into drops and the generation of foam at the wave crest, and the liquid entrainment by bubbles. As for the viscosity of the liquid, it affects both the local and global shape wave shapes, again the pressures are changed. Based on this parametric study, The effects of Weber number and Reynolds number are studied by comparing the results for two different scales 1:40 and 1:20, when the same fluids are used. Moreover, considering different cases of comparison with Reynolds number similarity and / or Weber number similarity, the results show that both similarities are essential to obtain a scaleindependent wave shape. However, the statistical distribution of pressure peaks also depends on other dimensionless numbers, namely the Mach number of the liquid and the Mach number of the gas.

The diploma thesis deals with the deposition of layers of 2-ethyl-2-oxazoline monomer in a dielectric barrier discharge in a nitrogen atmosphere. The theoretical part describes polyoxazolines, dielectric barrier discharge, plasma polymerization and diagnostic methods that describe the characteristics of the prepared layers. The experimental part deals with the deposition of layers and the determination of physicochemical properties. Finally, the results are compared with layers that were made of 2-methyl-2-oxazoline.