Resumo: Sistemas dinâmicos híbridos se diferenciam por exibir simultaneamente variados tipos de comportamento dinâmico (contínuo, discreto, eventos discretos) em diferentes partes do sistema. Neste trabalho foram estudados resultados de estabilidade no sentido de Lyapunov para sistemas dinâmicos híbridos gerais, que utilizam uma noção de tempo generalizado, definido em um espaço métrico totalmente ordenado. Mostrou-se que estes sistemas podem ser imersos em sistemas dinâmicos descontínuos definidos em R+, de forma que sejam preservadas suas propriedades qualitativas. Como foco principal, estudou-se resultados de estabilidade para sistemas dinâmicos descontínuos modelados por semigrupos de operadores, em que os estados do sistema pertencem à espaços de Banach. Neste caso, de forma alternativa à teoria clássica de estabilidade, os resultados não utilizam as usuais funções de Lyapunov, sendo portanto mais fáceis de se aplicar, tendo em vista a dificuldade em se encontrar tais funções para muitos sistemas. Além disso, os resultados foram aplicados à uma classe de equações diferenciais com retardo. / Abstract: Hybrid dynamical systems are characterized for showing simultaneously a variety of dynamic behaviors (continuous, discrete, discrete events) in different parts of the System. This work discusses stability results in the Lyapunov sense for general hybrid dynamical systems that use a generalized notion of time, defined in a completely ordered metric space. It has been shown that these systems may be immersed in discontinuous dynamical systems defined in R+, so that their quality properties are preserved. As the main focus, it is studied stability results for discontinuous dynamical systems modeled by semigroup operators, in which the states belong to Banach spaces. In this case, an alternative to the classical theory of stability, the results do not make use of the usual Lyapunov functions, and therefore are easier to apply, in view of the difficulty in finding such functions for many systems. Furthermore, the results were applied to a class of time-delay discontinuous differential equations. / Orientador: Geraldo Nunes Silva / Coorientador: Luís Antônio Fernandes de Oliveira / Banca: Carlos Alberto Raposo da Cunha / Banca: Waldemar Donizete Bastos / Mestre

Made available in DSpace on 2014-07-29T16:02:17Z (GMT). No. of bitstreams: 1 UBIRAJARA JOSe GAMA DE CASTRO.pdf: 516477 bytes, checksum: c0a8c62202b2da19be1a2dc69a29e416 (MD5) Previous issue date: 2011-02-24 / The main goal of this work is to study the behavior of Discontinuous Vector Fields in a neighborhood of a tipical singularity (tangency) using for this the regularization process developed by Teixeira and Sotomayor [9] and, using also, some technics of the Geometric Singular Perturbation Theory [2]. / O principal objetivo deste trabalho é estudar o comportamento numa vizinhança de uma singularidade típica (tangência) dos campos vetoriais descontínuos utilizando o processo de regularização desenvolvido por Teixeira e Sotomayor [9] e perturbações singulares [2].

The main emphasis of this thesis is to study a posteriori error analysis of discontinuous Galerkin (DG) methods for the elliptic variational inequalities. The DG methods have become very pop-ular in the last two decades due to its nature of handling complex geometries, allowing irregular meshes with hanging nodes and different degrees of polynomial approximation on different ele-ments. Moreover they are high order accurate and stable methods. Adaptive algorithms refine the mesh locally in the region where the solution exhibits irregular behaviour and a posteriori error estimates are the main ingredients to steer the adaptive mesh refinement. The solution of linear elliptic problem exhibits singularities due to change in boundary con-ditions, irregularity of coefficients and reentrant corners in the domain. Apart from this, the solu-tion of variational inequality exhibits additional irregular behaviour due to occurrence of the free boundary (the part of the domain which is a priori unknown and must be found as a component of the solution). In the lack of full elliptic regularity of the solution, uniform refinement is inefficient and it does not yield optimal convergence rate. But adaptive refinement, which is based on the residuals ( or a posteriori error estimator) of the problem, enhance the efficiency by refining the mesh locally and provides the optimal convergence. In this thesis, we derive a posteriori error estimates of the DG methods for the elliptic variational inequalities of the first kind and the second kind. This thesis contains seven chapters including an introductory chapter and a concluding chap-ter. In the introductory chapter, we review some fundamental preliminary results which will be used in the subsequent analysis. In Chapter 2, a posteriori error estimates for a class of DG meth-ods have been derived for the second order elliptic obstacle problem, which is a prototype for elliptic variational inequalities of the first kind. The analysis of Chapter 2 is carried out for the general obstacle function therefore the error estimator obtained therein involves the min/max func-tion and hence the computation of the error estimator becomes a bit complicated. With a mild assumption on the trace of the obstacle, we have derived a significantly simple and easily com-putable error estimator in Chapter 3. Numerical experiments illustrates that this error estimator indeed behaves better than the error estimator derived in Chapter 2. In Chapter 4, we have carried out a posteriori analysis of DG methods for the Signorini problem which arises from the study of the frictionless contact problems. A nonlinear smoothing map from the DG finite element space to conforming finite element space has been constructed and used extensively, in the analysis of Chapter 2, Chapter 3 and Chapter 4. Also, a common property shared by all DG methods allows us to carry out the analysis in unified setting. In Chapter 5, we study the C0 interior penalty method for the plate frictional contact problem, which is a fourth order variational inequality of the second kind. In this chapter, we have also established the medius analysis along with a posteriori analy-sis. Numerical results have been presented at the end of every chapter to illustrate the theoretical results derived in respective chapters. We discuss the possible extension and future proposal of the work presented in the Chapter 6. In the last chapter, we have documented the FEM codes used in the numerical experiments.

Cette thèse, comporte deux parties distinctes.Dans la première partie, on étudie l'existence et l'inexistence d'une inégalité qu'on a appelée l'inégalité de Hopf Uniforme (IHU), pour une équation linéaire de la forme Lv = f à coefficients bornés mesurables et sous les conditions de Dirichlet homogènes. L'IHU est une variante du principe de maximum, on l'a appliquée dans la preuve de la régularité W1;p 0 pour un problème semi-linéaire singulier : Lu = F(u) où les coefficients de L sont dans l'espace vmor (fonctions à oscillation moyenne évanescente) et F(u) est singulier en u = 0 F(0) = +∞. De plus, si les coefficients sont lipschitziens, on prouve que la régularité optimale du gradient de la solution u est bmor (fonctions à oscillation moyenne bornée i.e Grad u dans bmor).Dans la seconde partie, on s'intéresse à la régularité du système d'élasticité (équations stationnaires des ondes élastiques) avec une fonction source singulière au sens qu'elle n’est qu'intégrable par rapport à la fonction distance au bord du domaine. Via la dualité, nous montrons, selon ~f , que le problème admet une solution dite très faible dont le gradient n'est pas nécessairement intégrable sur tout le domaine mais uniquement localement. Nous déterminons aussi les fonctions vectorielles ~f pour lesquelles, ~u a son gradient intégrable sur tout l'espace de travail. / We discuss the existence and non existence of the so called Hopf uniform Inequality (variant of a maximum principle) for the linear equation Lv = f with measurable coefficients and under the homogeneous Dirichlet Boundary condition. Then we apply such inequality to prove the W1;p 0 -regularity of a semi linear problem Lu = F(u), singular at u = 0, with the coefficients of the main operator of L in the space of vanishing mean oscillation. Moreover, when those coefficients are Lipschitz, we show that the gradient of the solution is at most in the space of bounded mean oscillation : bmor. In the last part of this thesis, we are concerned with the linear easticity system (Stationnary equation of the waves elasticity). But, here the second terms varies with respect to the distance function until the boundary.Using the duality method, we study the regularity of the solution of the elasticity system for the data belonging to various weighted spaces.

L’objectif de cette thèse est de développer une méthode Galerkine discontinue d’ordre élevé capable de prendre en considération des simulations réalistes liées à la nanophotonique. Au cours des dernières décennies, l’évolution des techniques de lithographie a permis la création de structure géométriques de tailles nanométriques, révélant ainsi une large gamme de phénomènes nouveaux nés de l’interaction lumière-matière à ces échelles. Ces effets apparaissent généralement pour des objets de taille égale ou (très) inférieure à la longueur d’onde du champ incident. Ce travail repose sur le développement et l’implémentation de modèles de dispersion appropriés (principalement pour les métaux), ainsi que sur un large éventail de méthodes computationnelles classiques. Deux développements méthodologiques majeurs sont présentés et étudiés en détails: (i) les éléments courbes, et (ii) l’ordre d’approximation local. Ces études sont accompagnées de plusieurs cas-tests réalistes tirés de la nanophotonique. / The goal of this thesis is to develop a discontinuous Galerkin time-domain method to be able to handle realistic nanophotonics computations. During the last decades, the evolution of lithography techniques allowed the creation of geometrical structures at the nanometer scale, thus unveiling a variety of new phenomena arising from light-matter interactions at such levels. These effects usually occur when the device is of comparable size or (much) smaller than the wavelength of the incident field. This work relies on the development and implementation of appropriate models for dispersive materials (mostly metals), as well as on a large panel of classical computational techniques. Two major methodological developments are presented and studied in details: (i) curvilinear elements, and (ii) local order of approximation. This work is complemented with several physical studies of real-life nanophotonics applications.

Les pièces métalliques constituant le circuit primaire des installations nucléaires sont élaborées par forgeage à chaud. Pendant ce procédé, les transformations microstructurales induites par la déformation et les recuits déterminent une partie des propriétés mécaniques des produits finaux. L’orientation de la microstructure lors du processus de fabrication nécessite une connaissance précise des mécanismes physiques qui opèrent dans le matériau. Dans le cas de la déformation à chaud de l’acier austénitique 304L, ces modifications microstructurales dépendent de la recristallisation dynamique discontinue (DDRX) et de la recristallisation post-dynamique (PDRX). L’objet de ce projet est : (i) l’étude de la DDRX et de la PDRX dans les conditions de déformation du procédé de forgeage, (ii) l’étude de l’influence d’un ajout de niobium sur ces mécanismes, (iii) la modélisation de ces mécanismes afin de prédire les caractéristiques de la microstructure (moyenne et distribution de la taille des grains) à l’issue d’un procédé multipasses. Dans le cadre de l’étude, les conditions de déformation rencontrées lors du forgeage à chaud sont reproduites à l’aide d’essais de torsion sur des matériaux modèles contenant des teneurs en niobium différentes. La caractérisation et la modélisation des microstructures a permis de comprendre les effets respectifs de la température, de la vitesse de déformation ainsi que de l’ajout de niobium sur les mécanismes de la DDRX et de la PDRX. Dans cette étude, une nouvelle approche semi-topologique de l’hypothèse champs moyens est développée afin de permettre la prédiction de distributions de la taille de grain cohérentes avec les données expérimentales. / Cooling system of nuclear power plants is constituted of metallic parts obtained by hot forging. Thus during the manufacturing process, the microstructural trans- formations induced by the deformation and annealing process define partially the mechanical properties of the final products. A sharp knowledge of the physical mechanisms generated within the material is required to handle the microstructure. In the case of hot deformation of 304L austenitic stainless steel, the microstructural modifications depend on the discontinuous dynamic recrystallization (DDRX) and the post-dynamic recrystallization (PDRX). The aim of this project is: (i) the study of the DDRX and the PDRX under the conditions of deformation inherent in the forging process, (ii) the study of the influence of niobium addition on these mechanisms, (iii) the modeling of these me- chanisms in order to predict the microstructure characteristics (mean grain size and distribution) following a multipass process. As part of the research, the deformation conditions experienced during the hot forging process are replicated through torsion tests with model materials containing various niobium concentrations. Characterization and modeling of microstructures enable to understand the respective e ects of temperature, strain rate as well as niobium addittion on the DDRX and PDRX mechanisms. In this study, a new topological approach of mean-field hypothesis is developed in order to allow the prediction of realistic grain size distributions.

Les directives de l'Union Européenne conditionnent la R&D du secteur automobile concernant l'utilisation de matériaux plus légers ayant pour but de réduire la consommation de carburant et une diminution de l’émission de gaz d’échappement. L’objectif est de mettre au point des aciers allégés d’au moins 10% (ρmax ≈ 7g/cm3). Les alliages fer-aluminium possèdent des propriétés physiques et mécaniques prometteuses mais présentent des défauts de surface appelés roping, apparaissant après l’emboutissage à froid. Cette étude a consisté à mieux comprendre les conditions de recristallisation au cours du laminage à chaud afin de contrôler la microstructure et ainsi limiter ces défauts. Il est généralement admis, lors d’une déformation à chaud, que les alliages ferritiques, à haute énergie de défaut d’empilement, donnent lieu aux processus de recristallisation dynamique géométrique (RDG) et de recristallisation dynamique continue (RDC). Dans cette étude, l’existence d’une transition entre les mécanismes de RDC et de recristallisation dynamique discontinue (RDD) a été mise en évidence pour des températures comprises entre 900 et 1100°C et des vitesses de déformation comprises entre 0,1 et 50s 1. La recristallisation post dynamique a aussi été étudiée afin d’observer l’évolution de la microstructure lors de maintiens en température. Un modèle développé antérieurement pour la RDC de l’aluminium a ensuite été utilisé afin de simuler une passe de laminage. Bien que la comparaison des résultats expérimentaux et simulés fasse apparaître un certain nombre de différences, ce modèle permet de reproduire qualitativement les évolutions de la microstructure. / The instructions of the European Union pilot the R&D in the automotive industry regarding the use of lightweight materials which aims at reducing fuel consumption and emission of exhaust gases.The objective is to develop steels of density reduced by at least 10% (ρmax ≈ 7g/cm3). Iron aluminum alloys display promising physical and mechanical properties but they often exhibit surface defects, referred to as roping, appearing after the deep drawing process. This study was carried out to better understand the conditions of recrystallization during hot rolling to control the microstructure and thereby limit these defects.During hot deformation, it is generally agreed that geometric dynamic recrystallization (GDRX) and continuous dynamic recrystallization (CDRX) operate in ferritic alloys with high stacking fault energy. In this study, the existence of a transition between CDRX and the mechanism of discontinuous dynamic recrystallization (DDRX) has been brought into evidence in the temperature range 900 1100°C and strain rate range 0.1-50s-1. Post-dynamic recrystallization was also studied to observe the evolution of microstructure during holding temperatures.A model formerly developed for the CDRX of aluminum was then used to simulate a rolling pass. Comparison of computed and experimental results shows some differences but this model can reproduce microstructural changes qualitatively.

This thesis considers the problem of assigning a sentence its syntactic structure, which may be discontinuous. It proposes a class of models based on probabilistic grammars that are obtained by the automatic refinement of a given grammar. Different strategies for parsing with a refined grammar are developed. The induction, refinement, and application of two types of grammars (linear context-free rewriting systems and hybrid grammars) are evaluated empirically on two German and one Dutch corpus.

L'objectif principal de ce travail est d'améliorer la précision et l'efficacité des modèles LES au moyen des méthodes Galerkine discontinues (DG). Deux thématiques principales ont été étudiées: les stratégies d'adaptation spatiale et les modèles LES pour les méthodes d'ordre élevé.Concernant le premier thème, dans le cadre des méthodes DG la résolution spatiale peut être efficacement adaptée en modifiant localement soit le maillage (adaptation-h) soit le degré polynômial de la solution (adaptation-p). L'adaptation automatique de la résolution nécessite l'estimation des erreurs pour analyser la qualité de la solution locale et les exigences de résolution. L'efficacité de différentes stratégies de la littérature est comparée en effectuant des simulations h- et p-adaptatives.Sur la base de cette étude comparative, des algorithmes statiques et dynamiques p-adaptatifs pour la simulation des écoulements instationnaires sont ensuite développés et analysés. Les simulations numériques réalisées montrent que les algorithmes proposés peuvent réduire le coût de calcul des simulations des écoulements transitoires et statistiquement stationnaires.Un nouvel estimateur d'erreur est ensuite proposé. Il est local, car n'exige que des informations de l'élément et de ses voisins directs, et peut être calculé en cours de simulation pour un coût limité. Il est démontré que l'algorithme statique p-adaptatif basé sur cet estimateur d'erreur peut être utilisé pour améliorer la précision des simulations LES sur des écoulements turbulents statistiquement stationnaires.Concernant le second thème, une nouvelle méthode, consistante avec la discrétisation DG, est développée pour l'analyse a-priori des modèles DG-LES à partir des données DNS. Elle permet d'identifier le transfert d'énergie idéal entre les échelles résolues et non résolues. Cette méthode est appliquée à l'analyse de l'approche VMS (Variational Multiscale). Il est démontré que pour les résolutions fines, l'approche DG-VMS est capable de reproduire le transfert d'énergie idéal. Cependant, pour les résolutions grossières, typique de la LES à nombres de Reynolds élevés, un meilleur accord peut être obtenu en utilisant un modèle mixte Smagorinsky-VMS. / The main goal of this work is to improve the accuracy and computational efficiency of Large Eddy Simulations (LES) by means of discontinuous Galerkin (DG) methods. To this end, two main research topics have been investigated: resolution adaptation strategies and LES models for high-order methods.As regards the first topic, in the framework of DG methods the spatial resolution can be efficiently adapted by modifying either the local mesh size (h-adaptation) or the degree of the polynomial representation of the solution (p-adaptation).The automatic resolution adaptation requires the definition of an error estimation strategy to analyse the local solution quality and resolution requirements.The efficiency of several strategies derived from the literature are compared by performing p- and h-adaptive simulations. Based on this comparative study a suitable error indicator for the adaptive scale-resolving simulations is selected.Both static and dynamic p-adaptive algorithms for the simulation of unsteady flows are then developed and analysed. It is demonstrated by numerical simulations that the proposed algorithms can provide a reduction of the computational cost for the simulation of both transient and statistically steady flows.A novel error estimation strategy is then introduced. It is local, requiring only information from the element and direct neighbours, and can be computed at run-time with limited overhead. It is shown that the static p-adaptive algorithm based on this error estimator can be employed to improve the accuracy for LES of statistically steady turbulent flows.As regards the second topic, a novel framework consistent with the DG discretization is developed for the a-priori analysis of DG-LES models from DNS databases. It allows to identify the ideal energy transfer mechanism between resolved and unresolved scales.This approach is applied for the analysis of the DG Variational Multiscale (VMS) approach. It is shown that, for fine resolutions, the DG-VMS approach is able to replicate the ideal energy transfer mechanism.However, for coarse resolutions, typical of LES at high Reynolds numbers, a more accurate agreement is obtained by a mixed Smagorinsky-VMS model.

We present the analysis for the higher order continuous Galerkin−Petrov (cGP) time discretization schemes in combination with the one-level local projection stabilization in space applied to time-dependent convection-diffusion-reaction problems. Optimal a priori error estimates will be proved. Numerical studies support the theoretical results. Furthermore, a numerical comparison between continuous Galerkin−Petrov and discontinuous Galerkin time discretization schemes will be given.