Partial Differential Equations Commons^™

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 …

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 …

Efficient Local Comparison Of Images Using Krawtchouk Descriptors, Julian Deville

Online Theses and Dissertations

It is known that image comparison can prove cumbersome in both computational complexity and runtime, due to factors such as the rotation, scaling, and translation of the object in question. Due to the locality of Krawtchouk polynomials, relatively few descriptors are necessary to describe a given image, and this can be achieved with minimal memory usage. Using this method, not only can images be described efficiently as a whole, but specific regions of images can be described as well without cropping. Due to this property, queries can be found within a single large image, or collection of large images, which …

Physical Applications Of The Geometric Minimum Action Method, George L. Poppe Jr.

Dissertations, Theses, and Capstone Projects

This thesis extends the landscape of rare events problems solved on stochastic systems by means of the \textit{geometric minimum action method} (gMAM). These include partial differential equations (PDEs) such as the real Ginzburg-Landau equation (RGLE), the linear Schroedinger equation, along with various forms of the nonlinear Schroedinger equation (NLSE) including an application towards an ultra-short pulse mode-locked laser system (MLL).

Additionally we develop analytical tools that can be used alongside numerics to validate those solutions. This includes the use of instanton methods in deriving state transitions for the linear Schroedinger equation and the cubic diffusive NLSE.

These analytical solutions are …

Introduction To The Usu Library Of Solutions To The Einstein Field Equations, Ian M. Anderson, Charles G. Torre

Tutorials on... in 1 hour or less

This is a Maple worksheet providing an introduction to the USU Library of Solutions to the Einstein Field Equations. The library is part of the DifferentialGeometry software project and is a collection of symbolic data and metadata describing solutions to the Einstein equations.

Examining The Electrical Excitation, Calcium Signaling, And Mechanical Contraction Cycle In A Heart Cell, Kristen Deetz, Nygel Foster, Darius Leftwich, Chad Meyer, Shalin Patel, Carlos Barajas, Matthias K. Gobbert, Zana Coulibaly

Spora: A Journal of Biomathematics

As the leading cause of death in the United States, heart disease has become a principal concern in modern society. Cardiac arrhythmias can be caused by a dysregulation of calcium dynamics in cardiomyocytes. Calcium dysregulation, however, is not yet fully understood and is not easily predicted; this provides motivation for the subsequent research. Excitation-contraction coupling (ECC) is the process through which cardiomyocytes undergo contraction from an action potential. Calcium induced calcium release (CICR) is the mechanism through which electrical excitation is coupled with mechanical contraction through calcium signaling. The study of the interplay between electrical excitation, calcium signaling, and mechanical …

High Performance Techniques Applied In Partial Differential Equations Library, Shilei Lin

All College Thesis Program, 2016-2019

This thesis explores various Trilinos packages to determine a method for updating the deal.ii library, which specializes in solving partial differential equations by finite element methods. It begins with introducing related concepts and general goals, followed by exploring computational and mathematical methods which are analytical solutions of one-dimensional Boussinesq equations and developing newer prototypes for solvers in deal.ii based on Trilinos packages. After demonstrating the methods, it indicates the reducing solving time in newer prototypes. Based on results from the prototype, similar methods are applied to update the deal.ii library. In the end, a testing program is exploited to demonstrate …

A Physics-Based Approach To Modeling Wildland Fire Spread Through Porous Fuel Beds, Tingting Tang

Theses and Dissertations--Mechanical and Aerospace Engineering

Wildfires are becoming increasingly erratic nowadays at least in part because of climate change. CFD (computational fluid dynamics)-based models with the potential of simulating extreme behaviors are gaining increasing attention as a means to predict such behavior in order to aid firefighting efforts. This dissertation describes a wildfire model based on the current understanding of wildfire physics. The model includes physics of turbulence, inhomogeneous porous fuel beds, heat release, ignition, and firebrands. A discrete dynamical system for flow in porous media is derived and incorporated into the subgrid-scale model for synthetic-velocity large-eddy simulation (LES), and a general porosity-permeability model is …

Computational Fluid Dynamics In A Terminal Alveolated Bronchiole Duct With Expanding Walls: Proof-Of-Concept In Openfoam, Jeremy Myers

Theses and Dissertations

Mathematical Biology has found recent success applying Computational Fluid Dynamics (CFD) to model airflow in the human lung. Detailed modeling of flow patterns in the alveoli, where the oxygen-carbon dioxide gas exchange occurs, has provided data that is useful in treating illnesses and designing drug-delivery systems. Unfortunately, many CFD software packages have high licensing fees that are out of reach for independent researchers. This thesis uses three open-source software packages, Gmsh, OpenFOAM, and ParaView, to design a mesh, create a simulation, and visualize the results of an idealized terminal alveolar sac model. This model successfully demonstrates that OpenFOAM can be …

An Algorithm For The Machine Calculation Of Minimal Paths, Robert Whitinger

Electronic Theses and Dissertations

Problems involving the minimization of functionals date back to antiquity. The mathematics of the calculus of variations has provided a framework for the analytical solution of a limited class of such problems. This paper describes a numerical approximation technique for obtaining machine solutions to minimal path problems. It is shown that this technique is applicable not only to the common case of finding geodesics on parameterized surfaces in R³, but also to the general case of finding minimal functionals on hypersurfaces in Rⁿ associated with an arbitrary metric.

Stereographic Visualization Of Bose-Einstein Condensate Clouds To Measure The Gravitational Constant, Ed Wesley Wells

College of Graduate Studies: Theses & Dissertations

This thesis describes a set of tools that can be used for the rapid design of atom interferometer schemes suitable for measuring Newton's Universal Gravitation constant also known as "Big G". This tool set is especially applicable to Bose--Einstein--condensed systems present in NASA's Cold Atom Laboratory experiment to be deployed to the International Space Station in 2017. These tools include a method of approximating the solutions of the nonlinear Schrödinger or Gross--Pitaevskii equation (GPE) using the Lagrangian Variational Method. They also include a set of software tools for translating the approximate solutions of the GPE into images of the optical …

Integrability And Exact Solutions For A (2+1)-Dimensional Variable-Coefficient Kdv Equation, Zhang Yu, Xu Gui-Qiong

Applications and Applied Mathematics: An International Journal (AAM)

By using the WTC method and symbolic computation, we apply the Painlevé test for a (2+1)-dimensional variable-coefficient Kortweg-de Vries (KdV) equation, and the considered equation is found to possess the Painlevé property without any parametric constraints. The auto-Bǎcklund transformation and several types of exact solutions are obtained by using the Painlevé truncated expansion method. Finally, the Hirota’s bilinear form is presented and multi-soliton solutions are also constructed.

Stability Aware Delaunay Refinement, Bishal Acharya

UNLV Theses, Dissertations, Professional Papers, and Capstones

Good quality meshes are extensively used for finding approximate solutions for partial differential equations for fluid flow in two dimensional surfaces. We present an overview of existing algorithms for refinement and generation of triangular meshes. We introduce the concept of node stability in the refinement of Delaunay triangulation. We present two algorithms for generating stable refinement of Delaunay triangulation. We also present an experimental investigation of a triangulation refinement algorithm based on the location of the center of gravity and the location of the center of circumcircle. The results show that the center of gravity based refinement is more effective …

Numerical Studies For Solving Fractional Riccati Differential Equation, N. H. Sweilam, M. M. Khader, A. M. S. Mahdy

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, finite difference method (FDM) and Pade'-variational iteration method (Pade'- VIM) are successfully implemented for solving the nonlinear fractional Riccati differential equation. The fractional derivative is described in the Caputo sense. The existence and the uniqueness of the proposed problem are given. The resulting nonlinear system of algebraic equations from FDM is solved by using Newton iteration method; moreover the condition of convergence is verified. The convergence's domain of the solution is improved and enlarged by Pade'-VIM technique. The results obtained by using FDM is compared with Pade'-VIM. It should be noted that the Pade'-VIM is preferable because …

Introducing An Efficient Modification Of The Variational Iteration Method By Using Chebyshev Polynomials, M. M. Khader

Applications and Applied Mathematics: An International Journal (AAM)

In this article an efficient modification of the variational iteration method (VIM) is presented using Chebyshev polynomials. Special attention is given to study the convergence of the proposed method. The new modification is tested for some examples to demonstrate reliability and efficiency of the proposed method. A comparison of our numerical results those of the conventional numerical method, the fourth-order Runge-Kutta method (RK4) are given. The comparison shows that the solution using our modification is fast-convergent and is in excellent conformance with the exact solution. Finally, we conclude that the proposed method can be applied to a large class of …

Parallel And Distributed Simulation Of Parabolic And Telegraphic Equations., Ewedafe Simon Uzezi

Student Works (2010-2019)

In this thesis, a parallel implementation of explicit/implicit parallel algorithms such as the stationary iterative methods and the class of iterating alternating methods which includes: Alternating Direction Implicit (ADI), Iterative Alternating Direction Explicit (IADE), for D’Yakonov (IADE-DY), Double sweep Mitchell and Fairweather (MF-DS) and Alternating Group Explicit (AGE) method for solving 1-Dimensional (1-D), 2-Dimensional (2-D) Parabolic (special examples including 1-D, 2-D Bio-Heat Equation) and 1-D, 2-D and 3-D Telegraphic Equations on a distributed environment of Message Passing Interface (MPI) and Parallel Virtual Machine (PVM) platform is presented. To correlate the communication activity with computation, we counted events between significant MPI/PVM …

Parallel Newton-Krylov-Schwarz Solvers For The Full Potential Flow Equation, Jie Zhang

Computer Science Theses & Dissertations

Newton-Krylov-Schwarz methods are increasingly applied in Computational Fluid Dynamics (CFD). We develop a parallel analysis code based on this method for the full potential flow model. The full potential model consists of a single nonlinear second-order partial differential equation of mixed type (elliptic/hyperbolic), which we solve as a steady boundary-value problem.

We use a nine-point finite-difference stencil to discretize the equation. A Newtonlike linearization and correction method is used to solve the resulting set of nonlinear algebraic equations. To solve the inner linear equations, we employ a Krylov space method. Preconditioners are used to improve the convergence rate. In order …