Matching Items (28)
Filtering by
- Member of: ASU Electronic Theses and Dissertations
Description
Advances in experimental techniques have allowed for investigation of molecular dynamics at ever smaller temporal and spatial scales. There is currently a varied and growing body of literature which demonstrates the phenomenon of \emph{anomalous diffusion} in physics, engineering, and biology. In particular many diffusive type processes in the cell have been observed to follow a power law $\left \propto t^\alpha$ scaling of the mean square displacement of a particle. This contrasts with the expected linear behavior of particles undergoing normal diffusion. \emph{Anomalous sub-diffusion} ($\alpha<1$) has been attributed to factors such as cytoplasmic crowding of macromolecules, and trap-like structures in the subcellular environment non-linearly slowing the diffusion of molecules. Compared to normal diffusion, signaling molecules in these constrained spaces can be more concentrated at the source, and more diffuse at longer distances, potentially effecting the signalling dynamics. As diffusion at the cellular scale is a fundamental mechanism of cellular signaling and additionally is an implicit underlying mathematical assumption of many canonical models, a closer look at models of anomalous diffusion is warranted. Approaches in the literature include derivations of fractional differential diffusion equations (FDE) and continuous time random walks (CTRW). However these approaches are typically based on \emph{ad-hoc} assumptions on time- and space- jump distributions. We apply recent developments in asymptotic techniques on collisional kinetic equations to develop a FDE model of sub-diffusion due to trapping regions and investigate the nature of the space/time probability distributions assosiated with trapping regions. This approach both contrasts and compliments the stochastic CTRW approach by positing more physically realistic underlying assumptions on the motion of particles and their interactions with trapping regions, and additionally allowing varying assumptions to be applied individually to the traps and particle kinetics.
ContributorsHoleva, Thomas Matthew (Author) / Ringhofer, Christian (Thesis advisor) / Baer, Steve (Thesis advisor) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Taylor, Jesse (Committee member) / Arizona State University (Publisher)
Created2014
Description
Numerical simulations are very helpful in understanding the physics of the formation of structure and galaxies. However, it is sometimes difficult to interpret model data with respect to observations, partly due to the difficulties and background noise inherent to observation. The goal, here, is to attempt to bridge this gap between simulation and observation by rendering the model output in image format which is then processed by tools commonly used in observational astronomy. Images are synthesized in various filters by folding the output of cosmological simulations of gasdynamics with star-formation and dark matter with the Bruzual- Charlot stellar population synthesis models. A variation of the Virgo-Gadget numerical simulation code is used with the hybrid gas and stellar formation models of Springel and Hernquist (2003). Outputs taken at various redshifts are stacked to create a synthetic view of the simulated star clusters. Source Extractor (SExtractor) is used to find groupings of stellar populations which are considered as galaxies or galaxy building blocks and photometry used to estimate the rest frame luminosities and distribution functions. With further refinements, this is expected to provide support for missions such as JWST, as well as to probe what additional physics are needed to model the data. The results show good agreement in many respects with observed properties of the galaxy luminosity function (LF) over a wide range of high redshifts. In particular, the slope (alpha) when fitted to the standard Schechter function shows excellent agreement both in value and evolution with redshift, when compared with observation. Discrepancies of other properties with observation are seen to be a result of limitations of the simulation and additional feedback mechanisms which are needed.
ContributorsMorgan, Robert (Author) / Windhorst, Rogier A (Thesis advisor) / Scannapieco, Evan (Committee member) / Rhoads, James (Committee member) / Gardner, Carl (Committee member) / Belitsky, Andrei (Committee member) / Arizona State University (Publisher)
Created2012
Description
Diseases have been part of human life for generations and evolve within the population, sometimes dying out while other times becoming endemic or the cause of recurrent outbreaks. The long term influence of a disease stems from different dynamics within or between pathogen-host, that have been analyzed and studied by many researchers using mathematical models. Co-infection with different pathogens is common, yet little is known about how infection with one pathogen affects the host's immunological response to another. Moreover, no work has been found in the literature that considers the variability of the host immune health or that examines a disease at the population level and its corresponding interconnectedness with the host immune system. Knowing that the spread of the disease in the population starts at the individual level, this thesis explores how variability in immune system response within an endemic environment affects an individual's vulnerability, and how prone it is to co-infections. Immunology-based models of Malaria and Tuberculosis (TB) are constructed by extending and modifying existing mathematical models in the literature. The two are then combined to give a single nine-variable model of co-infection with Malaria and TB. Because these models are difficult to gain any insight analytically due to the large number of parameters, a phenomenological model of co-infection is proposed with subsystems corresponding to the individual immunology-based model of a single infection. Within this phenomenological model, the variability of the host immune health is also incorporated through three different pathogen response curves using nonlinear bounded Michaelis-Menten functions that describe the level or state of immune system (healthy, moderate and severely compromised). The immunology-based models of Malaria and TB give numerical results that agree with the biological observations. The Malaria--TB co-infection model gives reasonable results and these suggest that the order in which the two diseases are introduced have an impact on the behavior of both. The subsystems of the phenomenological models that correspond to a single infection (either of Malaria or TB) mimic much of the observed behavior of the immunology-based counterpart and can demonstrate different behavior depending on the chosen pathogen response curve. In addition, varying some of the parameters and initial conditions in the phenomenological model yields a range of topologically different mathematical behaviors, which suggests that this behavior may be able to be observed in the immunology-based models as well. The phenomenological models clearly replicate the qualitative behavior of primary and secondary infection as well as co-infection. The mathematical solutions of the models correspond to the fundamental states described by immunologists: virgin state, immune state and tolerance state. The phenomenological model of co-infection also demonstrates a range of parameter values and initial conditions in which the introduction of a second disease causes both diseases to grow without bound even though those same parameters and initial conditions did not yield unbounded growth in the corresponding subsystems. This results applies to all three states of the host immune system. In terms of the immunology-based system, this would suggest the following: there may be parameter values and initial conditions in which a person can clear Malaria or TB (separately) from their system but in which the presence of both can result in the person dying of one of the diseases. Finally, this thesis studies links between epidemiology (population level) and immunology in an effort to assess the impact of pathogen's spread within the population on the immune response of individuals. Models of Malaria and TB are proposed that incorporate the immune system of the host into a mathematical model of an epidemic at the population level.
ContributorsSoho, Edmé L (Author) / Wirkus, Stephen (Thesis advisor) / Castillo-Chavez, Carlos (Thesis advisor) / Chowell-Puente, Gerardo (Committee member) / Arizona State University (Publisher)
Created2011
Description
Factory production is stochastic in nature with time varying input and output processes that are non-stationary stochastic processes. Hence, the principle quantities of interest are random variables. Typical modeling of such behavior involves numerical simulation and statistical analysis. A deterministic closure model leading to a second order model for the product density and product speed has previously been proposed. The resulting partial differential equations (PDE) are compared to discrete event simulations (DES) that simulate factory production as a time dependent M/M/1 queuing system. Three fundamental scenarios for the time dependent influx are studied: An instant step up/down of the mean arrival rate; an exponential step up/down of the mean arrival rate; and periodic variation of the mean arrival rate. It is shown that the second order model, in general, yields significant improvement over current first order models. Specifically, the agreement between the DES and the PDE for the step up and for periodic forcing that is not too rapid is very good. Adding diffusion to the PDE further improves the agreement. The analysis also points to fundamental open issues regarding the deterministic modeling of low signal-to-noise ratio for some stochastic processes and the possibility of resonance in deterministic models that is not present in the original stochastic process.
ContributorsWienke, Matthew (Author) / Armbruster, Dieter (Thesis advisor) / Jones, Donald (Committee member) / Platte, Rodrigo (Committee member) / Gardner, Carl (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2015
Description
Swarms of animals, fish, birds, locusts etc. are a common occurrence but their coherence and method of organization poses a major question for mathematics and biology.The Vicsek and the Attraction-Repulsion are two models that have been proposed to explain the emergence of collective motion. A major issue for the Vicsek Model is that its particles are not attracted to each other, leaving the swarm with alignment in velocity but without spatial coherence. Restricting the particles to a bounded domain generates global spatial coherence of swarms while maintaining velocity alignment. While individual particles are specularly reflected at the boundary, the swarm as a whole is not. As a result, new dynamical swarming solutions are found.
The Attraction-Repulsion Model set with a long-range attraction and short-range repulsion interaction potential typically stabilizes to a well-studied flock steady state solution. The particles for a flock remain spatially coherent but have no spatial bound and explore all space. A bounded domain with specularly reflecting walls traps the particles within a specific region. A fundamental refraction law for a swarm impacting on a planar boundary is derived. The swarm reflection varies from specular for a swarm dominated by
kinetic energy to inelastic for a swarm dominated by potential energy. Inelastic collisions lead to alignment with the wall and to damped pulsating oscillations of the swarm. The fundamental refraction law provides a one-dimensional iterative map that allows for a prediction and analysis of the trajectory of the center of mass of a flock in a channel and a square domain.
The extension of the wall collisions to a scattering experiment is conducted by setting two identical flocks to collide. The two particle dynamics is studied analytically and shows a transition from scattering: diverging flocks to bound states in the form of oscillations or parallel motions. Numerical studies of collisions of flocks show the same transition where the bound states become either a single translating flock or a rotating (mill).
The Attraction-Repulsion Model set with a long-range attraction and short-range repulsion interaction potential typically stabilizes to a well-studied flock steady state solution. The particles for a flock remain spatially coherent but have no spatial bound and explore all space. A bounded domain with specularly reflecting walls traps the particles within a specific region. A fundamental refraction law for a swarm impacting on a planar boundary is derived. The swarm reflection varies from specular for a swarm dominated by
kinetic energy to inelastic for a swarm dominated by potential energy. Inelastic collisions lead to alignment with the wall and to damped pulsating oscillations of the swarm. The fundamental refraction law provides a one-dimensional iterative map that allows for a prediction and analysis of the trajectory of the center of mass of a flock in a channel and a square domain.
The extension of the wall collisions to a scattering experiment is conducted by setting two identical flocks to collide. The two particle dynamics is studied analytically and shows a transition from scattering: diverging flocks to bound states in the form of oscillations or parallel motions. Numerical studies of collisions of flocks show the same transition where the bound states become either a single translating flock or a rotating (mill).
ContributorsThatcher, Andrea (Author) / Armbruster, Hans (Thesis advisor) / Motsch, Sebastien (Committee member) / Ringhofer, Christian (Committee member) / Platte, Rodrigo (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2015
Description
Inverse problems model real world phenomena from data, where the data are often noisy and models contain errors. This leads to instabilities, multiple solution vectors and thus ill-posedness. To solve ill-posed inverse problems, regularization is typically used as a penalty function to induce stability and allow for the incorporation of a priori information about the desired solution. In this thesis, high order regularization techniques are developed for image and function reconstruction from noisy or misleading data. Specifically the incorporation of the Polynomial Annihilation operator allows for the accurate exploitation of the sparse representation of each function in the edge domain.
This dissertation tackles three main problems through the development of novel reconstruction techniques: (i) reconstructing one and two dimensional functions from multiple measurement vectors using variance based joint sparsity when a subset of the measurements contain false and/or misleading information, (ii) approximating discontinuous solutions to hyperbolic partial differential equations by enhancing typical solvers with l1 regularization, and (iii) reducing model assumptions in synthetic aperture radar image formation, specifically for the purpose of speckle reduction and phase error correction. While the common thread tying these problems together is the use of high order regularization, the defining characteristics of each of these problems create unique challenges.
Fast and robust numerical algorithms are also developed so that these problems can be solved efficiently without requiring fine tuning of parameters. Indeed, the numerical experiments presented in this dissertation strongly suggest that the new methodology provides more accurate and robust solutions to a variety of ill-posed inverse problems.
This dissertation tackles three main problems through the development of novel reconstruction techniques: (i) reconstructing one and two dimensional functions from multiple measurement vectors using variance based joint sparsity when a subset of the measurements contain false and/or misleading information, (ii) approximating discontinuous solutions to hyperbolic partial differential equations by enhancing typical solvers with l1 regularization, and (iii) reducing model assumptions in synthetic aperture radar image formation, specifically for the purpose of speckle reduction and phase error correction. While the common thread tying these problems together is the use of high order regularization, the defining characteristics of each of these problems create unique challenges.
Fast and robust numerical algorithms are also developed so that these problems can be solved efficiently without requiring fine tuning of parameters. Indeed, the numerical experiments presented in this dissertation strongly suggest that the new methodology provides more accurate and robust solutions to a variety of ill-posed inverse problems.
ContributorsScarnati, Theresa (Author) / Gelb, Anne (Thesis advisor) / Platte, Rodrigo (Thesis advisor) / Cochran, Douglas (Committee member) / Gardner, Carl (Committee member) / Sanders, Toby (Committee member) / Arizona State University (Publisher)
Created2018
Description
Predicting resistant prostate cancer is critical for lowering medical costs and improving the quality of life of advanced prostate cancer patients. I formulate, compare, and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). I accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). I demonstrate that the inverse problem of parameter estimation might be too complicated and simply relying on data fitting can give incorrect conclusions, since there is a large error in parameter values estimated and parameters might be unidentifiable. I provide confidence intervals to give estimate forecasts using data assimilation via an ensemble Kalman Filter. Using the ensemble Kalman Filter, I perform dual estimation of parameters and state variables to test the prediction accuracy of the models. Finally, I present a novel model with time delay and a delay-dependent parameter. I provide a geometric stability result to study the behavior of this model and show that the inclusion of time delay may improve the accuracy of predictions. Also, I demonstrate with clinical data that the inclusion of the delay-dependent parameter facilitates the identification and estimation of parameters.
ContributorsBaez, Javier (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric (Committee member) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Nagy, John (Committee member) / Arizona State University (Publisher)
Created2017
Description
Need-based transfers (NBTs) are a form of risk-pooling in which binary welfare exchanges
occur to preserve the viable participation of individuals in an economy, e.g. reciprocal gifting
of cattle among East African herders or food sharing among vampire bats. With the
broad goal of better understanding the mathematics of such binary welfare and risk pooling,
agent-based simulations are conducted to explore socially optimal transfer policies
and sharing network structures, kinetic exchange models that utilize tools from the kinetic
theory of gas dynamics are utilized to characterize the wealth distribution of an NBT economy,
and a variant of repeated prisoner’s dilemma is analyzed to determine whether and
why individuals would participate in such a system of reciprocal altruism.
From agent-based simulation and kinetic exchange models, it is found that regressive
NBT wealth redistribution acts as a cutting stock optimization heuristic that most efficiently
matches deficits to surpluses to improve short-term survival; however, progressive
redistribution leads to a wealth distribution that is more stable in volatile environments and
therefore is optimal for long-term survival. Homogeneous sharing networks with low variance
in degree are found to be ideal for maintaining community viability as the burden and
benefit of NBTs is equally shared. Also, phrasing NBTs as a survivor’s dilemma reveals
parameter regions where the repeated game becomes equivalent to a stag hunt or harmony
game, and thus where cooperation is evolutionarily stable.
occur to preserve the viable participation of individuals in an economy, e.g. reciprocal gifting
of cattle among East African herders or food sharing among vampire bats. With the
broad goal of better understanding the mathematics of such binary welfare and risk pooling,
agent-based simulations are conducted to explore socially optimal transfer policies
and sharing network structures, kinetic exchange models that utilize tools from the kinetic
theory of gas dynamics are utilized to characterize the wealth distribution of an NBT economy,
and a variant of repeated prisoner’s dilemma is analyzed to determine whether and
why individuals would participate in such a system of reciprocal altruism.
From agent-based simulation and kinetic exchange models, it is found that regressive
NBT wealth redistribution acts as a cutting stock optimization heuristic that most efficiently
matches deficits to surpluses to improve short-term survival; however, progressive
redistribution leads to a wealth distribution that is more stable in volatile environments and
therefore is optimal for long-term survival. Homogeneous sharing networks with low variance
in degree are found to be ideal for maintaining community viability as the burden and
benefit of NBTs is equally shared. Also, phrasing NBTs as a survivor’s dilemma reveals
parameter regions where the repeated game becomes equivalent to a stag hunt or harmony
game, and thus where cooperation is evolutionarily stable.
ContributorsKayser, Kirk (Author) / Armbruster, Dieter (Thesis advisor) / Lampert, Adam (Committee member) / Ringhofer, Christian (Committee member) / Motsch, Sebastien (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2018
Description
The most advanced social insects, the eusocial insects, form often large societies in which there is reproductive division of labor, queens and workers, have overlapping generations, and cooperative brood care where daughter workers remain in the nest with their queen mother and care for their siblings. The eusocial insects are composed of representative species of bees and wasps, and all species of ants and termites. Much is known about their organizational structure, but remains to be discovered.
The success of social insects is dependent upon cooperative behavior and adaptive strategies shaped by natural selection that respond to internal or external conditions. The objective of my research was to investigate specific mechanisms that have helped shaped the structure of division of labor observed in social insect colonies, including age polyethism and nutrition, and phenomena known to increase colony survival such as egg cannibalism. I developed various Ordinary Differential Equation (ODE) models in which I applied dynamical, bifurcation, and sensitivity analysis to carefully study and visualize biological outcomes in social organisms to answer questions regarding the conditions under which a colony can survive. First, I investigated how the population and evolutionary dynamics of egg cannibalism and division of labor can promote colony survival. I then introduced a model of social conflict behavior to study the inclusion of different response functions that explore the benefits of cannibalistic behavior and how it contributes to age polyethism, the change in behavior of workers as they age, and its biological relevance. Finally, I introduced a model to investigate the importance of pollen nutritional status in a honeybee colony, how it affects population growth and influences division of labor within the worker caste. My results first reveal that both cannibalism and division of labor are adaptive strategies that increase the size of the worker population, and therefore, the persistence of the colony. I show the importance of food collection, consumption, and processing rates to promote good colony nutrition leading to the coexistence of brood and adult workers. Lastly, I show how taking into account seasonality for pollen collection improves the prediction of long term consequences.
The success of social insects is dependent upon cooperative behavior and adaptive strategies shaped by natural selection that respond to internal or external conditions. The objective of my research was to investigate specific mechanisms that have helped shaped the structure of division of labor observed in social insect colonies, including age polyethism and nutrition, and phenomena known to increase colony survival such as egg cannibalism. I developed various Ordinary Differential Equation (ODE) models in which I applied dynamical, bifurcation, and sensitivity analysis to carefully study and visualize biological outcomes in social organisms to answer questions regarding the conditions under which a colony can survive. First, I investigated how the population and evolutionary dynamics of egg cannibalism and division of labor can promote colony survival. I then introduced a model of social conflict behavior to study the inclusion of different response functions that explore the benefits of cannibalistic behavior and how it contributes to age polyethism, the change in behavior of workers as they age, and its biological relevance. Finally, I introduced a model to investigate the importance of pollen nutritional status in a honeybee colony, how it affects population growth and influences division of labor within the worker caste. My results first reveal that both cannibalism and division of labor are adaptive strategies that increase the size of the worker population, and therefore, the persistence of the colony. I show the importance of food collection, consumption, and processing rates to promote good colony nutrition leading to the coexistence of brood and adult workers. Lastly, I show how taking into account seasonality for pollen collection improves the prediction of long term consequences.
ContributorsRodríguez Messan, Marisabel (Author) / Kang, Yun (Thesis advisor) / Castillo-Chavez, Carlos (Thesis advisor) / Kuang, Yang (Committee member) / Page Jr., Robert E (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2018
Description
Mathematical models are important tools for addressing problems that exceed experimental capabilities. In this work, I present ordinary and partial differential equation (ODE, PDE) models for two problems: Vicodin abuse and impact cratering.
The prescription opioid Vicodin is the nation's most widely prescribed pain reliever. The majority of Vicodin abusers are first introduced via prescription, distinguishing it from other drugs in which the most common path to abuse begins with experimentation. I develop and analyze two mathematical models of Vicodin use and abuse, considering only those patients with an initial Vicodin prescription. Through adjoint sensitivity analysis, I show that focusing efforts on prevention rather than treatment has greater success at reducing the total population of abusers. I prove that solutions to each model exist, are unique, and are non-negative. I also derive conditions for which these solutions are asymptotically stable.
Verification and Validation (V&V) are necessary processes to ensure accuracy of computational methods. Simulations are essential for addressing impact cratering problems, because these problems often exceed experimental capabilities. I show that the Free Lagrange (FLAG) hydrocode, developed and maintained by Los Alamos National Laboratory, can be used for impact cratering simulations by verifying FLAG against two analytical models of aluminum-on-aluminum impacts at different impact velocities and validating FLAG against a glass-into-water laboratory impact experiment. My verification results show good agreement with the theoretical maximum pressures, and my mesh resolution study shows that FLAG converges at resolutions low enough to reduce the required computation time from about 28 hours to about 25 minutes.
Asteroid 16 Psyche is the largest M-type (metallic) asteroid in the Main Asteroid Belt. Radar albedo data indicate Psyche's surface is rich in metallic content, but estimates for Psyche's composition vary widely. Psyche has two large impact structures in its Southern hemisphere, with estimated diameters from 50 km to 70 km and estimated depths up to 6.4 km. I use the FLAG hydrocode to model the formation of the largest of these impact structures. My results indicate an oblique angle of impact rather than a vertical impact. These results also support previous claims that Psyche is metallic and porous.
The prescription opioid Vicodin is the nation's most widely prescribed pain reliever. The majority of Vicodin abusers are first introduced via prescription, distinguishing it from other drugs in which the most common path to abuse begins with experimentation. I develop and analyze two mathematical models of Vicodin use and abuse, considering only those patients with an initial Vicodin prescription. Through adjoint sensitivity analysis, I show that focusing efforts on prevention rather than treatment has greater success at reducing the total population of abusers. I prove that solutions to each model exist, are unique, and are non-negative. I also derive conditions for which these solutions are asymptotically stable.
Verification and Validation (V&V) are necessary processes to ensure accuracy of computational methods. Simulations are essential for addressing impact cratering problems, because these problems often exceed experimental capabilities. I show that the Free Lagrange (FLAG) hydrocode, developed and maintained by Los Alamos National Laboratory, can be used for impact cratering simulations by verifying FLAG against two analytical models of aluminum-on-aluminum impacts at different impact velocities and validating FLAG against a glass-into-water laboratory impact experiment. My verification results show good agreement with the theoretical maximum pressures, and my mesh resolution study shows that FLAG converges at resolutions low enough to reduce the required computation time from about 28 hours to about 25 minutes.
Asteroid 16 Psyche is the largest M-type (metallic) asteroid in the Main Asteroid Belt. Radar albedo data indicate Psyche's surface is rich in metallic content, but estimates for Psyche's composition vary widely. Psyche has two large impact structures in its Southern hemisphere, with estimated diameters from 50 km to 70 km and estimated depths up to 6.4 km. I use the FLAG hydrocode to model the formation of the largest of these impact structures. My results indicate an oblique angle of impact rather than a vertical impact. These results also support previous claims that Psyche is metallic and porous.
ContributorsCaldwell, Wendy K (Author) / Wirkus, Stephen (Thesis advisor) / Asphaug, Erik (Committee member) / Camacho, Erika T (Committee member) / Crook, Sharon (Committee member) / Plesko, Catherine S (Committee member) / Smith, Hal (Committee member) / Arizona State University (Publisher)
Created2019