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.

The numerical modelling of the VKI-Longshot facility remains a challeng-ing task as it requires multi-physical numerical methods in order to simulate all the components. In the current dissertation, numerical tools were developed in order to study each component of the facility separately and a deep investigations of each stage of the shot were performed. This helped to better understand the different processes involved in the flow development inside this hypersonic wind tunnel. However the numerical computation of different regions of the facility treated as independent from each others remains an approximation at best.The accuracy of the rebuilding code for determining the free stream conditions and the total enthalpy in the VKI-Longshot facility was investigated by using a series of unsteady numerical computations of axisymmetric hypersonic flow over a heat flux probe. Good agreement was obtained between the numerical results and the measured data for both the stagnation pressure and the heat flux dur- ing the useful test time.The driver-driven part of the Longshot facility was modelled using the quasi one-dimensional Lagrangian solver L1d2. The three main conditions used for the experiments —low, medium and high Reynolds number —were considered.The chambrage effect due to the junction between the driver and the driven tubes in the VKI-Longshot facility was investigated. The computation showed great ben- efit of the chambrage in increasing the speed of the piston and thus the final compression ratio of the test gas.Two dimensional simulations of the flow in the driver and the driven tube were performed using Arbitrary Lagrangian Eulerian (ALE) solver in COOLFLuiD. A parallel multi-domain strategy was developed in order to integrate the moving piston within the computational domain.The computed pressure in the reservoir is compared to the one provided by the experiment and good agreement was obtained for both con- editions.Finally, an attempt was made to compute the starting process of the flow in the contoured nozzle. The transient computation of the flow showed how the primary shock initiates the flow in the nozzle before reaching the exit plan at time of 1.5 [ms] after the diaphragm rupture. The complex interactions of the reflected shocks in the throat raise the temperature above 9500 [K] which was not expected. Chemical dissociation of Nitrogen was not taken into account during this transient investigation which may play a key role considering the range of temperature reached near the throat. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Problems governed by partial differential equations (PDEs) in deformable domains, t Rd; d = 2; 3; are of fundamental importance in science and engineering. They are of particular relevance in the design of many engineering systems e.g., aircrafts and bridges as well as to the analysis of several biological phenomena e.g., blood ow in arteries. However, developing numerical scheme for such problems is still very challenging even when the deformation of the boundary of domain is prescribed a priori. Possibility of excessive mesh distortion is one of the major challenge when solving such problems with numerical methods using boundary tted meshes. The arbitrary Lagrangian- Eulerian (ALE) approach is a way to overcome this difficulty. Numerical simulations of convection-dominated problems have for long been the subject to many researchers. Galerkin formulations, which yield the best approximations for differential equations with high diffusivity, tend to induce spurious oscillations in the numerical solution of convection dominated equations. Though such spurious oscillations can be avoided by adaptive meshing, which is computationally very expensive on ne grids. Alternatively, stabilization methods can be used to suppress the spurious oscillations. In this work, the considered equation is designed within the framework of ALE formulation. In the first part, Streamline Upwind Petrov-Galerkin (SUPG) finite element method with conservative ALE formulation is proposed. Further, the first order backward Euler and the second order Crank-Nicolson methods are used for the temporal discretization. It is shown that the stability of the semi-discrete (continuous in time) ALE-SUPG equation is independent of the mesh velocity, whereas the stability of the fully discrete problem is unconditionally stable for implicit Euler method and is only conditionally stable for Crank-Nicolson time discretization. Numerical results are presented to support the stability estimates and to show the influence of the SUPG stabilization parameter in a time-dependent domain. In the second part of this work, SUPG stabilization method with non-conservative ALE formulation is proposed. The implicit Euler, Crank-Nicolson and backward difference methods are used for the temporal discretization. At the discrete level in time, the ALE map influences the stability of the corresponding discrete scheme with different time discretizations, and it leads to schemes where conservative and non-conservative formulations are no longer equivalent. The stability of the fully discrete scheme, irrespective of the temporal discretization, is only conditionally stable. It is observed from numerical results that the Crank-Nicolson scheme induces high oscillations in the numerical solution compare to the implicit Euler and the backward difference time discretiza-tions. Moreover, the backward difference scheme is more sensitive to the stabilization parameter k than the other time discretizations. Further, the difference between the solutions obtained with the conservative and non-conservative ALE forms is significant when the deformation of domain is large, whereas it is negligible in domains with small deformation. Finally, the local projection stabilization (LPS) and the higher order dG time stepping scheme are studied for convection dominated problems. The analysis is based on the quadrature formula for approximating the integrals in time. We considered the exact integration in time, which is impractical to implement and the Radau quadrature in time, which can be used in practice. The stability and error estimates are shown for the mathematical basis of considered numerical scheme with both time integration methods. The numerical analysis reveals that the proposed stabilized scheme with exact integration in time is unconditionally stable, whereas Radau quadrature in time is conditionally stable with time-step restriction depending on the ALE map. The theoretical estimates are illustrated with appropriate numerical examples with distinct features. The second order dG(1) time discretization is unconditionally stable while Crank-Nicolson gives the conditional stable estimates only. The convergence order for dG(1) is two which supports the error estimate.

This thesis details the development, verification and validation of an unsteady unstructured high order (≥ 3) h/p Discontinuous Galerkin - Fourier solver for the incompressible Navier-Stokes equations on static and rotating meshes in two and three dimensions. This general purpose solver is used to provide insight into cross-flow (wind or tidal) turbine physical phenomena. Simulation of this type of turbine for renewable energy generation needs to account for the rotational motion of the blades with respect to the fixed environment. This rotational motion implies azimuthal changes in blade aero/hydro-dynamics that result in complex flow phenomena such as stalled flows, vortex shedding and blade-vortex interactions. Simulation of these flow features necessitates the use of a high order code exhibiting low numerical errors. This thesis presents the development of such a high order solver, which has been conceived and implemented from scratch by the author during his doctoral work. To account for the relative mesh motion, the incompressible Navier-Stokes equations are written in arbitrary Lagrangian-Eulerian form and a non-conformal Discontinuous Galerkin (DG) formulation (i.e. Symmetric Interior Penalty Galerkin) is used for spatial discretisation. The DG method, together with a novel sliding mesh technique, allows direct linking of rotating and static meshes through the numerical fluxes. This technique shows spectral accuracy and no degradation of temporal convergence rates if rotational motion is applied to a region of the mesh. In addition, analytical mappings are introduced to account for curved external boundaries representing circular shapes and NACA foils. To simulate 3D flows, the 2D DG solver is parallelised and extended using Fourier series. This extension allows for laminar and turbulent regimes to be simulated through Direct Numerical Simulation and Large Eddy Simulation (LES) type approaches. Two LES methodologies are proposed. Various 2D and 3D cases are presented for laminar and turbulent regimes. Among others, solutions for: Stokes flows, the Taylor vortex problem, flows around square and circular cylinders, flows around static and rotating NACA foils and flows through rotating cross-flow turbines, are presented.