Feedback control is used in almost every aspect of modern life and is essential in almost all engineering systems. Since no mathematical model is perfect and disturbances occur frequently, feedback is required. The design of a feedback control has been widely investigated in finite-dimensional space. However, many systems of interest, such as fluid flow and large structural vibrations are described by nonlinear partial differential equations and their state evolves on an infinite-dimensional Hilbert space. Developing controller design methods for nonlinear infinite-dimensional systems is not trivial. The objectives of this thesis are divided into multiple tasks. First, the well-posedness of some classes of nonlinear partial differential equations defined on a Hilbert space are investigated. The following nonlinear affine system defined on the Hilbert space H is considered z ̇(t)=F(z(t))+Bu(t), t≥0 z (0) = z0, where z(t) ∈ H is the state vector and z0 is the initial condition. The vector u(t) ∈ U, where U is a Hilbert space, is a state-feedback control. The nonlinear operator F : D ⊂ H → H is densely defined in H and the linear operator B : U → H is a linear bounded operator. Conditions for the closed-loop system to have a unique solution in the Hilbert space H are given. Next, finding a single bounded state-feedback control for nonlinear partial differential equations is discussed. In particular, Lyapunov-indirect method is considered to control nonlinear infinite-dimensional systems and conditions on when this method achieves the goal of local asymptotic stabilization of the nonlinear infinite-dimensional system are given. The Kuramoto-Sivashinsky (KS) equation defined in the Hilbert space L2(−π,π) with periodic boundary conditions is considered. ∂z/∂t =−ν∂4z/∂x4 −∂2z/∂x2 −z∂z/∂x, t≥0 z (0) = z0 (x) , where the instability parameter ν > 0. The KS equation is a nonlinear partial differential equation that is first-order in time and fourth-order in space. It models reaction-diffusion systems and is related to various pattern formation phenomena where turbulence or chaos appear. For instance, it models long wave motions of a liquid film over a vertical plane. When the instability parameter ν < 1, this equation becomes unstable. This is shown by analyzing the stability of the linearized system and showing that the nonlinear C0- semigroup corresponding to the nonlinear KS equation is Fr ́echet differentiable. There are a number of papers establishing the stabilization of this equation via boundary control. In this thesis, we consider distributed control with a single bounded feedback control for the KS equation with periodic boundary conditions. First, it is shown that sta- bilizing the linearized KS equation implies local asymptotical stability of the nonlinear KS equation. This is done by establishing Fr ́echet differentiability of the associated nonlinear C0-semigroup and showing that it is equal to the linear C0-semigroup generated by the linearization of the equation. Next, a single state-feedback control that locally asymptot- ically stabilizes the KS equation is constructed. The same approach to stabilize the KS equation from one equilibrium point to another is used. Finally, the solution of the uncontrolled/state-feedback controlled KS equation is ap- proximated numerically. This is done using the Galerkin projection method to approximate infinite-dimensional systems. The numerical simulations indicate that the proposed Lyapunov-indirect method works in stabilizing the KS equation to a desired state. Moreover, the same approach can be used to stabilize the KS equation from one constant equilibrium state to another.

Feedback control is used in almost every aspect of modern life and is essential in almost all engineering systems. Since no mathematical model is perfect and disturbances occur frequently, feedback is required. The design of a feedback control has been widely investigated in finite-dimensional space. However, many systems of interest, such as fluid flow and large structural vibrations are described by nonlinear partial differential equations and their state evolves on an infinite-dimensional Hilbert space. Developing controller design methods for nonlinear infinite-dimensional systems is not trivial. The objectives of this thesis are divided into multiple tasks. First, the well-posedness of some classes of nonlinear partial differential equations defined on a Hilbert space are investigated. The following nonlinear affine system defined on the Hilbert space H is considered z ̇(t)=F(z(t))+Bu(t), t≥0 z (0) = z0, where z(t) ∈ H is the state vector and z0 is the initial condition. The vector u(t) ∈ U, where U is a Hilbert space, is a state-feedback control. The nonlinear operator F : D ⊂ H → H is densely defined in H and the linear operator B : U → H is a linear bounded operator. Conditions for the closed-loop system to have a unique solution in the Hilbert space H are given. Next, finding a single bounded state-feedback control for nonlinear partial differential equations is discussed. In particular, Lyapunov-indirect method is considered to control nonlinear infinite-dimensional systems and conditions on when this method achieves the goal of local asymptotic stabilization of the nonlinear infinite-dimensional system are given. The Kuramoto-Sivashinsky (KS) equation defined in the Hilbert space L2(−π,π) with periodic boundary conditions is considered. ∂z/∂t =−ν∂4z/∂x4 −∂2z/∂x2 −z∂z/∂x, t≥0 z (0) = z0 (x) , where the instability parameter ν > 0. The KS equation is a nonlinear partial differential equation that is first-order in time and fourth-order in space. It models reaction-diffusion systems and is related to various pattern formation phenomena where turbulence or chaos appear. For instance, it models long wave motions of a liquid film over a vertical plane. When the instability parameter ν < 1, this equation becomes unstable. This is shown by analyzing the stability of the linearized system and showing that the nonlinear C0- semigroup corresponding to the nonlinear KS equation is Fr ́echet differentiable. There are a number of papers establishing the stabilization of this equation via boundary control. In this thesis, we consider distributed control with a single bounded feedback control for the KS equation with periodic boundary conditions. First, it is shown that sta- bilizing the linearized KS equation implies local asymptotical stability of the nonlinear KS equation. This is done by establishing Fr ́echet differentiability of the associated nonlinear C0-semigroup and showing that it is equal to the linear C0-semigroup generated by the linearization of the equation. Next, a single state-feedback control that locally asymptot- ically stabilizes the KS equation is constructed. The same approach to stabilize the KS equation from one equilibrium point to another is used. Finally, the solution of the uncontrolled/state-feedback controlled KS equation is ap- proximated numerically. This is done using the Galerkin projection method to approximate infinite-dimensional systems. The numerical simulations indicate that the proposed Lyapunov-indirect method works in stabilizing the KS equation to a desired state. Moreover, the same approach can be used to stabilize the KS equation from one constant equilibrium state to another.

<p>We are interested in backward-in-time solution techniques for evolutionary PDE problems arising in fluid mechanics. In addition to their intrinsic interest, such techniques have applications in recently proposed retrograde data assimilation. As our model system we consider the terminal value problem for the Kuramoto-Sivashinsky equation in a l D periodic domain. The Kuramoto-Sivashinsky equation, proposed as a model for interfacial and combustion phenomena, is often also adopted as a toy model for hydrodynamic turbulence because of its multiscale and chaotic dynamics. Such backward problems are typical examples of ill-posed problems, where any disturbances are amplified exponentially during the backward march. Hence, regularization is required to solve such problems efficiently in practice. We consider regularization approaches in which the original ill-posed problem is approximated with a less ill-posed problem, which is achieved by adding a regularization term to the original equation. While such techniques are relatively well-understood for linear problems, it is still unclear what effect these techniques may have in the nonlinear setting. In addition to considering regularization terms with fixed magnitudes, we also explore a novel approach in which these magnitudes are adapted dynamically using simple concepts from the Control Theory.</p> / Thesis / Master of Science (MSc)

The problem of constructing data-based, predictive, reduced models for the Kuramoto–Sivashinsky equation is considered, under circumstances where one has observation data only for a small subset of the dynamical variables. Accurate prediction is achieved by developing a discrete-time stochastic reduced system, based on a NARMAX (Nonlinear Autoregressive Moving Average with eXogenous input) representation. The practical issue, with the NARMAX representation as with any other, is to identify an efficient structure, i.e., one with a small number of terms and coefficients. This is accomplished here by estimating coefficients for an approximate inertial form. The broader significance of the results is discussed.

This dissertation focuses on the further development of creating accurate numerical models of complex dynamical systems using the holistic discretisation technique [Roberts, Appl. Num. Model., 37:371-396, 2001]. I extend the application from second to fourth order systems and from only one spatial dimension in all previous work to two dimensions (2D). We see that the holistic technique provides useful and accurate numerical discretisations on coarse grids. We explore techniques to model the evolution of spatial patterns governed by pdes such as the Kuramoto-Sivashinsky equation and the real-valued Ginzburg-Landau equation. We aim towards the simulation of fluid flow and convection in three spatial dimensions. I show that significant steps have been taken in this dissertation towards achieving this aim. Holistic discretisation is based upon centre manifold theory [Carr, Applications of centre manifold theory, 1981] so we are assured that the numerical discretisation accurately models the dynamical system and may be constructed systematically. To apply centre manifold theory the domain is divided into elements and using a homotopy in the coupling parameter, subgrid scale fields are constructed consisting of actual solutions of the governing partial differential equation(pde). These subgrid scale fields interact through the introduction of artificial internal boundary conditions. View the centre manifold (macroscale) as the union of all states of the collection of subgrid fields (microscale) over the physical domain. Here we explore how to extend holistic discretisation to the fourth order Kuramoto-Sivashinsky pde. I show that the holistic models give impressive accuracy for reproducing the steady states and time dependent phenomena of the Kuramoto-Sivashinsky equation on coarse grids. The holistic method based on local dynamics compares favourably to the global methods of approximate inertial manifolds. The excellent performance of the holistic models shown here is strong evidence in support of the holistic discretisation technique. For shear dispersion in a 2D channel a one-dimensional numerical approximation is generated directly from the two-dimensional advection-diffusion dynamics. We find that a low order holistic model contains the shear dispersion term of the Taylor model [Taylor, IMA J. Appl. Math., 225:473-477, 1954]. This new approach does not require the assumption of large x scales, formerly absolutely crucial in deriving the Taylor model. I develop holistic discretisation for two spatial dimensions by applying the technique to the real-valued Ginzburg-Landau equation as a representative example of second order pdes. The techniques will apply quite generally to second order reaction-diffusion equations in 2D. This is the first study implementing holistic discretisation in more than one spatial dimension. The previous applications of holistic discretisation have developed algebraic forms of the subgrid field and its evolution. I develop an algorithm for numerical construction of the subgrid field and its evolution for 1D and 2D pdes and explore various alternatives. This new development greatly extends the class of problems that may be discretised by the holistic technique. This is a vital step for the application of the holistic technique to higher spatial dimensions and towards discretising the Navier-Stokes equations.

The PhD thesis considers four topics in dynamical systems and is based on one paper and three manuscripts. In Paper I we apply methods of interval analysis in order to compute the rigorous enclosure of rotation number. The described algorithm is supplemented with a method of proving the existence of periodic points which is used to check rationality of the rotation number. In Manuscript II we provide a numerical algorithm for computing critical points of the multiplier map for the quadratic family (i.e., points where the derivative of the multiplier with respect to the complex parameter vanishes). Manuscript III concerns continued fractions of quadratic irrationals. We show that the generating function corresponding to the sequence of denominators of the best rational approximants of a quadratic irrational is a rational function with integer coefficients. As a corollary we can compute the Lévy constant of any quadratic irrational explicitly in terms of its partial quotients. Finally, in Manuscript IV we develop a method for finding rigorous enclosures of all odd periodic solutions of the stationary Kuramoto-Sivashinsky equation. The problem is reduced to a bounded, finite-dimensional constraint satisfaction problem whose solution gives the desired information about the original problem. Developed approach allows us to exclude the regions in L2, where no solution can exist.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência. / A numerical approach is presented for amodel describing the pattern formation by ion beam sputtering on a material surface. This process is responsible for the appearance of unexpectedly organized patterns, such as ripples, nanodots, and hexagonal arrays of nanoholes. A numerical analysis of preexisting patterns is proposed to investigate surface dynamics, based on a model resumed in an anisotropic damped Kuramoto-Sivashinsky equation, in a two dimensional surface with periodic boundary conditions. While deterministic, its highly nonlinear character gives a rich range of results, making it possible to describe accurately different patterns. A finite-difference semi-implicit time splitting scheme is employed on the discretization of the governing equation. Simulations were conducted with realistic coefficients related to physical parameters (anisotropies, beam orientation, diffusion). The stability of the numerical scheme is analyzed with time step and grid spacing tests for the pattern evolution, and the Method ofManufactured Solutions has been used to verify the scheme. Ripples and hexagonal patterns were obtained from amonomodal initial condition for certain values of the damping coefficient, while spatiotemporal chaos appeared for lower values. The anisotropy effects on pattern formation were studied, varying the angle of incidence.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência. / A numerical approach is presented for amodel describing the pattern formation by ion beam sputtering on a material surface. This process is responsible for the appearance of unexpectedly organized patterns, such as ripples, nanodots, and hexagonal arrays of nanoholes. A numerical analysis of preexisting patterns is proposed to investigate surface dynamics, based on a model resumed in an anisotropic damped Kuramoto-Sivashinsky equation, in a two dimensional surface with periodic boundary conditions. While deterministic, its highly nonlinear character gives a rich range of results, making it possible to describe accurately different patterns. A finite-difference semi-implicit time splitting scheme is employed on the discretization of the governing equation. Simulations were conducted with realistic coefficients related to physical parameters (anisotropies, beam orientation, diffusion). The stability of the numerical scheme is analyzed with time step and grid spacing tests for the pattern evolution, and the Method ofManufactured Solutions has been used to verify the scheme. Ripples and hexagonal patterns were obtained from amonomodal initial condition for certain values of the damping coefficient, while spatiotemporal chaos appeared for lower values. The anisotropy effects on pattern formation were studied, varying the angle of incidence.