Ordinary Differential Equations and Applied Dynamics Commons^™

Type Ii Diabetes Treatment Comparison Via Compartment Modeling, Abigail M. Collins

Theses and Dissertations

Type II diabetes mellitus affects one in ten adults worldwide, yet the effects of treatment type and adherence level on developing complications and quality of life have not been well characterized at the population level, and mathematical modeling offers a structured way to examine these dynamics. This thesis adapts the Boutayeb et al. (2004) model to incorporate dynamic treatment types and levels of adherence, producing nine scenarios in which complication development rate and complication recovery rate differed, to compare peak complications and quality of life across treatment and adherence conditions. Using a system of ordinary differential equations and compartment modeling, …

Managing Multi-Drug Resistance: An Evolutionary Game Theory And Optimal Control Approach, Shukhrat Nasrulloev

Theses and Dissertations

Multi-drug resistance is an evolutionary process in which treatment eliminates sensitive cells, allowing resistant clones to dominate. This thesis investigates this process using a framework integrating population dynamics, evolutionary game theory, and optimal control theory. We develop a two-population logistic growth model describing competition between drug-sensitive and drug-resistant cells under treatment, construct dose-dependent payoff matrices and replicator dynamics to characterize evolutionary competition, and derive a critical drug level D_crit = (r_S - rR)/(dS - d_R) at which resistant cells gain a fitness advantage. An optimal control problem is formulated via Pontryagin's Maximum Principle to identify schedules …

Mathematical Modeling And Analysis Of Inflammation And Tissue Repair: Lung Inflammation And Wound Healing In Corals Under Stress, Quintessa Hay

Theses and Dissertations

A variety of insults, including tissue injury and/or exposure to pathogen, elicit an immune response in many organisms. An improperly regulated immune response can result in deleterious effects to the organism. Here we present models for lung injury in young and old mice and models for wound healing in coral reefs.

It is well known that the immune response becomes less effective in older individuals. This is of particular interest in pulmonary insults such as ventilator induced lung injury (VILI) or lung infection. We extended a mathematical model for the inflammatory response to VILI and used experimental data to select …

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 …

Mathematical Modeling Of Phage-Bacteria Population Dynamics, John Lawrence D. Palacios

Theses and Dissertations

Bacteriophages are viruses that infect and replicate within bacteria. Lytic phages cause the bacterial cell to burst, killing the bacteria. These types of phages can be used to treat patients with antibiotic-resistant bacterial infections. As a step in developing successful treatment protocols, we aim to understand the population dynamics of phages and bacteria using an in vitro model. We model the dynamics using the Campbell model, which consists of a delay differential equation (DDE), as a base model. We extended the model by including the emergence of phage resistance. We then compared the DDE model with a parallel ordinary differential …

Innovations In Drop Shape Analysis Using Deep Learning And Solving The Young-Laplace Equation For An Axisymmetric Pendant Drop, Andres P. Hyer

Theses and Dissertations

Axisymmetric Drop Shape Analysis (ADSA) is a technique commonly used to determine surface or interfacial tension. Applications of traditional ASDA methods to process analytical technologies are limited by computational speed and image quality. Here, we address these limitations using a novel machine learning approach to analysis. With a convolutional neural network (CNN), we were able to achieve an experimental fit precision of (+/-) 0.122 mN/m in predicting the surface tension of drop images at a rate of 1.5 ms^-1 versus 7.7 s^-1, which is more than 5,000 times faster than the traditional method. The results are validated on real images …

A Novel Chebyshev Wavelet Method For Solving Fractional-Order Optimal Control Problems, Ghodsieh Ghanbari

Theses and Dissertations

This thesis presents a numerical approach based on generalized fractional-order Chebyshev wavelets for solving fractional-order optimal control problems. The exact value of the Riemann– Liouville fractional integral operator of the generalized fractional-order Chebyshev wavelets is computed by applying the regularized beta function. We apply the given wavelets, the exact formula, and the collocation method to transform the studied problem into a new optimization problem. The convergence analysis of the proposed method is provided. The present method is extended for solving fractional-order, distributed-order, and variable-order optimal control problems. Illustrative examples are considered to show the advantage of this method in comparison …

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 …

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 …

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.

The Analysis Of Neural Heterogeneity Through Mathematical And Statistical Methods, Kyle Wendling

Theses and Dissertations

Diversity of intrinsic neural attributes and network connections is known to exist in many areas of the brain and is thought to significantly affect neural coding. Recent theoretical and experimental work has argued that in uncoupled networks, coding is most accurate at intermediate levels of heterogeneity. I explore this phenomenon through two distinct approaches: a theoretical mathematical modeling approach and a data-driven statistical modeling approach.

Through the mathematical approach, I examine firing rate heterogeneity in a feedforward network of stochastic neural oscillators utilizing a high-dimensional model. The firing rate heterogeneity stems from two sources: intrinsic (different individual cells) and network …

Delay Differential Equations And Their Application To Micro Electro Mechanical Systems, Asset Ospanov

Theses and Dissertations

Delay differential equations have a wide range of applications in engineering. This work is devoted to the analysis of delay Duffing equation, which plays a crucial role in modeling performance on demand Micro Electro Mechanical Systems (MEMS). We start with the stability analysis of a linear delay model. We also show that in certain cases the delay model can be efficiently approximated with a much simpler model without delay. We proceed with the analysis of a non-linear Duffing equation. This model is a significantly more complex mathematical model. For instance, the existence of a periodic solution for this equation is …

An Applied Mathematics Approach To Modeling Inflammation: Hematopoietic Bone Marrow Stem Cells, Systemic Estrogen And Wound Healing And Gas Exchange In The Lungs And Body, Racheal L. Cooper

Theses and Dissertations

Mathematical models apply to a multitude physiological processes and are used to make predictions and analyze outcomes of these processes. Specifically, in the medical field, a mathematical model uses a set of initial conditions that represents a physiological state as input and a set of parameter values are used to describe the interaction between variables being modeled. These models are used to analyze possible outcomes, and assist physicians in choosing the most appropriate treatment options for a particular situation. We aim to use mathematical modeling to analyze the dynamics of processes involved in the inflammatory process.

First, we create a …

A Comparison Of Obesity Interventions Using Energy Balance Models, Marcella Torres

Theses and Dissertations

An energy balance model of human metabolism developed by Hall et al. is extended to compare body composition outcomes among standard and proposed obesity interventions. Standard interventions include a drastic diet or a drastic diet with endurance training. Outcomes for these interventions are typically poor in clinical studies. Proposed interventions include a gradual diet and the addition of resistance training to preserve lean mass and metabolic rate. We see that resistance training, regardless of dietary strategy, achieves these goals. Finally, we observe that the optimal obesity intervention for continued maintenance of a healthy body composition following a diet includes a …

Development And Validation Of Reentry Simulation Using Matlab, Robert E. Jameson Jr.

Theses and Dissertations

This research effort develops a program using MATLAB to solve the equations of motion for atmospheric reentry and analyzes the validity of the program for use as a tool to expeditiously predict reentry profiles. The reentry vehicle is modeled as a point mass with constant aerodynamic properties as defined by the user. The equations of motion for reentry are based on the two-body problem. The atmosphere is modeled as a single layer exponentially decreasing in density. The MATLAB program has the ability to derive the initial trajectory conditions from the position and velocity relative to the rotating surface of the …

A Sampling And Transformation Approach To Solving Random Differential Equations, Roger A. Erich

Theses and Dissertations

This research explores an innovative sampling method used to conduct uncertainty analysis on a system with one random input. Given the distribution of the random input, X, we seek to find the distribution of the output random variable Y. When the functional form of the transformation Y=g(X) is not explicitly known, complicated procedures, such as stochastic projection or Monte Carlo simulation must be employed. The main focus of this research is determining the distribution of the random variable Y=g(X) where g(X) is the solution to an ordinary differential equation and X is a random parameter. Here, y=g(X) is approximated by …