Partial Differential Equations Commons^™

Detecting And Repairing Conflicting Constraints In Co-Trained Physics-Informed Neural Networks For Composite Curing Processes, Cooper J. Evans

Dissertations, Master's Theses and Master's Reports

Composite materials have become a critical component of modern manufacturing, especially in the automotive and aerospace industries. The curing process for these composites has been modeled using a variety of partial differential equations representing the heat transfer and composite curing kinetics. Optimizing the applied temperature profile is critical for maximizing the efficiency and capacity of composite part manufacturers. Constraints must be placed on the inputs and outputs of the model, including but not limited to, the applied temperature profile, part temperature, and final degree of cure. Conflicting sets of constraints are easy to unknowingly impose due to the highly coupled …

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 …

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 …

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.

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 …

An Introduction To The Time-Independent Schrödinger Equation And Methods To Solve It, Vu Giang, Alex Gnech

OUR Journal: ODU Undergraduate Research Journal

The Time-Independent Schrödinger Equation is a linear elliptic PDE that describes quantum-mechanical systems. Its significance in the science of submicroscopic phenomena, particularly quantum mechanics, is as central as Newton’s laws of motion are to classical mechanics. This study uses various methods, including novel neural networks and finite difference schemes, to solve the one-dimensional two-body equation.

Proof-Of-Concept For Converging Beam Small Animal Irradiator, Benjamin Insley

Dissertations and Theses (Open Access)

The Monte Carlo particle simulator TOPAS, the multiphysics solver COMSOL., and

several analytical radiation transport methods were employed to perform an in-depth proof-ofconcept

for a high dose rate, high precision converging beam small animal irradiation platform.

In the first aim of this work, a novel carbon nanotube-based compact X-ray tube optimized for

high output and high directionality was designed and characterized. In the second aim, an

optimization algorithm was developed to customize a collimator geometry for this unique Xray

source to simultaneously maximize the irradiator’s intensity and precision. Then, a full

converging beam irradiator apparatus was fit with a multitude …

Monolithic Multiphysics Simulation Of Hypersonic Aerothermoelasticity Using A Hybridized Discontinuous Galerkin Method, William Paul England

Theses and Dissertations

This work presents implementation of a hybridized discontinuous Galerkin (DG) method for robust simulation of the hypersonic aerothermoelastic multiphysics system. Simulation of hypersonic vehicles requires accurate resolution of complex multiphysics interactions including the effects of high-speed turbulent flow, extreme heating, and vehicle deformation due to considerable pressure loads and thermal stresses. However, the state-of-the-art procedures for hypersonic aerothermoelasticity are comprised of low-fidelity approaches and partitioned coupling schemes. These approaches preclude robust design and analysis of hypersonic vehicles for a number of reasons. First, low-fidelity approaches limit their application to simple geometries and lack the ability to capture small scale flow …

Mathematical Models Yield Insights Into Cnns: Applications In Natural Image Restoration And Population Genetics, Ryan Cecil

Electronic Theses and Dissertations

Due to a rise in computational power, machine learning (ML) methods have become the state-of-the-art in a variety of fields. Known to be black-box approaches, however, these methods are oftentimes not well understood. In this work, we utilize our understanding of model-based approaches to derive insights into Convolutional Neural Networks (CNNs). In the field of Natural Image Restoration, we focus on the image denoising problem. Recent work have demonstrated the potential of mathematically motivated CNN architectures that learn both `geometric' and nonlinear higher order features and corresponding regularizers. We extend this work by showing that not only can geometric features …

Ocean Wave Prediction And Characterization For Intelligent Maritime Transportation, Pujan Pokhrel

LSU New Orleans Theses and Dissertations

The national Earth System Prediction (ESPC) initiative aims to develop the predictions
for the next generation predictions of atmosphere, ocean, and sea-ice interactions in the scale of days to decades. This dissertation seeks to demonstrate the methods we can use to improve the ESPC models, especially the ocean prediction model. In the application side of the weather forecasts, this dissertation explores imitation learning with constraints to solve combinatorial optimization problems, focusing on the weather routing of surface vessels. Prediction of ocean waves is essential for various purposes, including vessel routing, ocean energy harvesting, agriculture, etc. Since the machine learning approaches …

A Bidirectional Formulation For Walk On Spheres, Yang Qi

Dartmouth College Master’s Theses

Poisson’s equations and Laplace’s equations are important linear partial differential equations (PDEs)
widely used in many applications. Conventional methods for solving PDEs numerically often need to
discretize the space first, making them less efficient for complex shapes. The random walk on spheres
method (WoS) is a grid-free Monte-Carlo method for solving PDEs that does not need to discrete the
space. We draw analogies between WoS and classical rendering algorithms, and find that the WoS
algorithm is conceptually identical to forward path tracing.
We show that solving the Poisson’s equation is equivalent to solving the Green’s function for every
pair of …

A Meshless Approach To Computational Pharmacokinetics, Anthony Matthew Khoury

Doctoral Dissertations and Master's Theses

The meshless method is an incredibly powerful technique for solving a variety of problems with unparalleled accuracy and efficiency. The pharmacokinetic problem of transdermal drug delivery (TDDD) is one such topic and is of significant complexity. The locally collocated meshless method (LCMM) is developed in solution to this topic. First, the meshless method is formulated to model this transport phenomenon and is then validated against an analytical solution of a pharmacokinetic problem set, to demonstrate this accuracy and efficiency. The analytical solution provides a locus by which convergence behavior are evaluated, demonstrating the super convergence of the locally collocated meshless …

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 …

Inference Of Surface Velocities From Oblique Time Lapse Photos And Terrestrial Based Lidar At The Helheim Glacier, Franklyn T. Dunbar Ii

Graduate Student Theses, Dissertations, & Professional Papers

Using time dependent observations derived from terrestrial LiDAR and oblique
time-lapse imagery, we demonstrate that a Bayesian approach to glacial motion es-
timation provides a concise way to incorporate multiple data products into a single
motion estimation procedure effectively producing surface velocity estimates with
an associated uncertainty. This approach brings both improved computational effi-
ciency, and greater scalability across observational time-frames when compared to
existing methods. To gauge efficacy, we apply these methods to a set of observa-
tions from the Helheim Glacier, a critical actor in contemporary mass loss trends
observed in the Greenland Ice Sheet. We find that …

Image Restoration Using Automatic Damaged Regions Detection And Machine Learning-Based Inpainting Technique, Chloe Martin-King

Computational and Data Sciences (PhD) Dissertations

In this dissertation we propose two novel image restoration schemes. The first pertains to automatic detection of damaged regions in old photographs and digital images of cracked paintings. In cases when inpainting mask generation cannot be completely automatic, our detection algorithm facilitates precise mask creation, particularly useful for images containing damage that is tedious to annotate or difficult to geometrically define. The main contribution of this dissertation is the development and utilization of a new inpainting technique, region hiding, to repair a single image by training a convolutional neural network on various transformations of that image. Region hiding is also …

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 …

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 …

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 …

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 …