Ordinary Differential Equations and Applied Dynamics Commons^™

Sir Models: Differential Equations That Support The Common Good, Lorelei Koss

CODEE Journal

This article surveys how SIR models have been extended beyond investigations of biologically infectious diseases to other topics that contribute to social inequality and environmental concerns. We present models that have been used to study sustainable agriculture, drug and alcohol use, the spread of violent ideologies on the internet, criminal activity, and health issues such as bulimia and obesity.

Newton's Law Of Cooling, Caleb J. Emmons

Journal of Humanistic Mathematics

A poem reflecting three different viewpoints on Newton's Law of Cooling.

The Global Stability Of The Solution To The Morse Potential In A Catastrophic Regime, Weerapat Pittayakanchit

HMC Senior Theses

Swarms of animals exhibit aggregations whose behavior is a challenge for mathematicians to understand. We analyze this behavior numerically and analytically by using the pairwise interaction model known as the Morse potential. Our goal is to prove the global stability of the candidate local minimizer in 1D found in A Primer of Swarm Equilibria. Using the calculus of variations and eigenvalues analysis, we conclude that the candidate local minimizer is a global minimum with respect to all solution smaller than its support. In addition, we manage to extend the global stability condition to any solutions whose support has a single …

Topological Data Analysis For Systems Of Coupled Oscillators, Alec Dunton

HMC Senior Theses

Coupled oscillators, such as groups of fireflies or clusters of neurons, are found throughout nature and are frequently modeled in the applied mathematics literature. Earlier work by Kuramoto, Strogatz, and others has led to a deep understanding of the emergent behavior of systems of such oscillators using traditional dynamical systems methods. In this project we outline the application of techniques from topological data analysis to understanding the dynamics of systems of coupled oscillators. This includes the examination of partitions, partial synchronization, and attractors. By looking for clustering in a data space consisting of the phase change of oscillators over a …

An Interactive Tool For The Computational Exploration Of Integrodifference Population Models, Kennedy Agwamba

HMC Senior Theses

Mathematical modeling of population dynamics can provide novel insight to the growth and dispersal patterns for a variety of species populations, and has become vital to the preservation of biodiversity on a global-scale. These growth and dispersal stages can be modeled using integrodifference equations that are discrete in time and continuous in space. Previous studies have identified metrics that can determine whether a given species will persist or go extinct under certain model parameters. However, a need for computational tools to compute these metrics has limited the scope and analysis within many of these studies. We aim to create computational …

A Mathematical Model Of The Effect Of Aspirin On Blood Clotting, Breeana J. Johng

Scripps Senior Theses

In this paper, we provide a mathematical model of the effect of aspirin on blood clotting. The model tracks the enzyme prostaglandin H synthase and an important blood clotting factor, thromboxane A₂, in the form of thromboxane B₂. Through model analysis, we determine conditions under which the reactions of prostaglandin H synthase are self-sustaining. Lastly, through numerical simulations, we demonstrate that the model accurately captures the steady-state chemical concentrations of interest in blood, both with and without aspirin treatment.

Arnold Diffusion, Florin Diacu

Journal of Humanistic Mathematics

No abstract provided.

Eradicating Malaria: Improving A Multiple-Timestep Optimization Model Of Malarial Intervention Policy, Taryn M. Ohashi

Scripps Senior Theses

Malaria is a preventable and treatable blood-borne disease whose complications can be fatal. Although many interventions exist in order to reduce the impacts of malaria, the optimal method of distributing these interventions in a geographical area with limited resources must be determined. This thesis refines a model that uses an integer linear program and a compartmental model of epidemiology called an SIR model of ordinary differential equations. The objective of the model is to find an intervention strategy over multiple time steps and multiple geographic regions that minimizes the number of days people spend infected with malaria. In this paper, …

A Comparison And Catalog Of Intrinsic Tumor Growth Models, Elizabeth A. Sarapata

HMC Senior Theses

Determining the dynamics and parameter values that drive tumor growth is of great interest to mathematical modelers, experimentalists and practitioners alike. We provide a basis on which to estimate the growth dynamics of ten different tumors by fitting growth parameters to at least five sets of published experimental data per type of tumor. These timescale tumor growth data are also used to determine which of the most common tumor growth models (exponential, power law, logistic, Gompertz, or von Bertalanffy) provides the best fit for each type of tumor. In order to compute the best-fit parameters, we implemented a hybrid local-global …

R₀ Analysis Of A Spatiotemporal Model For A Stream Population, H. W. Mckenzie, Y. Jin, Jon T. Jacobsen, M. A. Lewis

All HMC Faculty Publications and Research

Water resources worldwide require management to meet industrial, agricultural, and urban consumption needs. Management actions change the natural flow regime, which impacts the river ecosystem. Water managers are tasked with meeting water needs while mitigating ecosystem impacts. We develop process-oriented advection-diffusion-reaction equations that couple hydraulic flow to population growth, and we analyze them to assess the effect of water flow on population persistence. We present a new mathematical framework, based on the net reproductive rate R₀ for advection-diffusion-reaction equations and on related measures. We apply the measures to population persistence in rivers under various flow regimes. This work lays …

Analytic And Numerical Studies Of A Simple Model Of Attractive-Repulsive Swarms, Andrew S. Ronan

HMC Senior Theses

We study the equilibrium solutions of an integrodifferential equation used to model one-dimensional biological swarms. We assume that the motion of the swarm is governed by pairwise interactions, or a convolution in the continuous setting, and derive a continuous model from conservation laws. The steady-state solution found for the model is compactly supported and is shown to be an attractive equilibrium solution via linear perturbation theory. Numerical simulations support that the steady-state solution is attractive for all initial swarm distributions. Some initial results for the model in higher dimensions are also presented.

A Semilinear Wave Equation With Smooth Data And No Resonance Having No Continuous Solution, Jose F. Caicedo, Alfonso Castro

All HMC Faculty Publications and Research

We prove that a boundary value problem for a semilinear wave equation with smooth nonlinearity, smooth forcing, and no resonance cannot have continuous solutions. Our proof shows that this is due to the non-monotonicity of the nonlinearity.

Traveling Waves And Shocks In A Viscoelastic Generalization Of Burgers' Equation, Victor Camacho '07, Robert D. Guy, Jon T. Jacobsen

All HMC Faculty Publications and Research

We consider traveling wave phenomena for a viscoelastic generalization of Burgers' equation. For asymptotically constant velocity profiles we find three classes of solutions corresponding to smooth traveling waves, piecewise smooth waves, and piecewise constant (shock) solutions. Each solution type is possible for a given pair of asymptotic limits, and we characterize the dynamics in terms of the relaxation time and viscosity.

As Flat As Possible, Jon T. Jacobsen

All HMC Faculty Publications and Research

How does one determine a surface which is as flat as possible, such as those created by soap film surfaces? What does it mean to be as flat as possible? In this paper we address this question from two distinct points of view, one local and one global in nature. Continuing with this theme, we put a temporal twist on the question and ask how to evolve a surface so as to flatten it as efficiently as possible. This elementary discussion provides a platform to introduce a wide range of advanced topics in partial differential equations and helps students …

Approximations Of Continuous Newton's Method: An Extension Of Cayley's Problem, Jon T. Jacobsen, Owen Lewis '05, Bradley Tennis '06

All HMC Faculty Publications and Research

Continuous Newton's Method refers to a certain dynamical system whose associated flow generically tends to the roots of a given polynomial. An Euler approximation of this system, with step size h=1, yields the discrete Newton's method algorithm for finding roots. In this note we contrast Euler approximations with several different approximations of the continuous ODE system and, using computer experiments, consider their impact on the associated fractal basin boundaries of the roots

Strings, Chains, And Ropes, Darryl H. Yong

All HMC Faculty Publications and Research

Following Antman [Amer. Math. Mon., 87 (1980), pp. 359–370], we advocate a more physically realistic and systematic derivation of the wave equation suitable for a typical undergraduate course in partial differential equations. To demonstrate the utility of this derivation, three applications that follow naturally are described: strings, hanging chains, and jump ropes.

Optimal Therapy Regimens For Treatment-Resistant Mutations Of Hiv, Weiqing Gu, Helen Moore

All HMC Faculty Publications and Research

In this paper, we use control theory to determine optimal treatment regimens for HIV patients, taking into account treatment-resistant mutations of the virus. We perform optimal control analysis on a model developed previously for the dynamics of HIV with strains of various resistance to treatment (Moore and Gu, 2005). This model incorporates three types of resistance to treatments: strains that are not responsive to protease inhibitors, strains not responsive to reverse transcriptase inhibitors, and strains not responsive to either of these treatments. We solve for the optimal treatment regimens analytically and numerically. We find parameter regimes for which optimal dosing …

A Mathematical Model For Treatment-Resistant Mutations Of Hiv, Helen Moore, Weiqing Gu

All HMC Faculty Publications and Research

In this paper, we propose and analyze a mathematical model, in the form of a system of ordinary differential equations, governing mutated strains of human immunodeficiency virus (HIV) and their interactions with the immune system and treatments. Our model incorporates two types of resistant mutations: strains that are not responsive to protease inhibitors, and strains that are not responsive to reverse transcriptase inhibitors. It also includes strains that do not have either of these two types of resistance (wild-type virus) and strains that have both types. We perform our analysis by changing the system of ordinary differential equations (ODEs) to …

A Liouville-Gelfand Equation For K-Hessian Operators, Jon T. Jacobsen

All HMC Faculty Publications and Research

In this paper we establish existence and multiplicity results for a class of fully nonlinear elliptic equations of k-Hessian type with exponential nonlinearity. In particular, we characterize the precise dependence of the multiplicity of solutions with respect to both the space dimension and the value of k. The choice of exponential nonlinearity is motivated by the classical Liouville-Gelfand problem from combustible gas dynamics and prescribed curvature problems.

Examples Of Cayley 4-Manifolds, Weiqing Gu, Christopher Pries '03

All HMC Faculty Publications and Research

We determine several families of so-called Cayley 4-dimensional manifolds in the real Euclidean 8-space. Such manifolds are of interest because Cayley 4-manifolds are supersymmetric cycles that are candidates for representations of fundamental particles in String Theory. Moreover, some of the examples of Cayley manifolds discovered in this paper may be modified to construct explicit examples in our current search for new holomorphic invariants for Calabi-Yau 4-folds and for the further development of mirror symmetry.

We apply the classic results of Harvey and Lawson to find Cayley manifolds which are graphs of functions from the set of quaternions to itself. We …

Semilinear Equations With Discrete Spectrum, Alfonso Castro

All HMC Faculty Publications and Research

This is an overview of the solvability of semilinear equations where the linear part has discrete spectrum. Semilinear elliptic and hyperbolic equations, as well as Hammerstein integral equations, are used as motivating examples. The presentation is intended to be accessible to non experts.

The Liouville-Bratu-Gelfand Problem For Radial Operators, Jon T. Jacobsen, Klaus Schmitt

All HMC Faculty Publications and Research

We determine precise existence and multiplicity results for radial solutions of the Liouville–Bratu–Gelfand problem associated with a class of quasilinear radial operators, which includes perturbations of k-Hessian and p-Laplace operators.

Monotone Solutions Of A Nonautonomous Differential Equation For A Sedimenting Sphere, Andrew Belmonte, Jon T. Jacobsen, Anandhan Jayaraman

All HMC Faculty Publications and Research

We study a class of integrodifferential equations and related ordinary differential equations for the initial value problem of a rigid sphere falling through an infinite fluid medium. We prove that for creeping Newtonian flow, the motion of the sphere is monotone in its approach to the steady state solution given by the Stokes drag. We discuss this property in terms of a general nonautonomous second order differential equation, focusing on a decaying nonautonomous term motivated by the sedimenting sphere problem

A Bifurcation Theorem And Applications, Alfonso Castro, Jorge Cossio

All HMC Faculty Publications and Research

In this paper we give a sufficient condition on the nonlinear operator N for a point (λ, u) to be a local bifurcation point of equations of the form u + λL^-1(N(u)) = 0, where L is a linear operator in a real Hilbert space, L has compact inverse, and λ ∈ R is a parameter. Our result does not depend on the variational structure of the equation or the multiplicity of the eigenvalue of the linear operator L. Applications are made to systems of differential equations and to the existence of periodic solutions of nonlinear second order …

Nonnegative Solutions For A Class Of Radially Symmetric Nonpositone Problems, Alfonso Castro, Ratnasingham Shivaji

All HMC Faculty Publications and Research

We consider the existence of radially symmetric non-negative solutions for the boundary value problem

\begin{displaymath}\begin{array}{*{20}{c}} { - \Delta u(x) = \lambda f(u(x))\qua... ...\\ {u(x) = 0\quad \left\Vert x \right\Vert = 1} \\ \end{array} \end{displaymath}

where $ \lambda > 0,f(0) < 0$ (non-positone), $ f' \geq 0$ and $ f$ is superlinear. We establish existence of non-negative solutions for $ \lambda $ small which extends some work of our previous paper on non-positone problems, where we considered the case $ N = 1$. Our work also proves a recent conjecture by Joel Smoller and Arthur Wasserman.

Infinitely Many Radially Symmetric Solutions To A Superlinear Dirichlet Problem In A Ball, Alfonso Castro, Alexandra Kurepa

All HMC Faculty Publications and Research

In this paper we show that a radially symmetric superlinear Dirichlet problem in a ball has infinitely many solutions. This result is obtained even in cases of rapidly growing nonlinearities, that is, when the growth of the nonlinearity surpasses the critical exponent of the Sobolev embedding theorem. Our methods rely on the energy analysis and the phase-plane angle analysis of the solutions for the associated singular ordinary differential equation.

Existence And Uniqueness For A Variational Hyperbolic System Without Resonance, Peter W. Bates, Alfonso Castro

All HMC Faculty Publications and Research

In this paper, we study the existence of weak solutions of the problem

□u + ∇G(u) = f(t,x) ; (t,x) є Ω ≡ (0,π)x(0,π)

u(t,x) = 0 ; (t,x) є ∂Ω

where □ is the wave operator ∂²/∂t² - ∂²/∂x², G: Rⁿ→R is a function of class C² such that ∇G(0) = 0 and f:Ώ→R^n is a continuous function having first derivative with respect to t in (L₂,(Ω))ⁿ and satisfying

f(0,x) = f(π,x) = 0

for all x є [0,π].