Title: Weak solutions for a class of nonlinear integrodifferential equations Author: Ivan Soukup Department: Department of mathematical analysis Supervisor: RNDr. Tomáš Bárta, Ph.D. Supervisor's e-mail address: tomas.barta@mff.cuni.cz Abstract: The work investigates a system of evolutionary nonlinear partial integrodifferential equations in three dimensional space. In particular it stud- ies an existence of a solution to the system introduced in [1] with Dirichlet boundary condition and initial condition u0. We adopt the scheme of the proof from [9] and try to avoid the complications rising from the integral term. The procedure consists of an approximation of the convective term and an ap- proximation of the potentials of both nonlinearities using a quadratic function, proving the existence of the approximative solution and then returning to the original problem via regularity of the approximative solution and properties of the nonlinearities. The aim is to improve the results of the paper [1]. 1

This thesis deals with developing Galerkin based solution strategies for several important classes of differential equations involving derivatives and integrals of various fractional orders. Fractional order calculus finds use in several areas of science and engineering. The use of fractional derivatives may arise purely from the mathematical viewpoint, as in controller design, or it may arise from the underlying physics of the material, as in the damping behavior of viscoelastic materials. The physical origins of the fractional damping motivated us to study viscoelastic behavior of disordered materials at three levels. At the first level, we review two first principles models of rubber viscoelasticity. This leads us to study, at the next two levels, two simple disordered systems. The study of these two simplified systems prompted us towards an infinite dimensional system which is mathematically equivalent to a fractional order derivative or integral. This infinite dimensional system forms the starting point for our Galerkin projection based approximation scheme. In a simplified study of disordered viscoelastic materials, we show that the networks of springs and dash-pots can lead to fractional power law relaxation if the damping coefficients of the dash-pots follow a certain type of random distribution. Similar results are obtained when we consider a more simplified model, which involves a random system coefficient matrix. Fractional order derivatives and integrals are infinite dimensional operators and non-local in time: the history of the state variable is needed to evaluate such operators. This non-local nature leads to expensive long-time computations (O(t2) computations for solution up to time t). A finite dimensional approximation of the fractional order derivative can alleviate this problem. We present one such approximation using a Galerkin projection. The original infinite dimensional system is replaced with an equivalent infinite dimensional system involving a partial differential equation (PDE). The Galerkin projection reduces the PDE to a finite system of ODEs. These ODEs can be solved cheaply (O(t) computations). The shape functions used for the Galerkin projection are important, and given attention. Calculations with both global shape functions as well as finite elements are presented. The discretization strategy is improved in a few steps until, finally, very good performance is obtained over a user-specifiable frequency range (not including zero). In particular, numerical examples are presented showing good performance for frequencies varying over more than 7 orders of magnitude. For any discretization held fixed, however, errors will be significant at sufficiently low or high frequencies. We discuss why such asymptotics may not significantly impact the engineering utility of the method. Following this, we identify eight important classes of fractional differential equations (FDEs) and fractional integrodifferential equations (FIEs), and develop separate Galerkin based solution strategies for each of them. Distinction between these classes arises from the fact that both Riemann-Liouville as well as Caputo type derivatives used in this work do not, in general, follow either the law of exponents or the commutative property. Criteria used to identify these classes include; the initial conditions used, order of the highest derivative, integer or fractional order highest derivative, single or multiterm fractional derivatives and integrals. A key feature of our approximation scheme is the development of differential algebraic equations (DAEs) when the highest order derivative is fractional or the equation involves fractional integrals only. To demonstrate the effectiveness of our approximation scheme, we compare the numerical results with analytical solutions, when available, or with suitably developed series solutions. Our approximation scheme matches analytical/series solutions very well for all classes considered.

Dans ce mémoire, nous considérons un modèle épidémiologique multigroupe dans une population hétérogène, pour décrire la situation de l’épidémie de la COVID-19 dans une résidence pour personnes âgées. L’hétérogénéité liée ici à l’âge reflète une transmission élevée dûe à des interactions accrues, et un taux de mortalité plus élevé chez les personnes âgées. Du point de vue mathématique, nous obtenons un modèle SEIR multigroupe d’équations intégro-différentielles dans lequel nous considérons une distribution générale de la période infectieuse. Nous utilisons la méthode des fonctions de Lyapunov et une approche de la théorie des graphes pour déterminer le rôle du nombre de reproduction de base \(\mathcal{R}_0\) : l’état d’équilibre sans maladie est globalement asymptotiquement stable et l’épidémie s’éteint dans les deux groupes lorsque \(\mathcal{R}_0 \leq 1\), par contre elle persiste et l’état d’équilibre endémique est globalement asymptotiquement stable lorsque \(\mathcal{R}_0>1\). Les simulations numériques illustrent l’impact des stratégies de contrôle de la santé publique. / In this thesis, we consider a multiple group epidemiological model in a heterogeneous population to describe COVID-19 outbreaks in an elderly residential population. Age-based heterogeneity reflects higher transmission with enhanced interactions, and higher fatality rates in the elderly. Mathematically, we analyse a SEIR model in the form of a system of integro-differential equations with general distribution function for the infectious period. Lyapunov functions and graph-theoretical methods are employed to establish the role played by the basic reproduction ratio \(\mathcal{R}_0\) : global asymptotic stability of the disease-free equilibrium and no sustained outbreak when \(\mathcal{R}_0 \leq 1\), as opposed to persistent outbreak and globally asymptotic endemic equilibrium when \(\mathcal{R}_0>1\). Numerical simulations are presented to illustrate public health control strategies.