Matching Items (2)
Filtering by
- All Subjects: Ordinary differential equations
- Creators: Wirkus, Stephen
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
Description
Patients suffering from Retinitis Pigmentosa (RP), the most common type of inherited retinal degeneration, experience irreversible vision loss due to photoreceptor degeneration. The preservation of cone photoreceptors has been deemed medically relevant as a therapy aimed at preventing blindness in patients with RP. Cones rely on aerobic glycolysis to supply the metabolites necessary for outer segment (OS) renewal and maintenance. The rod-derived cone viability factor (RdCVF), a protein secreted by the rod photoreceptors that preserves the cones, accelerates the flow of glucose into the cone cell stimulating aerobic glycolysis. This dissertation presents and analyzes ordinary differential equation (ODE) models of cellular and molecular level photoreceptor interactions in health and disease to examine mechanisms leading to blindness in patients with RP.
First, a mathematical model composed of four ODEs is formulated to investigate the progression of RP, accounting for the new understanding of RdCVF’s role in enhancing cone survival. A mathematical analysis is performed, and stability and bifurcation analyses are used to explore various pathways to blindness. Experimental data are used for parameter estimation and model validation. The numerical results are framed in terms of four stages in the progression of RP. Sensitivity analysis is used to determine mechanisms that have a significant affect on the cones at each stage of RP. Utilizing a non-dimensional form of the RP model, a numerical bifurcation analysis via MATCONT revealed the existence of stable limit cycles at two stages of RP.
Next, a novel eleven dimensional ODE model of molecular and cellular level interactions is described. The subsequent analysis is used to uncover mechanisms that affect cone photoreceptor functionality and vitality. Preliminary simulations show the existence of oscillatory behavior which is anticipated when all processes are functioning properly. Additional simulations are carried out to explore the impact of a reduction in the concentration of RdCVF coupled with disruption in the metabolism associated with cone OS shedding, and confirms cone-on-rod reliance. The simulation results are compared with experimental data. Finally, four cases are considered, and a sensitivity analysis is performed to reveal mechanisms that significantly impact the cones in each case.
First, a mathematical model composed of four ODEs is formulated to investigate the progression of RP, accounting for the new understanding of RdCVF’s role in enhancing cone survival. A mathematical analysis is performed, and stability and bifurcation analyses are used to explore various pathways to blindness. Experimental data are used for parameter estimation and model validation. The numerical results are framed in terms of four stages in the progression of RP. Sensitivity analysis is used to determine mechanisms that have a significant affect on the cones at each stage of RP. Utilizing a non-dimensional form of the RP model, a numerical bifurcation analysis via MATCONT revealed the existence of stable limit cycles at two stages of RP.
Next, a novel eleven dimensional ODE model of molecular and cellular level interactions is described. The subsequent analysis is used to uncover mechanisms that affect cone photoreceptor functionality and vitality. Preliminary simulations show the existence of oscillatory behavior which is anticipated when all processes are functioning properly. Additional simulations are carried out to explore the impact of a reduction in the concentration of RdCVF coupled with disruption in the metabolism associated with cone OS shedding, and confirms cone-on-rod reliance. The simulation results are compared with experimental data. Finally, four cases are considered, and a sensitivity analysis is performed to reveal mechanisms that significantly impact the cones in each case.
ContributorsBrager, Danielle Christine (Author) / Camacho, Erika (Thesis advisor) / Wirkus, Stephen (Thesis advisor) / Fricks, John (Committee member) / Gardner, Carl (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2020