Cette thèse vise à améliorer la précision numérique des simulations aéroacoustiques d’écoulements dans un contexte précis, celui du cadre industriel avec un partenariat Dassault Aviation. Pour répondre à cette problématique, la modélisation aux grandes échelles est utilisée afin de la rendre plus efficace et adaptée à la méthode numérique des éléments finis stabilisée par SUPG/GLS. Afin de préciser la méthode numérique, une première partie est consacrée à la formulation théorique et pratique du code AETHER utilisé. La précision des schémas numériques spatial et temporel est aussi présentée. L’idéologie principale issue de la famille des modèles Variational Multi-Scale a été retenue afin de construire le nouveau modèle de sous-maille. En effet, une précédente thèse avait démontré la pertinence de ce type d’approche pour les éléments finis. Même si le cadre est applicatif, cette thèse propose une réflexion générale sur le filtrage numérique en éléments finis ainsi qu’un nouveau procédé pour filtrer le plus efficacement l’écoulement calculé. Cette nouvelle approche de filtrage est particulièrement bien adaptée aux éléments finis et à la montée en ordre spatial. Un modèle hybride de gestion des parois est aussi développé afin de pouvoir utiliser le nouveau modèle de sous-maille dans des configurations complexes comprenant des surfaces solides. Le processus de filtrage est testé sur le cas académique des tourbillons de Taylor-Green et présente un réel gain. Enfin le modèle global est utilisé pour calculer une configuration industrielle de tri-corps hypersustenté nommée LEISA II. Grâce au nouveau modèle proposé et validé par les résultats expérimentaux, il a été possible de fournir des interprétations physiques pointues sur le comportement complexe de l’écoulement du bec et du bruit qu’il génère. Cette dernière partie est une illustration pertinente de l’utilisation des modèles aux grandes échelles pourtant coûteux, et cela même dans un contexte industriel. / The goal of this thesis is to improve the numerical accuracy for aeroacoustic flow simulations in a given scope, that is an industrial application for a partnership with the aircraft company Dassault Aviation. These works are then looking for a new large eddy simulation (LES) model which is efficient and well suited for the finite element formulation and the SUPG/GLS stabilisation method. In order to clarify the scientific environment and numerical tools, a first part is devoted to the theoretical and practical framework of the AETHER code. The spatial and temporal performances of its numerical schemes are assessed too. The philosophy of the Variational Multi-scale models has been selected to build an improvement for the new subgrid model. Indeed, a previous thesis had already demonstrated the relevance of this kind of models especially for the finite element method. Despite the industrial framework, a general reflection on the numerical filtering in finite elements is suggested and a new filtering process is developed in order to sort efficiently the scales of the simulated flow. This new filtering method is especially well fitted to finite element simulations and the high spatial order schemes. An hybrid model has been developed too in order to be able to use the new VMS model in complex configurations involving solid bodies. The filtering process is assessed on an academic case called Taylor-Green vortices and shows a real benefit compare to classical approaches. Finally the whole model is used to compute an industrial configuration, a three-element high-lift device called LEISA II. Thank to the validation of the new model with the experimental results, it has been possible to find accurate explanations about the complex flow behaviour of the slat and its noise generation. This last part is a relevant demonstration of the LES models use in the industrial world even if they are still costly in computation ressources.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção. / Convection dominated flows represent a challenge for finite element method simulation. Many methods have been developed to address this problem. In this work we compare the performance of two methods in the solution of the convectiondiffusion and Navier-Stokes equations on environmental flow problems: the Streamline Upwind Petrov Galerkin (SUPG) and the semi-Lagrangian method. In Galerkin finite element methods for fluid flows, the matrix associated with the convective term is non-symmetric, and as a result, the best approximation property is lost. In practice, solutions are often corrupted by espurious oscillations. In this work, we present a semi- Lagrangian method, which is implicitly an upwind method, therefore solving the spurious oscillations problem, and a comparison between this semi-Lagrangian method and the Streamline Upwind Petrov Galerkin (SUPG), an stabilizing method of recognized performance. The SUPG method takes the interpolation and the weighting functions in different spaces, creating an effect so that the spurious oscillations are drastically attenuated. The semi-Lagrangean method is a integration factor method, in which the factor is an operator that shifts to a coordinate system that moves with the fluid, but it resets the Lagrangian coordinate system after each time step. This provides stability and the possibility to take bigger time steps. There are many works in the literature analyzing stabilized methods, but they do not analyze the semi-Lagrangian method, which represents the main contribution of this work: to recognize the strengths and weaknesses of the semi-Lagrangian method in convection dominated flows.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção. / Convection dominated flows represent a challenge for finite element method simulation. Many methods have been developed to address this problem. In this work we compare the performance of two methods in the solution of the convectiondiffusion and Navier-Stokes equations on environmental flow problems: the Streamline Upwind Petrov Galerkin (SUPG) and the semi-Lagrangian method. In Galerkin finite element methods for fluid flows, the matrix associated with the convective term is non-symmetric, and as a result, the best approximation property is lost. In practice, solutions are often corrupted by espurious oscillations. In this work, we present a semi- Lagrangian method, which is implicitly an upwind method, therefore solving the spurious oscillations problem, and a comparison between this semi-Lagrangian method and the Streamline Upwind Petrov Galerkin (SUPG), an stabilizing method of recognized performance. The SUPG method takes the interpolation and the weighting functions in different spaces, creating an effect so that the spurious oscillations are drastically attenuated. The semi-Lagrangean method is a integration factor method, in which the factor is an operator that shifts to a coordinate system that moves with the fluid, but it resets the Lagrangian coordinate system after each time step. This provides stability and the possibility to take bigger time steps. There are many works in the literature analyzing stabilized methods, but they do not analyze the semi-Lagrangian method, which represents the main contribution of this work: to recognize the strengths and weaknesses of the semi-Lagrangian method in convection dominated flows.

In this work, we deal with the finite element method Streamline Upwind/Petrov-Galerkin (SUPG) and use it to solve boundary value problem for the stationary convection-diffusion equation with dominant convection with Dirichlet boundary condition on the whole boundary of bounded polyhedral computational domain of dimension 1 and 2, respectively. We consider a quadratic Lagrangian finite elements on the line segments and triangles, respectively. The core of the work is a proposition of choice of stabilizing parameter of SUPG method as an elementwise affine function in outflow boundary layer and as an elementwise constant function in the rest of the computational domain. We show that this choice gives a more accurate solution than the choice of the stabilization parameter as a constant in each element. 1

Concrete is a popular construction material for bridges, due to its high durability and energy efficiency. An important concern for concrete bridges is the possible occurrence of chloride- induced corrosion in prestressing strands and reinforcing bars, which may substantially impact the service life of such structures. Chloride- induced corrosion is a complicated electrochemical process which is affected by heat transfer, moisture flow and transport of chemical species through the concrete pore network. Reliable and robust analytical tools are required to allow multi-physics simulations of steel corrosion. This study has developed a nonlinear finite element analysis program, called VT-MultiPhys, to enable multi-physics simulations, including analyses of chloride-induced corrosion. The program includes constitutive laws, element formulations and global solution schemes to allow the analysis of steady-state (static) and time-dependent (dynamic) problems, involving multiple, coupled processes such as mechanical deformation, heat transfer, mass flow and chemical reactions combined with advective/diffusive transport of the various species. Special analysis schemes, based on the streamline-upwind Petrov-Galerkin (SUPG) method, have also been implemented to address the spatial instabilities which characterize analyses of advection-dominated transport. The finite element modeling scheme, constitutive laws and boundary conditions for analysis of chloride-induced corrosion are described in detail. The constitutive laws can be combined with inelastic material models to capture the damage (e.g., cracking) due to chloride-induced corrosion. A set of verification analyses is presented, to demonstrate the capabilities of VT-MultiPhys to conduct different types of simulations and reproduce the closed-form analytical solutions of simple cases. Validation analyses for heat conduction, moisture flow and chloride transport, using data from experimental tests in the literature, are also presented. / Master of Science / The deterioration of concrete structures and infrastructures due to the chloride-induced corrosion in prestressing strands and reinforcing bars may substantially impact the service life of such structures. Chloride-induced corrosion is a complicated electrochemical process which is initiated and proceeds due to the chloride attacks at the surfaces of concrete structures and ends in the volume expansion, cracking and spalling of concrete. Due to the lack of comprehensive modeling tool, which can simultaneously comprise the influential factors in chloride-induced corrosion, the realistic estimation of the service life of reinforced concrete structures is still challenging. Reliable and robust analytical tools are required to allow multi-physics simulations of steel corrosion. This study has developed a comprehensive finite element analysis program, called VT-MultiPhys, for calculating and monitoring the contribution of chloride ions to chloride-induced corrosion during service life of concrete structures. The present analysis program enables modeling of the coupled physical process including heat transfer, moisture flow and transport of chemical species through the concrete pore network. Also, by modeling the influence of flexural cracks on chloride transport in concrete, the analysis program is able to predict the rate of steel corrosion in cracked concrete structures. A set of verification analyses is presented, to demonstrate the capabilities of VT-MultiPhys to conduct different types of simulations of heat conduction, moisture flow and chloride transport and the comparison is found to be satisfactory. The element formulations and solution algorithms in VT-MultiPhys also allow the investigation of other long-term deterioration mechanisms, such as carbonation-induced corrosion, alkali-silika reaction (ASR) and sulfate attack. The present contribution will hopefully enable and facilitate future research in these topics, through the formulation and implementation of proper constitutive laws and chemical reaction equations.

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 deals with solving singularly perturbed convection- diffusion equations. Firstly, we construct a matched asymptotic expansion of the solution of the singularly perturbed convection-diffusion equation in 1D and derive a formula for the zeroth-order asymptotic expansion in several two- dimensional polygonal domains. Further, we present a set of stabilization meth- ods for solving singularly perturbed problems and prove the uniform convergence of the Il'in-Allen-Southwell scheme in 1D. Finally, we introduce a modification of the streamline upwind Petrov/Galerkin (SUPG) method on convection-oriented meshes. This new method enjoys several profitable properties such as the ful- filment of the discrete maximum principle. Besides the analysis of the method and derivation of a priori error estimates in respective energy norms we also carry out several numerical experiments verifying the theoretical results.

The evolution of numerical methods and computational facilities allow re- searchers to explore complex physical phenomenons such as multiphase flows. The specific regime of incompressible, turbulent, bubbly two-phase flow (where a car- rier fluid is infused with bubbles or particles) is also receiving increased attention due to it’s appearance in major industrial processes. The main challenges arise from coupling individual aspects of the physics into a unified model and to provide a robust numerical framework. The presented work aimed at to achieve the second part by employing the most frequently used dispersed two-phase flow model and another incompressible, turbulent single phase solver as a base flow provider for coupled Lagrangian or surface tracking tools. Among the numerical techniques, the finite element method is a powerful can- didate when the need arises for multiphysics simulations (for example coupling with an electrochemical module) where the counterpart has a node based ap- proach. Stabilization schemes such as PSPG/SUPG/BULK provide remedies for the pressure decoupling and the inherent instability of the central discretization when applied for convective flow problems. As an alternative to unsteady solvers based upon an explicit or a fully im- plicit nonlinear treatment of the convective terms, a semi-implicit scheme results in a method of second order accurate in both space and time, has absolute linear stability and requires only a single or two linear system solution per time step. The application of the skew symmetric approach to the convective term further stabilizes the solution procedure and in some cases it even prevents divergence. The Eulerian-Eulerian two-phase flow model poses various issues to be over- come. The major difficulty is the density ratio between the phases; for an ordinary engineering problem it is in the order of thousands or more. The seemingly minus- cule differences in the formulation of the stabilizations can cause very different end results and require careful analysis. Volume fraction boundedness is of concern as well, but it is treatable by solving for its logarithm. Since the equations allow jumps (even separation of the phases) in the volume fraction field, discontinuity capturing techniques are also needed. Besides the standard ’spatial’ stabilization temporal smoothing is also necessary, otherwise the limitation in time step size becomes too stringent. Designing a flow solver is one side of the adventure, but verification is equally important. Comparison against analytical solution (such as the single and two- phase Taylor-Green testcase) provides insight and confirmation about the mathe- matical and physical properties. Meanwhile comparing with real life experiments prove the industrialization and usability of a code, dealing with low quality meshes and effective utilization of computer clusters. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished