Physical Sciences and Mathematics Commons^™

Modeling An Infection Outbreak With Quarantine: The Sibkr Model, Mikenna Dew, Amanda Langosch, Theadora Baker-Wallerstein

Rose-Hulman Undergraduate Mathematics Journal

Influenza is a respiratory infection that places a substantial burden in the world population each year. In this project, we study and interpret a data set from a flu outbreak in a British boarding school in 1978 with mathematical modeling. First, we propose a generalization of the SIR model based on the quarantine measure in place and establish the long-time behavior of the model. By analyzing the model mathematically, we determine the analytic formulas of the basic reproduction number, the long-time limit of solutions, and the maximum number of infection population. Moreover, we estimate the parameters of the model based …

On Solutions Of First Order Pde With Two-Dimensional Dirac Delta Forcing Terms, Ian Robinson

Rose-Hulman Undergraduate Mathematics Journal

We provide solutions of a first order, linear partial differential equation of two variables where the nonhomogeneous term is a two-dimensional Dirac delta function. Our results are achieved by applying the unilateral Laplace Transform, solving the subsequently transformed PDE, and reverting back to the original space-time domain. A discussion of existence and uniqueness of solutions, a derivation of solutions of the PDE coupled with a boundary and initial condition, as well as a few worked examples are provided.

The Effect Of Habitat Fragmentation On Plant Communities In A Spatially-Implicit Grassland Model, Mika T. Cooney, Benjamin R. Hafner, Shelby E. Johnson, Sean Lee

Rose-Hulman Undergraduate Mathematics Journal

The spatially implicit Tilman-Levins ODE model helps to explain why so many plant species can coexist in grassland communities. This now-classic modeling framework assumes a trade-off between colonization and competition traits and predicts that habitat destruction can lead to long transient declines called ``extinction debts.'' Despite its strengths, the Tilman-Levins model does not explicitly account for landscape scale or the spatial configuration of viable habitat, two factors that may be decisive for population viability. We propose modifications to the model that explicitly capture habitat geometry and the spatial pattern of seed dispersal. The modified model retains implicit space and is …

Analysis Of A Quantum Attack On The Blum-Micali Pseudorandom Number Generator, Tingfei Feng

Mathematical Sciences Technical Reports (MSTR)

In 2012, Guedes, Assis, and Lula proposed a quantum attack on a pseudorandom number generator named the Blum-Micali Pseudorandom number generator. They claimed that the quantum attack can outperform classical attacks super-polynomially. However, this paper shows that the quantum attack cannot get the correct seed and provides another corrected algorithm that is in exponential time but still faster than the classical attack. Since the original classical attacks are in exponential time, the Blum-Micali pseudorandom number generator would be still quantum resistant.

Numerical Analysis Of A Model For The Growth Of Microorganisms, Alexander Craig Montgomery, Braden J. Carlson

Rose-Hulman Undergraduate Mathematics Journal

A system of first-order differential equations that arises in a model for the growth of microorganisms in a chemostat with Monod kinetics is studied. A new, semi-implicit numerical scheme is proposed to approximate solutions to the system. It is shown that the scheme is uniquely solvable and unconditionally stable, and further properties of the scheme are analyzed. The convergence rate of the numerical solution to the true solution of the system is given, and it is shown convergence of the numerical solutions to the true solutions is uniform over any interval [0, T ] for T > 0.

Using Differential Equations To Model A Cockatoo On A Spinning Wheel As Part Of The Scudem V Modeling Challenge, Miles Pophal, Chenming Zhen, Henry Bae

Rose-Hulman Undergraduate Mathematics Journal

For the SCUDEM V 2020 virtual challenge, we received an outstanding distinction for modeling a bird perched on a bicycle wheel utilizing the appropriate physical equations of rotational motion. Our model includes both theoretical calculations and numerical results from applying the Heaviside function for the swing motion of the bird. We provide a discussion on: our model and its numerical results, the overall limitations and future work of the model we constructed, and the experience we had participating in SCUDEM V 2020.

Implementation Of A Least Squares Method To A Navier-Stokes Solver, Jada P. Lytch, Taylor Boatwright, Ja'nya Breeden

Rose-Hulman Undergraduate Mathematics Journal

The Navier-Stokes equations are used to model fluid flow. Examples include fluid structure interactions in the heart, climate and weather modeling, and flow simulations in computer gaming and entertainment. The equations date back to the 1800s, but research and development of numerical approximation algorithms continues to be an active area. To numerically solve the Navier-Stokes equations we implement a least squares finite element algorithm based on work by Roland Glowinski and colleagues. We use the deal.II academic library , the C++ language, and the Linux operating system to implement the solver. We investigate convergence rates and apply the least squares …

The Primitive Root Problem: A Problem In Bqp, Shixin Wu

Mathematical Sciences Technical Reports (MSTR)

Shor’s algorithm proves that the discrete logarithm problem is in BQP. Based on his algorithm, we prove that the primitive root problem, a problem that verifies if some integer g is a primitive root modulo p where p is the largest prime number smaller than 2n for a given n, which is assumed to be harder than the discrete logarithm problem, is in BQP by using an oracle quantum Turing machine.

Structure Of Number Theoretic Graphs, Lee Trent

Mathematical Sciences Technical Reports (MSTR)

The tools of graph theory can be used to investigate the structure
imposed on the integers by various relations. Here we investigate two
kinds of graphs. The first, a square product graph, takes for its vertices
the integers 1 through n, and draws edges between numbers whose product
is a square. The second, a square product graph, has the same vertex set,
and draws edges between numbers whose sum is a square.
We investigate the structure of these graphs. For square product
graphs, we provide a rather complete characterization of their structure as
a union of disjoint complete graphs. For …

On The Consistency Of Alternative Finite Difference Schemes For The Heat Equation, Tran April

Rose-Hulman Undergraduate Mathematics Journal

While the well-researched Finite Difference Method (FDM) discretizes every independent variable into algebraic equations, Method of Lines discretizes all but one dimension, leaving an Ordinary Differential Equation (ODE) in the remaining dimension. That way, ODE's numerical methods can be applied to solve Partial Differential Equations (PDEs). In this project, Linear Multistep Methods and Method of Lines are used to numerically solve the heat equation. Specifically, the explicit Adams-Bashforth method and the implicit Backward Differentiation Formulas are implemented as Alternative Finite Difference Schemes. We also examine the consistency of these schemes.

Convergence Properties Of Solutions Of A Length-Structured Density-Dependent Model For Fish, Geigh Zollicoffer

Rose-Hulman Undergraduate Mathematics Journal

We numerically study solutions to a length-structured matrix model for fish populations in which the probability that a fish grows into the next length class is a decreasing nonlinear function of the total biomass of the population. We make conjectures about the convergence properties of solutions to this equation, and give numerical simulations which support these conjectures. We also study the distribution of biomass in the different age classes as a function of the total biomass.

Computer Program Simulation Of A Quantum Turing Machine With Circuit Model, Shixin Wu

Mathematical Sciences Technical Reports (MSTR)

Molina and Watrous present a variation of the method to simulate a quantum Turing machine employed in Yao’s 1995 publication “Quantum Circuit Complexity”. We use a computer program to implement their method with linear algebra and an additional unitary operator defined to complete the details. Their method is verified to be correct on a quantum Turing machine.

Compare And Contrast Maximum Likelihood Method And Inverse Probability Weighting Method In Missing Data Analysis, Scott Sun

Mathematical Sciences Technical Reports (MSTR)

Data can be lost for different reasons, but sometimes the missingness is a part of the data collection process. Unbiased and efficient estimation of the parameters governing the response mean model requires the missing data to be appropriately addressed. This paper compares and contrasts the Maximum Likelihood and Inverse Probability Weighting estimators in an Outcome-Dependendent Sampling design that deliberately generates incomplete observations. WE demonstrate the comparison through numerical simulations under varied conditions: different coefficient of determination, and whether or not the mean model is misspecified.

Numerical Integration Through Concavity Analysis, Daniel J. Pietz

Rose-Hulman Undergraduate Mathematics Journal

We introduce a relationship between the concavity of a C2 func- tion and the area bounded by its graph and secant line. We utilize this relationship to develop a method of numerical integration. We then bound the error of the approximation, and compare to known methods, finding an improvement in error bound over methods of comparable computational complexity.

On The Construction And Mathematical Analysis Of The Wavelet Transform And Its Matricial Properties, Diego Sejas Viscarra

Rose-Hulman Undergraduate Mathematics Journal

We study the properties of computational methods for the Wavelet Transform and its Inverse from the point of view of Linear Algebra. We present a characterization of such methods as matrix products, proving in particular that each iteration corresponds to the multiplication of an adequate unitary matrix. From that point we prove that some important properties of the Continuous Wavelet Transform, such as linearity, distributivity over matrix multiplication, isometry, etc., are inherited by these discrete methods.

This work is divided into four sections. The first section corresponds to the classical theoretical foundation of harmonic analysis with wavelets; it is used …

Dna Self-Assembly Design For Gear Graphs, Chiara Mattamira

Rose-Hulman Undergraduate Mathematics Journal

Application of graph theory to the well-known complementary properties of DNA strands has resulted in new insights about more efficient ways to form DNA nanostructures, which have been discovered as useful tools for drug delivery, biomolecular computing, and biosensors. The key concept underlying DNA nanotechnology is the formation of complete DNA complexes out of a given collection of branched junction molecules. These molecules can be modeled in the abstract as portions of graphs made up of vertices and half-edges, where complete edges are representations of double-stranded DNA pieces that have joined together. For efficiency, one aim is to minimize the …

Modeling Braids With Space-Varying And Time-Varying Stranded Cellular Automata, Brian Chan

Mathematical Sciences Technical Reports (MSTR)

Braids in a traditional sense and braids in a mathematical sense are wildly different outlooks on the same concept. Using cellular automata to represent and analyze braids is a way to bridge the gap between them. Joshua and Lana Holden and Hao Yang have previously worked on developing and expanding upon a Stranded Cellular Automata (SCA) model capable of representing many different braids and weaves. Continuing their work, we were able to devise a more user-friendly method for interacting with the model such that even those without a mathematical background can construct and analyze braids of their own. This paper …

The Game Of Life On The Hyperbolic Plane, Yuncong Gu

Mathematical Sciences Technical Reports (MSTR)

In this paper, we work on the Game of Life on the hyperbolic plane. We are interested in different tessellations on the hyperbolic plane and different Game of Life rules. First, we show the exponential growth of polygons on the pentagon tessellation. Moreover, we find that the Group of 3 can keep the boundary of a set not getting smaller. We generalize the existence of still lifes by computer simulations. Also, we will prove some propositions of still lifes and cycles. There exists a still life under rules B1, B2, and S3.

Using Differential Equations To Model Predator-Prey Relations As Part Of Scudem Modeling Challenge, Zachary Fralish, Bernard Tyson Iii, Anthony Stefan

Rose-Hulman Undergraduate Mathematics Journal

Differential equation modeling challenges provide students with an opportunity to improve their mathematical capabilities, critical thinking skills, and communication abilities through researching and presenting on a differential equations model. This article functions to display an archetype summary of an undergraduate student team’s response to a provided prompt. Specifically, the provided mathematical model estimates how certain stimuli from a predator are accumulated to trigger a neural response in a prey. Furthermore, it tracks the propagation of the resultant action potential and the physical flight of the prey from the predator through the analysis of larval zebrafish as a model organism. This …

Forward Selection Via Distance Correlation, Ty Adams

Mathematical Sciences Technical Reports (MSTR)

No abstract provided.

Periodicity And Invertibility Of Lattice Gas Cellular Automata, Jiawen Wang

Mathematical Sciences Technical Reports (MSTR)

A cellular automaton is a type of mathematical system that models the behavior of a set of cells with discrete values in progressing time steps. The often complicated behaviors of cellular automata are studied in computer science, mathematics, biology, and other science related fields. Lattice gas cellular automata are used to simulate the movements of particles. This thesis aims to discuss the properties of lattice gas models, including periodicity and invertibility, and to examine their accuracy in reflecting the physics of particles in real life. Analysis of elementary cellular automata is presented to introduce the concept of cellular automata and …

Algorithms To Approximate Solutions Of Poisson's Equation In Three Dimensions, Ray Dambrose

Rose-Hulman Undergraduate Mathematics Journal

The focus of this research was to develop numerical algorithms to approximate solutions of Poisson's equation in three dimensional rectangular prism domains. Numerical analysis of partial differential equations is vital to understanding and modeling these complex problems. Poisson's equation can be approximated with a finite difference approximation. A system of equations can be formed that gives solutions at internal points of the domain. A computer program was developed to solve this system with inputs such as boundary conditions and a nonhomogenous source function. Approximate solutions are compared with exact solutions to prove their accuracy. The program is tested with an …

Mathematical Modeling Of Quadcopter Dynamics, Qikai Huang (Bruce Wingo)

Mathematical Sciences Technical Reports (MSTR)

Recently, Google, Amazon and others are attempting to develop delivery drones for commercial use, in particular Amazon Prime Air promising 30 minute delivery. One type of commonly used drone proposed for such purposes is a quadcopter. Quadcopters have been around for some time with original development in the 1920’s. They are popular now because they are mechanically simple and provide a good vehicle for unmanned flight. By controlling the speed of the four propellers, a quadcopter can roll, change pitch, change yaw, and accelerate. This research will focus on the study of classical mechanics theories on rigid body motion using …

Mathematical Modeling Of Quadcopter Dynamics, Qikai Huang (Bruce Wingo)

Rose-Hulman Undergraduate Research Publications

Recently, Google, Amazon and others are attempting to develop delivery drones for commercial use, in particular Amazon Prime Air promising 30 minute delivery. One type of commonly used drone proposed for such purposes is a quadcopter. Quadcopters have been around for some time with original development in the 1920’s. They are popular now because they are mechanically simple and provide a good vehicle for unmanned flight. By controlling the speed of the four propellers, a quadcopter can roll, change pitch, change yaw, and accelerate. This research will focus on the study of classical mechanics theories on rigid body motion using …

Variance Of Clusterings On Graphs, Thomas Vlado Mulc

Mathematical Sciences Technical Reports (MSTR)

Graphs that represent data often have structures or characteristics that can represent some relationships in the data. One of these structures is clusters or community structures. Most clustering algorithms for graphs are deterministic, which means they will output the same clustering each time. We investigated a few stochastic algorithms, and look into the consistency of their clusterings.

Spontaneous Synchrony On Graphs And The Emergence Of Order From Disorder, Dylan Linville, Daniel Trugillo Martins Fontes

Mathematical Sciences Technical Reports (MSTR)

From pulsars to pedestrians and bacteria to brain cells, objects that exhibit cyclical behavior, called oscillators, are found in a variety of different settings. When oscillators adjust their behavior in response to nearby oscillators, they often achieve a state of synchrony, in which they all have the same phase and frequency. Here, we explore the Kuramoto model, a simple and general model which describes oscillators as dynamical systems on a graph and has been used to study synchronization in systems ranging from firefly swarms to the power grid. We discuss analytical and numerical methods used to investigate the governing system …

Spontaneous Synchrony On Graphs And The Emergence Of Order From Disorder, Dylan Linville, Daniel Trugillo Martins Fontes

Rose-Hulman Undergraduate Research Publications

From pulsars to pedestrians and bacteria to brain cells, objects that exhibit cyclical behavior, called oscillators, are found in a variety of different settings. When oscillators adjust their behavior in response to nearby oscillators, they often achieve a state of synchrony, in which they all have the same phase and frequency. Here, we explore the Kuramoto model, a simple and general model which describes oscillators as dynamical systems on a graph and has been used to study synchronization in systems ranging from firefly swarms to the power grid. We discuss analytical and numerical methods used to investigate the governing system …

Analytical Solution Of The Symmetric Circulant Tridiagonal Linear System, Sean A. Broughton, Jeffery J. Leader

Mathematical Sciences Technical Reports (MSTR)

A circulant tridiagonal system is a special type of Toeplitz system that appears in a variety of problems in scientific computation. In this paper we give a formula for the inverse of a symmetric circulant tridiagonal matrix as a product of a circulant matrix and its transpose, and discuss the utility of this approach for solving the associated system.

Discovering Exoplanets Through Hidden Markov Model Analysis, Jon Drobny

Rose-Hulman Undergraduate Research Publications

The goal for the project is to develop a Hidden Markov Model for the detection and characterization of extrasolar planets through the analysis of light curves.

Thermal Detection Of Inaccessible Corrosion, Matthew Charnley, Andrew Rzeznik

Mathematical Sciences Technical Reports (MSTR)

In this paper, we explore the mathematical inverse problem of detecting corroded material on the reverse side of a partially accessible metal plate. We will show how a linearization can be used to simplify the initial problem and explain a regularization method used to obtain acceptable results for the corrosion profile. We will also state and perform some calculations for the full three-dimensional problem for possible future work.