Applied Mathematics Commons^™

A 1d Symmetric Interior Penalty Discontinuous Galerkin Solver In Rust, William Aey

Williams Honors College, Honors Research Projects

This honors project will build a 1D Symmetric Interior Discontinuous Galerkin (SIPDG) solver in Rust for Stum-Liouville type problems such as the Poisson equation, with Robin, Dirichlet, and Neumann boundary conditions. The work will cover the full pipeline: starting from the strong form of the PDE, deriving the DG weak form, implementing element and interface operators, and assembling or apply the discrete operator. Rust's safety and concurrency (e.g, via Rayon) will be used to explore serial and parallel performance. A test-driven development approach will be used to maintain a strong suite of tests. The project will result in a documented …

Eigenvalue Spacing Distributions And The Weak Disorder Limit For Random Schrodinger Operators, Kyle E. Hammer

Theses and Dissertations--Mathematics

We study a collection of discrete Schrodinger Operators with random potentials through the lens of global and local eigenvalue spacings. We discuss the three models: the standard scaled disorder Anderson Model, the Anderson-Bernoulli Polymer Model, and the Discrete Fractional Laplacian Anderson Model. First, we discuss the scaled disorder case using the invariant measure and its application to the density of states in the weak disorder limit. We also prove the limit of the local and global eigenvalue spacings in the non random case, and demonstrate numerically how randomness affects the eigenvalue spacings. We then discuss a special family of random …

A Mathematical Frameworks For Singular, Nonlinear Phenomena: Applications To Nematocyst Firing And Inhomogeneous Nls With Coulomb Potential, Abdulrahman Alharbi

Theses and Dissertations

Nematocysts are specialized cellular organelles found in all cnidarians, including corals and jellyfish, as well as in some single-celled protists such as dinoflagellates. These organelles display remarkable diversity in morphology and function, playing roles in prey capture and defense. The firing of a nematocyst is one of the fastest accelerations in nature, yet the underlying physical mechanisms remain not fully understood. In this work, we address key questions: how sufficient force is generated to overcome the fluid boundary layer, whether fluid–structure interaction models can reproduce observed dynamics, and what mechanisms trigger discharge.

Our research investigates models based on osmotic pressure …

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.

Entropic Foundation Of Finance And Physics: Securities Price Dynamics And Quantum Theory, Mohammad Abedi

Electronic Theses & Dissertations (2024 - present)

In many scientific and financial contexts, we must reason and make predictions under conditions of incomplete information. This dissertation develops Entropic Dynamics (ED) as a unified framework for deriving dynamical laws directly from principles of inference. Within this approach, probability distributions represent states of knowledge, and their evolution is determined through entropy maximization subject to relevant constraints. This leads to a novel concept of entropic time and a formulation of dynamics as an inferential process. In this talk, I will present how ED provides a common foundation across multiple domains. In physics, quantum dynamics for particles and scalar fields in …

The Ellipsoidal Principal Semi-Axis Geometry Of The Solution To An Ivp For A Matrix Diffusion Pde, Brandon M. Fox

Electronic Theses & Dissertations (2024 - present)

We analyze the underlying geometry of the solution to an IVP for a matrix diffusion PDE. We first derive the fundamental solution to the PDE. We then determine the unique solution to the IVP. From there, we begin analyzing its underlying geometry. We first observe that the geometry exhibits an ellipsoidal nature. Furthermore, we observe that it is described by the principal semi-axis geometry of the ellipsoids associated with the solution. This conclusion follows from applying the Principal Axis Theorem to the associated ellipsoids to establish their principal semi-axis geometry, as governed by the eigenstructure of the matrix. This thesis …

Performance Of Numerical Methods Applied To The Black–Scholes Model, Scott Cameron Williams

UNF Graduate Theses and Dissertations

We compare five numerical approaches for approximating solutions to the Black–Scholes partial differential equation for pricing European call options: FTCS, BTCS, Crank– Nicolson, Monte Carlo simulation, and a physics–informed neural network (PINN). These methods span finite difference techniques, probabilistic simulation, and machine learning. Performance is evaluated based on computational efficiency and accuracy relative to the analytical Black–Scholes solution.

Among the methods, Crank–Nicolson and the PINN demonstrated the strongest overall performance. Crank–Nicolson achieved the highest accuracy but exhibited increased runtime as the number of underlying stock price grid points grew. In contrast, the PINN produced slightly less accurate results but with …

A Reconstruction Method For The Anisotropic Electrical Conductivity In A Hybrid Planar Inverse Problem, Athanasios Dimitriadis

Graduate Studies Theses and Dissertations 2026

Some single-physics medical imaging methods fail to reliably discriminate benign from malignant biological tissue. In partial response, a new class of inverse problems employ multi-physics phenomena, coupling methods of high resolution with methods of high contrast. In this thesis we are concerned with one such problem in which an anisotropic electrical conductivity is to be recovered from some internal knowledge of the current density field generated by maintaining a fixed boundary voltage. This internal data is obtained by the method of Current Density Impedance Imaging (CDII), which combines Maxwell's equations with Magnetic Resonance Imaging (MRI) measurements. The anisotropy class is …

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 the loop expansion perturbatively renormalizable.

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 include existence for all cases that …

Mathematical Model Of Graphene, Douglas M. Sanor

Williams Honors College, Honors Research Projects

Graphene, a single-atom-thick layer of carbon arranged in a hexagonal lattice, exhibits exceptional mechanical, electrical, and thermal properties that make it a promising material for a wide range of engineering applications. This paper presents a mathematical framework for modeling the mechanical behavior of graphene, with a focus on atomistic-to-continuum approaches. We begin with a onedimensional Frenkel-Kontorova model that represents graphene as a discrete chain of particles interacting with both their nearest neighbors through harmonic spring potentials and an underlying substrate through van der Waals forces. Numerical simulations of this discrete model demonstrate the commensurate-toincommensurate phase transition, revealing how geometric mismatch …

Optimal Control Of Stochastic Systems: A Numerical Study Of The Stochastic Linear Quadratic Regulator Framework, Gülşen Orucova Büyüköz, Yaprak Güldoğan Dericioğlu, Tuğçem Partal

Mathematical Modelling and Numerical Simulation with Applications

Optimal control of stochastic linear systems is fundamental in control theory, with applications in robotics, finance, and engineering. The Stochastic Linear Quadratic Regulator (SLQR) derives optimal feedback laws via the Riccati equation but requires numerical discretization of the resulting stochastic dynamics. Despite extensive studies on numerical methods for stochastic differential equations, their performance within the SLQR framework remains insufficiently explored. This study compares two predictor–corrector schemes of different orders: the Order 1.0 Predictor-Corrector (PC) method and the Order 2.0 Weak PC method. A one-dimensional linear quadratic problem with a closed-form solution enables precise error evaluation against the analytical trajectory. Convergence …

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 …

On Sobolev Spaces And The Existence Of Weak Solutions To Boundary Value Problems, Skye X. Paul

Master's Theses

Many boundary value problems that arise in mathematical models have close connections to second order elliptic partial differential equations. This thesis introduces the idea of weak derivatives and Sobolev Spaces to generalize possible solutions. Using functional analysis centered around the Lax-Milgram theorem, we show the existence of these generalized solutions to boundary value problems including Laplace's Equation, 2nd order linear ODEs, and ultimately a general second order elliptic PDE. The work cumulates with recovering a number of central theorems of functional analysis in the context of Sobolev Spaces, creating a new perspective on the solvability of these boundary value problems.

Computational Models For Pre-Lens Tear Film Drug Concentration Dynamics With Drug Supply From A Contact Lens And Drug Exchange During Blinking, Mazen A. Althobaiti

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Mathematical Model Of Ovarian Cancer Tumor Growth In Mice, Jessica A. Hoffman

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Can We End The Hiv Epidemic In The U.S.? Linking Clinical And National Surveillance Data Through Multiscale Modeling From Patients To Populations, Necibe Tuncer

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Quantification Of Parameters To Predict The Rupture Of Intracranial Saccular Aneurysms Using Physics Informed Neural Networks, Alonso Gabriel Ogueda, Padmanabhan Seshaiyer

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

(Si15-093) Dynamics Of The Generalized Nonlinear Schrödinger Equation With A Source Using An Analytical Method, Manish Raghav, Manoj .

Applications and Applied Mathematics: An International Journal (AAM)

In the present article, an analytical method is used to obtain hyperbolic, trigonometric, and rational solutions of the generalized nonlinear Schrödinger (GNLS) equation with a source. The ability of solitons to preserve their shapes during propagation makes them suitable for optical fiber communication. Solutions to the generalized nonlinear Schrödinger equation with a source can also describe solitons, and understanding their dynamics helps to design communication systems based on solitons. The analytical method used is compelling and effective for finding exact solutions to various nonlinear evolution equations (NLEEs). To further understand the phenomena, we create 3−D, contour, and 2−D graphs of …

(Si15-113) Augmenting Cryptographic Security Through Inventive Application Of The Kharrat-Toma Transform Algorithm, Prabakaran Raghavendran, Tharmalingam Gunasekar, K. Sakthivel, Kamalendra Kumar, Shalini Gupta

Applications and Applied Mathematics: An International Journal (AAM)

This paper introduces a cryptographic technique combining the Kharrat-Toma Transform and congruence modulo operators to improve the security of message encryption. The proposed model uses the mathematical properties of the Kharrat-Toma Transform and its inverse for direct scrambling and unscrambling processes while embedding sufficient complexity to resist modern cryptanalytic attacks. The model is subjected to experimental tests, including encryption quality analysis, Shannon entropy, and NIST randomness tests, in order to prove the strength of the model. Through encryption quality analysis, symbol frequencies in the ciphertext are masked heavily from having much correlation between plaintext and ciphertext. Entropy values indicate near-theoretical …

Meshless Discrete Velocity Boltzmann Model For Porous Media Flow, Amandine Maidenberg

Doctoral Dissertations and Master's Theses

This dissertation explores the combination of two sophisticated techniques for addressing computational fluid dynamics: the discrete velocity Boltzmann equation (DVBE) and the localized collocation meshless model with upwinding (U-LCMM). The DVBE is a high-level model that describes the foundations of transport phenomena by addressing the microscale motions of particles themselves and the effect of their aggregate behaviors on continuum principles. This equation integrates multiple scales of phenomena; while it can be used for fluid flow at Navier-Stokes scales, it can also resolve fine features that can only be described at the molecular level. This type of model is necessary for …

Stabilized Weak-Gradient Discontinuous Finite Elements With Optimal Error Estimates For Second-Order Elliptic Pdes, Aymen Laadhari

Mathematical Modelling and Numerical Simulation with Applications

This work introduces an accurate finite element approach employing a new stabilized discrete weak gradient, designed for second-order elliptic problems on arbitrary conforming meshes. We formulate the approach within a discontinuous Galerkin framework and derive a consistent and coercive bilinear form. Appropriate error analysis on a model problem confirms optimal convergence. Building on the core analysis, we extend the method to more challenging settings, including time-dependent heterogeneous scenarios and a biophysically realistic optimal-control model of photobleaching in the budding yeast cell. We further illustrate the versatility of the weak-gradient construction by applying it to an unsteady level-set equation relevant to …

Solution Of Fractional Order Diffusion Equations With Clique Neural Network, Merve Zeynep Kaya, Mesut Karabacak, Ercan Çelik

Mathematical Modelling and Numerical Simulation with Applications

In this paper, the clique artificial neural network method is used to solve the fractional diffusion equation, which is a subclass of partial differential equations. The clique neural network architecture is constructed using input, hidden, and output layers. Several degrees of clique polynomials were used as activation functions, and the output layer was obtained by multiplying them with weight coefficients. Subsequently, the optimization equation was derived, and the exact solution, numerical solution, and error function graphs were obtained using a specialized algorithm. Analysis of the results demonstrates that the clique artificial neural network method provides quicker and more accurate results …

Accurate Temporal Integration Schemes For Nonlinear Adsorption Problems, Evan D. Butterworth

All Dissertations

We consider a nonlinear transport problem to model the chromatography process of high-capacity multimodal membranes. Robust and efficient algorithms that simulate these bioseparation processes are critical to developing therapeutics for various chronic illnesses and infectious diseases. However, much of the current methodology focuses on stabilization and linearization techniques, often implementing low-order time-discretizations and linearized adsorption, resulting in inefficiencies and inaccuracies in the numerical solution. Utilizing Rothe's method, we develop various time-discretization schemes coupled with the finite element method to solve the fully implicit problems. Stability and solvability results are presented for several methods. Through multiple high-level software implementations paired with …

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.

A Modified Sir Model Used To Investigate The Relationship Between Congenital And Adult Syphilis, Kaleesta R. Waysman

Honors Thesis

A complex SIR model integrating the relationship between adult and congenital syphilis was developed. The goal of the project was to determine the specific population(s) or control strategies that should be enforced, altered, or removed to decrease the number of children experiencing clinical sequelae due to congenital syphilis. Early clinical sequelae include hydrops fetalis, preterm birth, central nervous system infection, hepatosplenomegaly, hyperbilirubinemia, cholestasis, hemolytic anemia, snuffles, osteochondritis, and lesions or rashes in the palms and soles. Late clinical sequelae include interstitial keratitis, hearing loss, Hutchinson teeth, saber shins, Clutton joints, mulberry molars, and saddle nose. After implementing real-world data into …

Robust And Efficient Solvers For Physics-Based Pde’S, Elizabeth Hawkins

All Dissertations

This work was partially supported by the U.S. Department of Energy under award DE- SC0025292, by NSF grant DMS 2152623, and by NSF grant DMS 2011490.

This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Advanced Scientific Computing Research, Mathematical Multifaceted Integrated Capability Centers (MMICCs) program, under Field Work Proposal 22-025291 (Multifaceted Math- ematics for Predictive Digital Twins (M2dt)), Field Work Proposal 23-020467, and Computing and Information Sciences (CIS) investment area in the Laboratory Directed Research and Development program at Sandia National Laboratories. This written work is authored by an employee …

Domain Decomposition For Coupled Systems Of Fluid-Structure Interaction And Numerical Modeling For Thin Film Polymers, Amy De Castro

All Dissertations

We consider two primary areas of physical application in this work: fluid interaction systems with either linear elastic structures or with poroelastic structures, and thin film polymers, where the majority of the work focuses on the fluid-structure interaction systems.

In the first chapter, we present a strongly coupled partitioned method for fluid structure interaction (FSI) problems based on a monolithic formulation of the system which employs a Lagrange multiplier (LM). We prove that both the semi-discrete and fully discrete formulations are well-posed. To derive the partitioned scheme, a Schur complement equation, which implicitly expresses the Lagrange multiplier and the fluid …

Globally Adaptive Exponential Integrators For Stiff Systems Of Odes, Anzhelika Vasilyeva

Honors Theses

This thesis introduces a novel method for solving systems of Ordinary Differential Equations (ODEs) resulting from the spatial discretization of Partial Differential Equations (PDEs). The proposed approach builds upon an existing technique that employs Krylov projection, which requires evaluating a matrix function at each timestep. The innovation of the new method lies in its reuse strategy, which shifts the perspective from direct matrix function evaluation to polynomial interpolation. Numerical experiments conducted on constant and variable coefficient heat equations, with both smooth and discontinuous initial data, demonstrate the computational time advantage of the new approach. The results indicate that this method …

Properties Of Eigenvalues Of The Fractal Laplacian, Eric Stachura, Andrew Chincea

Symposium of Student Scholars

We investigate the properties of the eigenvalues of the fractal Laplacian. We begin by defining the fractal Laplacian operator in one dimension and formulate the corresponding Dirichlet eigenvalue problem. Analytical solutions are obtained for specific fractal parameters, and computational results illustrate the structure of eigenvalues and their associated eigenfunctions. We extend our analysis to two dimensions using separation of variables. Our findings contribute to a deeper understanding of how fractal geometry affects the spectral characteristics of differential operators.

Unfitted Finite Element Methods For Shape Optimization And Liquid Crystals, Jeremy T. Shahan

LSU Doctoral Dissertations

We present an approach to shape optimization problems that uses an unfitted finite element method (FEM). The domain geometry is represented, and optimized, using a (dis- crete) level set function and we consider objective functionals that are defined over bulk domains. For a discrete objective functional, defined in the unfitted FEM framework, we show that the exact discrete shape derivative essentially matches the shape derivative at the continuous level. In other words, our approach has the benefits of both optimize-then- discretize and discretize-then-optimize approaches.

Specifically, we establish the shape Fréchet differentiability of discrete (unfitted) bulk shape functionals using both the …