Partial Differential Equations Commons^™

Physics-Informed Neural Network Solution Of The 2d Helmholtz Equation With A Gaussian Source, Theodoros Panagiotakopoulos, Chris Velissaris, Aristotelis Nikolaos Rapsomanikis

Faculty Scholarship and Creative Works

We present a physics-informed neural network (PINN) framework for solving the complex-valued two-dimensional Helmholtz equation with a localized Gaussian source and spatially varying permittivity. Starting from Maxwell’s equations, the frequency-domain scalar Helmholtz formulation under transverse electric (TE) polarization is derived and enforced directly within the neural network loss function. The model employs a sinusoidal representation network (SIREN) architecture to capture the oscillatory nature of wave solutions and incorporates the Sommerfeld radiation condition to impose open boundary conditions. Training is performed using a hybrid collocation strategy combined with a two-stage optimization procedure consisting of Adam followed by L-BFGS. Numerical experiments in …

The Pure Yang-Mills Field: I: Fixed Time Existence, James Glimm, James Glimm

Department of Applied Mathematics & Statistics Faculty Publications

The convergence of renormalized perturbation theory to all finite orders is defined and shown to be valid for the fixed time perturbation theory of a pure Yang-Mills field.

Two pure Yang-Mills quantum gauge field theories are constructed, one based on short distance asymptotics and the other based on long distance asymptotics.

The construction depends on an assumed principle of a maximum rate of entropy production.

The Boson Yang-Mills Field: The Loop Expansion, James Glimm

Department of Applied Mathematics & Statistics Faculty Publications

This paper demonstrates convergence of the loop expansion for Yang-Mills fields.

The loop construction of perturbation theory is based on the axial gauge, ghost states, the BRST framework and the Gribov extension of the Hamiltonian, with a loop expansion cutoff to all finite orders for the dynamics.

The construction is established by renormalized perturbation theory convergent to all finite orders.

Two distinct Yang-Mills theories are constructed, one based on short distance asymptotics and the other based on long distance asymptotics.

The construction depends on an assumed principle of a maximum rate of entropy production.

The paper has sufficient generality to …

Advanced Soliton Structures And Elliptic Wave Patterns In A Sixth-Order Nonlinear Schrödinger Equation Using Improved Modified Extended Tanh Function Method, Mina Fahim, Hamdy Mohamed Ahmed, Islam Samir, Kamal Dib

Basic Science Engineering

In this work, a sixth–order extension of the nonlinear Schrödinger equation (NLSE) within its integrable hierarchy is investigated to model higher–order nonlinear and dispersive effects relevant to optical fiber systems and nonlinear wave propagation. By employing the Improved Modified Extended Tanh Function Method, a comprehensive family of exact analytical solutions is derived, encompassing bright and dark solitons, singular soliton structures, and singular periodic solutions. In addition, solution families expressed in terms of Jacobi elliptic functions, Weierstrass doubly periodic elliptic functions, and exponential profiles are obtained. The novelty of this study lies in extending the analytical framework of the NLSE hierarchy …

Constructing The Soliton Wave Structure And Stability Analysis To Generalized Calogero–Bogoyavlenskii–Schiff Equation Using Improved Simple Equation Method, Mina Fahim, Hamdy Mohamed Ahmed, Islam Samir Soliman, Mohamed Elsaid, Kamal Hassan Eldib

Basic Science Engineering

In this work, we investigated the (3+1)-dimensional generalized Calogero–Bogoyavlenskii–Schiff equation, which models long wave propagation in shallow water and plays a significant role in fluid mechanics and plasma physics. Using the improved simple equations method, we obtained various solutions, including dark, bright, and singular solitons, and combinations of singular periodic solutions and exponential rational solutions. Additionally, we performed a linear stability analysis to examine the stability properties of these wave solutions. To further illustrate their characteristics during propagation, we provided 3D and contour plots for some opted wave solutions.

Mathematical Modelling Of Hybrid Photonic Structures For Holographic Sensors, Jack Lyons

Doctoral

This thesis outlines a mathematical framework for modelling the formation of holographic gratings in hybrid photopolymer based nanocomposites with the aim of optimising their holographic recording properties for optical sensing applications. Thus, the second aim of the work is to model the change in optical properties of the grating in response to exposure to a target analyte. This work has been a collaborative research project between the School of Mathematics & Statistics at Technological University Dublin and the Centre for Industrial and Engineering Optics that have done extensive experimental work with holographic gratings recorded in photopolymer materials.

In recent years, …

Ergodic Switching Control For Markov-Feller Processes Ii, Jose L. Menaldi, Maurice Robin

Mathematics Faculty Research Publications

This is the continuation of Part I [14], where we considered control problems with long term average (or ergodic) cost for Markov switching processes (zt , nt ), nt being a discrete component with values in a finite set N . The control acts only on this discrete component and consists of immediate switching actions. We solve the ergodic problem in several situations extending previous works, mainly when zt is a reflected diffusion with or without jumps and when the set of control values is strictly smaller than N . In this Part II, we conclude our theoretical analysis with …

Ergodic Switching Control For Markov-Feller Processes I, Jose L. Menaldi, Maurice Robin

Mathematics Faculty Research Publications

We consider control problems with long term average (or ergodic) cost for Markov switching processes (zt, nt ), nt being a discrete component with values in a finite set N . The control acts only on this discrete component and consists of immediate switching actions. We solve the ergodic problem in several situations extending previous works, mainly when zt is a reflected diffusion with or without jumps and when the set of control values is strictly smaller than N .

Modelling The Formation Of Unslanted Holographic Gratings In Hybrid Photopolymer Media, Jack Lyons, Dana Mackey, Izabela Naydenova

Articles

The theoretical modelling of holographic recording in photopolymers has been an important tool in their optimisation. More complex, hybrid organic/inorganic photopolymers have been developed in pursuit of materials with higher sensitivity, low shrinkage, high dynamic range and environmental stability. Recent attempts to augment the existing models for the redistribution of inorganic nanoparticles in holographic recording were successful but there is still a knowledge gap in regards to modelling optical losses, mutual cross-diffusion, the formation of slanted holographic gratings and polymerization induced shrinkage in hybrid photopolymer media. This paper will describe a novel approach to modelling the formation of unslanted holographic …

Solving Fractional Differential Equations On A Quantum Computer: A Variational Approach, Fong Yew Leong, Dax Enshan Koh, Jian Feng Kong, Siong Thye Goh, Jun Yong Khoo, Wei Bin Ewe, Hongying Li, Jayne Thompson, Dario Poletti

Research Collection School Of Computing and Information Systems

We introduce an efficient variational hybrid quantum-classical algorithm designed for solving Caputo time-fractional partial differential equations. Our method employs an iterable cost function incorporating a linear combination of overlap history states. The proposed algorithm is not only efficient in terms of time complexity but also has lower memory costs compared to classical methods. Our results indicate that solution fidelity is insensitive to the fractional index and that gradient evaluation costs scale economically with the number of time steps. As a proof of concept, we apply our algorithm to solve a range of fractional partial differential equations commonly encountered in engineering …

Resonant Solutions Of The Non-Linear Schrödinger Equation With Periodic Potential, Arein Duaibes, Yulia Karpeshina

Mathematics Faculty Publications

The goal is construction of stationary solutions close to non-trivial combinations of two plane waves at high energies for a periodic non-linear Schrödinger Equation in dimension two. The corresponding isoenergetic surface is described for any sufficiently large energy k2. It is shown that the isoenergetic surface corresponding to k2 is essentially different from that for the zero potential even for small potentials. We use a combination of the perturbative results obtained earlier for the linear case and a method of successive approximation.

Bioheat Equation Analysis, Johnathan Makar

Mathematics Student Work

In our research, we are investigating Pennes Bioheat equation, which is used for simulating the propagation of heat energy in human tissues. This equation was proposed by Pennes in 1948 based on his experiments of measuring the radial temperature distribution in the forearm of nine subjects. Pennes' equation provides the theoretical basis for studying heat transfer in perfused tissue and has been widely studied since then. However, Pennes' equation has been criticized for various reasons, including the fact that his experimental data did not seem to match the model. One of the objectives of our work is to find the …

An Augmented Matched Interface And Boundary (Amib) Method For Solving Problems On Irregular 2d Domains, Benjamin Pentecost

Mathematics Student Work

A new method called Augmented Matched Interface and Boundary (AMIB) has been developed to solve partial differential equation models, such as the heat equation, over irregular two-dimensional domains. The original AMIB method features unique numerical treatments to solve problems with various boundary conditions and shapes, resulting in highly accurate and efficient numerical solutions. However, recent numerical experiments have revealed that the original AMIB method can fail when dealing with sharply curved boundaries. To address this issue, new numerical techniques have been introduced in our latest work to enhance the robustness of the AMIB method. These techniques have been numerically verified …

Total Variation Flow In R^N Dimensions With Examples Relating To Perimeters Of Level Sets, Luis Schneegans, Victoria Shumakovich

Undergraduate Research Symposium

In this project, we explore radial solutions to the Total Variation Flow (TVF) equation with the help of the Sign Fast Diffusion Equation (SFDE) and prior results in the 1-dimensional case. Specifically for radial solutions, we derive equations and explicit solutions relating to the n-dimensional case. Lastly, we look at how level sets and (time) profiles change.

Hamiltonian Models For The Propagation Of Long Gravity Waves, Higher-Order Kdv-Type Equations And Integrability, Rossen Ivanov

Book chapter/book

A single incompressible, inviscid, irrotational fluid medium bounded above by a free surface is considered. The Hamiltonian of the system is expressed in terms of the so-called Dirichlet-Neumann operators. The equations for the surface waves are presented in Hamiltonian form. Specific scaling of the variables is selected which leads to a KdV approximation with higher order nonlinearities and dispersion (higher-order KdV-type equation, or HKdV). The HKdV is related to the known integrable PDEs with an explicit nonlinear and nonlocal transformation.

A Robust Model Reduction For The Optimal Boundary Feedback Stabilization Of Magnetizable Piezoelectric Beams, Md Rafi As Sadeq Ibn Emran

Masters Theses & Specialist Projects

The one-dimensional Partial Differential Equation (PDE) model of the wave equation, which describes wave dynamics of a wide-range of controlled mechanical systems, with a state feedback controller at its boundary is known to have exponentially stable solutions. Moreover, it is also reported that the several model reductions of the wave equation by the standard Finite Differences and Finite Elements approximations suffer from the lack of exponential stability (and exact observability without a state feedback controller) uniformly as the discretization parameter tends to zero. This is due to the loss of a uniform gap among the high-frequency eigenvalues as the discretization …

U-No: U-Shaped Neural Operators, Md Ashiqur Rahman, Zachary E Ross, Kamyar Azizzadenesheli

Department of Computer Science Faculty Publications

Neural operators generalize classical neural networks to maps between infinite-dimensional spaces, e.g., function spaces. Prior works on neural operators proposed a series of novel methods to learn such maps and demonstrated unprecedented success in learning solution operators of partial differential equations. Due to their close proximity to fully connected architectures, these models mainly suffer from high memory usage and are generally limited to shallow deep learning models. In this paper, we propose U-shaped Neural Operator (U-NO), a U-shaped memory enhanced architecture that allows for deeper neural operators. U-NOs exploit the problem structures in function predictions and demonstrate fast training, data …

Effects Of Topography On Tornado Paths Using Navier Stokes Equations, Kayleigh Smith

Mathematics Senior Capstone Papers

Tornadoes are relatively common in Louisiana with an average of 55 tornadoes per year. Predicting tornado paths has been extremely challenging due to the many factors that play into the formation of a tornado. According to NASA/ADS topography can have significant influence on tornado direction. The main goal of this research is to analyze tornado patterns and to determine if local topography has an effect on tornadic activity. Navier-Stokes partial differential equations will be used to model the data that is collected from the national weather service and a finite difference method will be used to solve the equations. Data …

Integrable Systems On Symmetric Spaces From A Quadratic Pencil Of Lax Operators, Rossen Ivanov

Conference papers

The article surveys the recent results on integrable systems arising from quadratic pencil of Lax operator L, with values in a Hermitian symmetric space. The counterpart operator M in the Lax pair defines positive, negative and rational flows. The results are illustrated with examples from the A.III symmetric space. The modeling aspect of the arising higher order nonlinear Schrödinger equations is briefly discussed.

The Lagrangian Formulation For Wave Motion With A Shear Current And Surface Tension, Conor Curtin, Rossen Ivanov

Articles

The Lagrangian formulation for the irrotational wave motion is straightforward and follows from a Lagrangian functional which is the difference between the kinetic and the potential energy of the system. In the case of fluid with constant vorticity, which arises for example when a shear current is present, the separation of the energy into kinetic and potential is not at all obvious and neither is the Lagrangian formulation of the problem. Nevertheless, we use the known Hamiltonian formulation of the problem in this case to obtain the Lagrangian density function, and utilising the Euler-Lagrange equations we proceed to derive some …

Modeling Of Population Dynamics In Spatially Continuous And Heterogeneous Environments Using Partial Differential Equations, Ryan T. St. Clair

Masters Theses & Specialist Projects

This thesis provides an expansion of existing reaction diffusion population ecology models into patchy landscapes with greater spatial heterogeneity. The Sturm-Loiville eigenvalue problem for a heterogeneous 3 patch landscape is solved implicitly and the population growth rate and migration dynamics on such a landscape are thoroughly discussed. In addition a system of interface equations is found whose greatest solution is the dominant eigenvalue of an N-patch landscape. The results of this model agree well with previously published results and show the mathematical equivalence of two published approaches for incorporating organism behavior at habitat interfaces. Increased modeling capabilities in heterogeneous environments …

The Gelfand Problem For The Infinity Laplacian, Fernando Charro, Byungjae Son, Peiyong Wang

Mathematics Faculty Research Publications

We study the asymptotic behavior as p → ∞ of the Gelfand problem

−Δ_pu = λe^u in Ω ⊂ Rⁿ, u = 0 on ∂Ω.

Under an appropriate rescaling on u and λ, we prove uniform convergence of solutions of the Gelfand problem to solutions of

min{|∇_u|−Λe^u, −Δ_∞u} = 0 in Ω, u = 0 on ∂Ω.

We discuss existence, non-existence, and multiplicity of solutions of the limit problem in terms of Λ.

Asymptotic Mean-Value Formulas For Solutions Of General Second-Order Elliptic Equations, Pablo Blanc, Fernando Charro, Juan J. Manfredi, Julio D. Rossi

Mathematics Faculty Research Publications

We obtain asymptotic mean-value formulas for solutions of second-order elliptic equations. Our approach is very flexible and allows us to consider several families of operators obtained as an infimum, a supremum, or a combination of both infimum and supremum, of linear operators. The families of equations that we consider include well-known operators such as Pucci, Issacs, and k-Hessian operators.

Existence And Uniqueness Of Minimizers For A Nonlocal Variational Problem, Michael Pieper

Honors Program: Senior Projects (Public)

Nonlocal modeling is a rapidly growing field, with a vast array of applications and connections to questions in pure math. One goal of this work is to present an approachable introduction to the field and an invitation to the reader to explore it more deeply. In particular, we explore connections between nonlocal operators and classical problems in the calculus of variations. Using a well-known approach, known simply as The Direct Method, we establish well-posedness for a class of variational problems involving a nonlocal first-order differential operator. Some simple numerical experiments demonstrate the behavior of these problems for specific choices of …

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

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 …

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 …

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 …