No description available.

No description available.

Ill-posed inverse problems have wide applications in many fields such as oceanography, signal processing, machine learning, biomedical imaging, remote sensing, geophysics, and others. In this dissertation, we address the problem of solving unstable operator equations with iteratively regularized Newton-type algorithms. Important practical questions such as selection of regularization parameters, construction of generating (filtering) functions based on a priori information available for different models, algorithms for stopping rules and error estimates are investigated with equal attention given to theoretical study and numerical experiments.

No description available.

Adaptive methods in finite element analysis are essential tools in the efficient computation and error control of problems that may exhibit singularities. In this paper, we consider solving a boundary value problem which exhibits a singularity at the origin due to both the structure of the domain and the regularity of the exact solution. We introduce a hybrid mixed finite element method using Lagrange Multipliers to initially solve the partial differential equation for the both the flux and displacement. An a posteriori error estimate is then applied both locally and globally to approximate the error in the computed flux with that of the exact flux. Local estimation is the key tool in identifying where the mesh should be refined so that the error in the computed flux is controlled while maintaining efficiency in computation. Finally, we introduce a simple refinement process in order to improve the accuracy in the computed solutions. Numerical experiments are conducted to support the advantages of mesh refinement over a fixed uniform mesh.

NR 20140805

O trabalho aborda a questão das estimativas a posteriori dos erros de discretização e particularmente a recuperação dos gradientes de soluções numéricas obtidas com o método dos elementos finitos (MEF) e com o método dos elementos finitos generalizados (MEFG). Inicialmente, apresenta-se, em relação ao MEF, um resumido estado da arte e conceitos fundamentais sobre este tema. Em seguida, descrevem-se os estimadores propostos para o MEF denominados Estimador Z e \"Superconvergent Patch Recovery\" (SPR). No âmbito do MEF propõe-se de modo original a incorporação do \"Singular Value Decomposition\" (SVD) ao SPR aqui mencionada como SPR Modificado. Já no contexto do MEFG, apresenta-se um novo estimador do erro intitulado EPMEFG, estendendo-se para aquele método as idéias do SPR Modificado. No EPMEFG, a função polinomial local que permite recuperar os valores nodais dos gradientes da solução tem por suporte nuvens (conjunto de elementos finitos que dividem um nó comum) e resulta da aplicação de um critério de aproximação por mínimos quadrados em relação aos pontos de superconvergência. O número destes pontos é definido a partir de uma análise em cada elemento que compõe a nuvem, considerando-se o grau da aproximação local do campo de deslocamentos enriquecidos. Exemplos numéricos elaborados com elementos lineares triangulares e quadrilaterais são resolvidos com o Estimador Z, o SPR Modificado e o EPMEFG para avaliar a eficiência de cada estimador. Essa avaliação é realizada mediante o cálculo dos índices de efetividade. / The paper addresses the issue of a posteriori estimates of discretization errors and particularly the recovery of gradients of numerical solutions obtained with the finite element method (FEM) and the generalized finite element method (GFEM). Initially, it is presented, for the MEF, a brief state of the art and fundamental concepts about this topic. Next, it is described the proposed estimators for the FEM called Z-Estimator and Superconvergent Patch Recovery (SPR). It is proposed, originally, in the ambit of the FEM, the incorporation of the \"Singular Value Decomposition (SVD) to SPR mentioned here as Modified SPR. On the other hand, in the context of GFEM, it is presented a new error estimator entitled EPMEFG in order to expand the ideas of Modified SPR to that method. In EPMEFG, the local polynomial function that allows to recover the nodal values of the gradients of the solution has for support clouds (set of finite elements that share a common node) and results from the applying of a criterion of least squares approximation in relation to the superconvergent points. The number of these points is defined from an analysis of each cloud\'s element, considering the degree of local approximation of the displacement field enriched. Numerical examples elaborated with linear triangular and quadrilateral elements are solved with the Z-Estimator, the Modified SPR and the EPMEFG to evaluate the efficiency of each estimator. This evaluation is done calculating the effectivity indexes.

L'objectif de la thèse est l'étude mathématique et l'analyse numérique du problème non linéaire de Sobolev. Un premier chapitre est consacré à l'analyse a priori pour le problème de Sobolev où on utilise des méthodes de semi-discrétisation explicite en temps. Des estimations d'erreurs ont été obtenues assurant que les schémas numériques utilisés convergent lorsque le pas de discrétisation en temps et le pas de discrétisation en espace tendent vers zéro. Dans le second chapitre, on s'intéresse au problème de Sobolev singulièrement perturbé. En vue de la stabilité des schémas numériques, on utilise dans cette partie des méthodes numériques implicites (la méthode d'Euler et la méthode de Crank- Nicolson) pour discrétiser le problème par rapport au temps. Dans le troisième chapitre, on présente des applications et des illustrations où on utilise le logiciel "FreeFem++". Dans le dernier chapitre, on considère une équation de type Sobolev et on s'intéresse à la dérivation d'estimations d'erreur a posteriori pour la discrétisation de cette équation par la méthode des éléments finis conforme en espace et un schéma d'Euler implicite en temps. La borne supérieure est globale en espace et en temps et permet le contrôle effectif de l'erreur globale. A la fin du chapitre, on propose un algorithme adaptatif qui permet d'atteindre une précision relative fixée par l'utilisateur en raffinant les maillages adaptativement et en équilibrant les contributions en espace et en temps de l'erreur. On présente également des essais numériques. / The purpose of this work is the mathematical study and the numerical analysis of the nonlinear Sobolev problem. A first chapter is devoted to the a priori analysis for the Sobolev problem, where we use an explicit semidiscretization in time. A priori error estimates were obtained ensuring that the used numerical schemes converge when the time step discretization and the spatial step discretization tend to zero. In a second chapter, we are interested in the singularly perturbed Sobolev problem. For the stability of numerical schemes, we used in this part implicit semidiscretizations in time (the Euler method and the Crank-Nicolson method). Our estimates of Chapters 1 and 2 are confirmed in the third chapter by some numerical experiments. In the last chapter, we consider a Sobolev equation and we derive a posteriori error estimates for the discretization of this equation by a conforming finite element method in space and an implicit Euler scheme in time. The upper bound is global in space and time and allows effective control of the global error. At the end of the chapter, we propose an adaptive algorithm which ensures the control of the total error with respect to a user-defined relative precision by refining the meshes adaptively, equilibrating the time and space contributions of the error. We also present numerical experiments.

O trabalho aborda a questão das estimativas a posteriori dos erros de discretização e particularmente a recuperação dos gradientes de soluções numéricas obtidas com o método dos elementos finitos (MEF) e com o método dos elementos finitos generalizados (MEFG). Inicialmente, apresenta-se, em relação ao MEF, um resumido estado da arte e conceitos fundamentais sobre este tema. Em seguida, descrevem-se os estimadores propostos para o MEF denominados Estimador Z e \"Superconvergent Patch Recovery\" (SPR). No âmbito do MEF propõe-se de modo original a incorporação do \"Singular Value Decomposition\" (SVD) ao SPR aqui mencionada como SPR Modificado. Já no contexto do MEFG, apresenta-se um novo estimador do erro intitulado EPMEFG, estendendo-se para aquele método as idéias do SPR Modificado. No EPMEFG, a função polinomial local que permite recuperar os valores nodais dos gradientes da solução tem por suporte nuvens (conjunto de elementos finitos que dividem um nó comum) e resulta da aplicação de um critério de aproximação por mínimos quadrados em relação aos pontos de superconvergência. O número destes pontos é definido a partir de uma análise em cada elemento que compõe a nuvem, considerando-se o grau da aproximação local do campo de deslocamentos enriquecidos. Exemplos numéricos elaborados com elementos lineares triangulares e quadrilaterais são resolvidos com o Estimador Z, o SPR Modificado e o EPMEFG para avaliar a eficiência de cada estimador. Essa avaliação é realizada mediante o cálculo dos índices de efetividade. / The paper addresses the issue of a posteriori estimates of discretization errors and particularly the recovery of gradients of numerical solutions obtained with the finite element method (FEM) and the generalized finite element method (GFEM). Initially, it is presented, for the MEF, a brief state of the art and fundamental concepts about this topic. Next, it is described the proposed estimators for the FEM called Z-Estimator and Superconvergent Patch Recovery (SPR). It is proposed, originally, in the ambit of the FEM, the incorporation of the \"Singular Value Decomposition (SVD) to SPR mentioned here as Modified SPR. On the other hand, in the context of GFEM, it is presented a new error estimator entitled EPMEFG in order to expand the ideas of Modified SPR to that method. In EPMEFG, the local polynomial function that allows to recover the nodal values of the gradients of the solution has for support clouds (set of finite elements that share a common node) and results from the applying of a criterion of least squares approximation in relation to the superconvergent points. The number of these points is defined from an analysis of each cloud\'s element, considering the degree of local approximation of the displacement field enriched. Numerical examples elaborated with linear triangular and quadrilateral elements are solved with the Z-Estimator, the Modified SPR and the EPMEFG to evaluate the efficiency of each estimator. This evaluation is done calculating the effectivity indexes.

The presented work deals with the discontinuous Galerkin method with the anisotropic mesh adaptation for stationary convection-diffusion problems. Basic definitions are included in an introduction where we also present the used method. The following parts describe various methods for evaluating a Riemann metric, which is necessary for anisotropic mesh adaptation. The most important part of work follows - numerical experiments carried out with ADGFEM and ANGENER software packages. In these experiments, we compare different approaches for the definition of Riemann metrics and compare their efficiency. The main output of this thesis are subroutines for evaluation of the Riemann metric including its source code.