Physical Sciences and Mathematics Commons^™

Solving The Cable Equation, A Second-Order Time Dependent Pde For Non-Ideal Cables With Action Potentials In The Mammalian Brain Using Kss Methods, Nirmohi Charbe

Master's Theses

In this thesis we shall perform the comparisons of a Krylov Subspace Spectral method with Forward Euler, Backward Euler and Crank-Nicolson to solve the Cable Equation. The Cable Equation measures action potentials in axons in a mammalian brain treated as an ideal cable in the first part of the study. We shall subject this problem to the further assumption of a non-ideal cable. Assume a non-uniform cross section area along the longitudinal axis. At the present time, the effects of torsion, curvature and material capacitance are ignored. There is particular interest to generalize the application of the PDEs including and …

Learning The Game: Implementations Of Convolutional Networks In Automated Strategy Identification, Cameron Klig

Master's Theses

Games can be used to represent a wide variety of real world problems, giving rise to many applications of game theory. Various computational methods have been proposed for identifying game strategies, including optimized tree search algorithms, game-specific heuristics, and artificial intelligence. In the last decade, systems like AlphaGo and AlphaZero have significantly exceeded the performance of the best human players in Chess, Go, and other games. The most effective game engines to date employ convolutional neural networks (CNNs) to evaluate game boards, extract features, and predict the optimal next move. These engines are trained on billions of simulated games, wherein …

Distributed Control Of Servicing Satellite Fleet Using Horizon Simulation Framework, Scott Plantenga

Master's Theses

On-orbit satellite servicing is critical to maximizing space utilization and sustainability and is of growing interest for commercial, civil, and defense applications. Reliance on astronauts or anchored robotic arms for the servicing of next-generation large, complex space structures operating beyond Low Earth Orbit is impractical. Substantial literature has investigated the mission design and analysis of robotic servicing missions that utilize a single servicing satellite to approach and service a single target satellite. This motivates the present research to investigate a fleet of servicing satellites performing several operations for a large, central space structure.

This research leverages a distributed control approach, …

Unique Signed Minimal Wiring Diagrams And The Stanley-Reisner Correspondence, Vanessa Newsome-Slade

Master's Theses

Biological systems are commonly represented using networks consisting of interactions between various elements in the system. Reverse engineering, a method of mathematical modeling, is used to recover how the elements in the biological network are connected. These connections are encoded using wiring diagrams, which are directed graphs that describe how elements in a network affect one another. A signed wiring diagram provides additional information about the interactions between elements relating to activation and inhibition. Due to cost concerns, it is optimal to gain insight into biological networks with as few experiments and data as possible. Minimal wiring diagrams identify the …

A Network Analysis Of Covid-19 In The United States, Joseph C. Mcguire

Master's Theses

Through methods in network theory and time-series analysis, we will analyze the spread of COVID-19 in the United States by determining trends in state-by-state daily cases through a network construction. Previous researchers have found frameworks for approximating the spread of the COVID-19 pandemic and identifying potential rises in cases by a network construction based on correlation of cases between regions [1]. Applying this network construction we determine how this network and its structure act as a predictor for overall COVID-19 cases in the United States by preforming a trend analysis on a variety of network statistics and US COVID-19 cases.

Large Scale Disease Modeling, Walker Mattox

Master's Theses

In this we study large scale disease modeling. After understanding the mechanics behind the SIR disease model in an ODE sense, we will apply this knowledge to model disease spread in more and more increasing advanced cellular automata. Eventually, some of our cellular automata will include long distance travel. From this discrete data, we can then build an SIR model in the PDE sense to display large scale disease spread.

An Optimization Model For Minimization Of Systemic Risk In Financial Portfolios, Zachary Alexander Gelber

Master's Theses

In this thesis, we study how sovereign credit default swaps are able to measure systemic risk as well as how they can be used to construct optimal portfolios to minimize risk. We define the clustering coefficient as a proxy for systemic risk and design an optimization problem with the goal of minimizing the mean absolute deviation of the clustering coefficient on a group of nine European countries. Additionally, we define a metric we call the diversity score that measures the diversification of any given portfolio. We solve this problem for a baseline set of parameters, then spend the remainder of …

Multi-Valued Solutions For The Equation Of Motion, Darcy-Jordan Model, As A Cauchy Problem: A Shocking Event, Chandler Shimp

Master's Theses

Shocks are physical phenomenon that occur quite often around us. In this thesis we examine the occurrence of shocks in finite amplitude acoustic waves from a numerical perspective. These waves, or jump discontinuities, yield ill-behaved solutions when solved numerically. This study takes on the challenge of finding both single- and multi-valued solutions.

The previously unsolved problem in this study is the representation of the Equation of Motion (EoM) in the form of the Darcy-Jordan model (DJM) and expressed as a dimensionless IVP Cauchy problem. Prior attempts to solve have resulted only in implicit solutions or explicit solutions with certain initial …

On A Stochastic Model Of Epidemics, Rachel Prather

Master's Theses

This thesis examines a stochastic model of epidemics initially proposed and studied by Norman T.J. Bailey [1]. We discuss some issues with Bailey's stochastic model and argue that it may not be a viable theoretical platform for a more general epidemic model. A possible alternative approach to the solution of Bailey's stochastic model and stochastic modeling is proposed as well. Regrettably, any further study on those proposals will have to be discussed elsewhere due to a time constraint.

A Modified Preconditioned Conjugate Gradient Method For Approximating The Scattering Amplitude, Samson Ayo

Master's Theses

In this thesis, we look at an iterative method for approximating the scattering amplitude that involves solving two linear systems: a forward system Ax=b and an adjoint system A^Ty=g. Once these two systems are solved, the scattering amplitude, defined by g^Tx=y^Tb is easily obtained.

We derive a conjugate gradient-like iteration for a nonsymmetric saddle point matrix that is constructed to have a real positive spectrum. We investigate the use of Schur Complement preconditioners with block-diagonal factorization to speed up the convergence of our method and compare …

Rapid Implicit Diagonalization Of Variable-Coefficient Differential Operators Using The Uncertainty Principle, Carley Walker

Master's Theses

We propose to create a new numerical method for a class of time-dependent PDEs (second-order, one space dimension, Dirichlet boundary conditions) that can be used to obtain more accurate and reliable solutions than traditional methods. Previously, it was shown that conventional time-stepping methods could be avoided for time-dependent mathematical models featuring a finite number of homogeneous materials, thus assuming general piecewise constant coefficients. This proposed method will avoid the modeling shortcuts that are traditionally taken, and it will generalize the piecewise constant case of energy diffusion and wave propagation to work for an infinite number of smaller pieces, or a …

Diagonalization Of 1-D Schrodinger Operators With Piecewise Constant Potentials, Sarah Wright

Master's Theses

In today's world our lives are very layered. My research is meant to adapt current inefficient numerical methods to more accurately model the complex situations we encounter. This project focuses on a specific equation that is used to model sound speed in the ocean. As depth increases, the sound speed changes. This means the variable related to the sound speed is not constant. We will modify this variable so that it is piecewise constant. The specific operator in this equation also makes current time-stepping methods not practical. The method used here will apply an eigenfunction expansion technique used in previous …

Stability Analysis Of Krylov Subspace Spectral Methods For The 1-D Wave Equation In Inhomogeneous Media, Bailey Rester

Master's Theses

Krylov subspace spectral (KSS) methods are high-order accurate, explicit time-stepping methods for partial differential equations (PDEs) that also possess the stability characteristic of implicit methods. Unlike other time-stepping approaches, KSS methods compute each Fourier coefficient of the solution from an individualized approximation of the solution operator of the PDE. As a result, KSS methods scale effectively to higher spatial resolution. This thesis will present a stability analysis of a first-order KSS method applied to the wave equation in inhomogeneous media.

An Adaptive Approach To Gibbs’ Phenomenon, Jannatul Ferdous Chhoa

Master's Theses

Gibbs’ Phenomenon, an unusual behavior of functions with sharp jumps, is encountered while applying the Fourier Transform on them. The resulting reconstructions have high frequency oscillations near the jumps making the reconstructions far from being accurate. To get rid of the unwanted oscillations, we used the Lanczos sigma factor to adjust the Fourier series and we came across three cases. Out of the three, two of them failed to give us the right reconstructions because either it was removing the oscillations partially but not entirely or it was completely removing them but smoothing out the jumps a little too much. …

An Application Of The Unscented Kalman Filter For Spacecraft Attitude Estimation On Real And Simulated Light Curve Data, Kent A. Rush

Master's Theses

In the past, analyses of lightcurve data have been applied to asteroids in order to determine their axis of rotation, rotation rate and other parameters. In recent decades, these analyses have begun to be applied in the domain of Earth orbiting spacecraft. Due to the complex geometry of spacecraft and the wide variety of parameters that can influence the way in which they reflect light, these analyses require more complex assumptions and a greater knowledge about the object being studied. Previous investigations have shown success in extracting attitude parameters from unresolved spacecraft using simulated data. This paper presents a focused …

A Study Of The Design Of Adaptive Camber Winglets, Justin J. Rosescu

Master's Theses

A numerical study was conducted to determine the effect of changing the camber of a winglet on the efficiency of a wing in two distinct flight conditions. Camber was altered via a simple plain flap deflection in the winglet, which produced a constant camber change over the winglet span. Hinge points were located at 20%, 50% and 80% of the chord and the trailing edge was deflected between -5° and +5°. Analysis was performed using a combination of three-dimensional vortex lattice method and two-dimensional panel method to obtain aerodynamic forces for the entire wing, based on different winglet camber configurations. …

Using Modern Portfolio Theory To Analyze Virgil's Aeneid (Or Any Other Poem), David Patterson

Master's Theses

This paper demonstrates that it is possible to use mathematics to study literature as it has been used to study the social sciences. By focusing on mathematically defining economic and literary terms, it can be shown that the underlying mathematical structure behind key concepts in economics and literature are analogous. This opens the possibility of applying economic models in literature. Specifically, it is demonstrated that the economic mathematical model of modern portfolio theory can answer long standing questions around the Roman epic Aeneid by Virgil. The poet died before completing his poem. The relative completeness of the books of the …

Multi-Point Flux Approximations Via The O-Method, Christen Leggett

Master's Theses

When an oil refining company is drilling for oil, much of the oil gets left behind after the first drilling. Enhanced oil recovery techniques can be used to recover more of that oil, but these methods are quite expensive. When a company is deciding if it is worth their time and money to use enhanced oil recovery methods, simulations can be used to model oil flow, showing the behavior and location of the oil. While methods do exist to model this flow, these methods are often very slow and inaccurate due to a large domain and wide variance in coefficients. …

Krylov Subspace Spectral Methods With Non-Homogenous Boundary Conditions, Abbie Hendley

Master's Theses

For this thesis, Krylov Subspace Spectral (KSS) methods, developed by Dr. James Lambers, will be used to solve a one-dimensional, heat equation with non-homogenous boundary conditions. While current methods such as Finite Difference are able to carry out these computations efficiently, their accuracy and scalability can be improved. We will solve the heat equation in one-dimension with two cases to observe the behaviors of the errors using KSS methods. The first case will implement KSS methods with trigonometric initial conditions, then another case where the initial conditions are polynomial functions. We will also look at both the time-independent and time-dependent …

The Martingale Approach To Financial Mathematics, Jordan M. Rowley

Master's Theses

In this thesis, we will develop the fundamental properties of financial mathematics, with a focus on establishing meaningful connections between martingale theory, stochastic calculus, and measure-theoretic probability. We first consider a simple binomial model in discrete time, and assume the impossibility of earning a riskless profit, known as arbitrage. Under this no-arbitrage assumption alone, we stumble upon a strange new probability measure Q, according to which every risky asset is expected to grow as though it were a bond. As it turns out, this measure Q also gives the arbitrage-free pricing formula for every asset on our market. In …

The Long-Run Effects Of Tropical Cyclones On Infant Mortality, Isabel Miranda

Master's Theses

In the United States alone, each tropical cyclone causes an average of $14.6 billion worth of damages. In addition to the destruction of physical infrastructure, natural disasters also negatively impact human capital formation. These losses are often more difficult to observe, and therefore, are over looked when quantifying the true costs of natural disasters. One particular effect is an increase in infant mortality rates, an important indicator of a country’s general socioeconomic level. This paper utilizes a model created by Anttila-Hughes and Hsiang, that takes advantage of annual variation in tropical cyclones using annual spatial average maximum wind speeds and …

Automatic Construction Of Scalable Time-Stepping Methods For Stiff Pdes, Vivian Ashley Montiforte

Master's Theses

Krylov Subspace Spectral (KSS) Methods have been demonstrated to be highly scalable time-stepping methods for stiff nonlinear PDEs. However, ensuring this scalability requires analytic computation of frequency-dependent quadrature nodes from the coefficients of the spatial differential operator. This thesis describes how this process can be automated for various classes of differential operators to facilitate public-domain software implementation.

Rotordynamic Analysis Of Theoretical Models And Experimental Systems, Cameron R. Naugle

Master's Theses

This thesis is intended to provide fundamental information for the construction and

analysis of rotordynamic theoretical models, and their comparison the experimental

systems. Finite Element Method (FEM) is used to construct models using Timoshenko

beam elements with viscous and hysteretic internal damping. Eigenvalues

and eigenvectors of state space equations are used to perform stability analysis, produce

critical speed maps, and visualize mode shapes. Frequency domain analysis

of theoretical models is used to provide Bode diagrams and in experimental data

full spectrum cascade plots. Experimental and theoretical model analyses are used

to optimize the control algorithm for an Active Magnetic Bearing …

Eignefunctions For Partial Differential Equations On Two-Dimensional Domains With Piecewise Constant Coefficients, Abdullah Muheel Momit Aurko

Master's Theses

In this thesis, we develop a highly accurate and efficient algorithm for computing the solution of a partial differential equation defined on a two-dimensional domain with discontinuous coefficients. An example of such a problem is for modeling the diffusion of heat energy in two space dimensions, in the case where the spatial domain represents a medium consisting of two different but homogeneous materials, with periodic boundary conditions.

Since diffusivity changes based on the material, it will be represented using a piecewise constant function, and this results in the formation of a complicated mathematical model. Such a model is impossible to …

Modern Fair-Weather And Storm Sediment Transport Around Ship Island, Mississippi: Implications For Coastal Habitats And Restoration Efforts, Eve Rettew Eisemann

Master's Theses

The Mississippi – Alabama barrier island chain is experiencing accelerated sea level rise, decreased sediment supply, and frequent hurricane impacts. These three factors drive unprecedented rates of morphology change and ecosystem reduction. All islands in the chain have experienced land loss on the order of hectares per year since records began in the 1840s. In 1969, Hurricane Camille impacted as a Category 5, breaching Ship Island, and significantly reduced viable seagrass habitat. Hurricane Katrina impacted as a Category 3 in 2005, further widening Camille Cut. To better understand the sustainability of these important islands and the ecosystems they support, sediment …

Krylov Subspace Spectral Method With Multigrid For A Time-Dependent, Variable-Coefficient Partial Differential Equation, Haley Renee Dozier

Master's Theses

Krylov Subspace Spectral (KSS) methods are traditionally used to solve time-dependent, variable-coefficient PDEs. They are high-order accurate, component-wise methods that are efficient with variable input sizes.

This thesis will demonstrate how one can make KSS methods even more efficient by using a Multigrid-like approach for low-frequency components. The essential ingredients of Multigrid, such as restriction, residual correction, and prolongation, are adapted to the timedependent case. Then a comparison of KSS, KSS with Multigrid, KSS-EPI and standard Krylov projection methods will be demonstrated.

Approximation Of The Scattering Amplitude Using Nonsymmetric Saddle Point Matrices, Amber Sumner Robertson

Master's Theses

In this thesis we look at iterative methods for solving the primal (Ax = b) and dual (A^T y = g) systems of linear equations to approximate the scattering amplitude defined by g^Tx =y^Tb. We use a conjugate gradient-like iteration for a unsymmetric saddle point matrix that is contructed so as to have a real positive spectrum. We find that this method is more consistent than known methods for computing the scattering amplitude such as GLSQR or QMR. Then, we use techniques from "matrices, moments, and quadrature" to compute the scattering amplitude …

Smarticles: A Method For Identifying And Correcting Instability And Error Caused By Explicit Integration Techniques In Physically Based Simulations, Susan Aileen Marano

Master's Theses

Using an explicit integration method in physically based animations has many advantages including conceptual and computational simplicity, however, it re- quires small time steps to ensure low numerical instability. Simulations with large numbers of individually interacting components such as cloth, hair, and fluid models, are limited by the sections of particles most susceptible to error. This results in the need for smaller time steps than required for the majority of the system. These sections can be diverse and dynamic, quickly changing in size and location based on forces in the system. Identifying and handling these trou- blesome sections could allow …

An Examination Of The Yang-Baxter Equation, Alexandru Cibotarica

Master's Theses

The Yang-Baxter equation has been extensively studied due to its application in numerous fields of mathematics and physics. This thesis sets out to analyze the equation from the viewpoint of the algebraic product of matrices, i.e., the composition of linear maps, with the intent of characterizing the solutions of the Yang-Baxter equation.

We begin by examining the simple case of 22 matrices where it is possible to fully characterize the solutions. We connect the Yang-Baxter equation to the Cecioni-Frobenius Theorem and focus on obtaining solutions to the Yang-Baxter equation for special matrices where solutions are more easily found. Finally, …

To Shift Or Not To Shift: Maximizing The Efficiency Of A Wireless Sensor Network Using Multiple Sinks, Chabli Boler

Master's Theses

Energy efficiency is a vital part of wireless sensor networks (WSN). Using correct techniques, we can improve the overall data throughput of the network and, at the same time, cut down the cost of running it. The efficient utilization of the life of the sensor node in a WSN is a critical problem. This is a problem since the sensor nodes closest to the sink node expend more energy because they send more messages, mostly for routing, than any other node on the network.

This decreases the life of those sensor nodes and will lead to dead spots in the …