Matching Items (11)
Filtering by
- All Subjects: Applied Mathematics
- Creators: Kostelich, Eric
- Member of: ASU Electronic Theses and Dissertations
Description
In 1968, phycologist M.R. Droop published his famous discovery on the functional relationship between growth rate and internal nutrient status of algae in chemostat culture. The simple notion that growth is directly dependent on intracellular nutrient concentration is useful for understanding the dynamics in many ecological systems. The cell quota in particular lends itself to ecological stoichiometry, which is a powerful framework for mathematical ecology. Three models are developed based on the cell quota principal in order to demonstrate its applications beyond chemostat culture.
First, a data-driven model is derived for neutral lipid synthesis in green microalgae with respect to nitrogen limitation. This model synthesizes several established frameworks in phycology and ecological stoichiometry. The model demonstrates how the cell quota is a useful abstraction for understanding the metabolic shift to neutral lipid production that is observed in certain oleaginous species.
Next a producer-grazer model is developed based on the cell quota model and nutrient recycling. The model incorporates a novel feedback loop to account for animal toxicity due to accumulation of nitrogen waste. The model exhibits rich, complex dynamics which leave several open mathematical questions.
Lastly, disease dynamics in vivo are in many ways analogous to those of an ecosystem, giving natural extensions of the cell quota concept to disease modeling. Prostate cancer can be modeled within this framework, with androgen the limiting nutrient and the prostate and cancer cells as competing species. Here the cell quota model provides a useful abstraction for the dependence of cellular proliferation and apoptosis on androgen and the androgen receptor. Androgen ablation therapy is often used for patients in biochemical recurrence or late-stage disease progression and is in general initially effective. However, for many patients the cancer eventually develops resistance months to years after treatment begins. Understanding how and predicting when hormone therapy facilitates evolution of resistant phenotypes has immediate implications for treatment. Cell quota models for prostate cancer can be useful tools for this purpose and motivate applications to other diseases.
First, a data-driven model is derived for neutral lipid synthesis in green microalgae with respect to nitrogen limitation. This model synthesizes several established frameworks in phycology and ecological stoichiometry. The model demonstrates how the cell quota is a useful abstraction for understanding the metabolic shift to neutral lipid production that is observed in certain oleaginous species.
Next a producer-grazer model is developed based on the cell quota model and nutrient recycling. The model incorporates a novel feedback loop to account for animal toxicity due to accumulation of nitrogen waste. The model exhibits rich, complex dynamics which leave several open mathematical questions.
Lastly, disease dynamics in vivo are in many ways analogous to those of an ecosystem, giving natural extensions of the cell quota concept to disease modeling. Prostate cancer can be modeled within this framework, with androgen the limiting nutrient and the prostate and cancer cells as competing species. Here the cell quota model provides a useful abstraction for the dependence of cellular proliferation and apoptosis on androgen and the androgen receptor. Androgen ablation therapy is often used for patients in biochemical recurrence or late-stage disease progression and is in general initially effective. However, for many patients the cancer eventually develops resistance months to years after treatment begins. Understanding how and predicting when hormone therapy facilitates evolution of resistant phenotypes has immediate implications for treatment. Cell quota models for prostate cancer can be useful tools for this purpose and motivate applications to other diseases.
ContributorsPacker, Aaron (Author) / Kuang, Yang (Thesis advisor) / Nagy, John (Committee member) / Smith, Hal (Committee member) / Kostelich, Eric (Committee member) / Kang, Yun (Committee member) / Arizona State University (Publisher)
Created2014
Description
The tools developed for the use of investigating dynamical systems have provided critical understanding to a wide range of physical phenomena. Here these tools are used to gain further insight into scalar transport, and how it is affected by mixing. The aim of this research is to investigate the efficiency of several different partitioning methods which demarcate flow fields into dynamically distinct regions, and the correlation of finite-time statistics from the advection-diffusion equation to these regions.
For autonomous systems, invariant manifold theory can be used to separate the system into dynamically distinct regions. Despite there being no equivalent method for nonautonomous systems, a similar analysis can be done. Systems with general time dependencies must resort to using finite-time transport barriers for partitioning; these barriers are the edges of Lagrangian coherent structures (LCS), the analog to the stable and unstable manifolds of invariant manifold theory. Using the coherent structures of a flow to analyze the statistics of trapping, flight, and residence times, the signature of anomalous diffusion are obtained.
This research also investigates the use of linear models for approximating the elements of the covariance matrix of nonlinear flows, and then applying the covariance matrix approximation over coherent regions. The first and second-order moments can be used to fully describe an ensemble evolution in linear systems, however there is no direct method for nonlinear systems. The problem is only compounded by the fact that the moments for nonlinear flows typically don't have analytic representations, therefore direct numerical simulations would be needed to obtain the moments throughout the domain. To circumvent these many computations, the nonlinear system is approximated as many linear systems for which analytic expressions for the moments exist. The parameters introduced in the linear models are obtained locally from the nonlinear deformation tensor.
For autonomous systems, invariant manifold theory can be used to separate the system into dynamically distinct regions. Despite there being no equivalent method for nonautonomous systems, a similar analysis can be done. Systems with general time dependencies must resort to using finite-time transport barriers for partitioning; these barriers are the edges of Lagrangian coherent structures (LCS), the analog to the stable and unstable manifolds of invariant manifold theory. Using the coherent structures of a flow to analyze the statistics of trapping, flight, and residence times, the signature of anomalous diffusion are obtained.
This research also investigates the use of linear models for approximating the elements of the covariance matrix of nonlinear flows, and then applying the covariance matrix approximation over coherent regions. The first and second-order moments can be used to fully describe an ensemble evolution in linear systems, however there is no direct method for nonlinear systems. The problem is only compounded by the fact that the moments for nonlinear flows typically don't have analytic representations, therefore direct numerical simulations would be needed to obtain the moments throughout the domain. To circumvent these many computations, the nonlinear system is approximated as many linear systems for which analytic expressions for the moments exist. The parameters introduced in the linear models are obtained locally from the nonlinear deformation tensor.
ContributorsWalker, Phillip (Author) / Tang, Wenbo (Thesis advisor) / Kostelich, Eric (Committee member) / Mahalov, Alex (Committee member) / Moustaoui, Mohamed (Committee member) / Platte, Rodrigo (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
Two urban flows are analyzed, one concerned with pollutant transport in a Phoenix, Arizona neighborhood and the other with windshear detection at the Hong Kong International Airport (HKIA).
Lagrangian measures, identified with finite-time Lyapunov exponents, are first used to characterize transport patterns of inertial pollutant particles. Motivated by actual events the focus is on flows in realistic urban geometry. Both deterministic and stochastic transport patterns are identified, as inertial Lagrangian coherent structures. For the deterministic case, the organizing structures are well defined and are extracted at different hours of a day to reveal the variability of coherent patterns. For the stochastic case, a random displacement model for fluid particles is formulated, and used to derive the governing equations for inertial particles to examine the change in organizing structures due to ``zeroth-order'' random noise. It is found that, (1) the Langevin equation for inertial particles can be reduced to a random displacement model; (2) using random noise based on inhomogeneous turbulence, whose diffusivity is derived from $k$-$\epsilon$ models, major coherent structures survive to organize local flow patterns and weaker structures are smoothed out due to random motion.
A study of three-dimensional Lagrangian coherent structures (LCS) near HKIA is then presented and related to previous developments of two-dimensional (2D) LCS analyses in detecting windshear experienced by landing aircraft. The LCS are contrasted among three independent models and against 2D coherent Doppler light detection and ranging (LIDAR) data. Addition of the velocity information perpendicular to the lidar scanning cone helps solidify flow structures inferred from previous studies; contrast among models reveals the intramodel variability; and comparison with flight data evaluates the performance among models in terms of Lagrangian analyses. It is found that, while the three models and the LIDAR do recover similar features of the windshear experienced by a landing aircraft (along the landing trajectory), their Lagrangian signatures over the entire domain are quite different - a portion of each numerical model captures certain features resembling those LCS extracted from independent 2D LIDAR analyses based on observations. Overall, it was found that the Weather Research and Forecast (WRF) model provides the best agreement with the LIDAR data.
Finally, the three-dimensional variational (3DVAR) data assimilation scheme in WRF is used to incorporate the LIDAR line of sight velocity observations into the WRF model forecast at HKIA. Using two different days as test cases, it is found that the LIDAR data can be successfully and consistently assimilated into WRF. Using the updated model forecast LCS are extracted along the LIDAR scanning cone and compare to onboard flight data. It is found that the LCS generated from the updated WRF forecasts are generally better correlated with the windshear experienced by landing aircraft as compared to the LIDAR extracted LCS alone, which suggests that such a data assimilation scheme could be used for the prediction of windshear events.
Lagrangian measures, identified with finite-time Lyapunov exponents, are first used to characterize transport patterns of inertial pollutant particles. Motivated by actual events the focus is on flows in realistic urban geometry. Both deterministic and stochastic transport patterns are identified, as inertial Lagrangian coherent structures. For the deterministic case, the organizing structures are well defined and are extracted at different hours of a day to reveal the variability of coherent patterns. For the stochastic case, a random displacement model for fluid particles is formulated, and used to derive the governing equations for inertial particles to examine the change in organizing structures due to ``zeroth-order'' random noise. It is found that, (1) the Langevin equation for inertial particles can be reduced to a random displacement model; (2) using random noise based on inhomogeneous turbulence, whose diffusivity is derived from $k$-$\epsilon$ models, major coherent structures survive to organize local flow patterns and weaker structures are smoothed out due to random motion.
A study of three-dimensional Lagrangian coherent structures (LCS) near HKIA is then presented and related to previous developments of two-dimensional (2D) LCS analyses in detecting windshear experienced by landing aircraft. The LCS are contrasted among three independent models and against 2D coherent Doppler light detection and ranging (LIDAR) data. Addition of the velocity information perpendicular to the lidar scanning cone helps solidify flow structures inferred from previous studies; contrast among models reveals the intramodel variability; and comparison with flight data evaluates the performance among models in terms of Lagrangian analyses. It is found that, while the three models and the LIDAR do recover similar features of the windshear experienced by a landing aircraft (along the landing trajectory), their Lagrangian signatures over the entire domain are quite different - a portion of each numerical model captures certain features resembling those LCS extracted from independent 2D LIDAR analyses based on observations. Overall, it was found that the Weather Research and Forecast (WRF) model provides the best agreement with the LIDAR data.
Finally, the three-dimensional variational (3DVAR) data assimilation scheme in WRF is used to incorporate the LIDAR line of sight velocity observations into the WRF model forecast at HKIA. Using two different days as test cases, it is found that the LIDAR data can be successfully and consistently assimilated into WRF. Using the updated model forecast LCS are extracted along the LIDAR scanning cone and compare to onboard flight data. It is found that the LCS generated from the updated WRF forecasts are generally better correlated with the windshear experienced by landing aircraft as compared to the LIDAR extracted LCS alone, which suggests that such a data assimilation scheme could be used for the prediction of windshear events.
ContributorsKnutson, Brent (Author) / Tang, Wenbo (Thesis advisor) / Calhoun, Ronald (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Mahalov, Alex (Committee member) / Arizona State University (Publisher)
Created2018
Description
I investigate two models interacting agent systems: the first is motivated by the flocking and swarming behaviors in biological systems, while the second models opinion formation in social networks. In each setting, I define natural notions of convergence (to a ``flock" and to a ``consensus'', respectively), and study the convergence properties of each in the limit as $t \rightarrow \infty$. Specifically, I provide sufficient conditions for the convergence of both of the models, and conduct numerical experiments to study the resulting solutions.
ContributorsTheisen, Ryan (Author) / Motsch, Sebastien (Thesis advisor) / Lanchier, Nicholas (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2018
Description
The increased number of novel pathogens that potentially threaten the human population has motivated the development of mathematical and computational modeling approaches for forecasting epidemic impact and understanding key environmental characteristics that influence the spread of diseases. Yet, in the case that substantial uncertainty surrounds the transmission process during a rapidly developing infectious disease outbreak, complex mechanistic models may be too difficult to be calibrated quick enough for policy makers to make informed decisions. Simple phenomenological models that rely on a small number of parameters can provide an initial platform for assessing the epidemic trajectory, estimating the reproduction number and quantifying the disease burden from the early epidemic phase.
Chapter 1 provides background information and motivation for infectious disease forecasting and outlines the rest of the thesis.
In chapter 2, logistic patch models are used to assess and forecast the 2013-2015 West Africa Zaire ebolavirus epidemic. In particular, this chapter is concerned with comparing and contrasting the effects that spatial heterogeneity has on the forecasting performance of the cumulative infected case counts reported during the epidemic.
In chapter 3, two simple phenomenological models inspired from population biology are used to assess the Research and Policy for Infectious Disease Dynamics (RAPIDD) Ebola Challenge; a simulated epidemic that generated 4 infectious disease scenarios. Because of the nature of the synthetically generated data, model predictions are compared to exact epidemiological quantities used in the simulation.
In chapter 4, these models are applied to the 1904 Plague epidemic that occurred in Bombay. This chapter provides evidence that these simple models may be applicable to infectious diseases no matter the disease transmission mechanism.
Chapter 5, uses the patch models from chapter 2 to explore how migration in the 1904 Plague epidemic changes the final epidemic size.
The final chapter is an interdisciplinary project concerning within-host dynamics of cereal yellow dwarf virus-RPV, a plant pathogen from a virus group that infects over 150 grass species. Motivated by environmental nutrient enrichment due to anthropological activities, mathematical models are employed to investigate the relevance of resource competition to pathogen and host dynamics.
Chapter 1 provides background information and motivation for infectious disease forecasting and outlines the rest of the thesis.
In chapter 2, logistic patch models are used to assess and forecast the 2013-2015 West Africa Zaire ebolavirus epidemic. In particular, this chapter is concerned with comparing and contrasting the effects that spatial heterogeneity has on the forecasting performance of the cumulative infected case counts reported during the epidemic.
In chapter 3, two simple phenomenological models inspired from population biology are used to assess the Research and Policy for Infectious Disease Dynamics (RAPIDD) Ebola Challenge; a simulated epidemic that generated 4 infectious disease scenarios. Because of the nature of the synthetically generated data, model predictions are compared to exact epidemiological quantities used in the simulation.
In chapter 4, these models are applied to the 1904 Plague epidemic that occurred in Bombay. This chapter provides evidence that these simple models may be applicable to infectious diseases no matter the disease transmission mechanism.
Chapter 5, uses the patch models from chapter 2 to explore how migration in the 1904 Plague epidemic changes the final epidemic size.
The final chapter is an interdisciplinary project concerning within-host dynamics of cereal yellow dwarf virus-RPV, a plant pathogen from a virus group that infects over 150 grass species. Motivated by environmental nutrient enrichment due to anthropological activities, mathematical models are employed to investigate the relevance of resource competition to pathogen and host dynamics.
ContributorsPell, Bruce (Author) / Kuang, Yang (Thesis advisor) / Chowell-Puente, Gerardo (Committee member) / Nagy, John (Committee member) / Kostelich, Eric (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2016
Description
A general continuum model for simulating the flow of ions in the salt baths that surround and fill excitable neurons is developed and presented. The ion densities and electric potential are computed using the drift-diffusion equations. In addition, a detailed model is given for handling the electrical dynamics on interior membrane boundaries, including a model for ion channels in the membranes that facilitate the transfer of ions in and out of cells. The model is applied to the triad synapse found in the outer plexiform layer of the retina in most species. Experimental evidence suggests the existence of a negative feedback pathway between horizontal cells and cone photoreceptors that modulates the flow of calcium ions into the synaptic terminals of cones. However, the underlying mechanism for this feedback is controversial and there are currently three competing hypotheses: the ephaptic hypothesis, the pH hypothesis and the GABA hypothesis. The goal of this work is to test some features of the ephaptic hypothesis using detailed simulations that employ rigorous numerical methods. The model is first applied in a simple rectangular geometry to demonstrate the effects of feedback for different extracellular gap widths. The model is then applied to a more complex and realistic geometry to demonstrate the existence of strictly electrical feedback, as predicted by the ephaptic hypothesis. Lastly, the effects of electrical feedback in regards to the behavior of the bipolar cell membrane potential is explored. Figures for the ion densities and electric potential are presented to verify key features of the model. The computed steady state IV curves for several cases are presented, which can be compared to experimental data. The results provide convincing evidence in favor of the ephaptic hypothesis since the existence of feedback that is strictly electrical in nature is shown, without any dependence on pH effects or chemical transmitters.
ContributorsJones, Jeremiah (Author) / Gardner, Carl (Committee member) / Baer, Steven (Committee member) / Crook, Sharon (Committee member) / Kostelich, Eric (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2013
Description
High-dimensional systems are difficult to model and predict. The underlying mechanisms of such systems are too complex to be fully understood with limited theoretical knowledge and/or physical measurements. Nevertheless, redcued-order models have been widely used to study high-dimensional systems, because they are practical and efficient to develop and implement. Although model errors (biases) are inevitable for reduced-order models, these models can still be proven useful to develop real-world applications. Evaluation and validation for idealized models are indispensable to serve the mission of developing useful applications. Data assimilation and uncertainty quantification can provide a way to assess the performance of a reduced-order model. Real data and a dynamical model are combined together in a data assimilation framework to generate corrected model forecasts of a system. Uncertainties in model forecasts and observations are also quantified in a data assimilation cycle to provide optimal updates that are representative of the real dynamics. In this research, data assimilation is applied to assess the performance of two reduced-order models. The first model is developed for predicting prostate cancer treatment response under intermittent androgen suppression therapy. A sequential data assimilation scheme, the ensemble Kalman filter (EnKF), is used to quantify uncertainties in model predictions using clinical data of individual patients provided by Vancouver Prostate Center. The second model is developed to study what causes the changes of the state of stratospheric polar vortex. Two data assimilation schemes: EnKF and ES-MDA (ensemble smoother with multiple data assimilation), are used to validate the qualitative properties of the model using ECMWF (European Center for Medium-Range Weather Forecasts) reanalysis data. In both studies, the reduced-order model is able to reproduce the data patterns and provide insights to understand the underlying mechanism. However, significant model errors are also diagnosed for both models from the results of data assimilation schemes, which suggests specific improvements of the reduced-order models.
ContributorsWu, Zhimin (Author) / Kostelich, Eric (Thesis advisor) / Moustaoui, Mohamed (Thesis advisor) / Jones, Chris (Committee member) / Espanol, Malena (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2021
Description
A pneumonia-like illness emerged late in 2019 (coined COVID-19), caused by SARSCoV-2, causing a devastating global pandemic on a scale never before seen sincethe 1918/1919 influenza pandemic. This dissertation contributes in providing deeper
qualitative insights into the transmission dynamics and control of the disease in the
United States. A basic mathematical model, which incorporates the key pertinent
epidemiological features of SARS-CoV-2 and fitted using observed COVID-19 data,
was designed and used to assess the population-level impacts of vaccination and face
mask usage in mitigating the burden of the pandemic in the United States. Conditions
for the existence and asymptotic stability of the various equilibria of the model were
derived. The model was shown to undergo a vaccine-induced backward bifurcation
when the associated reproduction number is less than one. Conditions for achieving
vaccine-derived herd immunity were derived for three of the four FDA-approved vaccines (namely Pfizer, Moderna and Johnson & Johnson vaccine), and the vaccination
coverage level needed to achieve it decreases with increasing coverage of moderately and highly-effective face masks. It was also shown that using face masks as a singular
intervention strategy could lead to the elimination of the pandemic if moderate or
highly-effective masks are prioritized and pandemic elimination prospects are greatly
enhanced if the vaccination program is combined with a face mask use strategy that
emphasizes the use of moderate to highly-effective masks with at least moderate coverage. The model was extended in Chapter 3 to allow for the assessment of the
impacts of waning and boosting of vaccine-derived and natural immunity against
the BA.1 Omicron variant of SARS-CoV-2. It was shown that vaccine-derived herd
immunity can be achieved in the United States via a vaccination-boosting strategy
which entails fully vaccinating at least 72% of the susceptible populace. Boosting
of vaccine-derived immunity was shown to be more beneficial than boosting of natural immunity. Overall, this study showed that the prospects of the elimination of
the pandemic in the United States were highly promising using the two intervention
measures.
ContributorsSafdar, Salman (Author) / Gumel, Abba (Thesis advisor) / Kostelich, Eric (Committee member) / Kang, Yun (Committee member) / Fricks, John (Committee member) / Espanol, Malena (Committee member) / Arizona State University (Publisher)
Created2023
Description
This dissertation develops a second order accurate approximation to the magnetic resonance (MR) signal model used in the PARSE (Parameter Assessment by Retrieval from Single Encoding) method to recover information about the reciprocal of the spin-spin relaxation time function (R2*) and frequency offset function (w) in addition to the typical steady-state transverse magnetization (M) from single-shot magnetic resonance imaging (MRI) scans. Sparse regularization on an approximation to the edge map is used to solve the associated inverse problem. Several studies are carried out for both one- and two-dimensional test problems, including comparisons to the first order approximation method, as well as the first order approximation method with joint sparsity across multiple time windows enforced. The second order accurate model provides increased accuracy while reducing the amount of data required to reconstruct an image when compared to piecewise constant in time models. A key component of the proposed technique is the use of fast transforms for the forward evaluation. It is determined that the second order model is capable of providing accurate single-shot MRI reconstructions, but requires an adequate coverage of k-space to do so. Alternative data sampling schemes are investigated in an attempt to improve reconstruction with single-shot data, as current trajectories do not provide ideal k-space coverage for the proposed method.
ContributorsJesse, Aaron Mitchel (Author) / Platte, Rodrigo (Thesis advisor) / Gelb, Anne (Committee member) / Kostelich, Eric (Committee member) / Mittelmann, Hans (Committee member) / Moustaoui, Mohamed (Committee member) / Arizona State University (Publisher)
Created2019