Partial Differential Equations Commons^™

Delta-Shaped Approximation Based Homotopy Analysis Method For Nonlinear Poisson-Type Partial Differential Equations, Cyril Ocloo

Dissertations

This research aims to solve nonlinear Poisson-type partial differential equations (PDEs) by the approach of the homotopy analysis method (HAM) incorporated with approximate particular solutions (APS) using Delta-shaped basis (DSB) approximations.

With the inclusion of the h auxiliary parameters, we tackle nonlinear problems by studying the mathematical characteristics of the h curve. This is to ensure the numerical convergence of the HAM.

In the solution process, we use the homotopy analysis method to convert a nonlinear PDE into linear inhomogeneous PDEs, which are solved using the method of approximate particular solutions with DSB.

A proper value of the h is …

Scalable Solution Of Time-Dependent Pdes Through Component-Wise Exponential Integrator, Chelsea Drum

Dissertations

Exponential integrators, such as exponential Runge-Kutta or Rosenbrock methods, are designed specifically for the time integration of stiff systems of ordinary differential equations (ODEs) and allow the use of larger time steps than other general-purpose ODE solvers. However, these methods rely on computing matrix function-vector products that are traditionally computed using a Krylov projection, such as Lanczos or Arnoldi iteration, that involves substantial computational expense at high spatial resolution. Krylov Subspace Spectral (KSS) methods' frequency-dependent approach, designed to circumvent stiffness in linear problems, computes these products with greater scalability. We propose the combination of such KSS methods with exponential integrators …

A Tale Of Two Diagonalizations: Methods To Diagonalize A 1-D Piecewise Constant Indefinite Schrödinger Operator, Sarah Wright

Dissertations

We present two numerical methods for computing the solution of a partial differential equation (PDE) for modeling acoustic pressure, known as an extra-wide angle parabolic equation, that features the square root of a differential operator. The differential operator is the negative of an indefinite Schrödinger operator with a piecewise constant potential. This work primarily deals with the 3-piece case; however, a generalization is made to the case of an arbitrary number of pieces. In the first method, the Rayleigh-Secant Method, through restriction to a judiciously chosen lower-dimensional subspace, approximate eigenfunctions are used to obtain estimates for the eigenvalues of the …

Stochastic Modeling Of Flows In Membrane Pore Networks, Binan Gu

Dissertations

Membrane filters provide immediate solutions to many urgent problems such as water purification, and effective remedies to pressing environmental concerns such as waste and air treatment. The ubiquity of applications gives rise to a significant amount of research in membrane material selection and structural design to optimize filter efficiency. As physical experiments tend to be costly, numerical simulation and analysis of fluid flow, foulant transport and geometric evolution due to foulant deposition in complex geometries become particularly relevant. In this dissertation, several mathematical modeling and analytical aspects of the industrial membrane filtration process are investigated. A first-principles mathematical model for …

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.

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. …

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 …

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 …

Solution Of Pdes For First-Order Photobleaching Kinetics Using Krylov Subspace Spectral Methods, Somayyeh Sheikholeslami

Dissertations

We solve the first order reaction-diffusion equations which describe binding-diffusion kinetics using a photobleaching scanning profile of a confocal laser scanning microscope approximated by a Gaussian laser profile. We show how to solve these equations with prebleach steady-state initial conditions using a time-domain method known as a Krylov Subspace Spectral (KSS) method. KSS methods are explicit methods for solving time- dependent variable-coefficient partial differential equations (PDEs). KSS methods are advantageous compared to other methods because of their stability and their superior scalability. These advantages are obtained by applying Gaussian quadrature rules in the spectral domain developed by Golub and Meurant. …

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 …

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 …

Time Integration Methods Of Fundamental Solutions And Approximate Fundamental Solutions For Nonlinear Elliptic Partial Differential Equations, Corey Leon Jones

Dissertations

A time-dependent method is coupled with the Method of Approximate Particular Solutions (MAPS) of Delta-shaped basis functions, the Method of Fundamental Solutions (MFS), and the Method of Approximate Fundamental Solutions (MAFS) to solve a second order nonlinear elliptic partial differential equation (PDE) on regular and irregular shaped domains. The nonlinear PDE boundary value problem is first transformed into a time-dependent quasilinear problem by introducing a fictitious time. Forward Euler integration is then used to ultimately convert the problem into a sequence of time-dependent linear nonhomogeneous modified Helmholtz boundary value problems on which the superposition principle is applied to split the …