Matching Items (23)
Filtering by
Description
The main goal of this project is to study approximations of functions on circular and spherical domains using the cubed sphere discretization. On each subdomain, the function is approximated by windowed Fourier expansions. Of particular interest is the dependence of accuracy on the different choices of windows and the size of the overlapping regions. We use Matlab to manipulate each of the variables involved in these computations as well as the overall error, thus enabling us to decide which specific values produce the most accurate results. This work is motivated by problems arising in atmospheric research.
ContributorsSopa, Megan Grace (Author) / Platte, Rodrigo (Thesis director) / Kostelich, Eric (Committee member) / Department of Information Systems (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2018-05
Description
A specific species of the genus Geobacter exhibits useful electrical properties when processing a molecule often found in waste water. A team at ASU including Dr Cèsar Torres and Dr Sudeep Popat used that species to create a special type of solid oxide fuel cell we refer to as a microbial fuel cell. Identification of possible chemical processes and properties of the reactions used by the Geobacter are investigated indirectly by taking measurements using Electrochemical Impedance Spectroscopy of the electrode-electrolyte interface of the microbial fuel cell to obtain the value of the fuel cell's complex impedance at specific frequencies. Investigation of the multiple polarization processes which give rise to measured impedance values is difficult to do directly and so examination of the distribution function of relaxation times (DRT) is considered instead. The DRT is related to the measured complex impedance values using a general, non-physical equivalent circuit model. That model is originally given in terms of a Fredholm integral equation with a non-square integrable kernel which makes the inverse problem of determining the DRT given the impedance measurements an ill-posed problem. The original integral equation is rewritten in terms of new variables into an equation relating the complex impedance to the convolution of a function based upon the original integral kernel and a related but separate distribution function which we call the convolutional distribution function. This new convolutional equation is solved by reducing the convolution to a pointwise product using the Fourier transform and then solving the inverse problem by pointwise division and application of a filter function (equivalent to regularization). The inverse Fourier transform is then taken to get the convolutional distribution function. In the literature the convolutional distribution function is then examined and certain values of a specific, less general equivalent circuit model are calculated from which aspects of the original chemical processes are derived. We attempted to instead directly determine the original DRT from the calculated convolutional distribution function. This method proved to be practically less useful due to certain values determined at the time of experiment which meant the original DRT could only be recovered in a window which would not normally contain the desired information for the original DRT. This limits any attempt to extend the solution for the convolutional distribution function to the original DRT. Further research may determine a method for interpreting the convolutional distribution function without an equivalent circuit model as is done with the regularization method used to solve directly for the original DRT.
ContributorsBaker, Robert Simpson (Author) / Renaut, Rosemary (Thesis director) / Kostelich, Eric (Committee member) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
A semi-implicit, fourth-order time-filtered leapfrog numerical scheme is investigated for accuracy and stability, and applied to several test cases, including one-dimensional advection and diffusion, the anelastic equations to simulate the Kelvin-Helmholtz instability, and the global shallow water spectral model to simulate the nonlinear evolution of twin tropical cyclones. The leapfrog scheme leads to computational modes in the solutions to highly nonlinear systems, and time-filters are often used to damp these modes. The proposed filter damps the computational modes without appreciably degrading the physical mode. Its performance in these metrics is superior to the second-order time-filtered leapfrog scheme developed by Robert and Asselin.
ContributorsBurke, Lee Matthew Moore (Author) / Moustaoui, Mohamed (Thesis director) / Kostelich, Eric (Committee member) / School of Mathematical and Statistical Sciences (Contributor) / Mechanical and Aerospace Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
Glioblastoma multiforme (GBM) is a malignant, aggressive and infiltrative cancer of the central nervous system with a median survival of 14.6 months with standard care. Diagnosis of GBM is made using medical imaging such as magnetic resonance imaging (MRI) or computed tomography (CT). Treatment is informed by medical images and includes chemotherapy, radiation therapy, and surgical removal if the tumor is surgically accessible. Treatment seldom results in a significant increase in longevity, partly due to the lack of precise information regarding tumor size and location. This lack of information arises from the physical limitations of MR and CT imaging coupled with the diffusive nature of glioblastoma tumors. GBM tumor cells can migrate far beyond the visible boundaries of the tumor and will result in a recurring tumor if not killed or removed. Since medical images are the only readily available information about the tumor, we aim to improve mathematical models of tumor growth to better estimate the missing information. Particularly, we investigate the effect of random variation in tumor cell behavior (anisotropy) using stochastic parameterizations of an established proliferation-diffusion model of tumor growth. To evaluate the performance of our mathematical model, we use MR images from an animal model consisting of Murine GL261 tumors implanted in immunocompetent mice, which provides consistency in tumor initiation and location, immune response, genetic variation, and treatment. Compared to non-stochastic simulations, stochastic simulations showed improved volume accuracy when proliferation variability was high, but diffusion variability was found to only marginally affect tumor volume estimates. Neither proliferation nor diffusion variability significantly affected the spatial distribution accuracy of the simulations. While certain cases of stochastic parameterizations improved volume accuracy, they failed to significantly improve simulation accuracy overall. Both the non-stochastic and stochastic simulations failed to achieve over 75% spatial distribution accuracy, suggesting that the underlying structure of the model fails to capture one or more biological processes that affect tumor growth. Two biological features that are candidates for further investigation are angiogenesis and anisotropy resulting from differences between white and gray matter. Time-dependent proliferation and diffusion terms could be introduced to model angiogenesis, and diffusion weighed imaging (DTI) could be used to differentiate between white and gray matter, which might allow for improved estimates brain anisotropy.
ContributorsAnderies, Barrett James (Author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / Stepien, Tracy (Committee member) / Harrington Bioengineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
Chebfun is a collection of algorithms and an open-source software system in object-oriented Matlab that extends familiar powerful methods of numerical computation involving numbers to continuous or piecewise-continuous functions. The success of this strategy is based on the mathematical fact that smooth functions can be represented very efficiently by polynomial interpolation at Chebyshev points or by trigonometric interpolation at equispaced points for periodic functions. More recently, the system has been extended to handle bivariate functions and vector fields. These two new classes of objects are called Chebfun2 and Chebfun2v, respectively. We will show that Chebfun2 and Chebfun2v, and can be used to accurately and efficiently perform various computations on parametric surfaces in two or three dimensions, including path trajectories and mean and Gaussian curvatures. More advanced surface computations such as mean curvature flows are also explored. This is also the first work to use the newly implemented trigonometric representation, namely Trigfun, for computations on surfaces.
ContributorsPage-Bottorff, Courtney Michelle (Author) / Platte, Rodrigo (Thesis director) / Kostelich, Eric (Committee member) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
Magnetic resonance imaging (MRI) data of metastatic brain cancer patients at the Barrow Neurological Institute sparked interest in the radiology department due to the possibility that tumor size distributions might mimic a power law or an exponential distribution. In order to consider the question regarding the growth trends of metastatic brain tumors, this thesis analyzes the volume measurements of the tumor sizes from the BNI data and attempts to explain such size distributions through mathematical models. More specifically, a basic stochastic cellular automaton model is used and has three-dimensional results that show similar size distributions of those of the BNI data. Results of the models are investigated using the likelihood ratio test suggesting that, when the tumor volumes are measured based on assuming tumor sphericity, the tumor size distributions significantly mimic the power law over an exponential distribution.
ContributorsFreed, Rebecca (Co-author) / Snopko, Morgan (Co-author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / WPC Graduate Programs (Contributor) / School of Accountancy (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2018-12
Description
Coherent vortices are ubiquitous structures in natural flows that affect mixing and transport of substances and momentum/energy. Being able to detect these coherent structures is important for pollutant mitigation, ecological conservation and many other aspects. In recent years, mathematical criteria and algorithms have been developed to extract these coherent structures in turbulent flows. In this study, we will apply these tools to extract important coherent structures and analyze their statistical properties as well as their implications on kinematics and dynamics of the flow. Such information will aide representation of small-scale nonlinear processes that large-scale models of natural processes may not be able to resolve.
ContributorsCass, Brentlee Jerry (Author) / Tang, Wenbo (Thesis director) / Kostelich, Eric (Committee member) / Department of Information Systems (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2018-05
Description
Using weather data from the Weather Research and Forecasting model (WRF), we analyze the transport of inertial particles in Hurricane Katrina in order to identify coherent patterns of motion. For our analysis, we choose a Lagrangian approach instead of an Eulerian approach because the Lagrangian approach is objective and frame-independent, and gives results which are better defined. In particular, we locate Lagrangian Coherent Structures (LCS), which are smooth sets of fluid particles which are locally most hyperbolic (either attracting or repelling). We implement a variational method for locating LCS and compare the results to previous computation of LCS using Finite-Time Lyapunov Exponents (FTLE) to identify regions of high stretching in the fluid flow.
ContributorsDeibel, Angelica Rae (Author) / Tang, Wenbo (Thesis director) / Moustaoui, Mohamed (Committee member) / Kostelich, Eric (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Department of Physics (Contributor)
Created2013-05
Description
In this research we consider stochastic models of Glioblastoma Multiforme brain tumors. We first look at a model by K. Swanson et al., which describes the dynamics as random diffusion plus deterministic logistic growth. We introduce a stochastic component in the logistic growth in the form of a random growth rate defined by a Poisson process. We show that this stochastic logistic growth model leads to a more accurate evaluation of the tumor growth compared its deterministic counterpart. We also discuss future plans to incorporate individual patient geometry, extend the model to three dimensions and to incorporate effects of different treatments into our model, in collaboration with a local hospital.
ContributorsManning, Michael Clare (Author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / Gardner, Carl (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / School of Letters and Sciences (Contributor) / School of Human Evolution and Social Change (Contributor)
Created2013-12
Description
Tumor associated microglia-and-macrophages (TAMS) may constitute up to 30% of the composition of glioblastoma. Through mechanisms not well understood, TAMS are thought to aid the progression and invasiveness of glioblastoma. In an effort to investigate properties of TAMS in the context of glioblastoma, I utilized data from a PDGF-driven rat model of glioma that highly resembles human glioblastoma. Data was collected from time-lapse microscopy of slice cultures that differentially labels glioma cells and also microglia cells within and outside the tumor microenvironment. Here I show that microglia localize in the tumor and move with greater speed and migration than microglia outside the tumor environment. Following previous studies that show microglia can be characterized by certain movement distributions based on environmental influences, in this study, the majority of microglia movement was characterized by a power law distribution with a characteristic power law exponent lower than outside the tumor region. This indicates that microglia travel at greater distances within the tumor region than outside of it.
ContributorsJuliano, Joseph Dominic (Author) / Kostelich, Eric (Thesis director) / Nagy, John (Committee member) / Swanson, Kristin (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Department of Chemistry and Biochemistry (Contributor)
Created2013-12