Partial Differential Equations Commons^™

Finite Element Approximation Of Solutions Of The Equations Of Electroporoelasticity, Yu Hu

Mathematics Theses and Dissertations

In this thesis we consider the solution of the equations of electroporoelasticity, which are a combination of Maxwell's equations and the poroelasticity equations. Included is a description of suitable initial and boundary conditions, weak formulation of the equations, and the error estimate for a general numerical method.

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.

Smooth Global Approximation For Continuous Data Assimilation, Kenneth R. Brown

Theses and Dissertations

This thesis develops the finite element method, constructs local approximation operators, and bounds their error. Global approximation operators are then constructed with a partition of unity. Finally, an application of these operators to data assimilation of the two-dimensional Navier-Stokes equations is presented, showing convergence of an algorithm in all Sobolev topologies.

The Exact Factorization Equations For One- And Two-Level Systems, Bart Rosenzweig

Theses and Dissertations

Exact Factorization is a framework for studying quantum many-body problems. This decomposes the wavefunctions of such systems into conditional and marginal components. We derive corresponding evolution equations for molecular systems whose conditional electronic subsystems are described by one or two Born-Oppenheimer levels and develop a program for their mathematical study.

Qualitative Analysis Of A Modified Leslie-Gower Predator-Prey Model With Weak Allee Effect Ii, Manoj K. Singh, B. S. Bhadauria

Applications and Applied Mathematics: An International Journal (AAM)

The article aims to study a modified Leslie-Gower predator-prey model with Allee effect II, affecting the functional response with the assumption that the extent to which the environment provides protection to both predator and prey is the same. The model has been studied analytically as well as numerically, including stability and bifurcation analysis. Compared with the predator-prey model without Allee effect, it is found that the weak Allee effect II can bring rich and complicated dynamics, such as the model undergoes to a series of bifurcations (Homoclinic, Hopf, Saddle-node and Bogdanov-Takens). The existence of Hopf bifurcation has been shown for …

A High Order Finite Difference Method To Solve The Steady State Navier-Stokes Equations, Nihal J. Siriwardana, Saroj P. Pradhan

Applications and Applied Mathematics: An International Journal (AAM)

In this article, we develop a fourth order finite difference method to solve the system of steady state Navier-Stokes equations and apply it to the benchmark problem known as the square cavity flow problem. The numerical results of 𝑢-velocity components and 𝑣-velocity components obtained at the center of the cavity are compared with the results obtained by the method developed by Greenspan and Casulli to solve the time dependent system of Navier-Stokes equations. The method described in this article is easy to implement and it has been shown to be more efficient and stable than the method by Greenspan and …

Hybrid Algorithm For Singularly Perturbed Delay Parabolic Partial Differential Equations, Imiru T. Daba, Gemechis F. Duressa

Applications and Applied Mathematics: An International Journal (AAM)

This study aims at constructing a numerical scheme for solving singularly perturbed parabolic delay differential equations. Taylor’s series expansion is applied to approximate the shift term. The obtained result is approximated by using the implicit Euler method in the temporal discretization on a uniform step size with the hybrid numerical scheme consisting of the midpoint upwind method in the outer layer region and the cubic spline method in the inner layer region on a piecewise uniform Shishkin mesh in the spatial discretization. The constructed scheme is an ε−uniformly convergent accuracy of order one. Some test examples are considered to testify …

An Examination Of Fontan Circulation Using Differential Equation Models And Numerical Methods, Vanessa Maybruck

Honors Student Research

Certain congenital heart defects can lead to the development of only a single pumping chamber, or ventricle, in the heart instead of the usual two ventricles. Individuals with this defect undergo a corrective, three-part surgery, the third step of which is the Fontan procedure, but as the patients age, their cardiovascular health will likely deteriorate. Using computational fluid dynamics and differential equations, Fontan circulation can be modeled to investigate why the procedure fails and how Fontan failure can be maximally prevented. Borrowing from well-established literature on RC circuits, the differential equation models simulate systemic blood flow in a piecewise, switch-like …

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.

The Effect Of Initial Conditions On The Weather Research And Forecasting Model, Aaron D. Baker

Electronic Theses and Dissertations

Modeling our atmosphere and determining forecasts using numerical methods has been a challenge since the early 20th Century. Most models use a complex dynamical system of equations that prove difficult to solve by hand as they are chaotic by nature. When computer systems became more widely adopted and available, approximating the solution of these equations, numerically, became easier as computational power increased. This advancement in computing has caused numerous weather models to be created and implemented across the world. However a challenge of approximating these solutions accurately still exists as each model have varying set of equations and variables to …

Lecture 08: Partial Eigen Decomposition Of Large Symmetric Matrices Via Thick-Restart Lanczos With Explicit External Deflation And Its Communication-Avoiding Variant, Zhaojun Bai

Mathematical Sciences Spring Lecture Series

There are continual and compelling needs for computing many eigenpairs of very large Hermitian matrix in physical simulations and data analysis. Though the Lanczos method is effective for computing a few eigenvalues, it can be expensive for computing a large number of eigenvalues. To improve the performance of the Lanczos method, in this talk, we will present a combination of explicit external deflation (EED) with an s-step variant of thick-restart Lanczos (s-step TRLan). The s-step Lanczos method can achieve an order of s reduction in data movement while the EED enables to compute eigenpairs in batches along with a number …

Lecture 07: Nonlinear Preconditioning Methods And Applications, Xiao-Chuan Cai

Mathematical Sciences Spring Lecture Series

We consider solving system of nonlinear algebraic equations arising from the discretization of partial differential equations. Inexact Newton is a popular technique for such problems. When the nonlinearities in the system are well-balanced, Newton's method works well, but when a small number of nonlinear functions in the system are much more nonlinear than the others, Newton may converge slowly or even stagnate. In such a situation, we introduce some nonlinear preconditioners to balance the nonlinearities in the system. The preconditioners are often constructed using a combination of some domain decomposition methods and nonlinear elimination methods. For the nonlinearly preconditioned problem, …

Lecture 10: Preconditioned Iterative Methods For Linear Systems, Edmond Chow

Mathematical Sciences Spring Lecture Series

Iterative methods for the solution of linear systems of equations – such as stationary, semi-iterative, and Krylov subspace methods – are classical methods taught in numerical analysis courses, but adapting these methods to run efficiently at large-scale on high-performance computers is challenging and a constantly evolving topic. Preconditioners – necessary to aid the convergence of iterative methods – come in many forms, from algebraic to physics-based, are regularly being developed for linear systems from different classes of problems, and similarly are evolving with high-performance computers. This lecture will cover the background and some recent developments on iterative methods and preconditioning …

A Direct Method For Modeling And Simulations Of Elliptic And Parabolic Interface Problems, Kumudu Gamage, Yan Peng

College of Sciences Posters

Interface problems have many applications in fluid dynamics, molecular biology, electromagnetism, material science, heat distribution in engines, and hyperthermia treatment of cancer. Mathematically, interface problems commonly lead to partial differential equations (PDE) whose in- put data are discontinuous or singular across the interfaces in the solution domain. Many standard numerical methods designed for smooth solutions poorly work for interface problems as solutions of the interface problems are mostly non-smoothness or discontinuous. Moving interface problems depends on the accuracy of the gradient of the solution at the interface. Therefore, it became essential to derive a method for interface problems that gives …

Analysis Of Boundary Observability Of Strongly Coupled One-Dimensional Wave Equations With Mixed Boundary Conditions, Wilson Dennis Horner

Masters Theses & Specialist Projects

*see note below

In control theory, the time it takes to receive a signal after it is sent is referred to as the observation time. For certain types of materials, the observation time to receive a wave signal differs depending on a variety of factors, such as material density, flexibility, speed of the wave propagation, etc. Suppose we have a strongly coupled system of two wave equations describing the longitudinal vibrations on a piezoelectric beam of length L. These two wave equations have non-identical wave propagation speeds c1 and c2. First, we prove the exact observability inequality with the optimal …

Discontinuous Galerkin Method Applied To Navier-Stokes Equations, Ian Deruiter, Mahboub Baccouch

UNO Student Research and Creative Activity Fair

Discontinuous Galerkin (DG) finite element methods are becoming important techniques for the computational solution of many real-world problems describe by differential equations. They combine many attractive features of the finite element and the finite volume methods. These methods have been successfully applied to many important PDEs arising from a wide range of applications. DG methods are highly accurate numerical methods and have considerable advantages over the classical numerical methods available in the literature. DG methods can easily handle meshes with hanging nodes, elements of various types and shapes, and local spaces of different orders. Furthermore, DG methods provide accurate and …

Stochastic Navier-Stokes Equations With Markov Switching, Po-Han Hsu

LSU Doctoral Dissertations

This dissertation is devoted to the study of three-dimensional (regularized) stochastic Navier-Stokes equations with Markov switching. A Markov chain is introduced into the noise term to capture the transitions from laminar to turbulent flow, and vice versa. The existence of the weak solution (in the sense of stochastic analysis) is shown by studying the martingale problem posed by it. This together with the pathwise uniqueness yields existence of the unique strong solution (in the sense of stochastic analysis). The existence and uniqueness of a stationary measure is established when the noise terms are additive and autonomous. Certain exit time estimates …

The Pencil Code, A Modular Mpi Code For Partial Differential Equations And Particles: Multipurpose And Multiuser-Maintained, The Pencil Code Collaboration, Chao-Chin Yang

Physics & Astronomy Faculty Research

The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) (A. Brandenburg & Dobler, 2010) to meteorological cloud microphysics (Li et al., 2017) and engineering applications in combustion (Babkovskaia et al., 2011). Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to …

Principles For Determining The Motion Of Blood Through Arteries, Sylvio R. Bistafa

Euleriana

Translation of Principia pro motu sanguinis per arterias determinando (E855). This work of 1775 by L. Euler is considered to be the first mathematical treatment of circulatory physiology and hemodynamics.

Multicomponent Fokas-Lenells Equations On Hermitian Symmetric Spaces, Vladimir Gerdjikov, Rossen Ivanov

Articles

Multi-component integrable generalizations of the Fokas-Lenells equation, associated with each irreducible Hermitian symmetric space are formulated. Description of the underlying structures associated to the integrability, such as the Lax representation and the bi-Hamiltonian formulation of the equations is provided. Two reductions are considered as well, one of which leads to a nonlocal integrable model. Examples with Hermitian symmetric spaces of all classical series of types A.III, BD.I, C.I and D.III are presented in details, as well as possibilities for further reductions in a general form.