Matching Items (48)
141463-Thumbnail Image.png
Description

Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of

Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm3 to 62 mm3, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.

ContributorsRutter, Erica (Author) / Stepien, Tracy (Author) / Anderies, Barrett (Author) / Plasencia, Jonathan (Author) / Woolf, Eric C. (Author) / Scheck, Adrienne C. (Author) / Turner, Gregory H. (Author) / Liu, Qingwei (Author) / Frakes, David (Author) / Kodibagkar, Vikram (Author) / Kuang, Yang (Author) / Preul, Mark C. (Author) / Kostelich, Eric (Author) / College of Liberal Arts and Sciences (Contributor)
Created2017-05-31
141494-Thumbnail Image.png
Description

Background:
Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and

Background:
Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.

Results:
We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.

Conclusions:
The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.

ContributorsKostelich, Eric (Author) / Kuang, Yang (Author) / McDaniel, Joshua (Author) / Moore, Nina Z. (Author) / Martirosyan, Nikolay L. (Author) / Preul, Mark C. (Author) / College of Liberal Arts and Sciences (Contributor)
Created2011-12-21
141503-Thumbnail Image.png
Description

This paper studies the effect of targeted observations on state and parameter estimates determined with Kalman filter data assimilation (DA) techniques. We first provide an analytical result demonstrating that targeting observations within the Kalman filter for a linear model can significantly reduce state estimation error as opposed to fixed or

This paper studies the effect of targeted observations on state and parameter estimates determined with Kalman filter data assimilation (DA) techniques. We first provide an analytical result demonstrating that targeting observations within the Kalman filter for a linear model can significantly reduce state estimation error as opposed to fixed or randomly located observations. We next conduct observing system simulation experiments for a chaotic model of meteorological interest, where we demonstrate that the local ensemble transform Kalman filter (LETKF) with targeted observations based on largest ensemble variance is skillful in providing more accurate state estimates than the LETKF with randomly located observations. Additionally, we find that a hybrid ensemble Kalman filter parameter estimation method accurately updates model parameters within the targeted observation context to further improve state estimation.

ContributorsBellsky, Thomas (Author) / Kostelich, Eric (Author) / Mahalov, Alex (Author) / College of Liberal Arts and Sciences (Contributor)
Created2014-06-01
Description

In the last two decades, fantasy sports have grown massively in popularity. Fantasy football in particular is the most popular fantasy sport in the United States. People spend hours upon hours every year building, researching, and perfecting their teams to compete with others for money or bragging rights. One problem,

In the last two decades, fantasy sports have grown massively in popularity. Fantasy football in particular is the most popular fantasy sport in the United States. People spend hours upon hours every year building, researching, and perfecting their teams to compete with others for money or bragging rights. One problem, however, is that National Football League (NFL) players are human and will not perform the same as they did last week or last season. Because of this, there is a need to create a machine learning model to help predict when players will have a tough game or when they can perform above average. This report discusses the history and science of fantasy football, gathering large amounts of player data, manipulating the information to create more insightful data points, creating a machine learning model, and how to use this tool in a real-world situation. The initial model created significantly accurate predictions for quarterbacks and running backs but not receivers and tight ends. Improvements significantly increased the accuracy by reducing the mean average error to below one for all positions, resulting in a successful model for all four positions.

ContributorsCase, Spencer (Author) / Johnson, Jarod (Co-author) / Kostelich, Eric (Thesis director) / Zhuang, Houlong (Committee member) / Barrett, The Honors College (Contributor) / Department of Psychology (Contributor) / Mechanical and Aerospace Engineering Program (Contributor)
Created2023-05
Description

The main purpose of this project is to create a method for determining the absolute position of an accelerometer. Acceleration and angular speed were obtained from an accelerometer attached to a vehicle as it moves around. As the vehicle moves to collect information the orientation of the accelerometer changes, so

The main purpose of this project is to create a method for determining the absolute position of an accelerometer. Acceleration and angular speed were obtained from an accelerometer attached to a vehicle as it moves around. As the vehicle moves to collect information the orientation of the accelerometer changes, so a rotation matrix is applied to the data based on the angular change at each time. The angular change and distance are obtained by using the trapezoidal approximation of the integrals. This method was first validated by using simple sets of "true" data which are explicitly known sets of data to compare the results to. Then, an analysis of how different time steps and levels of noise affect the error of the results was performed to determine the optimal time step of 0.1 sec that was then used for the actual tests. The tests that were performed were: a stationary test for uses of calibration, a straight line test to verify a simple test, and a closed loop test to test the accuracy. The graphs for these tests give no indication of the actual paths, so the final results can only show that the data from the accelerometer is too noisy and inaccurate for this method to be used by this sensor. The future work would be to test different ways to get more accurate data and then use it to verify this methods. These ways could include using more sensors to interpolate the data, reducing noise by using a different sensor, or adding a filter. Then, if this method is considered accurate enough, it could be implemented into control systems.

ContributorsHorner, Devon (Author) / Kostelich, Eric (Thesis director) / Crook, Sharon (Committee member) / Barrett, The Honors College (Contributor) / Mechanical and Aerospace Engineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor)
Created2023-05
168448-Thumbnail Image.png
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

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
187669-Thumbnail Image.png
Description
Advancements to a dual scale Large Eddy Simulation (LES) modeling approach for immiscible turbulent phase interfaces are presented. In the dual scale LES approach, a high resolution auxiliary grid, used to capture a fully resolved interface geometry realization, is linked to an LES grid that solves the filtered Navier-Stokes equations.

Advancements to a dual scale Large Eddy Simulation (LES) modeling approach for immiscible turbulent phase interfaces are presented. In the dual scale LES approach, a high resolution auxiliary grid, used to capture a fully resolved interface geometry realization, is linked to an LES grid that solves the filtered Navier-Stokes equations. Exact closure of the sub-filter interface terms is provided by explicitly filtering the fully resolved quantities from the auxiliary grid. Reconstructing a fully resolved velocity field to advance the phase interface requires modeling several sub-filter effects, including shear and accelerational instabilities and phase change. Two sub-filter models were developed to generate these sub-filter hydrodynamic instabilities: an Orr-Sommerfeld model and a Volume-of-Fluid (VoF) vortex sheet method. The Orr-Sommerfeld sub-filter model was found to be incompatible with the dual scale approach, since it is unable to generate interface rollup and a process to separate filtered and sub-filter scales could not be established. A novel VoF vortex sheet method was therefore proposed, since prior vortex methods have demonstrated interface rollup and following the LES methodology, the vortex sheet strength could be decomposed into its filtered and sub-filter components. In the development of the VoF vortex sheet method, it was tested with a variety of classical hydrodynamic instability problems, compared against prior work and linear theory, and verified using Direct Numerical Simulations (DNS). An LES consistent approach to coupling the VoF vortex sheet with the LES filtered equations is presented and compared against DNS. Finally, a sub-filter phase change model is proposed and assessed in the dual scale LES framework with an evaporating interface subjected to decaying homogeneous isotropic turbulence. Results are compared against DNS and the interplay between surface tension forces and evaporation are discussed.
ContributorsGoodrich, Austin Chase (Author) / Herrmann, Marcus (Thesis advisor) / Dahm, Werner (Committee member) / Kim, Jeonglae (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2023
187415-Thumbnail Image.png
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,

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
171849-Thumbnail Image.png
Description
This thesis focuses on the turbulent bluff body wakes in incompressible and compressible flows. An incompressible wake flow past an axisymmetric body of revolution at a diameter-based Reynolds number Re=5000 is investigated via a direct numerical simulation. It is followed by the development of a compressible solver using a split-form

This thesis focuses on the turbulent bluff body wakes in incompressible and compressible flows. An incompressible wake flow past an axisymmetric body of revolution at a diameter-based Reynolds number Re=5000 is investigated via a direct numerical simulation. It is followed by the development of a compressible solver using a split-form discontinuous Galerkin spectral element method framework with shock capturing. In the study on incompressible wake flows, three dominant coherent vortical motions are identified in the wake: the vortex shedding motion with the frequency of St=0.27, the bubble pumping motion with St=0.02, and the very-low-frequency (VLF) motion originated in the very near wake of the body with the frequencies St=0.002 and 0.005. The very-low-frequency motion is associated with a slow precession of the wake barycenter. The vortex shedding pattern is demonstrated to follow a reflectional symmetry breaking mode, with the detachment location rotating continuously and making a full circle over one vortex shedding period. The VLF radial motion with St=0.005 originates as m = 1 mode, but later transitions into m = 2 mode in the intermediate wake. Proper orthogonaldecomposition (POD) and dynamic mode decomposition (DMD) are further performed to analyze the spatial structure associated with the dominant coherent motions. Results of the POD and DMD analysis are consistent with the results of the azimuthal Fourier analysis. To extend the current incompressible code to be able to solve compressible flows, a computational methodology is developed using a high-order approximation for the compressible Navier-Stokes equations with discontinuities. The methodology is based on a split discretization framework with a summation-by-part operator. An entropy viscosity method and a subcell finite volume method are implemented to capture discontinuities. The developed high-order split-form with shock-capturing methodology is subject to a series of evaluation on cases from subsonic to hypersonic, from one-dimensional to three dimensional. The Taylor-Green vortex case and the supersonic sphere wake case show the capability to handle three-dimensional turbulent flows without and with the presence of shocks. It is also shown that higher-order approximations yield smaller errors than lower-order approximations, for the same number of total degrees of freedom.
ContributorsZhang, Fengrui (Author) / Peet, Yulia (Thesis advisor) / Kostelich, Eric (Committee member) / Kim, Jeonglae (Committee member) / Hermann, Marcus (Committee member) / Adrian, Ronald (Committee member) / Arizona State University (Publisher)
Created2022
Description
Glioblastoma Multiforme is a prevalent and aggressive brain tumor. It has an average 5-year survival rate of 6% and average survival time of 14 months. Using patient-specific MRI data from the Barrow Neurological Institute, this thesis investigates the impact of parameter manipulation on reaction-diffusion models for predicting and simulating glioblastoma

Glioblastoma Multiforme is a prevalent and aggressive brain tumor. It has an average 5-year survival rate of 6% and average survival time of 14 months. Using patient-specific MRI data from the Barrow Neurological Institute, this thesis investigates the impact of parameter manipulation on reaction-diffusion models for predicting and simulating glioblastoma growth. The study aims to explore key factors influencing tumor morphology and to contribute to enhancing prediction techniques for treatment.
ContributorsShayegan, Tara (Author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / Barrett, The Honors College (Contributor) / School of Human Evolution & Social Change (Contributor)
Created2024-05