We explore the behavior of an integrodifferential equation used to model one-dimensional biological swarms. In this model, we assume the motion of the swarm is determined by pairwise interactions, which in a continuous setting corresponds to a convolution of the swarm density with a pairwise interaction kernel. For a large class of interaction kernels, we derive conditions that lead to solutions which spread, blow up, or reach a steady state. For a smaller class of interaction kernels, we are able to make more quantitative predictions. In the spreading case, we predict the approximate shape and scaling of a similarity profile, as well as the approximate behavior at the endpoints of the swarm (via solutions to a traveling wave problem). In the blow up case, we derive an upper bound for the time to blow up. In the steady state case, we use previous results to predict the equilibrium swarm density. We support our predictions with numerical simulations. We also consider an extension of the original model which incorporates external forces. By analyzing and simulating particular cases, we determine that the addition of an external force can qualitatively change the behavior of the system.

The onset of flowering time in a plant is extremely significant when evaluating population success. Floral growth, seed production, and dispersal are all dependent upon flowering time. Flowering early (and hence longer) increases the prospect of pollination but typically reduces vegetative growth and yields fewer/smaller flowers. Flowering late (and hence shorter) guarantees more/bigger flowers but carries the risk of insufficient pollination. This fundamental trade-off between growth and flowering time suggests that there may be an optimal time to initiate flowering. In this thesis, we consider a deterministic hybrid integrodifferential model where we represent the growing season in continuous time and the time between seasons as a discrete map. We track the evolution of flowering time, as a phenotype, by explicitly considering it as a variable in our model. The model is analyzed from two different viewpoints: (1) by mutual invasion analysis in the sense of adaptive dynamics; and (2) by deriving equations for the mean trait value and total population density when flowering time is considered to be Gamma-distributed. In both cases evolution to an intermediary flowering time was observed.

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

In this study, we consider numerical schemes for a class of singular integrodifferential equations with a nonatomic difference operator. Our interest in this particular class has been motivated by an application in aeroelasticity. By applying nonconforming finite element methods, we are able to establish convergence for a semi-discrete scheme. We use an ordinary differential equation solver for the semi-discrete scheme and then improve the result by using a fully discretized scheme. We report our numerical findings and comment on the rates of convergence. / Ph. D.

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.

The thesis deals with qualitative properties of solutions of singular initial value problems for ordinary differential and integrodifferential equations which occur in the theory of linear and nonlinear electrical circuits and the theory of therminionic currents. The research is concentrated especially on questions of existence and uniqueness of solutions, asymptotic estimates of solutions and modications of Adomian decomposition method for singular initial problems. Solution algoritms are derived for scalar differential equations of Lane-Emden type using Taylor series and modication of the Adomian decomposition method. For certain classes of nonlinear of integrodifferential equations asymptotic expansions of solutions are constructed in a neighbourhood of a singular point. By means of the combination of Wazewski's topological method and Schauder xed-point theorem there are proved asymptotic estimates of solutions in a region which is homeomorphic to a cone having vertex coinciding with the initial point. Using Banach xed-point theorem the uniqueness of a solution of the singular initial value problem is proved for systems of integrodifferential equations of Volterra and Fredholm type including implicit systems. Moreover, conditions of continuous dependence of a solution on a parameter are determined. Obtained results are presented in illustrative examples.