Estuda-se o problema de equilíbrio sem força de corpo de uma esfera anisotrópica sob compressão radial uniformemente distribuída sobre o seu contorno no contexto da teoria da elasticidade linear clássica. A solução deste problema prediz o fenômeno inaceitável da auto-intersecção em uma região próxima ao centro da esfera para uma dada faixa de parâmetros materiais. Sob o contexto de uma teoria de minimização do funcional de energia potencial total da elasticidade linear clássica com a restrição de que o determinante do gradiente da função mudança de configuração seja injetivo, este fenômeno é eliminado. Aplicam-se duas formulações do Método dos Elementos Finitos de Galerkin Descontínuo (MEFGD) para obter soluções aproximadas para o problema de equilíbrio da esfera sem restrição. A primeira formulação do MEFGD aproxima diretamente os campos de deslocamento e deformação infinitesimal. A consideração do campo adicional de deformação na formulação do MEFGD aumenta o número de graus de liberdade associados aos nós da malha de elementos finitos e, consequentemente, o custo computacional. Com o objetivo de reduzir o número de graus de liberdade, introduz-se neste trabalho uma formulação alternativa do MEFGD. Nesta formulação, o campo de deformação infinitesimal não é obtido diretamente da inversão do sistema de equações resultante, mas sim por pós-processamento, a partir do campo de deslocamento aproximado. As soluções aproximadas obtidas com ambas as formulações do MEFGD são comparadas com a solução exata do problema sem restrição e com soluções aproximadas obtidas com o Método dos Elementos Finitos de Galerkin Clássico (MEFGC). Ambas as formulações do MEFGD fornecem melhores aproximações para a solução exata do que as aproximações obtidas com o MEFGC. Os erros entre a solução exata e as soluções aproximadas obtidas com a formulação alternativa do MEFGD são um pouco maiores do que os erros correspondentes obtidos com a formulação original do MEFGD. Este aumento nos erros é compensado pelo menor esforço computacional exigido pela formulação alternativa. Este trabalho serve de base para o estudo de problemas com restrição de injetividade utilizando o método de Galerkin descontínuo. / The equilibrium problem without body force of an anisotropic sphere under radial compression that is uniformly distributed on the sphere\'s boundary is investigated in the context of the classical linear elasticity theory. The solution of this problem predicts the unacceptable phenomenon of self-intersection in a vicinity of the center of the sphere for a given range of material parameters. This phenomenon can be eliminated in the context of a theory that minimizes the total potential energy of classical linear elasticity subjected to the restriction that the deformation field be injective. Two formulations of the Finite Element Method using Discontinuous Galerkin (MEFGD) are used to obtain approximate solutions for the unconstrained problem. The first formulation of the MEFGD approximates both the displacement and the strain fields. The consideration of the strain as an additional field in the formulation of the MEFGD increases the number of degrees of freedom associated to the finite elements and, therefore, the computational cost. With the objective of reducing the number of degrees of freedom, an alternative formulation of the MEFGD is introduced in this work. In this formulation, the strain field is not obtained directly from the inversion of the resulting linear system of equations, but from a post-processing calculation using the approximate displacement field. The approximate solutions obtained with both formulations of the MEFGD are compared with the exact solution of the problem without restriction and with approximate solutions obtained with the Finite Element Method using Classical Galerkin (MEFGC). Both formulations of the MEFGD yield better approximations for the exact solution than the approximations obtained with the MEFGC. The errors between the exact solution and the approximate solutions obtained with the alternative formulation of the MEFGD are slightly higher than the corresponding errors obtained with the original formulation of the MEFGD. These errors are compensated by the fact that the alternative formulation requires less computational effort than the computational effort required by the original formulation. This work serves as a basis for the study of problems with the injectivity restriction using the discontinuous Galerkin method.

Orientador: João Frederico da Costa Azevedo Meyer / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica / Made available in DSpace on 2018-08-21T04:55:31Z (GMT). No. of bitstreams: 1 CarmonaTabares_PauloCesar_D.pdf: 24565019 bytes, checksum: 8ebe9aed1d258a0712f49e9711f8d107 (MD5) Previous issue date: 2012 / Resumo: Neste trabalho, construímos um modelo matemático para avaliar as conjecturas existentes acerca do impacto que tem o material inorgânico particulado (sedimento) nas populações bentônicas predominantes na Enseada Potter. Na construção do modelo são utilizadas informações do fenômeno, proporcionadas pelas pesquisas permanentes na região de estudo. Como resultado, logramos comprovar mediante simulações numéricas, o efeito que produz o sedimento na distribuição e abundância das espécies do substrato marinho, constatando neste ecossistema particular as consequências do aquecimento global nessa parte da região antártica. A modelagem é feita com um sistema de equações diferenciais parciais não- lineares sobre um domínio bidimensional irregular (descritiva da região original), o qual é discretizado nas variáveis espaciais por elementos finitos de primeira ordem e na variável temporal pelo Método de Crank-Nicolson. A resolução do sistema não-linear resultante é aproximada através de um método preditor-corretor cuja solução aproximada é visualizada e valorada qualitativamente usando gráficos evolutivos obtidos por simulações em ambiente MATLAB / Abstract: In this work, we built a mathematical model to evaluate existing conjectures about the impact that inorganic particulate material (sediment) has upon predominating benthic populations in Potter Cove. For the mathematical model, phenomena information was that provided by permanent researches in the study area. As a result, by means of numerical simulations, we were able to confirm the effect of sediment over distribution and abundance for species of marine substrate, verifying in this particular ecosystem, the effects of global warming in this specific Antarctic region. Modeling is done with a system of nonlinear partial differential equations over an irregular two-dimensional domain (descriptive of the original region), which is discretized in the spatial variables by first order finite elements and in the time variable by Crank-Nicolson. The resolution of the resulting nonlinear system is approximated by a predictor-corrector method and the solution is displayed and qualitatively valorized using evolutive graphics, obtain in a MATLAB environment / Doutorado / Matematica Aplicada / Doutor em Matemática Aplicada

Dans ce manuscrit, nous nous inéressons à l'étude du mouvement des vagues soumises uniquement à leur poids par le biais d'équations asymptotiques. Nous commençons par rappeler la dérivation des principaux modèles généralement utilisés (Boussinesq, Green-Naghdi,...). Nous introduisons également un nouveau modèle exprimé en amplitude-flux qui correspond à une variante des équations de Nwogu. Dans le second chapitre, nous démontrons un résultat d'existence en temps long pour ces nouvelles équations et nous étudions l'existence d'ondes solitaires pour les équations de Boussinesq. Ce travail permet notamment de calculer avec une grande précision ces solutions exactes. Le troisième chapitre détaille les différences non linéaires que l'on retrouve entre les différentes équations de Boussinesq (modèles en flux-amplitude comparés aux modèles en vitesse-amplitude). Enfin, les deux derniers chapitres introduisent un nouveau paradigme pour trouver des schémas numériques adaptés aux modèles asymptotiques. L'idée est d'appliquer une analyse asymptotique aux équations d'Euler discrétisées. Ce nouveau paradigme est appliqué aux équations de Peregrine, de Nwogu et de Green-Naghdi. Plusieurs cas tests sont proposés dans ces deux chapitres. / In this work, we are interested in the evolution of water waves under the gravity force using asymptotics models. We start by recalling the derivation of most used models (Boussinesq, Green-Naghdi,...) and we introduce a new model expressed amplitude-flux, which is an alternative version of the Nwogu equations. In the second chapter, we prove a long time existence result for the new model and we investigate the existence of solitary waves for the Boussinesq models. This work allow us to compute these solutions with a good precision. The third chapter highlights the nonlinear differences between the Boussinesq equations (amplitude-flux models versus amplitude-velocity models). Finally, the two last chapter introduce a new paradigm in order to find numerical schemes adapted to asymptotics models. The idea is to apply an asymptotic analysis to a discretized Euler system. This new paradigm is applied to Peregrine equations, Nwogu equations and Green-Naghdi equations. Test cases are presented in these two chapters

This thesis is concerned with the analysis of the finite element method and the discontinuous Galerkin method for the numerical solution of an elliptic boundary value problem with a nonlinear Newton boundary condition in a two-dimensional polygonal domain. The weak solution loses regularity in a neighbourhood of boundary singularities, which may be at corners or at roots of the weak solution on edges. The main attention is paid to the study of error estimates. It turns out that the order of convergence is not dampened by the nonlinearity, if the weak solution is nonzero on a large part of the boundary. If the weak solution is zero on the whole boundary, the nonlinearity only slows down the convergence of the function values but not the convergence of the gradient. The same analysis is carried out for approximate solutions obtained with numerical integration. The theoretical results are verified by numerical experiments. 1

La résolution analytique du problème de la stabilité d’une suspension solide (particules rigides de forme sphérique) dans le système de Taylor-Couette cylindrique a été menée. On s’est basé sur les travaux de Ali and Lueptow (2002) pour formuler les équations régissant la stabilité de l’écoulement dans le cadre d’une théorie linéaire. Ces équations sont valables dans l’approximation du faible espace annulaire et ont pour but la prévision de l’instabilité primaire. A cet effet, nous avons utilisé une méthode variationnelle telle que la méthode de Galerkin pour résoudre le problème aux valeurs propres conduisant à établir le diagramme de stabilité liée au nombre d’onde au voisinage de l’état critique du développement de la première instabilité. Tout d’abord, on a cherché à mettre au point les calculs dans le cas de l’écoulement de Taylor-Couette classique en se référant aux travaux de Ali and Lueptow (2002). Ensuite on a procédé à la résolution systématique des équations du mouvement et l’on cherche à déterminer le critère d’apparition des instabilités en présence de particules en suspension et l’on détermine simultanément les paramètres de couplage entre forces d’interaction liquide-solide. L’ensemble des travaux ainsi réalisés permettront de lever la contradiction fondamentale entre la théorie et l’expérience. L’étude expérimentale a permis d’analyser les effets de la concentration des particules en suspension (disques) et du rapport d’aspect radial ’ sur l’apparition des instabilités dans le système de Taylor-Dean. Le dispositif expérimental utilisé est constitué d’un cylindre intérieur tournant et le cylindre extérieur est maintenu fixe. Le comportement rhéologique du fluide utilisé est viscoplastique obéissant au modèle de Herschel Bulkley. L’apparition des instabilités est examinée en utilisant une technique de visualisation. Pour une concentration donnée dans l’intervalle étudié, la nature des structures apparaissant dans le système d’écoulement dépend ’, alors que pour une valeur donnée de ’ dans l’intervalle étudié, la valeur du nombre de Taylor critique Tac dépend de la concentration des particules. Nous avons obtenu que le nombre de Taylor critique Tac correspondant au déclenchement de la première instabilité évolue non linéairement en fonction de ’. De plus, nous avons examiné expérimentalement les effets de limitation axiale (effet de bords) sur le déclenchement des instabilités dans le système de Taylor-Dean. Les résultats obtenus montrent que les bords tournants n’affectent pas le type de structures qui apparaissent dans le système d’écoulement de Taylor-Dean. Cependant, ils influencent le seuil critique d’apparition des instabilités qui est marquée par des valeurs élevées du nombre de Taylor critique pour des bords tournants ce qui indique un effet stabilisant des bords mobiles. / The analytical solution of the stability problem of a solid suspension (rigid spherical particles) in the system of cylindrical Taylor-Couette was conducted. We are based on the work of Ali and Lueptow (2002) to formulate the equations governing the stability of the flow in a linear theory. These equations are valid in the approximation of small gap configuration and aim to predict the primary instability. For this purpose, we used a variational method such as the Galerkin method to solve the eigenvalue problem leading to establish the stability diagram related to the wave number in the vicinity of the critical state of development of the first instability. First, we develop the calculations in the case of Taylor-Couette flow with reference to classic work of Ali and Lueptow (2002). Then, we carried out a systematic solution of the equations of motion and we search to determine the criterion of onset of instabilities in the presence of suspended particles and coupling parameters are simultaneously determined from liquid-solid interaction force. All work carried out and will remove the fundamental contradiction between theory and experiment. The experimental study has analyzed the effect of the concentration of suspended particles (disks) and radial aspect ratio ' on the occurrence of instabilities in the Taylor-Dean flow system. The experimental device used consists of a rotating inner cylinder and the outer cylinder is stationary. The rheological behavior of the fluid is viscoplastic obeying to Herschel Bulkley model. The onset of instability is examined using a visualization technique. For a given concentration in the range studied, the nature of the structure appearing in the flow system depends on ', while for a given value of ' in the range studied, the value of the critical Taylor number Tac depends on the particle concentration. We obtain that the critical Taylor number Tac corresponding to the onset of the first instability evolves nonlinearly versus '. In addition, we examined experimentally the effect of axial limitation (endwall effects) on the onset of instabilities in the Taylor-Dean flow system. The results show that the rotating endwalls do not affect the type of structures that appear in the Taylor-Dean flow system. However, they influence the threshold of appearance of instabilities which is characterized by high values of the critical Taylor number for rotating endwalls indicating a stabilizing effect of the rotating endwalls.

Le contexte scientifique de cette thèse est l'imagerie sismique dont le but est de reconstituer la structure du sous-sol de la Terre. Comme le forage a un coût assez élevé, l'industrie pétrolière s'intéresse à des méthodes capables de reconstituer les images de la structure terrestre interne avant de le faire. La technique d'imagerie sismique la plus utilisée est la technique de sismique-réflexion qui est basée sur le modèle de l'équation d'ondes. L'imagerie sismique est un problème inverse qui requiert de résoudre un grand nombre de problèmes directs. Dans ce contexte, nous nous intéressons dans cette thèse à la résolution du problème direct en régime harmonique, soit à la résolution des équations d'Helmholtz. L'objectif principal est de proposer et de développer un nouveau type de solveur élément fini (EF) caractérisé par un opérateur discret de taille réduite (comparée à la taille des solveurs déjà existants) sans pour autant altérer la précision de la solution numérique. Nous considérons les méthodes de Galerkine discontinues (DG). Comme les méthodes DG classiques sont plus coûteuses que les méthodes EF continues si l'on considère un même problème à cause d'un grand nombre de degrés de liberté couplés, résultat des approximations discontinues, nous développons une nouvelle classe de méthode DG réduisant ce problème : la méthode DG hybride (HDG). Pour valider l'efficacité de la méthode HDG proposée, nous comparons les résultats obtenus avec ceux obtenus avec une méthode DG basée sur des flux décentrés en 2D. Comme l'industrie pétrolière s'intéresse au traitement de données réelles, nous développons ensuite la méthode HDG pour les équations élastiques d'Helmholtz 3D. / The scientific context of this thesis is seismic imaging which aims at recovering the structure of the earth. As the drilling is expensive, the petroleum industry is interested by methods able to reconstruct images of the internal structures of the earth before the drilling. The most used seismic imaging method in petroleum industry is the seismic-reflection technique which uses a wave equation model. Seismic imaging is an inverse problem which requires to solve a large number of forward problems. In this context, we are interested in this thesis in the modeling part, i.e. the resolution of the forward problem, assuming a time-harmonic regime, leading to the so-called Helmholtz equations. The main objective is to propose and develop a new finite element (FE) type solver characterized by a reduced-size discrete operator (as compared to existing such solvers) without hampering the accuracy of the numerical solution. We consider the family of discontinuous Galerkin (DG) methods. However, as classical DG methods are much more expensive than continuous FE methods when considering steady-like problems, because of an increased number of coupled degrees of freedom as a result of the discontinuity of the approximation, we develop a new form of DG method that specifically address this issue: the hybridizable DG (HDG) method. To validate the efficiency of the proposed HDG method, we compare the results that we obtain with those of a classical upwind flux-based DG method in a 2D framework. Then, as petroleum industry is interested in the treatment of real data, we develop the HDG method for the 3D elastic Helmholtz equations.

Voluntary coughs are used as a diagnostic tool to detect lung diseases. Understanding the mechanics of a cough is therefore crucial to accurately interpreting the test results. A cough is characterised by a dynamic compression of the airways, resulting in large flow velocities and producing transient peak expiratory flows. Existing models for pulmonary flow have one or more of the following limitations: 1) they assume quasi-steady flows, 2) they assume low speed flows, 3) they assume a symmetrical branching airway system. The main objective of this thesis is to develop a model for a cough in the branching pulmonary airway system. First, the time-dependent one-dimensional equations for flow in a compliant tube is used to simulate a cough in a single airway. Using anatomical and physiological data, the tube law coupling the fluid and airway mechanics is constructed to accurately mimic the airway behaviour in its inflated and collapsed states. Next, a novel model for air flow in an airway bifurcation is constructed. The model is the first to capture successfully subcritical and supercritical flows across the bifurcation and allows for free time evolution from one case to another. The model is investigated by simulating a cough in both symmetric and asymmetric airway bifurcations. Finally, a cough model for the complete branching airway system is developed. The model takes into account the key factors involved in a cough; namely, the compliance of the lungs and the airways, the coughing effort and the sudden opening of the glottis. The reliability of the model is assessed by comparing the model predictions with previous experimental results. The model captures the main characteristics of forced expiatory flows; namely, the flow limitation phenomenon (the flow out of the lungs becomes independent of the applied expiratory effort) and the negative effort dependence phenomenon (the flow out of the lungs decreases with increasing expiratory effort). The model also gives a good qualitative agreement with the measured values of airway resistance. The location of the collapsed airway segment during forced expiration is, however, inconsistent with previous experimental results. The effect of changing the model parameters on the model predictions is therefore discussed.

Thesis deals with solving the problems of continuum mechanics by method of Isogeometric analysis. This relatively young method combines the advantages of precise NURBS geometry and robustness of the classical finite element method. The method is described on procedure of solving a plane Poissons boundary value problem. Solver is implemented in MatLab and algorithms are attached to the text.

This work is concerned with the applicability of several approximation methods (finite section method, Galerkin and collocation methods with maximum defect splines for uniform and non uniform meshes) to operators belonging to the closed subalgebra generated by operators of multiplication bz piecewise continuous functions and convolution operators also with piecewise continuous generating function.

The compressible Euler equations describe the motion of compressible inviscid fluids. They are used in many areas ranging from aerospace, automotive, and nuclear engineering to chemistry, ecology, climatology, and others. Mathematically, the compressible Euler equations represent a hyperbolic system consisting of several nonlinear partial differential equations (conservation laws). These equations are solved most frequently by means of Finite Volume Methods (FVM) and low-order Finite Element Methods (FEM). However, both these approaches are lacking higher order accuracy and moreover, it is well known that conforming FEM is not the optimal tool for the discretization of first-order equations. The most promissing approach to the approximate solution of the compressible Euler equations is the discontinuous Galerkin method that combines the stability of FVM, with excellent approximation properties of higher-order FEM. The objective of this Master Thesis was to develop, implement and test new adaptive algorithms for the nonstationary compressible Euler equations based on higher-order discontinuous Galerkin (hp-DG) methods. The basis for the new methods were the discontinuous Galerkin methods and space-time adaptive hp-FEM algorithms on dynamical meshes for nonstationary second-order problems. The new algorithms...