Applied Mathematics Commons^™

Social Aggregation In Pea Aphids: Experiment And Random Walk Modeling, Christa Nilsen, John Paige, Olivia Warner, Benjamin Mayhew, Ryan Sutley, Matthew Lam '15, Andrew J. Bernoff, Chad M. Topaz

All HMC Faculty Publications and Research

From bird flocks to fish schools and ungulate herds to insect swarms, social biological aggregations are found across the natural world. An ongoing challenge in the mathematical modeling of aggregations is to strengthen the connection between models and biological data by quantifying the rules that individuals follow. We model aggregation of the pea aphid, Acyrthosiphon pisum. Specifically, we conduct experiments to track the motion of aphids walking in a featureless circular arena in order to deduce individual-level rules. We observe that each aphid transitions stochastically between a moving and a stationary state. Moving aphids follow a correlated random walk. …

Nonlocal Aggregation Models: A Primer Of Swarm Equilibria, Andrew J. Bernoff, Chad M. Topaz

All HMC Faculty Publications and Research

Biological aggregations such as fish schools, bird flocks, bacterial colonies, and insect swarms have characteristic morphologies governed by the group members' intrinsic social interactions with each other and by their interactions with the external environment. Starting from a simple discrete model treating individual organisms as point particles, we derive a nonlocal partial differential equation describing the evolving population density of a continuum aggregation. To study equilibria and their stability, we use tools from the calculus of variations. In one spatial dimension, and for several choices of social forces, external forces, and domains, we find exact analytical expressions for the equilibria. …

Existence And Qualitative Properties Of Solutions For Nonlinear Dirichlet Problems, Alfonso Castro, Jorge Cossio, Carlos Vélez

All HMC Faculty Publications and Research

Sign-changing solutions to semilinear elliptic problems in connection with their Morse indices. To this end, we first establish a priori bounds for one-sign solutions. Secondly, using abstract saddle point principles we find large augmented Morse index solutions. In this part, extensive use is made of critical groups, Morse index arguments, Lyapunov-Schmidt reduction, and Leray-Schauder degree. Finally, we provide conditions under which these solutions necessarily change sign and we comment about further qualitative properties.

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 …

Squaring, Cubing, And Cube Rooting, Arthur T. Benjamin

All HMC Faculty Publications and Research

We present mentally efficient algorithms for mentally squaring and cubing 2-digit and 3-digit numbers and for finding cube roots of numbers with 2-digit or 3-digit answers.

Existence Of Solutions For A Semilinear Wave Equation With Non-Monotone Nonlinearity, Alfonso Castro, Benjamin Preskill '09

All HMC Faculty Publications and Research

For double-periodic and Dirichlet-periodic boundary conditions, we prove the existence of solutions to a forced semilinear wave equation with asymptotically linear nonlinearity, no resonance, and non-monotone nonlinearity when the forcing term is not flat on characteristics. The solutions are in L^∞ when the forcing term is in L^∞ and continous when the forcing term is continuous. This is in contrast with the results in [4], where the non-enxistence of continuous solutions is established even when forcing term is of class C^∞ but is flat on a characteristic.

Voting In Agreeable Societies, Deborah E. Berg '06, Serguei Norine, Francis E. Su, Robin Thomas, Paul Wollan

All HMC Faculty Publications and Research

No abstract provided in this article.

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

Turing Patterns On Growing Spheres: The Exponential Case, Julijana Gjorgjieva, Jon T. Jacobsen

All HMC Faculty Publications and Research

We consider Turing patterns for reaction-diffusion systems on the surface of a growing sphere. In particular, we are interested in the effect of dynamic growth on the pattern formation. We consider exponential isotropic growth of the sphere and perform a linear stability analysis and compare the results with numerical simulations.

Communicating Applied Mathematics: Four Examples, Daniel E. Finkel, Christopher Kuster, Matthew Lasater, Rachel Levy, Jill P. Reese, Ilse C. F. Ipsen

All HMC Faculty Publications and Research

Communicating Applied Mathematics is a writing- and speaking-intensive graduate course at North Carolina State University. The purpose of this article is to provide a brief description of the course objectives and the assignments. Parts A–D of of this article represent the class projects and illustrate the outcome of the course:

• The Evolution of an Optimization Test Problem: From Motivation to Implementation, by Daniel E. Finkel and Jill P. Reese

• Finding the Volume of a Powder from a Single Surface Height Measurement, by Christopher Kuster

• Finding Oscillations in Resonant Tunneling Diodes, by Matthew Lasater

• …

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 …

Lower Bounds For Simplicial Covers And Triangulations Of Cubes, Adam Bliss '03, Francis E. Su

All HMC Faculty Publications and Research

We show that the size of a minimal simplicial cover of a polytope P is a lower bound for the size of a minimal triangulation of P, including ones with extra vertices. We then use this fact to study minimal triangulations of cubes, and we improve lower bounds for covers and triangulations in dimensions 4 through at least 12 (and possibly more dimensions as well). Important ingredients are an analysis of the number of exterior faces that a simplex in the cube can have of a specified dimension and volume, and a characterization of corner simplices in terms of their …

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.

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.

An Existence Result For A Class Of Sublinear Semipositone Systems, Alfonso Castro, C. Maya, Ratnasingham Shivaji

All HMC Faculty Publications and Research

We consider the existence of positive solutions for the system

-Δu_i = λ[f_i(u₁,u₂,...,u_m) - h_i]; Ω

u_i = 0; ∂Ω

where λ > 0 is a parameter, Δ is the Laplacian operator, Ω is a bounded domain in Rⁿ; n ≥ 1 with a smooth boundary ∂Ω, f_i are C¹ functions satisfying f₁(0,0,...,0) = 0, lim z→∞ f_i(z,z,...,z) = ∞ and lim z→∞ f_i(z,z,...,z)/z = 0, and h_i are nonnegative continuous functions in Ω for i = 1,2,...,m. …

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 …

The Effect Of The Domain Topology On The Number Of Minimal Nodal Solutions Of An Elliptic Equation At Critical Growth In A Symmetric Domain, Alfonso Castro, Mónica Clapp

All HMC Faculty Publications and Research

We consider the Dirichlet problem Δu + λu + |u|^2*−2u = 0 in Ω, u = 0 on ∂Ω where Ω is a bounded smooth domain in R^N, N≥4, and 2* = 2N/(N−2) is the critical Sobolev exponent. We show that if Ω is invariant under an orthogonal involution then, for λ>0 sufficiently small, there is an effect of the equivariant topology of Ω on the number of solutions which change sign exactly once.

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.

Analysis Of The N-Card Version Of The Game Le Her, Arthur T. Benjamin, Alan J. Goldman

All HMC Faculty Publications and Research

We present a complete solution to a card game with historical origins. Our analysis exploits the convexity properties in the payoff matrix, allowing this discrete game to be resolved by continuous methods.

Modeling Control Of Hiv Infection Through Structured Treatment Interruptions With Recommendations For Experimental Protocol, Shannon Kubiak, Heather Lehr, Rachel Levy, Todd Moeller, Albert Parker, Edward Swim

All HMC Faculty Publications and Research

Highly Active Anti-Retroviral Therapy (HAART) of HIV infection has significantly reduced morbidity and mortality in developed countries. However, since these treatments can cause side effects and require strict adherence to treatment protocol, questions about whether or not treatment can be interrupted or discontinued with control of infection maintained by the host immune system remain to be answered. We present sensitivity analysis of a compartmental model for HIV infection that allows for treatment interruptions, including the sensitivity of the compartments themselves to our parameters as well as the sensitivity of the cost function used in parameter estimation. Recommendations are made about …

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

Stability Of Self-Similar Solutions For Van Der Waals Driven Thin Film Rupture, Thomas P. Witelski, Andrew J. Bernoff

All HMC Faculty Publications and Research

Recent studies of pinch-off of filaments and rupture in thin films have found infinite sets of first-type similarity solutions. Of these, the dynamically stable similarity solutions produce observable rupture behavior as localized, finite-time singularities in the models of the flow. In this letter we describe a systematic technique for calculating such solutions and determining their linear stability. For the problem of axisymmetric van der Waals driven rupture (recently studied by Zhang and Lister), we identify the unique stable similarity solution for point rupture of a thin film and an alternative mode of singularity formation corresponding to annular “ring rupture.”

Positive Solutions For A Concave Semipositone Dirichlet Problem, Alfonso Castro, Ratnasingham Shivaji

All HMC Faculty Publications and Research

No abstract provided for this article

On Multiple Solutions Of A Nonlinear Dirichlet Problem, Alfonso Castro, Jorge Cossio, John M. Neuberger

All HMC Faculty Publications and Research

We prove that a semilinear elliptic boundary value problem has five solutions when the range of the derivative of the nonlinearity includes at least the first two eigenvalues. We also prove that if the region is a ball the semilinear elliptic problem has two solutions that change sign and are nonradial.

Quasi-Steady Monopole And Tripole Attractors In Relaxing Vortices, Louis F. Rossi, Joseph F. Lingevitch, Andrew J. Bernoff

All HMC Faculty Publications and Research

Using fully nonlinear simulations of the two-dimensional Navier–Stokes equations at large Reynolds number (Re), we bracket a threshold amplitude above which a perturbed Gaussian monopole will relax to a quasi-steady, rotating tripole, and below which will relax to an axisymmetric monopole. The resulting quasi-steady structures are robust to small perturbations. We propose a means of measuring the decay rate of disturbances to asymptotic vortical structures wherein streamlines and lines of constant vorticity correspond in some rotating or translating frame. These experiments support the hypothesis that small or moderate deviations from asymptotic structures decay through inviscid and viscous mixing.

Positive Solution Curves Of Semipositone Problems With Concave Nonlinearities, Alfonso Castro, Sudhasree Gadam, Ratnasingham Shivaji

All HMC Faculty Publications and Research

We consider the positive solutions to the semilinear equation:

-Δu(x) = λf(u(x)) for x ∈ Ω

u(x) = 0 for x ∈ ∂Ω

where Ω denotes a smooth bounded region in R^N (N > 1) and λ > 0. Here f :[0, ∞)→R is assumed to be monotonically increasing, concave and such that f(0) < 0 (semipositone). Assuming that f'(∞) ≡ lim t→∞ f'(t) > 0, we establish the stability and uniqueness of large positive solutions in terms of (f(t)/t)'. When Ω is a ball, we determine the exact number of positive solutions for each λ > 0. We also obtain the geometry of the branches of positive solutions completely and establish how …