Numerical Analysis and Scientific Computing Commons^™

Muon Lifetime: Theory, Experiment, And Simulation, Soren Agustin Munoz

Physics

This senior project investigates the muon lifetime through three complementary approaches: theoretical calculation, laboratory measurement, and computational simulation. The theoretical component develops the necessary background from relativistic field theory to the effective weak interaction, culminating in the leading-order Fermi-theory prediction of τ ≈ 2.2 μs, which explains why the muon lifetime lies on the microsecond scale.

The experimental component measures the lifetime of stopped cosmic-ray muons using a plastic scintillator, photomultiplier tube, and analog timing electronics. A binned Poisson likelihood fit to the primary 15-day acquisition run τ = 2.17^+0.03_-0.09 μs, consistent with the accepted value within the …

Effects Of Coulomb Collisions On The Structures Of Fast Magnetosonic Shocks, Bo Farnell

Physics and Astronomy Undergraduate Senior Theses

We present fully kinetic particle-in-cell simulations demonstrating qualitative effects of binary Coulomb collisions on fast magnetosonic shocks. We find that with a sufficiently high collisional frequency, shock rippling and reformation can be inhibited, creating stationary, laminar shocks. We see that collisions rapidly bring the reflected population into equilibrium with the upstream population, adding credibility to existing theories that the separation of these two populations create these dynamical behaviors around the shock transition region. We also see that Ohmic heating from Coulomb collisions can suppress instabilities.

Towards Physics-Informed Neural Networks For Simulating Multiphase Geothermal Convection​*, Daniel C. Patton, Andrew Harrison Eno

Campus Research Month

Water and steam flow through porous rock, transferring heat via conduction and buoyancy-driven convection caused by density differences. Traditional numerical methods (finite-volume/finite-element) model this well but can become memory-intensive and unstable for long, high-detail simulations. This work demonstrates a Physics-Informed Neural Network (PINN) using a finite-difference approach within the NVIDIA PhysicsNeMo framework to simulate magma chambers in 2D. Tested on the Rio Pisco pluton in Peru, results are compared with the USGS HYDROTHERM model. PINNs learn from physical laws, offering accurate, flexible solutions with less data and development effort.

Developing A Multi-Band Gnss Remote Sensing Station To Measure Precipitable Water Vapor For Flash Flood Nowcasting, Madison Gleydura

Doctoral Dissertations and Master's Theses

Flash flood nowcasting in Central and Southern Appalachia is particularly challenging due to steep terrain, narrow valleys, highly localized rainfall patterns, and limited measurement coverage. Traditional remote sensing methods, such as Doppler radar and microwave radiometry, suffer from reduced resolution at long range and signal blockage by mountains. GNSS-meteorology offers an established alternative for measuring precipitable water vapor and is currently integrated into several numerical weather models. Recent research demonstrates that commercial-grade GNSS receivers can produce tropospheric products comparable to those from geodetic-grade equipment. The gaps in mountain coverage can be addressed by developing a low-cost, self-contained embedded system that …

Basis Design For Electronic Structure And Beyond, Weishi Wang

Dartmouth College Ph.D Dissertations

At the intersection of quantum physics, quantum chemistry, and materials science, electronic structure is the study of electrons in solid-state and molecular systems. Electronic-structure computation relies on discretizing the many-electron Hamiltonian with a finite single-particle basis set. However, basis-set construction is conventionally treated as an ad hoc preprocessing step. This thesis develops an expressive and flexible framework for active, system-oriented basis-set design and numerical modeling strategies that treat basis functions as tunable representations to encode electronic ground-state information.

We first introduce a multi-layered, differentiable basis-construction framework that embeds a set of primitive parameters into mixed-contracted Gaussian-type orbitals. We then develop …

Utilizing Computer Modeling To Optimize Electric Fields Within Xenon Time Projection Chambers, Miles Meloni

Honors Theses

XENONnT is a physics experiment designed with the goal of detecting dark matter particles. The detector is a time projection chamber; a series of charged electrodes creates an electric field, surrounding a central body filled with liquid and gaseous xenon. Photomultiplier tubes (PMTs), positioned on either end of the chamber, serve to detect light signals. We seek to minimize the root mean square of the electric field norms experienced by the PMTs. This quantity corresponds to the variance in the electric field observed by the PMTs. Establishing a consistent electric field is important to maintaining these sensitive components. The electric …

Stability Analysis Of Thermohaline Convection With A Time-Varying Shear Flow Using The Lyapunov Method, Kalin Kochnev

Honors Scholar Theses

This work applies the Lyapunov method to identify instabilities and compute the growth rate of a linear time-varying system. The linear system studied describes cold fresh water on top of hot salty water with a periodically time-varying background shear flow. A time-dependent weighting matrix is employed to construct a Lyapunov function candidate. The resulting linear matrix inequalities are discretized in time using the forward Euler method. As the number of temporal discretization points increases, the growth rate predicted by the Lyapunov method or Floquet theory, used for comparison, will converge to the same value obtained from numerical simulations. Furthermore, the …

Establishing Convergence Thresholds For Pre-Trajectory Sampling With Batched Execution Across Random Quantum Circuits, Taylor L. Eskew, Jerome F. Gonthier, Taylor L. Patti, Andrew N. Jordan

Student Scholar Symposium Abstracts and Posters

A crucial aspect of validating quantum protocols is understanding the noise produced by quantum computing devices. Using simulations that can replicate this noise allows for a lower-cost alternative to hardware experiments. Stochastic, so-called "trajectory" methods are often used as a quadratically reduced approximation to density matrix simulations, but traditional implementations have limited sampling capacity and provide no error-based metadata. The Pre-Trajectory Sampling with Batched Execution (PTSBE) [Patti et al., 2025] algorithm provides a solution by combining fine-tuned, well-documented noise sampling with computational intermediate caching.

While the original work is effective on quantum error correction circuits, its performance on general circuits …

Modeling Synaptic Dysfunction As Neural Contagion: A Graph-Based Sedr Framework For Simulating Signal Spread, Michelle Marfo, Dr. Padmanabhan Seshaiyer, Alonso Ogueda-Oliva

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 …

Design And Development Of A Standalone Digital Holographic Microscope Employing Phase-Driven Reconstruction And Classification For Biomedical Imaging And Optical Diagnostics, Charlotte Kyeremah

Graduate Doctoral Dissertations

Access to advanced biomedical imaging technologies remains a significant challenge in resource-limited settings, especially for early disease detection and monitoring of diseases such as malaria, HIV, and other blood-borne diseases. Although point-of-care (POC) devices have gained popularity in global health, many rely on antibody-based tests, lateral flow strips, or optical readouts that often lack quantitative capabilities, sensitivity to early infections, or versatility in different diagnostic targets. In addition, these systems are typically dependent on disposable reagents or manual interpretation, which limits their effectiveness in remote areas. Digital Holographic Microscopy (DHM) presents a promising alternative as a label-free imaging method capable …

Simulating Chemiluminescence Spectra Of Oh Radical Using Quantum Dynamical Simulations, Triet M. Cao

Research from the Berry Summer Thesis Institute, 2025

Chemiluminescence is a phenomenon of light emission that occurs when molecules transition from an excited state to the ground state. This phenomenon is widely utilized in applications ranging from quantifying product yields in chemical reactions to biological imaging and atmospheric chemistry studies. It has also been used in many crucial applications in the scarce conditions of heat and electricity, such as light sources for astronomy discovery or deep-sea diving. A notable example of such reactions is combustion, where highly reactive hydroxyl radicals (OH) are products. These radicals exhibit a chemiluminescence spectrum because of the interaction between their ground state and …

Algebraic Multigrid Methods For Nonsymmetric And Indefinite Problems: Theory And Applications, Ahsan Ali

Mathematics & Statistics ETDs

Algebraic multigrid (AMG) is a well-established and highly efficient solver for symmetric positive definite (SPD) systems arising from elliptic and parabolic PDEs, while nonsymmetric systems from hyperbolic PDEs remain a significant challenge. This dissertation develops AMG methods and theory for nonsymmetric problems. First, we develop a novel approach combining mode constraints from energy-minimization AMG with local approximations of ideal restriction in $\ell$AIR, resulting in constrained $\ell$AIR (C$\ell$AIR), which demonstrates scalable convergence across advective and diffusive problems. Second, we extend optimal AMG theory by deriving spectral radius estimates for the two-grid error transfer operator using matrix-induced orthogonality, enabling convergence predictions for …

Physics Embedded Neural Network: A Novel Data-Free Numerical Method For Solving Computational Physics Problems, Pawan Gaire

Dissertations and Doctoral Documents, University of Nebraska-Lincoln, 2023–

A novel approach for solving partial differential equations (PDEs) using neural networks for scientific computing is introduced. The proposed approach, referred to as physics-embedded neural network (PENN), features a unique architecture that incorporates the PDE and boundary conditions information directly within the final fully-connected layer of the feed-forward neural network (NN). The key aspect of PENN is the parallel numerical embedding of a differential equation associated with physical problems within the activation function of the network’s final layer. This integration leads to a new class of computational solvers competitive with classical methods like the Finite Element Method (FEM) and capable …

Robust Spacecraft Autonomy For Deep Space Exploration In Special Euclidean Group Se(3), Matthew Wittal

Doctoral Dissertations and Master's Theses

Over the past half-century, humanity has gained extensive experience conducting manned spaceflight near Earth. Arguably, "near Earth" could even include the Moon — the most distant destination humans have reached. However, "near" in this work primarily refers low Earth orbit (LEO). One could argue that we have not truly left Earth since the Apollo, as spacecraft in some LEOs remain subject to atmospheric drag thus emphasizing their continued connection to Earth's immediate environment. Reflecting on this, it becomes clear that humanity has largely remained bound to Earth’s immediate vicinity since the Apollo missions reached the Moon. However, that is set …

Scalable Parallel-In-Time Integration For Equations Of Motion, Nathan W. Chapman

All Master's Theses

Physical simulations always need to balance accuracy and run-time. This work implements the Parareal Algorithm using graphics processing units across a distributed system to accurately simulate time-dependent physics while attempting to minimize runtime. Data-transfer latency is identified as the primary bottleneck, for which mitigation methods are provided. Benchmarks comparing single-GPU and distributed implementations on a spectrum of coarse and fine discretizations are analyzed.

Development And Application Of Computational Tools For Data-Driven Materials Science., Logan L. Lang

Graduate Theses, Dissertations, and Problem Reports (ETD)

Modern materials science generates vast amounts of data from computational simulations and experiments, creating significant challenges for data processing and analysis. This thesis addresses these challenges through the development and application of computational tools within the framework of Material Data Science (MDS). Contributions span the four pillars of MDS: Material/Molecular Data, Algorithms, Databases, and High-Throughput Processes—with a primary focus on the Algorithm, Data, Database pillars.

For the Algorithm pillar, two Python libraries were developed to streamline common analysis tasks. PyProcar simplifies the post-processing and visualization of electronic structure data (band structures, density of states, Fermi surfaces) obtained from various Density …

Modeling Relativistic Fluids In Dynamical Spacetimes, Terrence Alphonse Pierre Jacques

Graduate Theses, Dissertations, and Problem Reports (ETD)

Multi-messenger astrophysics opens a new era in our understanding of the most dynamic and energetic systems in the Universe. Correlating gravitational-wave and electromagnetic signals in space and time enables stringent tests of models for core-collapse supernovae, merging supermassive black-hole binaries with accretion disks and jets, and mergers of compact object binaries such as binary neutron stars (BNS) and white dwarfs. Comparisons between models and multi-messenger observations may be used to constrain the neutron-star equation of state (EOS), formation channels for compact-object binaries, and emission mechanisms behind short gamma-ray bursts.

In modeling such astrophysical systems, great success has been achieved by …

Personalized Physics Learning Through Ai: Insights From Problem Generation, Chatbot Dialogues, And Intelligent Tutoring Systems, Atharva Dange

Physics Dissertations - Archive

Artificial intelligence (AI) is poised to transform science education, yet questions remain on how best to integrate these technologies into teaching and learning. This dissertation investigates the use of AI-driven tools in university physics courses through three complementary studies. In the first study, a generative language model (ChatGPT) was used to create novel physics homework problems aligned with course objectives. Analysis showed that, after expert vetting, AI-generated questions can foster higher-order problem-solving and reduce student reliance on solution memorization, though careful instructor oversight is required to ensure accuracy. The second study embedded an AI chatbot as a learning aid in …

The Presence Of Outer Giant Planets And Their Role In Inner Planet Formation With And Without Their Influence, Mateo E. Guerra Toro

Graduate Theses/Dissertations

We performed dynamical simulations of the giant impact phase of planet formation to investigate the formation of inner terrestrial planets under the influence of 4 solar system-like outer giant planets. We developed a new code using the N-body simulation suite REBOUND and REBOUNDx (Rein et al. (2019) and Tamayo et al. (2019)) to simulate 2 stages of planetary formation: a residual gaseous protoplanetary disk phase and subsequent dynamical evolution after the disk photoevaporates. The initial conditions for the inner planetary embryos were taken by Morrison et al. (2020) based on a range of solid surface densities that produced Super-Earth terrestrial …

Linking Empirical Data And Numerical Simulation To Characterize Dynamic Fire Behavior Associated With Interacting Firelines, Marta Sergeevna Jerebets

Graduate Student Theses, Dissertations, & Professional Papers

Understanding fuel pattern-fire process relationships is key for predicting fire behavior and effects with follow-on benefits to proactive fire management and model validation. To characterize dynamic fire behavior, this thesis leverages empirical data and numerical simulation through two complementary studies.

In the first study, longwave thermal sensors aboard unmanned aerial systems (UAS) were used to capture fine-scale fire behavior in two experimental grass burns. A novel paired design was used to quantify the effects of fuel arrangement on fire behavior with 3.66 m diameter treatments cut to a height of 0.15 m. The treatments ephemerally reduced fire rate of spread …

A Universal Implementation Of Radiative Effects In Neutrino Event Generators, Júlia Tena-Vidal, Adi Ashkkenazi, Lawrence B. Weinstein, Peter Blunden, Steven Dytman, Noah Steinberg

Physics Faculty Publications

Due to the similarities between electron-nucleus (eA) and neutrino-nucleus scattering (νA), eA data can contribute key information to improve cross-section modeling in eA and hence in νA event generators. However, to compare data and generated events, either the data must be radiatively corrected or radiative effects need to be included in the event generators. We implemented a universal radiative corrections program that can be used with all reaction mechanisms and any eA event generator. Our program includes real photon radiation by the incident and scattered electrons, and virtual photon exchange and photon vacuum polarization diagrams. It …

State Preparation Of Lattice Field Theories Using Quantum Optimal Control, Jack Y. Araz, Siddhanth Bhowmick, Matt Grau, Thomas J. Mcentire, Felix Ringer

Physics Faculty Publications

We explore the application of quantum optimal control (QOC) techniques to state preparation of lattice field theories on quantum computers. As a first example, we focus on the Schwinger model, quantum electrodynamics in 1 + 1 dimensions. We demonstrate that QOC can significantly speed up the ground state preparation compared to gate-based methods, even for models with long-range interactions. Using classical simulations, we explore the dependence on the interqubit coupling strength and the device connectivity, and we study the optimization in the presence of noise. While our simulations indicate potential speedups, the results strongly depend on the device specifications. In …

Geometric Gnns For Charged Particle Tracking At Gluex, Ahmed Hossam Mohammed, Kishansingh Rajput, Simon Taylor, Denis Furletov, Sergey Furletov, Malachi Schram

Computer Science Faculty Publications

Nuclear physics experiments are aimed at uncovering the fundamental building blocks of matter. The experiments involve high-energy collisions that produce complex events with many particle trajectories. Tracking charged particles resulting from collisions in the presence of a strong magnetic field is critical to enable the reconstruction of particle trajectories and precise determination of interactions. It is traditionally achieved through combinatorial approaches that scale worse than linearly as the number of hits grows. Since particle hit data naturally form a point cloud and can be structured as graphs, graph neural networks (GNNs) emerge as an intuitive and effective choice for this …

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.

Numerical Techniques For Improving Simulations Of Tropical Cyclones, Yassine Tissaoui

Dissertations

The increasing frequency and intensity of tropical cyclones (TCs) due to climate change pose significant challenges for forecasting and mitigating their impacts. Despite advancements, accurately predicting TC rapid intensification (RI) remains a challenge. Large eddy simulation (LES) allows for explicitly resolving the large eddies involved in TC turbulence, thus providing an avenue for studying the mechanisms behind their intensification and RI. LES of a full tropical cyclone is very computationally expensive and its accuracy will depend on both explicit and implicit dissipation within an atmospheric model. This dissertation presents two novel numerical methodologies with the potential to improve the efficiency …

Using Cache Files To Improve The Efficiency Of The Ucgretina Simulation Code, Blake Mcnulty

Physics and Astronomy Summer Fellows

UCGretina is a program that is used to simulate the Gretina gamma-ray tracking array used in experiments at the Facility for Rare Isotope Beams (FRIB) to collect gamma rays produced in reactions such as proton scattering. This is done by generating a simulation of the beam particles used and sending them into a simulated recreation that mimics the experimental environment. My job this summer was to tackle a major issue that the program had which is the speed or efficiency at which it ran. The main issue was that it would take days for certain simulation sets to run which …

Simulating Cross-Scale Solid-Fluid Interaction Phenomena, Jinyuan Liu

Dartmouth College Ph.D Dissertations

Solid-fluid interactions are ubiquitous in nature, and accurate simulation methods are essential for realistic animation, industrial design, and engineering analysis. Com- pared to large-scale coupling phenomena, simulating fine-scale interactions poses extra challenges due to factors such as surface tension, material wettability, and geometric complexity. In this thesis, we pursue novel methodologies to accurately model in- terfacial dynamics between surface-tension fluids and codimensional solids, involving capillary interactions, controllable wettability, and robust contact behaviors. Our ini- tial approach involves developing a novel three-way coupling method, which utilizes a thin liquid membrane, modelled as a simplicial mesh, to facilitate accurate momen- tum transfer, …

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 …

Development Of On-The-Fly Quasi-Steady State Approximation For Chemical Kinetics In Cfd, Abhinav Balamurugan

Doctoral Dissertations and Master's Theses

This study analyzes the feasibility of On-The-Fly Quasi-Steady-State Approximation (OTF-QSSA) application for solving chemical kinetics within Computational Fluid Dynamics (CFD) simulations, aiming to reduce the computational demand of detailed mechanisms. An algorithm that dynamically identifies and designates Quasi-Steady-State (QSS) species at specific grid locations and instances during the simulation was developed. With this information, our method pseudo-delays the advancement of concentrations for these QSS species—effectively setting their rate of concentration change to zero for a set number iteration before updating using the detailed mechanism and thereby omitting the computationally intensive processes typically required for their calculation during those skipped iteration. …