Numerical Analysis and Computation Commons^™

Mathematical Models: The Lanchester Equations And The Zombie Apocalypse, Hailey Bauer

Undergraduate Theses and Capstone Projects

This research study used mathematical models to analyze and depicted specific battle situations and the outcomes of the zombie apocalypse. The original models that predicted warfare were the Lanchester models, while the zombie apocalypse models were fictional expansions upon mathematical models used to examine infectious diseases. In this paper, I analyzed and compared different mathematical models by examining each model’s set of assumptions and the impact of the change in variables on the population classes. The purpose of this study was to understand the basics of the discrete dynamical systems and to determine the similarities between imaginary and realistic models. …

Reducing Memory Access Latencies Using Data Compression In Sparse, Iterative Linear Solvers, Neil Lindquist

All College Thesis Program, 2016-2019

Solving large, sparse systems of linear equations plays a significant role in certain scientific computations, such as approximating the solutions of partial differential equations. However, solvers for these types of problems usually spend most of their time fetching data from main memory. In an effort to improve the performance of these solvers, this work explores using data compression to reduce the amount of data that needs to be fetched from main memory. Some compression methods were found that improve the performance of the solver and problem found in the HPCG benchmark, with an increase in floating point operations per second …

Predicting Win Rates In Competitive Overwatchtm, Andrea Sibley

Mathematics Senior Capstone Papers

OverwatchTM is a video game published by Blizzard Entertainment R where two teams comprised of six people each compete against one another to accomplish a specific goal. The goal of each game is dependent on which map is being played. The maps are divided into four categories: Assault, Escort, Control, and Hybrid. A data set comprised of 3000 games of competitive OverwatchTM is used to determine how likely a team is to win their match. The factors used to determine the likelihood of winning are the map type and the skill ranking for each team. The data set is pre-processed …

Parameter Estimation And Optimal Design Techniques To Analyze A Mathematical Model In Wound Healing, Nigar Karimli

Masters Theses & Specialist Projects

For this project, we use a modified version of a previously developed mathematical model, which describes the relationships among matrix metalloproteinases (MMPs), their tissue inhibitors (TIMPs), and extracellular matrix (ECM). Our ultimate goal is to quantify and understand differences in parameter estimates between patients in order to predict future responses and individualize treatment for each patient. By analyzing parameter confidence intervals and confidence and prediction intervals for the state variables, we develop a parameter space reduction algorithm that results in better future response predictions for each individual patient. Moreover, use of another subset selection method, namely Structured Covariance Analysis, that …

Partially Isometric Matrices: A Brief And Selective Survey, Stephan Ramon Garcia, Matthew Okubo Patterson, William T. Ross

Department of Math & Statistics Faculty Publications

We survey a variety of results about partially isometric matrices. We focus primarily on results that are distinctly finite-dimensional. For example, we cover a recent solution to the similarity problem for partial isometries. We also discuss the unitary similarity problem and several other results.

Modeling The Distribution Of Lightning Strike Distances Outside A Preexisting Lightning Area, Dawn L. Sanderson

Theses and Dissertations

Air Force Instruction 91-203 (AFI 91-203) directs that a lightning warning be issued when lightning is occurring or imminent within a 5 nautical mile (NM) radius of a predetermined location or activity. The 45 Weather Squadron (WS), located on the central eastern coast of Florida, balances the safety of personnel and space launch vehicles with lost productivity of taking shelter from lightning. The primary objective of this study investigates if this 5 NM safety radius can be reduced while maintaining a desired level of safety. The research uses processed Lightning Detection and Ranging (LDAR) data to map the movement of …

Algorithms To Approximate Solutions Of Poisson's Equation In Three Dimensions, Ray Dambrose

Rose-Hulman Undergraduate Mathematics Journal

The focus of this research was to develop numerical algorithms to approximate solutions of Poisson's equation in three dimensional rectangular prism domains. Numerical analysis of partial differential equations is vital to understanding and modeling these complex problems. Poisson's equation can be approximated with a finite difference approximation. A system of equations can be formed that gives solutions at internal points of the domain. A computer program was developed to solve this system with inputs such as boundary conditions and a nonhomogenous source function. Approximate solutions are compared with exact solutions to prove their accuracy. The program is tested with an …

Numerical Investigation Of Inclination On The Thermal Performance Of Porous Fin Heatsink Using Pseudospectral Collocation Method, George Oguntala, Gbeminiyi Sobamowo, Raed Abd-Alhameed, James Noras

Karbala International Journal of Modern Science

Numerical investigation of inclination effect on the thermal performance of a porous fin heat sink is presented. The developed thermal model is solved using pseudo-spectral collocation method (PSCM). Parametric studies are carried out using PSCM, and the thermal characterization of heat sink with the inclined porous fin of rectangular geometry is presented. Results show that heat sink of inclined porous fin exhibits higher thermal performance than heat sink of vertical porous fin operating under the same thermal conditions with the same geometrical configurations. Performance of inclined or tilted fin increases with decrease in length-thickness aspect ratio. However, increase in the …

Studies Of Boundary Layer Transition From Laminar To Turbulent, Elizabeth Spaulding

Honors Program: Senior Projects (Public)

In engineering applications, there is a strong desire to reduce energy losses due to turbulent energy production. However, the theoretical understanding of turbulent and transition flow is still lacking due to the absence of an exact mathematical solution for turbulent flow. In the current project, transition and turbulent behavior in wall-bounded flow is studied, with an emphasis on concepts from boundary layer theory. Organized structures, such as streaks, waves, and vortices, which occur in such flows, are analyzed. Especially for near-wall structures, control strategies for suppressing such structures are discussed. Using mathematical equations which govern fluid flow, models for different …

Mhd Boundary Layer Slip Flow Over A Flat Plate With Soret And Dufour Effects, B. Shashidar Reddy, K. Saritha

Applications and Applied Mathematics: An International Journal (AAM)

The present paper studies the effects of Soret and Dufour on MHD boundary layer slip flow over a flat plate. The governing partial differential equations are converted to a set of nonlinear ordinary differential equations by using similarity transformations. Then, these equations are solved numerically by implicit Finite Difference Scheme. The numerical solutions for Velocity, Temperature and Concentration profiles for the related essential physical parameters are visualized through graphs and discussed. Results show that the velocity rises whereas the temperature and concentration reduces with the respective slip parameters. The increase in Soret number or decrease in Dufour number reduces the …

Spectral Discretization Errors In Filtered Subspace Iteration, Jay Gopalakrishnan, Luka Grubišić, Jeffrey S. Ovall

Mathematics and Statistics Faculty Publications and Presentations

We consider filtered subspace iteration for approximating a cluster of eigenvalues (and its associated eigenspace) of a (possibly unbounded) selfadjoint operator in a Hilbert space. The algorithm is motivated by a quadrature approximation of an operator-valued contour integral of the resolvent. Resolvents on infinite dimensional spaces are discretized in computable finite-dimensional spaces before the algorithm is applied. This study focuses on how such discretizations result in errors in the eigenspace approximations computed by the algorithm. The computed eigenspace is then used to obtain approximations of the eigenvalue cluster. Bounds for the Hausdorff distance between the computed and exact eigenvalue clusters …

Parametric Models For Signature Prediction And Feature Extraction, Julie Ann Jackson

Faculty Publications

This paper compares and contrasts numerical electromagnetic (EM) prediction methods and parametric scattering models for radar signature prediction and feature extraction.

Universal Quantum Computation, Junya Kasahara

Theses, Dissertations and Capstones

We study quantum computers and their impact on computability. First, we summarize the history of computer science. Only a few articles have determined the direction of computer science and industry despite the fact that many works have been dedicated to the present success. We choose articles by A. M. Turing and D. Deutsch, because A. M. Turing proposed the basic architecture of modern computers while D. Deutsch proposed an architecture for the next generation of computers called quantum computers. Second, we study the architecture of modern computers using Turing machines. The Turing machine has the basic design of modern computers …

Estimation Of Multivariate Asset Models With Jumps, Angela Loregian, Laura Ballotta, Gianluca Gianluca Fusai, Marcos Fabricio Perez

Business Faculty Publications

We propose a consistent and computationally efficient two-step methodology for the estimation of multidimensional non-Gaussian asset models built using Levy processes. The proposed framework allows for dependence between assets and different tail behaviors and jump structures for each asset. Our procedure can be applied to portfolios with a large number of assets as it is immune to estimation dimensionality problems. Simulations show good finite sample properties and significant efficiency gains. This method is especially relevant for risk management purposes such as, for example, the computation of portfolio Value at Risk and intra-horizon Value at Risk, as we show in detail …

Stair Climbing Hand Truck, James Mcpherson

All Undergraduate Projects

Abstract

Getting a heavy object up a flight of stairs usually requires a team of two or more people. Even with a team of people, the task is often still difficult, dangerous, and possibly insurmountable by one person. This problem is especially prevalent in for those who are moving into apartment complexes. Most apartment complexes have many buildings with two or more floors of living quarters, and elevators are often missing. This project sought to offer a solution to this problem. The solution in question; a motorized hand-truck with 2, trigonal planar pinwheels in place of the stock wheels. The …

Pawnee Dam Inflow Design Flood (Idf) Update And Stage-Frequency Curve Development Using Rmcrfa, Jennifer P. Christensen, Joshua J. Melliger

United States Geological Survey: Water Reports and Publications

Pawnee Dam is one of the ten Salt Creek Dams designed and built in the 1960s to mitigate flooding in Lincoln, Nebraska. This short paper illustrates the update of the Pawnee Dam inflow design flood (IDF) through calibration to recent high flow events and the development of its stage-frequency or hydrologic loading curve with the U.S. Army Corps of Engineers’ Risk Management Center Reservoir Frequency Analysis (RMC-RFA) model. The IDF update follows Engineering Regulation 1110-8-2, Inflow Design Flood for Dams and Reservoirs, including unit hydrograph peaking and two antecedent pool elevations. Background information on the original design of the dam …

On The Role Of Ill-Conditioning: Biharmonic Eigenvalue Problem And Multigrid Algorithms, Kasey Bray

Theses and Dissertations--Mathematics

Very fine discretizations of differential operators often lead to large, sparse matrices A, where the condition number of A is large. Such ill-conditioning has well known effects on both solving linear systems and eigenvalue computations, and, in general, computing solutions with relative accuracy independent of the condition number is highly desirable. This dissertation is divided into two parts.

In the first part, we discuss a method of preconditioning, developed by Ye, which allows solutions of Ax=b to be computed accurately. This, in turn, allows for accurate eigenvalue computations. We then use this method to develop discretizations that yield accurate computations …

Discontinuous Galerkin Methods For Convection-Diffusion Equations And Applications In Petroleum Engineering, Nattaporn Chuenjarern

Dissertations, Master's Theses and Master's Reports

This dissertation contains research in discontinuous Galerkin (DG) methods applying to convection-diffusion equations. It contains both theoretical analysis and applications. Initially, we develop a conservative local discontinuous Galerkin (LDG) method for the coupled system of compressible miscible displacement problem in two space dimensions. The main difficulty is how to deal with the discontinuity of approximations of velocity, u, in the convection term across the cell interfaces. To overcome the problems, we apply the idea of LDG with IMEX time marching using the diffusion term to control the convection term. Optimal error estimates in L^infinity(0, T; L^{2 …}

Approximation Of The Generalized Singular Value Expansion, Matthew Jacob Roberts

Dissertations, Master's Theses and Master's Reports

Let $X$, $Y$, and $Z$ be real separable Hilbert spaces, let $T:X \to Y$ be a compact operator, and let $L:D(L) \to Z$ be a closed and densely defined linear operator. Then the generalized singular value expansion (GSVE) is an expansion that expresses $T$ and $L$ in terms of a common orthonormal basis. Under certain hypotheses on discretization, the GSVE of an approximate operator pair $(T_j,L_j)$, where $T_j:X_j \to Y_j$ and $L_j:X_j \to Z_j$, converges to the GSVE of $(T,L)$. Error estimates establish a rate of convergence that is consistent with numerical experiments in the case of discretization using piecewise …

High Order Bound-Preserving Discontinuous Galerkin Methods And Their Applications In Petroleum Engineering, Ziyao Xu

Dissertations, Master's Theses and Master's Reports

This report contains researches in the theory of high-order bound-preserving (BP) discontinuous Galerkin (DG) method and their applications in petroleum engineering. It contains both theoretical analysis and numerical experiments. The compressible miscible displacements and wormhole propagation problem, arising in petroleum engineering, is used to describe the evolution of the pressure and concentrations of different components of fluid in porous media. The important physical features of concentration and porosity include their boundedness between 0 and 1, as well as the monotone increasing for porosity in wormhole propagation model. How to keep these properties in the simulation is crucial to the robustness …