Matching Items (18)
Description
The theme for this work is the development of fast numerical algorithms for sparse optimization as well as their applications in medical imaging and source localization using sensor array processing. Due to the recently proposed theory of Compressive Sensing (CS), the $\ell_1$ minimization problem attracts more attention for its ability to exploit sparsity. Traditional interior point methods encounter difficulties in computation for solving the CS applications. In the first part of this work, a fast algorithm based on the augmented Lagrangian method for solving the large-scale TV-$\ell_1$ regularized inverse problem is proposed. Specifically, by taking advantage of the separable structure, the original problem can be approximated via the sum of a series of simple functions with closed form solutions. A preconditioner for solving the block Toeplitz with Toeplitz block (BTTB) linear system is proposed to accelerate the computation. An in-depth discussion on the rate of convergence and the optimal parameter selection criteria is given. Numerical experiments are used to test the performance and the robustness of the proposed algorithm to a wide range of parameter values. Applications of the algorithm in magnetic resonance (MR) imaging and a comparison with other existing methods are included. The second part of this work is the application of the TV-$\ell_1$ model in source localization using sensor arrays. The array output is reformulated into a sparse waveform via an over-complete basis and study the $\ell_p$-norm properties in detecting the sparsity. An algorithm is proposed for minimizing a non-convex problem. According to the results of numerical experiments, the proposed algorithm with the aid of the $\ell_p$-norm can resolve closely distributed sources with higher accuracy than other existing methods.
ContributorsShen, Wei (Author) / Mittlemann, Hans D (Thesis advisor) / Renaut, Rosemary A. (Committee member) / Jackiewicz, Zdzislaw (Committee member) / Gelb, Anne (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2011
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
In vertebrate outer retina, changes in the membrane potential of horizontal cells affect the calcium influx and glutamate release of cone photoreceptors via a negative feedback. This feedback has a number of important physiological consequences. One is called background-induced flicker enhancement (BIFE) in which the onset of dim background enhances the center flicker response of horizontal cells. The underlying mechanism for the feedback is still unclear but competing hypotheses have been proposed. One is the GABA hypothesis, which states that the feedback is mediated by gamma-aminobutyric acid (GABA), an inhibitory neurotransmitter released from horizontal cells. Another is the ephaptic hypothesis, which contends that the feedback is non-GABAergic and is achieved through the modulation of electrical potential in the intersynaptic cleft between cones and horizontal cells. In this study, a continuum spine model of the cone-horizontal cell synaptic circuitry is formulated. This model, a partial differential equation system, incorporates both the GABA and ephaptic feedback mechanisms. Simulation results, in comparison with experiments, indicate that the ephaptic mechanism is necessary in order for the model to capture the major spatial and temporal dynamics of the BIFE effect. In addition, simulations indicate that the GABA mechanism may play some minor modulation role.
ContributorsChang, Shaojie (Author) / Baer, Steven M. (Thesis advisor) / Gardner, Carl L (Thesis advisor) / Crook, Sharon M (Committee member) / Kuang, Yang (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2012
Description
This project attempts to create an accurate numerical simulation of the eastern limb of the HH 901 jet in the Mystic Mountain formation located in the Carina Nebula. Using a 3rd order accurate WENO numerical scheme in space, and a 3rd order accurate RK method in time, the temperature, density, radiative cooling, length, and average jet velocity of this astrophysical phenomenon were simulated based on observations made by Hubble Space Telescope and the work of Reiter and Smith (2013) and (2014). The results of this simulation are displayed in three figures, one each for temperature, radiative cooling, and density, which show a jet displaying morphology consistent with that of the HH 901 eastern limb without adjustment for stellar wind. Also discussed are the effects of different jet speeds, initial conditions, and pulse parameters on the shape and behavior of the simulated jets, as well as continuing work to be done on the simulation to enhance its accuracy and usefulness.
ContributorsKreitzer, Kyle (Author) / Gardner, Carl (Thesis director) / Jones, Jeremiah (Committee member) / School of Mathematical and Statistical Sciences (Contributor) / Mechanical and Aerospace Engineering Program (Contributor, Contributor) / Barrett, The Honors College (Contributor)
Created2019-05
Description
The purpose of this thesis is to accurately simulate the surface brightness in various spectral emission lines of the HH 901 jets in the Mystic Mountain Formation of the Carina Nebula. To accomplish this goal, we gathered relevant spectral emission line data for [Fe II] 12660 Å, Hα 6563 Å, and [S II] 6720 Å to compare with Hubble Space Telescope observations of the HH 901 jets presented in Reiter et al. (2016). We derived the emissivities for these lines from the spectral synthesis code Cloudy by Ferland et al. (2017). In addition, we used WENO simulations of density, temperature, and radiative cooling to model the jet. We found that the computed surface brightness values agreed with most of the observational surface brightness values. Thus, the 3D cylindrically symmetric simulations of surface brightness using the WENO code and Cloudy spectral emission models are accurate for jets like HH 901. After detailing these agreements, we discuss the next steps for the project, like adding an external ambient wind and performing the simulations in full 3D.
ContributorsMohan, Arun (Author) / Gardner, Carl (Thesis director) / Jones, Jeremiah (Committee member) / Computer Science and Engineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2020-05
Description
The retina is the lining in the back of the eye responsible for vision. When light photons hits the retina, the photoreceptors within the retina respond by sending impulses to the optic nerve, which connects to the brain. If there is injury to the eye or heredity retinal problems, this part can become detached. Detachment leads to loss of nutrients, such as oxygen and glucose, to the cells in the eye and causes cell death. Sometimes the retina is able to be surgically reattached. If the photoreceptor cells have not died and the reattachment is successful, then these cells are able to regenerate their outer segments (OS) which are essential for their functionality and vitality. In this work we will explore how the regrowth of the photoreceptor cells in a healthy eye after retinal detachment can lead to a deeper understanding of how eye cells take up nutrients and regenerate. This work uses a mathematical model for a healthy eye in conjunction with data for photoreceptors' regrowth and decay. The parameters for the healthy eye model are estimated from the data and the ranges of these parameter values are centered +/- 10\% away from these values are used for sensitivity analysis. Using parameter estimation and sensitivity analysis we can better understand how certain processes represented by these parameters change within the model as a result of retinal detachment. Having a deeper understanding for any sort of photoreceptor death and growth can be used by the greater scientific community to help with these currently irreversible conditions that lead to blindness, such as retinal detachment. The analysis in this work shows that maximizing the carrying capacity of the trophic pool and the rate of RDCVF, as well as minimizing nutrient withdrawal of the rods and the cones from the trophic pool results in both the most regrowth and least cell death in retinal detachment.
ContributorsGoldman, Miriam Ayla (Author) / Camacho, Erikia (Thesis director) / Wirkus, Stephen (Committee member) / School of Mathematical and Natural Sciences (Contributor, Contributor) / Barrett, The Honors College (Contributor)
Created2017-05
Description
The purpose of this thesis is to accurately simulate in 3D the HH901 jet in the Mystic Mountain Formation of the Carina Nebula. Astronomers present a narrow-band Wide Field Camera image of Carina and the morphology of some astrophysical jets, including HH901. The simulation attempts to replicate features of the jet, among which are pulses, bow shock, terminal Mach disk, and Kelvin-Helmholtz rollup. We use the gas dynamical equations to solve for density, velocity, and temperature. The numerical methods used to solve the equations are third-order WENO (weighted essentially non-oscillatory) and third-order Runge-Kutta. Graphs of density and radiative cooling demonstrate the effect of adding wind (nonzero ambient velocity). The paper discusses the altering of the ambient velocity and final time to fit the shape of the jet in the Hubble image. The suggested next steps are simulating the other HH901 jet and comparing the jets’ atomic makeups to advance understanding of astrophysical jets.
ContributorsBuyer, Michael (Author) / Gardner, Carl (Thesis director) / Jones, Jeremiah (Committee member) / Barrett, The Honors College (Contributor) / Computer Science and Engineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor)
Created2022-05
Description
Silicon photonic technology continues to dominate the solar industry driven by steady improvement in device and module efficiencies. Currently, the world record conversion efficiency (~26.6%) for single junction silicon solar cell technologies is held by silicon heterojunction (SHJ) solar cells based on hydrogenated amorphous silicon (a-Si:H) and crystalline silicon (c-Si). These solar cells utilize the concept of carrier selective contacts to improve device efficiencies. A carrier selective contact is designed to optimize the collection of majority carriers while blocking the collection of minority carriers. In the case of SHJ cells, a thin intrinsic a-Si:H layer provides crucial passivation between doped a-Si:H and the c-Si absorber that is required to create a high efficiency cell. There has been much debate regarding the role of the intrinsic a-Si:H passivation layer on the transport of photogenerated carriers, and its role in optimizing device performance. In this work, a multiscale model is presented which utilizes different simulation methodologies to study interfacial transport across the intrinsic a-Si:H/c-Si heterointerface and through the a-Si:H passivation layer. In particular, an ensemble Monte Carlo simulator was developed to study high field behavior of photogenerated carriers at the intrinsic a-Si:H/c-Si heterointerface, a kinetic Monte Carlo program was used to study transport of photogenerated carriers across the intrinsic a-Si:H passivation layer, and a drift-diffusion model was developed to model the behavior in the quasi-neutral regions of the solar cell. This work reports de-coupled and self-consistent simulations to fully understand the role and effect of transport across the a-Si:H passivation layer in silicon heterojunction solar cells, and relates this to overall solar cell device performance.
ContributorsMuralidharan, Pradyumna (Author) / Goodnick, Stephen M (Thesis advisor) / Vasileska, Dragica (Thesis advisor) / Honsberg, Christiana (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2019
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
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