Ordinary Differential Equations and Applied Dynamics Commons^™

Analytical And Numerical Analysis Of The Sirs Model, Catherine Nguyen

Student Research Submissions

Mathematical models in epidemiology describe how diseases affect and spread within a population. By understanding the trends of a disease, more effective public health policies can be made. In this paper, the Susceptible-Infected-Recovered-Susceptible (SIRS) Model was examined analytically and numerically to compare with the data for Coronavirus Disease 2019 (COVID-19). Since the SIRS model is a complex model, analytical techniques were used to solve simplified versions of the SIRS model in order to understand general trends that occur. Then by Euler's Method, the Runge-Kutta Method, and the Predictor-Corrector Method, computational approximations were obtained to solve and plot the SIRS model. …

Proof-Of-Concept For Converging Beam Small Animal Irradiator, Benjamin Insley

Dissertations & Theses (Open Access)

The Monte Carlo particle simulator TOPAS, the multiphysics solver COMSOL., and

several analytical radiation transport methods were employed to perform an in-depth proof-ofconcept

for a high dose rate, high precision converging beam small animal irradiation platform.

In the first aim of this work, a novel carbon nanotube-based compact X-ray tube optimized for

high output and high directionality was designed and characterized. In the second aim, an

optimization algorithm was developed to customize a collimator geometry for this unique Xray

source to simultaneously maximize the irradiator’s intensity and precision. Then, a full

converging beam irradiator apparatus was fit with a multitude …

Mathematical Analysis Of Eukaryotic Pericentromere, Puranjan Ghimire

Theses and Dissertations

The centromere is crucial for chromosomal stability and their proper segregation during cell division in eukaryotes. Surrounding the centromere are pericentromeres, made of repetitive DNA elements called pericentromeric repeats, varying from 10 in fission yeast to thousands in humans. These repeats form densely packed heterochromatin, where genes are usually silenced. The silencing mechanism across different pericentromeric repeats remains unclear.

Despite variations in sequence and length, pericentromeric repeats are conserved across eukaryotes, indicating their functional importance. This dissertation presents mathematical models to quantify gene silencing in fission yeast and humans. In fission yeast, my model predicts that silencing occurs only with …

Multiscale Modelling Of Brain Networks And The Analysis Of Dynamic Processes In Neurodegenerative Disorders, Hina Shaheen

Theses and Dissertations (Comprehensive)

The complex nature of the human brain, with its intricate organic structure and multiscale spatio-temporal characteristics ranging from synapses to the entire brain, presents a major obstacle in brain modelling. Capturing this complexity poses a significant challenge for researchers. The complex interplay of coupled multiphysics and biochemical activities within this intricate system shapes the brain's capacity, functioning within a structure-function relationship that necessitates a specific mathematical framework. Advanced mathematical modelling approaches that incorporate the coupling of brain networks and the analysis of dynamic processes are essential for advancing therapeutic strategies aimed at treating neurodegenerative diseases (NDDs), which afflict millions of …

Controlled Manipulation And Transport By Microswimmers In Stokes Flows, Jake Buzhardt

All Dissertations

Remotely actuated microscale swimming robots have the potential to revolutionize many aspects of biomedicine. However, for the longterm goals of this field of research to be achievable, it is necessary to develop modelling, simulation, and control strategies which effectively and efficiently account for not only the motion of individual swimmers, but also the complex interactions of such swimmers with their environment including other nearby swimmers, boundaries, other cargo and passive particles, and the fluid medium itself. The aim of this thesis is to study these problems in simulation from the perspective of controls and dynamical systems, with a particular focus …

Rigid Body Constrained Motion Optimization And Control On Lie Groups And Their Tangent Bundles, Brennan S. Mccann

Doctoral Dissertations and Master's Theses

Rigid body motion requires formulations where rotational and translational motion are accounted for appropriately. Two Lie groups, the special orthogonal group SO(3) and the space of quaternions H, are commonly used to represent attitude. When considering rigid body pose, that is spacecraft position and attitude, the special Euclidean group SE(3) and the space of dual quaternions DH are frequently utilized. All these groups are Lie groups and Riemannian manifolds, and these identifications have profound implications for dynamics and controls. The trajectory optimization and optimal control problem on Riemannian manifolds presents significant opportunities for theoretical development. Riemannian optimization is an attractive …

Dynamical Aspects In (4+1)-Body Problems, Ryan Gauthier

Theses and Dissertations (Comprehensive)

The n-body problem models a system of n-point masses that attract each other via some binary interaction. The (n + 1)-body problem assumes that one of the masses is located at the origin of the coordinate system. For example, an (n+1)-body problem is an ideal model for Saturn, seen as the central mass, and one of its outer rings. A relative equilibrium (RE) is a special solution of the (n+1)-body problem where the non-central bodies rotate rigidly about the centre of mass. In rotating coordinates, these solutions become equilibria.

In this thesis we study dynamical aspects of planar (4 + …

A Novel Method For Sensitivity Analysis Of Time-Averaged Chaotic System Solutions, Christian A. Spencer-Coker

Theses and Dissertations

The direct and adjoint methods are to linearize the time-averaged solution of bounded dynamical systems about one or more design parameters. Hence, such methods are one way to obtain the gradient necessary in locally optimizing a dynamical system’s time-averaged behavior over those design parameters. However, when analyzing nonlinear systems whose solutions exhibit chaos, standard direct and adjoint sensitivity methods yield meaningless results due to time-local instability of the system. The present work proposes a new method of solving the direct and adjoint linear systems in time, then tests that method’s ability to solve instances of the Lorenz system that exhibit …

A Weak Fractional Calculus Theory And Numerical Methods For Fractional Differential Equations, Mitchell D. Sutton

Doctoral Dissertations

This dissertation is comprised of four integral parts. The first part comprises a self-contained new theory of weak fractional differential calculus in one-dimension. The crux of this new theory is the introduction of a weak fractional derivative notion which is a natural generalization of integer order weak derivatives; it also helps to unify multiple existing fractional derivative definitions.

The second part of this work presents three new families of fractional Sobolev spaces and their accompanying theory in one-dimension. The new construction and theory are based on a newly developed notion of weak fractional derivatives, which are natural generalizations of the …

Sensitivity Analysis Of Basins Of Attraction For Nelder-Mead, Sonia K. Shah

Honors Projects

The Nelder-Mead optimization method is a numerical method used to find the minimum of an objective function in a multidimensional space. In this paper, we use this method to study functions - specifically functions with three-dimensional graphs - and create images of the basin of attraction of the function. Three different methods are used to create these images named the systematic point method, randomized centroid method, and systemized centroid method. This paper applies these methods to different functions. The first function has two minima with an equivalent function value. The second function has one global minimum and one local minimum. …

Role Of Inhibition And Spiking Variability In Ortho- And Retronasal Olfactory Processing, Michelle F. Craft

Theses and Dissertations

Odor perception is the impetus for important animal behaviors, most pertinently for feeding, but also for mating and communication. There are two predominate modes of odor processing: odors pass through the front of nose (ortho) while inhaling and sniffing, or through the rear (retro) during exhalation and while eating and drinking. Despite the importance of olfaction for an animal’s well-being and specifically that ortho and retro naturally occur, it is unknown whether the modality (ortho versus retro) is transmitted to cortical brain regions, which could significantly instruct how odors are processed. Prior imaging studies show different …

Stroke Clustering And Fitting In Vector Art, Khandokar Shakib

Senior Independent Study Theses

Vectorization of art involves turning free-hand drawings into vector graphics that can be further scaled and manipulated. In this paper, we explore the concept of vectorization of line drawings and study multiple approaches that attempt to achieve this in the most accurate way possible. We utilize a software called StrokeStrip to discuss the different mathematics behind the parameterization and fitting involved in the drawings.

Numerical Study Of The Seiqr Model For Covid-19, Caitlin Holt

Student Research Submissions

In this research project, we used numerical methods to investigate trends in the susceptible, exposed, infectious, quarantined, recovered, closed cases and insusceptible populations for the COVID-19 pandemic in 2021. We used the SEIQR model containing seven ordinary differential equations, based on the SIR model for epidemics. An analytical solution was derived from a simplified version of the model, created by making various assumptions about the original model. Numerical solutions were generated for the first 100 days of 2021 using algorithms based on Euler's Method, Runge-Kutta Method, and Multistep Methods. Our goal is to show that numerical methods can help us …

Dynamic Parameter Estimation From Partial Observations Of The Lorenz System, Eunice Ng

Theses and Dissertations

Recent numerical work of Carlson-Hudson-Larios leverages a nudging-based algorithm for data assimilation to asymptotically recover viscosity in the 2D Navier-Stokes equations as partial observations on the velocity are received continuously-in-time. This "on-the-fly" algorithm is studied both analytically and numerically for the Lorenz equations in this thesis.

An Examination Of Fontan Circulation Using Differential Equation Models And Numerical Methods, Vanessa Maybruck

Honors Student Research

Certain congenital heart defects can lead to the development of only a single pumping chamber, or ventricle, in the heart instead of the usual two ventricles. Individuals with this defect undergo a corrective, three-part surgery, the third step of which is the Fontan procedure, but as the patients age, their cardiovascular health will likely deteriorate. Using computational fluid dynamics and differential equations, Fontan circulation can be modeled to investigate why the procedure fails and how Fontan failure can be maximally prevented. Borrowing from well-established literature on RC circuits, the differential equation models simulate systemic blood flow in a piecewise, switch-like …

High-Order Flexible Multirate Integrators For Multiphysics Applications, Rujeko Chinomona

Mathematics Theses and Dissertations

Traditionally, time integration methods within multiphysics simulations have been chosen to cater to the most restrictive dynamics, sometimes at a great computational cost. Multirate integrators accurately and efficiently solve systems of ordinary differential equations that exhibit different time scales using two or more time steps. In this thesis, we explore three classes of time integrators that can be classified as one-step multi-stage multirate methods for which the slow dynamics are evolved using a traditional one step scheme and the fast dynamics are solved through a sequence of modified initial value problems. Practically, the fast dynamics are subcycled using a small …

Lexicographic Sensitivity Functions For Nonsmooth Models In Mathematical Biology, Matthew D. Ackley

Electronic Theses and Dissertations

Systems of ordinary differential equations (ODEs) may be used to model a wide variety of real-world phenomena in biology and engineering. Classical sensitivity theory is well-established and concerns itself with quantifying the responsiveness of such models to changes in parameter values. By performing a sensitivity analysis, a variety of insights can be gained into a model (and hence, the real-world system that it represents); in particular, the information gained can uncover a system's most important aspects, for use in design, control or optimization of the system. However, while the results of such analysis are desirable, the approach that classical theory …

Buckling Loads Of A Graphene Layer Interacting With Rigid Substrates, Bradley Beckwith

Williams Honors College, Honors Research Projects

The goal of this project is to formulate a model that can predict the buckling of a graphene layer between two rigid substrates. The model will predict the buckling of the graphene layer when it is parallel to the substrates and an edge load is applied to the ends of the layer. Our main focus is to use the model to predict buckling loads given different assumptions for interaction forces between the graphene layer and the substrates. For this project continuum modeling will be used to create a model for the graphene buckling problem. This modeling leads to a total …

Modelling The Transition From Homogeneous To Columnar States In Locust Hopper Bands, Miguel Velez

HMC Senior Theses

Many biological systems form structured swarms, for instance in locusts, whose swarms are known as hopper bands. There is growing interest in applying mathematical models to understand the emergence and dynamics of these biological and social systems. We model the locusts of a hopper band as point particles interacting through repulsive and attractive social "forces" on a one dimensional periodic domain. The primary goal of this work is to modify this well studied modelling framework to be more biological by restricting repulsion to act locally between near neighbors, while attraction acts globally between all individuals. This is a biologically motivated …

Sum Of Cubes Of The First N Integers, Obiamaka L. Agu

Electronic Theses, Projects, and Dissertations

In Calculus we learned that 􏰅Sum^{n}_{k=1} k = [n(n+1)]/2 , that Sum^{􏰅n}_{k=1} k^2 = [n(n+1)(2n+1)]/6 , and that Sum^{n}_{k=1} k^{3} = (n(n+1)/2)^{2}. These formulas are useful when solving for the area below quadratic or cubic function over an interval [a, b]. This tedious process, solving for areas under a quadratic or a cubic, served as motivation for the introduction of Riemman integrals. For the overzealous math student, these steps were replaced by a simpler method of evaluating antiderivatives at the endpoints a and b. From my recollection, a former instructor informed us to do the value of memorizing these formulas. …

An Epidemiological Model With Simultaneous Recoveries, Ariel B. Farber

Electronic Theses and Dissertations

Epidemiological models are an essential tool in understanding how infection spreads throughout a population. Exploring the effects of varying parameters provides insight into the driving forces of an outbreak. In this thesis, an SIS (susceptible-infectious-susceptible) model is built partnering simulation methods, differential equations, and transition matrices with the intent to describe how simultaneous recoveries influence the spread of a disease in a well-mixed population. Individuals in the model transition between only two states; an individual is either susceptible — able to be infected, or infectious — able to infect others. Events in this model (infections and recoveries) occur by way …

Mathematical Models: The Lanchester Equations And The Zombie Apocalypse, Hailey Bauer

Undergraduate Theses and Capstone Projects

This research study used mathematical models to analyze and depicted specific battle situations and the outcomes of the zombie apocalypse. The original models that predicted warfare were the Lanchester models, while the zombie apocalypse models were fictional expansions upon mathematical models used to examine infectious diseases. In this paper, I analyzed and compared different mathematical models by examining each model’s set of assumptions and the impact of the change in variables on the population classes. The purpose of this study was to understand the basics of the discrete dynamical systems and to determine the similarities between imaginary and realistic models. …

Risk Assessment Of Dropped Cylindrical Objects In Offshore Operations, Adelina Steven

University of New Orleans Theses and Dissertations

Dropped object are defined as any object that fall under its own weight from a previously static position or fell due to an applied force from equipment or a moving object. It is among the top ten causes of injuries and fatality in oil and gas industry. To solve this problem, several in-house tools and guidelines is developed over time to assess the risk of dropped objects on the sub-sea structures. This thesis focuses on compiling and comparing those methods in hope to improve the recommended practices available in the market. A simple modification is done on the in-house tools …

Rotordynamic Analysis Of Theoretical Models And Experimental Systems, Cameron R. Naugle

Master's Theses

This thesis is intended to provide fundamental information for the construction and

analysis of rotordynamic theoretical models, and their comparison the experimental

systems. Finite Element Method (FEM) is used to construct models using Timoshenko

beam elements with viscous and hysteretic internal damping. Eigenvalues

and eigenvectors of state space equations are used to perform stability analysis, produce

critical speed maps, and visualize mode shapes. Frequency domain analysis

of theoretical models is used to provide Bode diagrams and in experimental data

full spectrum cascade plots. Experimental and theoretical model analyses are used

to optimize the control algorithm for an Active Magnetic Bearing …

Theoretical Analysis Of Nonlinear Differential Equations, Emily Jean Weymier

Electronic Theses and Dissertations

Nonlinear differential equations arise as mathematical models of various phenomena. Here, various methods of solving and approximating linear and nonlinear differential equations are examined. Since analytical solutions to nonlinear differential equations are rare and difficult to determine, approximation methods have been developed. Initial and boundary value problems will be discussed. Several linear and nonlinear techniques to approximate or solve the linear or nonlinear problems are demonstrated. Regular and singular perturbation theory and Magnus expansions are our particular focus. Each section offers several examples to show how each technique is implemented along with the use of visuals to demonstrate the accuracy, …

A Model To Predict Concentrations And Uncertainty For Mercury Species In Lakes, Ashley Hendricks

Dissertations, Master's Theses and Master's Reports

To increase understanding of mercury cycling, a seasonal mass balance model was developed to predict mercury concentrations in lakes and fish. Results indicate that seasonality in mercury cycling is significant and is important for a northern latitude lake. Models, when validated, have the potential to be used as an alternative to measurements; models are relatively inexpensive and are not as time intensive. Previously published mercury models have neglected to perform a thorough validation. Model validation allows for regulators to be able to make more informed, confident decisions when using models in water quality management. It is critical to quantify uncertainty; …

Understanding The Nature Of Nanoscale Wetting Through All-Atom Simulations, Oliver Evans

Williams Honors College, Honors Research Projects

The spreading behavior of spherical and cylindrical water droplets between 30Å and 100Å in radius on a sapphire surface is investigated using all-atom molecular dynamics simulations for durations on the order of tens of nanoseconds. A monolayer film develops rapidly and wets the surface, while the bulk of the droplet spreads on top of the monolayer, maintaining the shape of a spherical cap. Unlike previous simulations in the literature, the bulk radius is found to increase to a maximum value and receed as the monolayer continues to expand. Simple time and droplet size dependence is observed for monolayer radius and …

High-Order Relaxed Multirate Infinitesimal Step Methods For Multiphysics Applications, Jean Sexton

Mathematics Theses and Dissertations

In this work, we consider numerical methods for integrating multirate ordinary differential equations. We are interested in the development of new multirate methods with good stability properties and improved efficiency over existing methods. We discuss the development of multirate methods, particularly focusing on those that are based on Runge-Kutta theory. We introduce the theory of Generalized Additive Runge-Kutta methods proposed by Sandu and Günther. We also introduce the theory of Recursive Flux Splitting Multirate Methods with Sub-cycling described by Schlegel, as well as the Multirate Infinitesimal Step methods this work is based on. We propose a generic structure called Flexible …

On The Ramberg-Osgood Stress-Strain Model And Large Deformations Of Cantilever Beams, Ronald J. Giardina Jr

University of New Orleans Theses and Dissertations

In this thesis the Ramberg-Osgood nonlinear model for describing the behavior of many different materials is investigated. A brief overview of the model as it is currently used in the literature is undertaken and several misunderstandings and possible pitfalls in its application is pointed out, especially as it pertains to more recent approaches to finding solutions involving the model. There is an investigation of the displacement of a cantilever beam under a combined loading consisting of a distributed load across the entire length of the beam and a point load at its end and new solutions to this problem are …

Application Of Symplectic Integration On A Dynamical System, William Frazier

Electronic Theses and Dissertations

Molecular Dynamics (MD) is the numerical simulation of a large system of interacting molecules, and one of the key components of a MD simulation is the numerical estimation of the solutions to a system of nonlinear differential equations. Such systems are very sensitive to discretization and round-off error, and correspondingly, standard techniques such as Runge-Kutta methods can lead to poor results. However, MD systems are conservative, which means that we can use Hamiltonian mechanics and symplectic transformations (also known as canonical transformations) in analyzing and approximating solutions. This is standard in MD applications, leading to numerical techniques known as symplectic …