A singularly perturbed convection-diffusion problem in two and three space dimensions is discretized using the streamline upwind Petrov Galerkin (SUPG) variant of the finite element method. The dominant convection frequently gives rise to solutions with layers; hence anisotropic finite elements can be applied advantageously. The main focus is on a posteriori energy norm error estimation that is robust in the perturbation parameter and with respect to the mesh anisotropy. A residual error estimator and a local problem error estimator are proposed and investigated. The analysis reveals that the upper error bound depends on the alignment of the anisotropies of the mesh and of the solution. Hence reliable error estimation is possible for suitable anisotropic meshes. The lower error bound depends on the problem data via a local mesh Peclet number. Thus efficient error estimation is achieved for small mesh Peclet numbers. Altogether, error estimation approaches for isotropic meshes are successfully extended to anisotropic elements. Several numerical experiments support the analysis.

O Método dos Elementos Finitos Generalizados, MEFG, compartilha importantes características dos métodos sem malha. As funções de aproximação do MEFG, atreladas aos pontos nodais, são enriquecidas de modo análogo ao refinamento p realizado no Método das Nuvens hp. Por outro lado, por empregar uma malha de elementos para construir as funções partição da unidade, ele também pode ser entendido como uma forma não convencional do Método dos Elementos Finitos. Neste trabalho, ambas as interpretações são consideradas. Os métodos sem malha, particularmente o Método de Galerkin Livre de Elementos e o Método das Nuvens hp, são introduzidos com o propósito de estabelecer os conceitos fundamentais para a descrição do MEFG. Na seqüência, apresentam-se aplicações numéricas em análise linear e evidenciam-se características que tornam o MEFG interessante para a simulação da propagação de descontinuidades. Após discutir os modelos de dano adotados para representar o comportamento não-linear do material, são introduzidos exemplos de aplicação, inicialmente do Método das Nuvens hp e depois do MEFG, na análise de estruturas de concreto. Os resultados obtidos servem de argumento para a implementação de um procedimento p-adaptativo, particularmente com o MEFG. Propõe-se, então a adaptação do Método dos Resíduos em Elementos Equilibrados à formulação do MEFG. Com vistas ao seu emprego em problemas não-lineares, algumas modificações são introduzidas à formulação do estimador. Mostra-se que a medida obtida para representar o erro, apesar de fundamentada em diversas hipóteses nem sempre possíveis de serem satisfeitas, ainda assim viabiliza a análise não-linear p-adaptativa. Ao final, são enumeradas propostas para a aplicação do MEFG em problemas caracterizados pela propagação de defeitos / The Generalized Finite Element Method, GFEM, shares several features with the so called meshless methods. The approximation functions used in the GFEM are associated with nodal points like in meshless methods. In addition, the enrichment of the approximation spaces can be done in the same fashion as in the meshless hp-Cloud method. On the other hand, the partition of unity used in the GFEM is provided by Lagrangian finite element shape functions. Therefore, this method can also be understood as a variation of the Finite Element Method. Indeed, both interpretations of the GFEM are valid and give unique insights into the method. The meshless character of the GFEM justified the investigation of meshless methods in this work. Among them, the Element Free Galerkin Method and the hp-Cloud Method are described aiming to introduce key concepts of the GFEM formulation. Following that, several linear problems are solved using these three methods. Such linear analysis demonstrates several features of the GFEM and its suitability to simulate propagating discontinuities. Next, damage models employed to model the nonlinear behavior of concrete structures are discussed and numerical analysis using the hp-Cloud Method and the GFEM are presented. The results motivate the implementation of a p-adaptive procedure tailored to the GFEM. The technique adopted is the Equilibrated Element Residual Method. The estimator is modified to take into account nonlinear peculiarities of the problems considered. The hypotheses assumed in the definition of the error measure are sometimes violated. Nonetheless, it is shown that the proposed error indicator is effective for the class of p-adaptive nonlinear analysis investigated. Finally, several suggestions are enumerated considering future applications of the GFEM, specially for the simulation of damage and crack propagation

The vector finite element method has gained great attention since overcoming the deficiencies incurred by the scalar basis functions for the vector Helmholtz equation. Most implementations of vector FEM have been non-adaptive, where a mesh of the domain is generated entirely in advance and used with a constant degree polynomial basis to assign the degrees of freedom. To reduce the dependency on the users' expertise in analyzing problems with complicated boundary structures and material characteristics, and to speed up the FEM tool, the demand for adaptive FEM grows high. For efficient adaptive FEM, error estimators play an important role in assigning additional degrees of freedom. In this proposal study, hierarchical vector basis functions and four error estimators for p-refinement are investigated for electromagnetic applications.

O Método dos Elementos Finitos Generalizados, MEFG, compartilha importantes características dos métodos sem malha. As funções de aproximação do MEFG, atreladas aos pontos nodais, são enriquecidas de modo análogo ao refinamento p realizado no Método das Nuvens hp. Por outro lado, por empregar uma malha de elementos para construir as funções partição da unidade, ele também pode ser entendido como uma forma não convencional do Método dos Elementos Finitos. Neste trabalho, ambas as interpretações são consideradas. Os métodos sem malha, particularmente o Método de Galerkin Livre de Elementos e o Método das Nuvens hp, são introduzidos com o propósito de estabelecer os conceitos fundamentais para a descrição do MEFG. Na seqüência, apresentam-se aplicações numéricas em análise linear e evidenciam-se características que tornam o MEFG interessante para a simulação da propagação de descontinuidades. Após discutir os modelos de dano adotados para representar o comportamento não-linear do material, são introduzidos exemplos de aplicação, inicialmente do Método das Nuvens hp e depois do MEFG, na análise de estruturas de concreto. Os resultados obtidos servem de argumento para a implementação de um procedimento p-adaptativo, particularmente com o MEFG. Propõe-se, então a adaptação do Método dos Resíduos em Elementos Equilibrados à formulação do MEFG. Com vistas ao seu emprego em problemas não-lineares, algumas modificações são introduzidas à formulação do estimador. Mostra-se que a medida obtida para representar o erro, apesar de fundamentada em diversas hipóteses nem sempre possíveis de serem satisfeitas, ainda assim viabiliza a análise não-linear p-adaptativa. Ao final, são enumeradas propostas para a aplicação do MEFG em problemas caracterizados pela propagação de defeitos / The Generalized Finite Element Method, GFEM, shares several features with the so called meshless methods. The approximation functions used in the GFEM are associated with nodal points like in meshless methods. In addition, the enrichment of the approximation spaces can be done in the same fashion as in the meshless hp-Cloud method. On the other hand, the partition of unity used in the GFEM is provided by Lagrangian finite element shape functions. Therefore, this method can also be understood as a variation of the Finite Element Method. Indeed, both interpretations of the GFEM are valid and give unique insights into the method. The meshless character of the GFEM justified the investigation of meshless methods in this work. Among them, the Element Free Galerkin Method and the hp-Cloud Method are described aiming to introduce key concepts of the GFEM formulation. Following that, several linear problems are solved using these three methods. Such linear analysis demonstrates several features of the GFEM and its suitability to simulate propagating discontinuities. Next, damage models employed to model the nonlinear behavior of concrete structures are discussed and numerical analysis using the hp-Cloud Method and the GFEM are presented. The results motivate the implementation of a p-adaptive procedure tailored to the GFEM. The technique adopted is the Equilibrated Element Residual Method. The estimator is modified to take into account nonlinear peculiarities of the problems considered. The hypotheses assumed in the definition of the error measure are sometimes violated. Nonetheless, it is shown that the proposed error indicator is effective for the class of p-adaptive nonlinear analysis investigated. Finally, several suggestions are enumerated considering future applications of the GFEM, specially for the simulation of damage and crack propagation

Some boundary value problems yield anisotropic solutions, e.g. solutions with boundary layers. If such problems are to be solved with the finite element method (FEM), anisotropically refined meshes can be advantageous. In order to construct these meshes or to control the error one aims at reliable error estimators. For isotropic meshes such estimators are known but they fail when applied to anisotropic meshes. Rectangular (or cuboidal) anisotropic meshes were already investigated. In this paper an error estimator is presented for tetrahedral or triangular meshes which offer a much greater geometrical flexibility.

We consider a singularly perturbed reaction-diffusion problem and derive and rigorously analyse an a posteriori residual error estimator that can be applied to anisotropic finite element meshes. The quotient of the upper and lower error bounds is the so-called matching function which depends on the anisotropy (of the mesh and the solution) but not on the small perturbation parameter. This matching function measures how well the anisotropic finite element mesh corresponds to the anisotropic problem. Provided this correspondence is sufficiently good, the matching function is O(1). Hence one obtains tight error bounds, i.e. the error estimator is reliable and efficient as well as robust with respect to the small perturbation parameter. A numerical example supports the anisotropic error analysis.

We consider a posteriori error estimators that can be applied to anisotropic tetrahedral finite element meshes, i.e. meshes where the aspect ratio of the elements can be arbitrarily large. Two kinds of Zienkiewicz-Zhu (ZZ) type error estimators are derived which are both based on some recovered gradient. Two different, rigorous analytical approaches yield the equivalence of both ZZ error estimators to a known residual error estimator. Thus reliability and efficiency of the ZZ error estimation is obtained. Particular attention is paid to the requirements on the anisotropic mesh. The analysis is complemented and confirmed by several numerical examples.

Le but de cette étude est d'améliorer les performances, en termes d'espace mémoire et de temps de calcul, des simulations actuelles de l'Interaction mécanique Pastille-Gaine (IPG), phénomène complexe pouvant avoir lieu lors de fortes montées en puissance dans les réacteurs à eau sous pression. Parmi les méthodes de raffinement de maillage, méthodes permettant de simuler efficacement des singularités locales, une approche multi-grille locale a été choisie car elle présente l'intérêt de pouvoir utiliser le solveur en boîte noire tout en ayant un faible nombre de degrés de liberté à traiter par niveau. La méthode Local Defect Correction (LDC), adaptée à une discrétisation de type éléments finis, a tout d'abord été analysée et vérifiée en élasticité linéaire, sur des configurations issues de l'IPG, car son utilisation en mécanique des solides est peu répandue. Différentes stratégies concernant la mise en oeuvre pratique de l'algorithme multi-niveaux ont également été comparées. La combinaison de la méthode LDC et de l'estimateur d'erreur a posteriori de Zienkiewicz-Zhu, permettant d'automatiser la détection des zones à raffiner, a ensuite été testée. Les performances obtenues sur des cas bidimensionnels et tridimensionnels sont très satisfaisantes, l'algorithme proposé se montrant plus performant que des méthodes de raffinement h-adaptatives. Enfin, l'algorithme a été étendu à des problèmes mécaniques non linéaires. Les questions d'un raffinement espace/temps mais aussi de la transmission des conditions initiales lors du remaillage ont entre autres été abordées. Les premiers résultats obtenus sont encourageants et démontrent l'intérêt de la méthode LDC pour des calculs d'IPG. / The aim of this study is to improve the performances, in terms of memory space and computational time, of the current modelling of the Pellet-Cladding mechanical Interaction (PCI),complex phenomenon which may occurs during high power rises in pressurised water reactors. Among the mesh refinement methods - methods dedicated to efficiently treat local singularities - a local multi-grid approach was selected because it enables the use of a black-box solver while dealing few degrees of freedom at each level. The Local Defect Correction (LDC) method, well suited to a finite element discretisation, was first analysed and checked in linear elasticity, on configurations resulting from the PCI, since its use in solid mechanics is little widespread. Various strategies concerning the implementation of the multilevel algorithm were also compared. Coupling the LDC method with the Zienkiewicz-Zhu a posteriori error estimator in orderto automatically detect the zones to be refined, was then tested. Performances obtained on two-dimensional and three-dimensional cases are very satisfactory, since the algorithm proposed is more efficient than h-adaptive refinement methods. Lastly, the LDC algorithm was extended to nonlinear mechanics. Space/time refinement as well as transmission of the initial conditions during the remeshing step were looked at. The first results obtained are encouraging and show the interest of using the LDC method for PCI modelling.

Vários modelos computacionais que implementam o transporte de soluto em meio poroso saturado surgem constantemente em publicações científicas devido à suma importância dada à compreensão e previsão do transporte de constituintes dissolvidos em água subterrânea. As soluções numéricas obtidas por esquemas computacionais não estão imunes aos erros de discretização. No entanto, a confiabilidade nos resultados obtidos das complexas operações provenientes da dinâmica de fluidos computacional pode ser aumentada através de estimadores de erro a posteriori que indicam a precisão da solução numérica de um modelo matemático que simula o fenômeno físico de interesse. Neste trabalho é apresentado um estimador residual para a equação parabólica que descreve os fenômenos de advecção-dispersão-reação (ADR) em meio poroso saturado, considerando o transporte em regime de pequena advecção. A solução numérica da equação ADR é obtida pelo método dos elementos finitos que emprega termos upwind para minimizar as inconvenientes oscilações espúrias. A implementação do código computacional para obter essa solução numérica e o seu correspondente erro a posteriori, é feita em linguagem JAVA na plataforma Eclipse seguindo o paradigma da Programação Orientada a Objetos (POO). A solução numérica da equação elíptica do fluxo subterrâneo e o seu estimador de erro com características de recuperação do gradiente, o estimador ZZ, também são implementados no código JAVA. Assim, a solução da equação do transporte é obtida em função da reusabilidade POO prevista na implementação da equação do fluxo. A comparação da solução numérica do modelo ADR 2D com a correspondente solução analítica disponível na literatura, demonstra que o estimador residual apresenta excelentes índices de eficiência. Os resultados numéricos obtidos mostraram que o estimador residual encontra-se limitado inferior e superiormente pelo erro real da solução em malha grosseira. O estimador ZZ mostrou-se inadequado para a análise do erro de aproximação das equações ADR. Os exemplos selecionados para verificação e aplicação do estimador residual abrangem, em diferentes escalas, modelos que descrevem reação de primeira ordem e modelos com fenômenos de sorção e retardamento na migração do contaminante em meio poroso saturado. Em conseqüência, o estimador residual proposto provou ser computável, eficiente e robusto no sentido de abranger uma grande variedade das aplicações dos fenômenos de transporte de contaminantes em meio poroso saturado e regime de pequena advecção. / Several computational models that implement the solute migration in saturated porous media constantly appear in scientific publications due to the great importance given to the understanding and forecast of the solute transport in groundwater. The numerical solutions obtained by computational schemes are not immune to errors related to the discretization process. However, the reliability of the results obtained by the complex operations of the computational fluids dynamics can be enhanced by a posteriori error estimates that indicate the accuracy of the numerical solution. In this work a residual error estimator is presented for the parabolic equation that describes the advection-dispersion-reaction phenomena (ADR) in saturated porous media, considering the transport in small advection regime. The numerical solution of the ADR equation is obtained by the finite element method using upwind terms to minimize the spurious oscillations. The computational code and the correspondent a posteriori error estimates are implemented in Java language following the Object Oriented Programming (OOP) paradigm in Eclipse platform. The numerical solution of the elliptic groundwater flow equation and the respective error estimates with gradient recovery characteristic, the ZZ-estimator, are also implemented in the JAVA code. The solution of the transport equation is obtained as a consequence of the OOP reusability intended in the implementation of the flow equation. The numerical solution of the ADR 2D simulation compared to the analytical solution available in the literature, demonstrate the excellent effectivity index presented by the residual error estimator. The obtained results indicate that the residual error estimator is lower and upper bounded by a solution in coarse mesh. The ZZ-estimator showed to be inadequate for the error analysis of the ADR equations. The examples selected for validation and application of the residual estimator include, in distinct scales, models that describe reaction of first order and models with sorption and retardation phenomena in the pollutant migration in saturated porous media. Therefore, the proposed residual error estimator proved to be computable, efficient and robust in the sense of solving a great variety of applications of transport phenomena in saturated porous media at small advection regime.

Les travaux de la présente thèse portent sur l’estimation d'erreur a posteriori pour les solutions numériques par éléments finis de l'équation des ondes élastiques dans les milieux hétérogènes. Deux types d’estimation ont été développés. Le premier considère directement l’équation élastodynamique et conduit à un nouvel estimateur d'erreur a posteriori explicite en norme L∞ en temps. Les principales caractéristiques de cet estimateur explicite sont l'utilisation de la méthode de résidus et le développement de reconstructions en temps et en espace selon les différentes régularités exigées par les différents termes contribuant à l’obtention d’une borne supérieure. L’analyse numérique de cet estimateur dans le cas des maillages uniformes montre qu’il assure bien une borne supérieure mais avec une propriété asymptotique qui reste à améliorer. Le deuxième type d’estimateur d’erreur est développé dans le contexte de la propagation des ondes à haute fréquence dans des milieux hétérogènes à l’échelle mésoscopique. Il s’agit d’une nouvelle erreur en résidus basée sur l'équation de transfert radiatif, qui est obtenue par un développement asymptotique multi-échelle de l'équation d'onde en utilisant la transformation de Wigner en espace-temps. Les résidus sont exprimés en termes de densités énergétiques calculés dans l’espace des phases pour les solutions d’onde numériques transitoires par éléments finis. L’analyse numérique de cette erreur appliquée aux milieux homogènes et hétérogènes en 1D a permis de valider notre approche. Les champs d’application visés sont la propagation des ondes sismiques dans les milieux géophysiques ou la propagation des ondes ultrasonores dans les milieux polycristallins. / This thesis work deals with a posteriori error estimates for finite element solutions of the elastic wave equation in heterogeneous media. Two different a posteriori estimation approaches are developed. The first one, in a classical way, considers directly the elastodynamic equation and results in a new explicit error estimator in a non-natural L∞ norm in time. Its key features are the use of the residual method and the development of space and time reconstructions with respect to regularities required by different residual operators contributing to the proposed error bound. Numerical applications of the error bound with different mesh sizes show that it gives rise to a fully computable upper bound. However, its effectivity index and its asymptotic accuracy remain to be improved. The second error estimator is derived for high frequency wave propagation problem in heterogeneous media in the weak coupling regime. It is a new residual-type error based on the radiative transfer equation, which is derived by a multi-scale asymptotic expansion of the wave equation in terms of the spatio-temporal Wigner transforms of wave fields. The residual errors are in terms of angularly resolved energy quantities of numerical solutions of waves by finite element method. Numerical calculations of the defined errors in 1D homogeneous and heterogeneous media allow validating the proposed error estimation approach. The application field of this work is the numerical modelling of the seismic wave propagation in geophysical media or the ultrasonic wave propagation in polycrystalline materials.