Partial Differential Equations Commons^™

On Sobolev Spaces And The Existence Of Weak Solutions To Boundary Value Problems, Skye X. Paul

Master's Theses

Many boundary value problems that arise in mathematical models have close connections to second order elliptic partial differential equations. This thesis introduces the idea of weak derivatives and Sobolev Spaces to generalize possible solutions. Using functional analysis centered around the Lax-Milgram theorem, we show the existence of these generalized solutions to boundary value problems including Laplace's Equation, 2nd order linear ODEs, and ultimately a general second order elliptic PDE. The work cumulates with recovering a number of central theorems of functional analysis in the context of Sobolev Spaces, creating a new perspective on the solvability of these boundary value problems.

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 …

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.

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 …

A Time Integration Method Of Approximate Particular Solutions For Nonlinear Ordinary Differential Equations, Cyril Ocloo

Master's Theses

We consider a time-dependent method which is coupled with the method of approximate particular solutions (MAPS) of Delta-shaped basis functions and the method of fundamental solutions (MFS) to solve nonlinear ordinary differential equations. Firstly, we convert a nonlinear problem into a sequence of time-dependent non-homogeneous boundary value problems through a fictitious time integration method. The superposition principle is applied to split the numerical solution at each time step into an approximate particular solution and a homogeneous solution. Delta-shaped basis functions are used to provide an approximation of the source function at each time step. The purpose of this is to …

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 …

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.

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 …

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.

High Order Finite Elements For Lagrangian Computational Fluid Dynamics, Truman Everett Ellis

Master's Theses

A general finite element method is presented to solve the Euler equations in a Lagrangian reference frame. This FEM framework allows for separate arbitrarily high order representation of kinematic and thermodynamic variables. An accompanying hydrodynamics code written in Matlab is presented as a test-bed to experiment with various basis function choices. A wide range of basis function pairs are postulated and a few choices are developed further, including the bi-quadratic Q2-Q1d and Q2-Q2d elements. These are compared with a corresponding pair of low order bi-linear elements, traditional Q1-Q0 and sub-zonal pressure Q1-Q1d. Several test problems are considered including static convergence …