This is a book on numerical methods for singular perturbation problems - in particular stationary convection-dominated convection-diffusion problems. More precisely it is devoted to the construction and analysis of layer-adapted meshes underlying these numerical methods. An early important contribution towards the optimization of numerical methods by means of special meshes was made by N.S. Bakhvalov in 1969. His paper spawned a lively discussion in the literature with a number of further meshes being proposed and applied to various singular perturbation problems. However, in the mid 1980s this development stalled, but was enlivend again by G.I. Shishkin's proposal of piecewise- equidistant meshes in the early 1990s. Because of their very simple structure they are often much easier to analyse than other meshes, although they give numerical approximations that are inferior to solutions on competing meshes. Shishkin meshes for numerous problems and numerical methods have been studied since and they are still very much in vogue. With this contribution we try to counter this development and lay the emphasis on more general meshes that - apart from performing better than piecewise-uniform meshes - provide a much deeper insight in the course of their analysis. In this monograph a classification and a survey are given of layer-adapted meshes for convection-diffusion problems. It tries to give a comprehensive review of state-of-the art techniques used in the convergence analysis for various numerical methods: finite differences, finite elements and finite volumes. While for finite difference schemes applied to one-dimensional problems a rather complete convergence theory for arbitrary meshes is developed, the theory is more fragmentary for other methods and problems and still requires the restriction to certain classes of meshes.

Réalisé en cotutelle avec Aix Marseille Université. / Cette thèse se découpe en trois parties. Les deux premières portent sur le développement de méthodes pour la construction de modèles géométriques moyens et pour la comparaison de modèles. Ces approches sont appliquées à la cornée humaine pour l’élaboration d’atlas et pour l’étude biométrique robuste. La troisième partie porte sur une méthode générique d'extraction d'informations dans un maillage en s'appuyant sur des propriétés différentielles discrètes afin de construire une structure par graphe permettant l'extraction de caractéristiques par une description sémantique. Les atlas anatomiques conventionnels (papier ou CD-ROM) sont limités par le fait qu'ils montrent généralement l'anatomie d'un seul individu qui ne représente pas nécessairement bien la population dont il est issu. Afin de remédier aux limitations des atlas conventionnels, nous proposons dans la première partie d’élaborer un atlas numérique 3D contenant les caractéristiques moyennes et les variabilités de la morphologie d'un organe, plus particulièrement de la cornée humaine. Plusieurs problématiques sont abordées, telles que la construction d'une cornée moyenne et la comparaison de cornées. Il existe à ce jour peu d'études ayant ces objectifs car la mise en correspondance de surfaces cornéennes est une problématique non triviale. En plus d'aider à développer une meilleure connaissance de l'anatomie cornéenne, la modélisation 3D de la cornée normale permet de détecter tout écart significatif par rapport à la "normale" permettant un diagnostic précoce de pathologies ou anomalies de la forme de la cornée. La seconde partie a pour objectif de développer une méthode pour reconnaître une surface parmi un groupe de surfaces à l’aide de leurs acquisitions 3D respectives, dans le cadre d’une application de biométrie sur la cornée. L’idée est de quantifier la différence entre chaque surface et une surface donnée, et de déterminer un seuil permettant la reconnaissance. Ce seuil est dépendant des variations normales au sein d’un même sujet, et du bruit inhérent à l’acquisition. Les surfaces sont rognées et trouées de façon imprévisible, de plus il n’y a pas de point de mise en correspondance commun aux surfaces. Deux méthodes complémentaires sont proposées. La première consiste à calculer le volume entre les surfaces après avoir effectué un recalage, et à utiliser ce volume comme un critère de similarité. La seconde approche s’appuie sur une décomposition en harmoniques sphériques en utilisant les coefficients comme des descripteurs de forme, qui permettront de comparer deux surfaces. Des résultats sont présentés pour chaque méthode en les comparant à la méthode la plus récemment décrite dans la littérature, les avantages et inconvénients de chacune sont détaillés. Une méthodologie en cascade utilisant ces deux méthodes afin de combiner les avantages de chacune est aussi proposée. La troisième et dernière partie porte sur une nouvelle méthode de décomposition en graphes de maillages 3D triangulés. Nous utilisons des cartes de courbures discrètes comme descripteur de forme afin de découper le maillage traité en huit différentes catégorie de carreaux (ou peak, ridge, saddle ridge, minimal, saddle valley, valley, pit et flat). Ensuite, un graphe d'adjacence est construit avec un nœud pour chaque carreau. Toutes les catégories de carreaux ne pouvant pas être adjacentes dans un contexte continu, des jonctions intermédiaires sont ajoutées afin d'assurer une cohérence continue entre les zones. Ces graphes sont utilisés pour extraire des caractéristiques géométriques décrites par des motifs (ou patterns), ce qui permet de détecter des régions spécifiques dans un modèle 3D, ou des motifs récurrents. Cette méthode de décomposition étant générique, elle peut être appliquée à de nombreux domaines où il est question d’analyser des modèles géométriques, en particulier dans le contexte de la cornée. / This thesis comprises three parts. The first two parts concern the development of methods for the construction of mean geometric models and for model comparison. These approaches are applied to the human cornea for the construction of atlases and a robust biometric study. The third part focuses on a generic method for the extraction of information in a mesh. This approach is based on discrete differential properties for building a graph structure to extract features using a semantic description. Conventional anatomical atlases (paper or CD-ROM) are limited by the fact they generally show the anatomy of a single individual who does not necessarily represent the population from which they originate. To address the limitations of conventional atlases, we propose in the first part of this thesis to construct a 3D digital atlas containing the average characteristics and variability of the morphology of an organ, especially that of the human cornea. Several issues are addressed, such as the construction of an average cornea and the comparison of corneas. Currently, there are few studies with these objectives because the matching of corneal surfaces is a non-trivial problem. In addition to help to develop a better understanding of the corneal anatomy, 3D models of normal corneas can be used to detect any significant deviation from the norm, thereby allowing for an early diagnosis of diseases or abnormalities using the shape of the cornea. The second part of this thesis aims to develop a method for recognizing a surface from a group of surfaces using their 3D acquisitions in a biometric application pertinent to the cornea. The concept behind this method is to quantify the difference between each surface and a given surface and to determine the threshold for recognition. This threshold depends on normal variations within the same subject and noise due to the acquisition system. The surfaces are randomly trimmed and pierced ; moreover, there is no common landmark on the surfaces. Two complementary methods are proposed. The first method consists of the computation of the volume between the surfaces after performing geometrical matching and the use of this volume as a criterion of similarity. The second approach is based on a decomposition of the surfaces into spherical harmonics using the coefficients as shape descriptors to compare the two surfaces. Each result of the proposed methods is compared to the most recent method described in the literature, with the benefits and disadvantages of each one described in detail. A cascading methodology using both methods to combine the advantages of each method is also proposed. The third and final part of this thesis focuses on a new method for decomposing 3D triangulated meshes into graphs. We use discrete curvature maps as the shape descriptor to split the mesh in eight different categories (peak, ridge, saddle ridge, minimal, saddle valley, valley, pit and flat). Next, an adjacency graph is built with a node for each patch. Because all categories of patches cannot be adjacent in a continuous context, intermediate junctions are added to ensure the continuous consistency between patches. These graphs are used to extract geometric characteristics described by patterns that allow for the detection of specific regions in a 3D model or recurrent characteristics. This decomposition method, being generic, can be used in many applications to analyze geometric models, especially in the context of the cornea.

Este trabalho apresenta uma estratégia numérica desenvolvida usando o método dos elementos finitos para simular o processo de falha de compósitos cimentícios reforçados com fibras de aço. O material é descrito como um compósito composto por três fases: matriz cimentícia (pasta, argamassa ou concreto), fibras descontínuas discretas, e interface fibra-matriz. Um novo esquema de acoplamento para malhas de elementos finitos não-conformes foi desenvolvido para acoplar as malhas geradas independentes, da matriz cimentícia e de uma nuvem de fibras de aço, baseado na utilização de novos elementos finitos desenvolvidos, denominados elementos finitos de acoplamento. Utilizando este esquema de acoplamento, um procedimento não-rígido é proposto para a modelagem do complexo comportamento não linear da interface fibra-matriz, utilizando um modelo constitutivo de dano apropriado para descrever a relação entre a tensão de cisalhamento (tensão de aderência) e deslizamento relativo entre a matriz e cada fibra de aço individualmente. Este esquema também foi adotado para considerar a presença de barras de aço para as análises de estruturas de concreto armado. As fibras de aço são modeladas usando elementos finitos lineares com dois nós (elementos de treliça) com modelo material elastoplástico. As fibras são posicionadas usando uma distribuição randômica uniforme isotrópica, considerando o efeito parede. Uma abordagem contínua e outra descontínua são investigadas para a modelagem do comportamento frágil da matriz cimentícia. Para a primeira, é utilizado um modelo de dano isotrópico com duas variáveis de dano para descrever o comportamento de dano à tração e à compressão. A segunda emprega uma técnica de fragmentação de malha que utiliza elementos finitos degenerados, posicionados entre todos os elementos finitos que formam a matriz cimentícia. Para esta técnica é proposto um modelo constitutivo à tração, compatível com a abordagem descontínua forte contínua, para prever a propagação de fissura. Para acelerar o cálculo e aumentar a robustez dos modelos de dano contínuos para simular o processamento de falhas, um esquema de integração implícito-explícito é utilizado. Exemplos numéricos são apresentados ao longo do desenvolvimento desta tese. Inicialmente, exemplos numéricos com um único reforço são apresentados para validar a técnica desenvolvida e para investigar à influência das propriedades geométricas 7 das fibras e sua posição em relação à superfície de falha. Posteriormente, exemplos mais complexos são considerados envolvendo uma nuvem de fibras. Nestes casos, atenção especial é dada à influência da distribuição das fibras no comportamento do compósito relacionado ao processo de fissuração. Comparações com resultados experimentais demonstram que a aplicação da ferramenta numérica para modelar o comportamento de compósitos cimentícios reforçados com fibras de aço é muito promissora e pode ser utilizada como uma importante ferramenta para melhor entender os efeitos dos diferentes aspectos envolvidos no processo de falha deste material. / This work presents a numerical strategy developed using the Finite Element Method (FEM) to simulate the failure process of Steel Fiber Reinforced Cementitious Composites (SFRCCs). The material is described as a composite made up by three phases: a cementitious matrix (paste, mortar or concrete), discrete discontinuous fibers, and a fiber-matrix interface. A novel coupling scheme for non-matching finite element meshes has been developed to couple the independent generated meshes of the bulk cementitious matrix and a cloud of discrete discontinuous fibers based on the use of special finite elements developed, termed Coupling Finite Elements (CFEs). Using this approach, a nonrigid coupling procedure is proposed for modeling the complex nonlinear behavior of the fiber-matrix interface by adopting an appropriate constitutive damage model to describe the relation between the shear stress (adherence stress) and the relative sliding between the matrix and each fiber individually. This scheme has also been adopted to account for the presence of regular reinforcing bars in the analysis of reinforced concrete structural elements. The steel fibers are modeled using two-node finite elements (truss elements) with a one-dimensional elastoplastic constitutive model. They are positioned using an isotropic uniform random distribution, considering the wall effect of the mold. Continuous and discontinuous approaches are developed to model the brittle behavior of the bulk cementitious matrix. For the former, an isotropic damage model including two independent scalar damage variables for describing the composite behavior under tension and compression is considered. The discontinuous approach is based on a mesh fragmentation technique that employs degenerated solid finite elements in between all regular (bulk) elements. In this case, a tensile damage constitutive model, compatible with the Continuum Strong Discontinuity Approach (CSDA), is proposed to predict crack propagation. To increase the computability and robustness of the continuum damage models used to simulate the failure processes in both of the strategies, an implicit-explicit integration scheme is used. Numerical analyses are performed throughout the presentation of the work. Initially, numerical examples with a single reinforcement are presented to validate the technique and to investigate the influence of the fibers geometrical properties and its position relative to the crack surface. Then, more complex examples involving a cloud of steel fibers are considered. In these cases, special attention is given to the analysis of the influence of the fiber distribution on the composite behavior relative to the cracking process. Comparisons with experimental results demonstrate that the application of the numerical tool for modeling the behavior of SFRCCs is very promising and may constitute an important tool for better understanding the effects of the different aspects involved in the failure process of this material.

Le calcul scientifique parallèle est un domaine en plein essor qui permet à la fois d’augmenter la vitesse des longs traitements, de traiter des problèmes de taille plus importante ou encore des problèmes plus précis. Ce domaine permet donc d’aller plus loin dans les calculs scientifiques, d’obtenir des résultats plus pertinents, car plus précis, ou d’étudier des problèmes plus volumineux qu’auparavant. Dans le monde plus particulier de la simulation numérique scientifique, la résolution d’équations aux dérivées partielles (EDP) est un calcul particulièrement demandeur de ressources parallèles. Si les ressources matérielles permettant le calcul parallèle sont de plus en plus présentes et disponibles pour les scientifiques, à l’inverse leur utilisation et la programmation parallèle se démocratisent difficilement. Pour cette raison, des modèles de programmation parallèle, des outils de développement et même des langages de programmation parallèle ont vu le jour et visent à simplifier l’utilisation de ces machines. Il est toutefois difficile, dans ce domaine dit du “parallélisme implicite”, de trouver le niveau d’abstraction idéal pour les scientifiques, tout en réduisant l’effort de programmation. Ce travail de thèse propose tout d’abord un modèle permettant de mettre en oeuvre des solutions de parallélisme implicite pour les simulations numériques et la résolution d’EDP. Ce modèle est appelé “Structured Implicit Parallelism for scientific SIMulations” (SIPSim), et propose une vision au croisement de plusieurs types d’abstraction, en tentant de conserver les avantages de chaque vision. Une première implémentation de ce modèle, sous la forme d’une librairie C++ appelée SkelGIS, est proposée pour les maillages cartésiens à deux dimensions. Par la suite, SkelGIS, et donc l’implémentation du modèle, est étendue à des simulations numériques sur les réseaux (permettant l’application de simulations représentant plusieurs phénomènes physiques). Les performances de ces deux implémentations sont évaluées et analysées sur des cas d’application réels et complexes et démontrent qu’il est possible d’obtenir de bonnes performances en implémentant le modèle SIPSim. / Parallel scientific computations is an expanding domain of computer science which increases the speed of calculations and offers a way to deal with heavier or more accurate calculations. Thus, the interest of scientific computations increases, with more precised results and bigger physical domains to study. In the particular case of scientific numerical simulations, solving partial differential equations (PDEs) is an especially heavy calculation and a perfect applicant to parallel computations. On one hand, it is more and more easy to get an access to very powerfull parallel machines and clusters, but on the other hand parallel programming is hard to democratize, and most scientists are not able to use these machines. As a result, high level programming models, framework, libraries, languages etc. have been proposed to hide technical details of parallel programming. However, in this “implicit parallelism” field, it is difficult to find the good abstraction level while keeping a low programming effort. This thesis proposes a model to write implicit parallelism solutions for numerical simulations such as mesh-based PDEs computations. This model is called “Structured Implicit Parallelism for scientific SIMulations” (SIPSim), and proposes an approach at the crossroads of existing solutions, taking advantage of each one. A first implementation of this model is proposed, as a C++ library called SkelGIS, for two dimensional Cartesian meshes. A second implementation of the model, and an extension of SkelGIS, proposes an implicit parallelism solution for network-simulations (which deals with simulations with multiple physical phenomenons), and is studied in details. A performance analysis of both these implementations is given on real case simulations, and it demonstrates that the SIPSim model can be implemented efficiently.

One of the characteristics of the high speed ows over blunt bodies is the detached shock formed in front of the body. The distance of the shock from the stagnation point measured along the stagnation streamline is termed as the shock stand o distance or the shock detachment distance. It is one of the most basic parameters in such ows. The need to know the shock stand o distance arises due to the high temperatures faced in these cases. The biggest challenge faced in high enthalpy ows is the high amounts of heat transfer to the body. The position of the shock is relevant in knowing the temperatures that the body being subjected to such ows will have to face and thus building an efficient system to reduce the heat transfer. Despite being a basic parameter, there is no theoretical means to determine the shock stand o distance which is accepted universally. Deduction of this quantity depends more or less on experimental or computational means until a successful theoretical model for its predictions is developed. The experimental data available in open literature for spherical bodies in high speed ows mostly lies beyond the 2 km/s regime. Experiments were conducted to determine the shock stand o distance in the velocity range of 1-2 km/s. Three different hemispherical bodies of radii 25, 40 and 50 mm were taken as test models. Since the shock stand o distance is known to depend on the density ratio across the shock and hence gamma (ratio of specific heats), two different test gases, air and carbon dioxide were used for the experiments here. Five different test cases were studied with air as the test gas; Mach 5.56 with Reynolds number of 5.71 million/m and enthalpy of 1.08 MJ/kg, Mach 5.39 with Reynolds number of 3.04 million/m and enthalpy of 1.42 MJ/kg Mach 8.42 with Reynolds number of 1.72 million/m and enthalpy of 1.21 MJ/kg, Mach 11.8 with Reynolds number of 1.09 million/m and enthalpy of 2.03 MJ/kg and Mach 11.25 with Reynolds number of 0.90 million/m and enthalpy of 2.88 MJ/kg. For the experiments conducted with carbon dioxide as test gas, typical freestream conditions were: Mach 6.66 with Reynolds number of 1.46 million/m and enthalpy of 1.23 MJ/kg. The shock stand o distance was determined from the images that were obtained through schlieren photography, the ow visualization technique employed here. The results obtained were found to follow the same trend as the existing experimental data in the higher velocity range. The experimental data obtained was compared with two different theoretical models given by Lobb and Olivier and was found to match. Simulations were carried out in HiFUN, an in-house CFD package for Euler and laminar own conditions for Mach 8 own over 50 mm body with air as the test gas. The computational data was found to match well with the experimental and theoretical data

La Physique des Hautes Densités d’Énergies (HEDP) est caractérisée par desécoulements multi-matériaux fortement compressibles. Le domaine contenant l’écoulementsubit de grandes variations de taille et est le siège d’ondes de chocs et dedétente intenses. La représentation Lagrangienne est bien adaptée à la descriptionde ce type d’écoulements. Elle permet en effet une très bonne description deschocs ainsi qu’un suivit naturel des interfaces multi-matériaux et des surfaces libres.En particulier, les schémas Volumes Finis centrés Lagrangiens GLACE (GodunovtypeLAgrangian scheme Conservative for total Energy) et EUCCLHYD (ExplicitUnstructured Cell-Centered Lagrangian HYDrodynamics) ont prouvé leur efficacitépour la modélisation des équations de la dynamique des gaz ainsi que de l’élastoplasticité.Le travail de cette thèse s’inscrit dans la continuité des travaux de Maireet Nkonga [JCP, 2009] pour la modélisation de l’hydrodynamique et des travauxde Kluth et Després [JCP, 2010] pour l’hyperelasticité. Plus précisément, cettethèse propose le développement de méthodes robustes et précises pour l’extension3D du schéma EUCCLHYD avec une extension d’ordre deux basée sur les méthodesMUSCL (Monotonic Upstream-centered Scheme for Conservation Laws) et GRP(Generalized Riemann Problem). Une attention particulière est portée sur la préservationdes symétries et la monotonie des solutions. La robustesse et la précision duschéma seront validées sur de nombreux cas tests Lagrangiens dont l’extension 3Dest particulièrement difficile. / High Energy Density Physics (HEDP) flows are multi-material flows characterizedby strong shock waves and large changes in the domain shape due to rarefactionwaves. Numerical schemes based on the Lagrangian formalism are good candidatesto model this kind of flows since the computational grid follows the fluid motion.This provides accurate results around the shocks as well as a natural tracking ofmulti-material interfaces and free-surfaces. In particular, cell-centered Finite VolumeLagrangian schemes such as GLACE (Godunov-type LAgrangian scheme Conservativefor total Energy) and EUCCLHYD (Explicit Unstructured Cell-CenteredLagrangian HYDrodynamics) provide good results on both the modeling of gas dynamicsand elastic-plastic equations. The work produced during this PhD thesisis in continuity with the work of Maire and Nkonga [JCP, 2009] for the hydrodynamicpart and the work of Kluth and Després [JCP, 2010] for the hyperelasticitypart. More precisely, the aim of this thesis is to develop robust and accurate methodsfor the 3D extension of the EUCCLHYD scheme with a second-order extensionbased on MUSCL (Monotonic Upstream-centered Scheme for Conservation Laws)and GRP (Generalized Riemann Problem) procedures. A particular care is taken onthe preservation of symmetries and the monotonicity of the solutions. The schemerobustness and accuracy are assessed on numerous Lagrangian test cases for whichthe 3D extensions are very challenging.

Este trabalho apresenta uma estratégia numérica desenvolvida usando o método dos elementos finitos para simular o processo de falha de compósitos cimentícios reforçados com fibras de aço. O material é descrito como um compósito composto por três fases: matriz cimentícia (pasta, argamassa ou concreto), fibras descontínuas discretas, e interface fibra-matriz. Um novo esquema de acoplamento para malhas de elementos finitos não-conformes foi desenvolvido para acoplar as malhas geradas independentes, da matriz cimentícia e de uma nuvem de fibras de aço, baseado na utilização de novos elementos finitos desenvolvidos, denominados elementos finitos de acoplamento. Utilizando este esquema de acoplamento, um procedimento não-rígido é proposto para a modelagem do complexo comportamento não linear da interface fibra-matriz, utilizando um modelo constitutivo de dano apropriado para descrever a relação entre a tensão de cisalhamento (tensão de aderência) e deslizamento relativo entre a matriz e cada fibra de aço individualmente. Este esquema também foi adotado para considerar a presença de barras de aço para as análises de estruturas de concreto armado. As fibras de aço são modeladas usando elementos finitos lineares com dois nós (elementos de treliça) com modelo material elastoplástico. As fibras são posicionadas usando uma distribuição randômica uniforme isotrópica, considerando o efeito parede. Uma abordagem contínua e outra descontínua são investigadas para a modelagem do comportamento frágil da matriz cimentícia. Para a primeira, é utilizado um modelo de dano isotrópico com duas variáveis de dano para descrever o comportamento de dano à tração e à compressão. A segunda emprega uma técnica de fragmentação de malha que utiliza elementos finitos degenerados, posicionados entre todos os elementos finitos que formam a matriz cimentícia. Para esta técnica é proposto um modelo constitutivo à tração, compatível com a abordagem descontínua forte contínua, para prever a propagação de fissura. Para acelerar o cálculo e aumentar a robustez dos modelos de dano contínuos para simular o processamento de falhas, um esquema de integração implícito-explícito é utilizado. Exemplos numéricos são apresentados ao longo do desenvolvimento desta tese. Inicialmente, exemplos numéricos com um único reforço são apresentados para validar a técnica desenvolvida e para investigar à influência das propriedades geométricas 7 das fibras e sua posição em relação à superfície de falha. Posteriormente, exemplos mais complexos são considerados envolvendo uma nuvem de fibras. Nestes casos, atenção especial é dada à influência da distribuição das fibras no comportamento do compósito relacionado ao processo de fissuração. Comparações com resultados experimentais demonstram que a aplicação da ferramenta numérica para modelar o comportamento de compósitos cimentícios reforçados com fibras de aço é muito promissora e pode ser utilizada como uma importante ferramenta para melhor entender os efeitos dos diferentes aspectos envolvidos no processo de falha deste material. / This work presents a numerical strategy developed using the Finite Element Method (FEM) to simulate the failure process of Steel Fiber Reinforced Cementitious Composites (SFRCCs). The material is described as a composite made up by three phases: a cementitious matrix (paste, mortar or concrete), discrete discontinuous fibers, and a fiber-matrix interface. A novel coupling scheme for non-matching finite element meshes has been developed to couple the independent generated meshes of the bulk cementitious matrix and a cloud of discrete discontinuous fibers based on the use of special finite elements developed, termed Coupling Finite Elements (CFEs). Using this approach, a nonrigid coupling procedure is proposed for modeling the complex nonlinear behavior of the fiber-matrix interface by adopting an appropriate constitutive damage model to describe the relation between the shear stress (adherence stress) and the relative sliding between the matrix and each fiber individually. This scheme has also been adopted to account for the presence of regular reinforcing bars in the analysis of reinforced concrete structural elements. The steel fibers are modeled using two-node finite elements (truss elements) with a one-dimensional elastoplastic constitutive model. They are positioned using an isotropic uniform random distribution, considering the wall effect of the mold. Continuous and discontinuous approaches are developed to model the brittle behavior of the bulk cementitious matrix. For the former, an isotropic damage model including two independent scalar damage variables for describing the composite behavior under tension and compression is considered. The discontinuous approach is based on a mesh fragmentation technique that employs degenerated solid finite elements in between all regular (bulk) elements. In this case, a tensile damage constitutive model, compatible with the Continuum Strong Discontinuity Approach (CSDA), is proposed to predict crack propagation. To increase the computability and robustness of the continuum damage models used to simulate the failure processes in both of the strategies, an implicit-explicit integration scheme is used. Numerical analyses are performed throughout the presentation of the work. Initially, numerical examples with a single reinforcement are presented to validate the technique and to investigate the influence of the fibers geometrical properties and its position relative to the crack surface. Then, more complex examples involving a cloud of steel fibers are considered. In these cases, special attention is given to the analysis of the influence of the fiber distribution on the composite behavior relative to the cracking process. Comparisons with experimental results demonstrate that the application of the numerical tool for modeling the behavior of SFRCCs is very promising and may constitute an important tool for better understanding the effects of the different aspects involved in the failure process of this material.

Cette thèse traite des équations de Maxwell en domaine temporel. Le principal objectif est de proposer des méthodes de type éléments finis d'ordre élevé pour les équations de Maxwell et des schémas d'intégration en temps efficaces sur des maillages localement raffinés. Nous considérons des méthodes GDDT (Galerkine Discontinues en Domaine Temporel) s'appuyant sur une interpolation polynomiale d'ordre arbitrairement élevé des composantes du champ électromagnétique. Les méthodes GDDT pour les équations de Maxwell s'appuient le plus souvent sur des schémas d'intégration en temps explicites dont la condition de stabilité peut être très restrictive pour des maillages raffinés. Pour surmonter cette limitation, nous considérons des schémas en temps qui consistent à appliquer un schéma implicite localement, dans les régions raffinées, tout en préservant un schéma explicite sur le reste du maillage. Nous présentons une étude théorique complète et une comparaison de deux méthodes GDDT localement implicites. Des expériences numériques en 2D et 3D illustrent l'utilité des schémas proposés. Le traitement numérique de milieux de propagation complexes est également l'un des objectifs. Nous considérons l'interaction des ondes électromagnétiques avec les tissus biologiques qui est au cœur de nombreuses applications dans le domaine biomédical. La modélisation numérique nécessite alors de résoudre le système de Maxwell avec des modèles appropriés de dispersion. Nous formulons une méthode GDDT localement implicite pour le modèle de Debye et proposons une analyse théorique et numérique complète du schéma. / This work deals with the time-domain formulation of Maxwell's equations. The main objective is to propose high-order finite element type methods for the discretization of Maxwell's equations and efficient time integration methods on locally refined meshes. We consider Discontinuous Galerkin Time-Domain (DGTD) methods relying on an arbitrary high-order polynomial interpolation of the components of the electromagnetic field. Existing DGTD methods for Maxwell's equations often rely on explicit time integration schemes and are constrained by a stability condition that can be very restrictive on highly refined meshes. To overcome this limitation, we consider time integration schemes that consist in applying an implicit scheme locally i.e. in the refined regions of the mesh, while preserving an explicit scheme in the complementary part. We present a full theoretical study and a comparison of two locally implicit DGTD methods. Numerical experiments for 2D and 3D problems illustrate the usefulness of the proposed time integration schemes. The numerical treatment of complex propagation media is also one of the objectives. We consider the interaction of electromagnetic waves with biological tissues that is of interest to applications in biomedical domain. Numerical modeling then requires to solve the system of Maxwell's equations coupled to appropriate models of physical dispersion. We derive a locally implicit DGTD method for the Debye model and we achieve a full theoretical and numerical analysis of the resulting scheme.

Cette thèse porte sur la modélisation et la discrétisation d’écoulements Darcy dans les milieux poreux fracturés. Nous suivons l’approche des modèles, dits à dimensions hybrides, qui représentent les réseaux de fractures comme des surfaces de codimension 1 immergées dans la matrice. Les modèles considérés prennent en compte les interactions entre matrice et fractures et permettent de traiter des fractures agissant comme conduites ou comme barrières, ce que nécessite de prendre en compte les sauts de pression aux interfaces matrice-fracture. Dans le cas des écoulements diphasiques, nous proposons des modèles, qui prennent en compte les sauts de saturations aux interfaces matrice-fracture, dû à la capillarité. L’analyse numérique est menée dans le cadre général de la méthode de discrétisations gradients, qui est étendue aux modèles considérés. Deux familles de schémas numériques, le schéma Vertex Approximate Gradient et le schéma Volumes Finis Hybrides sont adaptées aux modèles à dimensions hybrides. On prouve via des résultats de densité que ce sont des schémas gradients, pour lesquels la convergence est établie. En diphasique, l’existence d’une solution est obtenue en passant. Plusieurs cas tests sont présentés. En monophasique, on observe la convergence sur des différents types de mailles pour une famille de solutions dans un milieux fracturé hétérogène et anisotrope. En diphasique, nous présentons une série de cas tests afin de comparer les modèles à dimensions hybrides au modèle de référence, dans lequel les fractures ont la même dimension que la matrice. A part quantifier le gain en performance de calcul, ces tests montrent la qualité des différents modèles réduits. / This thesis investigates the modelling of Darcy flow through fractured porous media and its discretization on general polyhedral meshes. We follow the approach of hybrid dimensional models, invoking a complex network of planar fractures. The models account for matrix-fracture interactions and fractures acting either as drains or as barriers, i.e. we have to deal with pressure discontinuities at matrix-fracture interfaces. In the case of two phase flow, we present two models, which permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The numerical analysis is performed in the general framework of the Gradient Discretisation Method, which is extended to the models under consideration. Two families of schemes namely the Vertex Approximate Gradient scheme (VAG) and the Hybrid Finite Volume scheme (HFV) are detailed and shown to fit in the gradient scheme framework, which yields, in particular, convergence. For single phase flow, we obtain convergence of order 1 via density results. For two phase flow, the existence of a solution is obtained as a byproduct of the convergence analysis. Several test cases are presented. For single phase flow, we study the convergence on different types of meshes for a family of solutions. For two phase flow, we compare the hybrid-dimensional models to the reference equidimensional model, in which fractures have the same dimension as the matrix. This does not only provide quantitative evidence about computational gain, but also leads to deep insight about the quality of the proposed reduced models.

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.