Numerical Analysis and Computation 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 …

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 …

A Component-Wise Approach To Smooth Extension Embedding Methods, Vivian Montiforte

Dissertations

Krylov Subspace Spectral (KSS) Methods have demonstrated to be highly scalable methods for PDEs. However, a current limitation of these methods is the requirement of a rectangular or box-shaped domain. Smooth Extension Embedding Methods (SEEM) use fictitious domain methods to extend a general domain to a simple, rectangular or box-shaped domain. This dissertation describes how these methods can be combined to extend the applicability of KSS methods, while also providing a component-wise approach for solving the systems of equations produced with SEEM.

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 …

Variable Compact Multi-Point Upscaling Schemes For Anisotropic Diffusion Problems In Three-Dimensions, James Quinlan

Dissertations

Simulation is a useful tool to mitigate risk and uncertainty in subsurface flow models that contain geometrically complex features and in which the permeability field is highly heterogeneous. However, due to the level of detail in the underlying geocellular description, an upscaling procedure is needed to generate a coarsened model that is computationally feasible to perform simulations. These procedures require additional attention when coefficients in the system exhibit full-tensor anisotropy due to heterogeneity or not aligned with the computational grid. In this thesis, we generalize a multi-point finite volume scheme in several ways and benchmark it against the industry-standard routines. …

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 …

Enhancement Of Krylov Subspace Spectral Methods Through The Use Of The Residual, Haley Dozier

Dissertations

Depending on the type of equation, finding the solution of a time-dependent partial differential equation can be quite challenging. Although modern time-stepping methods for solving these equations have become more accurate for a small number of grid points, in a lot of cases the scalability of those methods leaves much to be desired. That is, unless the timestep is chosen to be sufficiently small, the computed solutions might exhibit unreasonable behavior with large input sizes. Therefore, to improve accuracy as the number of grid points increases, the time-steps must be chosen to be even smaller to reach a reasonable solution. …

Adaptive Meshfree Methods For Partial Differential Equations, Jaeyoun Oh

Dissertations

There are many types of adaptive methods that have been developed with different algorithm schemes and definitions for solving Partial Differential Equations (PDE). Adaptive methods have been developed in mesh-based methods, and in recent years, they have been extended by using meshfree methods, such as the Radial Basis Function (RBF) collocation method and the Method of Fundamental Solutions (MFS). The purpose of this dissertation is to introduce an adaptive algorithm with a residual type of error estimator which has not been found in the literature for the adaptive MFS. Some modifications have been made in developing the algorithm schemes depending …

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.

Radial Basis Function Differential Quadrature Method For The Numerical Solution Of Partial Differential Equations, Daniel Watson

Dissertations

In the numerical solution of partial differential equations (PDEs), there is a need for solving large scale problems. The Radial Basis Function Differential Quadrature (RBFDQ) method and local RBF-DQ method are applied for the solutions of boundary value problems in annular domains governed by the Poisson equation, inhomogeneous biharmonic equation, and the inhomogeneous Cauchy-Navier equations of elasticity. By choosing the collocation points properly, linear systems can be obtained so that the coefficient matrices have block circulant structures. The resulting systems can be efficiently solved using matrix decomposition algorithms (MDAs) and fast Fourier transforms (FFTs). For the local RBFDQ method, the …

Numerical Solution Of Partial Differential Equations Using Polynomial Particular Solutions, Thir R. Dangal

Dissertations

Polynomial particular solutions have been obtained for certain types of partial differential operators without convection terms. In this dissertation, a closed-form particular solution for more general partial differential operators with constant coefficients has been derived for polynomial basis functions. The newly derived particular solutions are further coupled with the method of particular solutions (MPS) for numerically solving a large class of elliptic partial differential equations. In contrast to the use of Chebyshev polynomial basis functions, the proposed approach is more flexible in selecting the collocation points inside the domain. Polynomial basis functions are well-known for yielding ill-conditioned systems when their …

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 …

Efficient Denoising And Sharpening Of Color Images Through Numerical Solution Of Nonlinear Diffusion Equations, Linh T. Duong

Honors Theses

The purpose of this project is to enhance color images through denoising and sharpening, two important branches of image processing, by mathematically modeling the images. Modifications are made to two existing nonlinear diffusion image processing models to adapt them to color images. This is done by treating the red, green, and blue (RGB) channels of color images independently, contrary to the conventional idea that the channels should not be treated independently. A new numerical method is needed to solve our models for high resolution images since current methods are impractical. To produce an efficient method, the solution is represented as …

Fast Method Of Particular Solutions For Solving Partial Differential Equations, Anup Raja Lamichhane

Dissertations

Method of particular solutions (MPS) has been implemented in many science and engineering problems but obtaining the closed-form particular solutions, the selection of the good shape parameter for various radial basis functions (RBFs) and simulation of the large-scale problems are some of the challenges which need to overcome. In this dissertation, we have used several techniques to overcome such challenges.

The closed-form particular solutions for the Matérn and Gaussian RBFs were not known yet. With the help of the symbolic computational tools, we have derived the closed-form particular solutions of the Matérn and Gaussian RBFs for the Laplace and biharmonic …

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.

On The Selection Of A Good Shape Parameter For Rbf Approximation And Its Application For Solving Pdes, Lei-Hsin Kuo

Dissertations

Meshless methods utilizing Radial Basis Functions~(RBFs) are a numerical method that require no mesh connections within the computational domain. They are useful for solving numerous real-world engineering problems. Over the past decades, after the 1970s, several RBFs have been developed and successfully applied to recover unknown functions and to solve Partial Differential Equations (PDEs).
However, some RBFs, such as Multiquadratic (MQ), Gaussian (GA), and Matern functions, contain a free variable, the shape parameter, c. Because c exerts a strong influence on the accuracy of numerical solutions, much effort has been devoted to developing methods for determining shape parameters which provide …