To understand the mechanisms of energy exchange between the tropics and the midlatitudes, it is necessary to develop simplified climate models. Motivated by linear wave theory, one such model is derived below. It captures the nonlinear interaction between barotropic and first baroclinic modes. In particular, it allows for the study of the barotropic response to a baroclinic forcing. Numerical methods for handling this nonlinear system are carefully developed and validated. The response generated by a physically realistic Kelvin wave forcing is studied and is found to consist mainly of one eastward propagating wave (phase-locked to the forcing) and two westward propagating (Rossby) waves. The Rossby waves are shown to be highly constrained by the initial parameters of the forcing and an explanation of this result is proposed.

There is often freedom in choosing the location of actuators on systems governed by partial differential equations. The actuator locations should be selected in order to optimize the performance criterion of interest. The main focus of this thesis is to consider H-∞-performance with state-feedback. That is, both the controller and the actuator locations are chosen to minimize the effect of disturbances on the output of a full-information plant. Optimal H-∞-disturbance attenuation as a function of actuator location is used as the cost function. It is shown that the corresponding actuator location problem is well-posed. In practice, approximations are used to determine the optimal actuator location. Conditions for the convergence of optimal performance and the corresponding actuator location to the exact performance and location are provided. Examples are provided to illustrate that convergence may fail when these conditions are not satisfied. Systems of large model order arise in a number of situations; including approximation of partial differential equation models and power systems. The system descriptions are sparse when given in descriptor form but not when converted to standard first-order form. Numerical calculation of H-∞-attenuation involves iteratively solving large H-∞-algebraic Riccati equations (H-∞-AREs) given in the descriptor form. An iterative algorithm that preserves the sparsity of the system description to calculate the solutions of large H-∞-AREs is proposed. It is shown that the performance of our proposed algorithm is similar to a Schur method in many cases. However, on several examples, our algorithm is both faster and more accurate than other methods. The calculation of H-∞-optimal actuator locations is an additional layer of optimization over the calculation of optimal attenuation. An optimization algorithm to calculate H-∞-optimal actuator locations using a derivative-free method is proposed. The results are illustrated using several examples motivated by partial differential equation models that arise in control of vibration and diffusion.

No description available.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In this thesis we present a proof of the Malgrange-Ehrenpreis theorem, which states that every operator with constant coefficients non identically zero has a fundamental solution. / Nessa dissertação apresentamos uma demonstração do Teorema de Malgrange-Ehrenpreis, que afirma que todo operador de coeficientes constantes não identicamente nulo tem uma solução fundamental.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / In this dissertation we present a proof of a Bernard Malgrange theorem, which establishes a necessary and sufficient condition for the global solvability of a linear operator with constant coefficients. / Nessa disserta¸c ao apresentamos uma demonstra¸c ao do Teorema de Bernard Malgrange, que estabelece condi¸c ao necess´aria e suficiente para que um operador linear com coeficientes constantes seja globalmente resol´uvel.

In the thesis we investigate two problems on Partial Differential Equations (PDEs) in differential geometry and fluid mechanics. First, we prove the weak L^p continuity of the Gauss-Codazzi-Ricci (GCR) equations, which serve as a compatibility condition for the isometric immersions of Riemannian and semi-Riemannian manifolds. Our arguments, based on the generalised compensated compactness theorems established via functional and micro-local analytic methods, are intrinsic and global. Second, we prove the vanishing viscosity limit of an incompressible fluid in three-dimensional smooth, curved domains, with the kinematic and Navier boundary conditions. It is shown that the strong solution of the Navier-Stokes equation in H^r+1 (r > 5/2) converges to the strong solution of the Euler equation with the kinematic boundary condition in H^r, as the viscosity tends to zero. For the proof, we derive energy estimates using the special geometric structure of the Navier boundary conditions; in particular, the second fundamental form of the fluid boundary and the vorticity thereon play a crucial role. In these projects we emphasise the linkages between the techniques in differential geometry and mathematical hydrodynamics.

Nous développons une nouvelle approche éléments finis pour des équations aux dérivées partielles elliptiques de type élasticité linéaire ou Stokes sur une surface fermée de R3. La surface considérée est décrite par le zéro d'une fonction de niveau assez régulière. Le problème se ramène à la minimisation d'une fonctionnelle énergie pour le champ de vitesse sous contraintes. Les contraintes sont de deux types : (i) la vitesse est tangentielle à la surface, (ii) la surface est inextensible. Cette deuxième contrainte équivaut à l'incompressibilité surfacique du champ de vitesse. Nous abordons ce problème de deux façons : la pénalisation et l'introduction de deux multiplicateurs de Lagrange. Cette dernière méthode a l'avantage de traiter le cas de la limite incompressible d'un écoulement en surface dont nous présentons pour la première fois l'analyse théorique et numérique. Nous montrons des estimations d'erreurs sur la solution discrète et les tests numériques confirment l'optimalité des ces estimations. Pour cela, nous proposons plusieurs approches pour le calcul numérique de la normale et la courbure de la surface. L'implémentation utilise la librairie libre d'éléments finis Rheolef. Nous présentons aussi des résultats de simulations numériques pour une application en biologie : la morphogenèse de l'embryon de la drosophile, durant laquelle des déformations tangentielles d'une monocouche de cellules avec une faible variation d'aire. Ce phénomène est connu sous le nom de l'extension de la bande germinale. / We develop a novel finite element approach for linear elasticity or Stokes-type PDEs set on a closed surface of $mathbb{R}^3$. The surface we consider is described as the zero of a sufficiently smooth level-set function. The problem can be written as the minimisation of an energy function over a constrained velocity field. Constraints areof two different types: (i) the velocity field is tangential to the surface, (ii) the surface is inextensible. This second constraint is equivalent to surface incompressibility of the velocity field. We address thisproblem in two different ways: a penalty method and a mixed method involving two Lagrange multipliers. This latter method allows us to solve the limiting case of incompressible surface flow, for which we present a novel theoretical and numerical analysis. Error estimates for the discrete solution are given andnumerical tests shows the optimality of the estimates. For this purpose, several approaches for the numerical computation of the normal and curvature of the surface are proposed. The implementation relies on the Rheolef open-source finite element library. We present numerical simulations for a biological application: the morphogenesis of Drosophila embryos, duringwhich tangential flows of a cell monolayer take place with a low surface-area variation. This phenomenon is known as germ-band extension.

Sejam M uma variedade de Hadamard com curvatura seccional KM ≤ −k2 < 0 e ∂ M sua fronteira assintótica. Dizemos que M satisfaz a condição de convexidade estrita se, dados x ∈ ∂∞M e W ⊂ ∂∞M aberto relativo contendo x, existe um aberto Ω ⊂ M de classe C2 tais que x ∈ Int (∂ Ω) ⊂ W e M \ Ω ´e convexo. Provamos que a condição de convexidade estrita implica que M éregular no infinito com relação ao operador Q[u] := div a(|∇u|) \ |∇u| ∇u definido no espa¸co de Sobolev W 1,p(M ), onde a ∈ C1([0, +∞)) satisfaz a(0) = 0, at(s) > 0 para todo s > 0, a(s) ≤ C (sp−1 + 1), ∀s ≥ 0, onde C > 0 é uma constante, e a(s) ≥ sq para algum q > 0 e para s ≈ 0 e supomos que é possível resolver problemas de Dirichlet em bolas (compactas) de M com dados contínuos no bordo. Segue disto que sob a condição de convexidade estrita, os problemas de Dirichlet para equação de hipersuperfície mínima e para o p-laplaciano, p > 1, são solúveis para qualquer dado contínuo prescrito no bordo assintótico. Também provamos que se M é rotacionalmente simétrica ou se inf BR+1 KM ≥ −e 2kR /R2+2 , R ≥ R∗, para certos R∗ e E > 0, então M satisfaz a condição de convexidade estrita. / Let M be Hadamard manifold with sectional curvature KM ≤ −k2, k > 0 and ∂∞M its asymptotic boundary. We say that M satisfies the strict convexity condition if, given x ∈ ∂∞M and a relatively open subset W ⊂ 2 ∂∞M containing x, there exists a C open subset Ω ⊂ M such that x ∈ Int (∂∞Ω) ⊂ W and M \ Ω is convex. We prove that the strict convexity condition implies that M is regular at infinity relative to the operator Q [u] := div a(|∇u|) \ |∇u| ∇u , defined on the Sobolev space W 1,p(M ), where a ∈ C 1 ([0, ∞)) satisfies a(0) = 0, at(s) > 0 for all s > 0, a(s) ≤ C (s p−1 + 1), ∀s ≥ 0, where C > 0 is a constant, and a(s) ≥ sq , for some q > 0 and for s ≈ 0 and we suppose that it is possible to solve Dirichlet problems on (compact) balls of M with continuous boundary data. It follows that under the strict convexity condition, the Dirichlet problems for the minimal hypersurface and the p-Laplacian, p > 1, equations are solvable for any prescribed continuous asymptotic boundary data. We also prove that if M is rotationally symmetric or if inf BR+1 KM ≥ −e2kR/R2+2 , R ≥ R∗, for some R∗ and E > 0, then M satisfies the SC condition.

This work is dedicated to the local and global well-possednes study of Cauchy s Problem associated to the nonlinear Schrödinger equation, to the initial data nonzero at infinity. / Conselho Nacional de Desenvolvimento Científico e Tecnológico / Este trabalho é dedicado ao estudo da boa colocação local e global do Problema de Cauchy associado à equação não linear de Schrödinger, com dado inicial não nulo no infinito.

Interval analysis is an essential tool in the construction of validated numerical solutions of Initial Value Problems (IVP) for Ordinary (ODE) and Partial (PDE) Differential Equations. A validated solution typically consists of guaranteed lower and upper bounds for the exact solution or set of exact solutions in the case of uncertain data, i.e. it is an interval function (enclosure) containing all solutions of the problem. IVP for ODE: The central point of discussion is the wrapping effect. A new concept of wrapping function is introduced and applied in studying this effect. It is proved that the wrapping function is the limit of the enclosures produced by any method of certain type (propagate and wrap type). Then, the wrapping effect can be quantified as the difference between the wrapping function and the optimal interval enclosure of the solution set (or some norm of it). The problems with no wrapping effect are characterized as problems for which the wrapping function equals the optimal interval enclosure. A sufficient condition for no wrapping effect is that there exist a linear transformation, preserving the intervals, which reduces the right-hand side of the system of ODE to a quasi-isotone function. This condition is also necessary for linear problems and "near" necessary in the general case. Hyperbolic PDE: The Initial Value Problem with periodic boundary conditions for the wave equation is considered. It is proved that under certain conditions the problem is an operator equation with an operator of monotone type. Using the established monotone properties, an interval (validated) method for numerical solution of the problem is proposed. The solution is obtained step by step in the time dimension as a Fourier series of the space variable and a polynomial of the time variable. The numerical implementation involves computations in Fourier and Taylor functoids. Propagation of discontinuo~swaves is a serious problem when a Fourier series is used (Gibbs phenomenon, etc.). We propose the combined use of periodic splines and Fourier series for representing discontinuous functions and a method for propagating discontinuous waves. The numerical implementation involves computations in a Fourier hyper functoid. / Mathematical Sciences / D. Phil. (Mathematics)