Partial Differential Equations Commons^™

(R1491) Numerical Solution Of The Time-Space Fractional Diffusion Equation With Caputo Derivative In Time By A-Polynomial Method, Saeid Abbasbandy, Jalal Hajishafieiha

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, a novel type of polynomial is defined which is equipped with an auxiliary parameter a. These polynomials are a combination of the Chebyshev polynomials of the second kind. The approximate solution of each equation is assumed as the sum of these polynomials and then, with the help of the collocation points, the unknown coefficients of each polynomial, as well as auxiliary parameter, is obtained optimally. Now, by placing the optimal value of a in polynomials, the polynomials are obtained without auxiliary parameter, which is the restarted step of the present method. The time discretization is performed …

(R1480) Heat Transfer In Peristaltic Motion Of Rabinowitsch Fluid In A Channel With Permeable Wall, Mahadev M. Channakote, Dilipkumar V. Kalse

Applications and Applied Mathematics: An International Journal (AAM)

This paper is intended to investigate the effect of heat transfer on the peristaltic flow of Rabinowitsch fluid in a channel lined with a porous material. The Navier -Stokes equation governs the channel's flow, and Darcy's law describes the permeable boundary. The Rabinowitsch fluid model's governing equations are solved by utilizing approximations of the long-wavelength and small number of Reynolds. The expressions for axial velocity, temperature distribution, pressure gradient, friction force, stream function are obtained. The influence on velocity, pressure gradient, friction force, and temperature on pumping action of different physical parameters is explored via graphs.

(R1494) Approximate Solutions Of The Telegraph Equation, Ilija Jegdić

Applications and Applied Mathematics: An International Journal (AAM)

In this paper the initial boundary value problems for the linear telegraph equation in one and two space dimensions are considered. To find approximate solutions, a recently proposed optimization-free approach that utilizes artificial neural networks with one hidden layer is used, in which the connecting weights from the input layer to the hidden layer are chosen randomly and the weights from the hidden layer to the output layer are found by solving a system of linear equations. One of the advantages of this method, in comparison to the usual discretization methods for the two-dimensional linear telegraph equation, is that this …

(R1496) Impact Of Electronic States Of Conical Shape Of Indium Arsenide/Gallium Arsenide Semiconductor Quantum Dots, Md. Fayz-Al-Asad, Md. Al-Rumman, Md. Nur Alam, Salma Parvin, Cemil Tunç

Applications and Applied Mathematics: An International Journal (AAM)

Semiconductor quantum dots (QDs) have unique atom-like properties. In this work, the electronic states of quantum dot grown on a GaAs substrate has been studied. The analytical expressions of electron wave function for cone-like quantum dot on the semiconductor surface has been obtained and the governing eigen value equation has been solved, thereby obtaining the dependence of ground state energy on radius and height of the cone-shaped -dots. In addition, the energy of eigenvalues is computed for various length and thickness of the wetting layer (WL). We discovered that the eigen functions and energies are nearly associated with the GaAs …

Survivals Of Two Cooperating Species Of Animals, Joon Hyuk Kang

Faculty Publications

The purpose of this paper is to give conditions for the existence and uniqueness of positive solution to a rather general type of elliptic system of the Dirichlet problem on a bounded domain Ω" role="presentation" style="box-sizing: border-box; margin: 0px; padding: 0px; display: inline-block; line-height: normal; font-size: 16.2px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; border: 0px; position: relative;">Ω in Rn" role="presentation" style="box-sizing: border-box; margin: 0px; padding: 0px; display: inline-block; line-height: normal; font-size: 16.2px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: …

Ecological Dynamics On Large Metapopulation Graphs, Daniel Cooney

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Mathematical Modelling And Simulation Using An Efficient Pinns Algorithm To Understand Spread Of Infection In Enclosed Spaces, Long Nguyen, Arkaprovo Ghosal, Rudra Nagalia, Padmanabhan Seshaiyer

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

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 …

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.

A Generalized Polar-Coordinate Integration Formula, Oscillatory Integral Techniques, And Applications To Convolution Powers Of Complex-Valued Functions On $\Mathbb{Z}^D$, Huan Q. Bui

Honors Theses

In this thesis, we consider a class of function on $\mathbb{R}^d$, called positive homogeneous functions, which interact well with certain continuous one-parameter groups of (generally anisotropic) dilations. Generalizing the Euclidean norm, positive homogeneous functions appear naturally in the study of convolution powers of complex-valued functions on $\mathbb{Z}^d$. As the spherical measure is a Radon measure on the unit sphere which is invariant under the symmetry group of the Euclidean norm, to each positive homogeneous function $P$, we construct a Radon measure $\sigma_P$ on $S=\{\eta \in \mathbb{R}^d:P(\eta)=1\}$ which is invariant under the symmetry group of $P$. With this measure, we prove …

Five-Wave Resonances In Deep Water Gravity Waves: Integrability, Numerical Simulations And Experiments, Dan Lucas, Marc Perlin, Dian-Yong Liu, Shane Walsh, Rossen Ivanov, Miguel D. Bustamante

Articles

In this work we consider the problem of finding the simplest arrangement of resonant deep water gravity waves in one-dimensional propagation, from three perspectives: Theoretical, numerical and experimental. Theoretically this requires using a normal-form Hamiltonian that focuses on 5-wave resonances. The simplest arrangement is based on a triad of wave vectors K₁ + K₂ = K₃ (satisfying specific ratios) along with their negatives, corresponding to a scenario of encountering wave packets, amenable to experiments and numerical simulations. The normal-form equations for these encountering waves in resonance are shown to be non-integrable, but they admit an integrable reduction …