Numerical Analysis and Computation Commons^™

Integrating Gis With Interim Payment Valuation In Road Construction Projects: A Conceptual Framework, Abdulrahman I. Iro, Juma M. Matindana, Julian Ijumulana

Tanzania Journal of Engineering and Technology (TJET)

Abstract

Interim Payment Valuation is a critical process in road construction contract administration, yet conventional valuation practices remain heavily dependent on manual measurements, fragmented documentation, spreadsheets, and professional judgement. These limitations can affect measurement accuracy, transparency, traceability, and the timeliness of payment certification. Although Geographic Information Systems have increasingly been applied to construction planning, quantity measurement, progress monitoring, infrastructure management, and decision support, their integration with contractual and financial processes for interim payment valuation remains insufficiently explored. This study therefore develops a conceptual framework for integrating GIS with Interim Payment Valuation in road construction projects. A PRISMA-guided structured literature review …

Voltage-Related Power Quality Issues And Impacts On Distribution Networks With Sensitive Loads - A Review, Godwin Elinazi Mnkeni, Jackson Justo, Aviti Thadei Mushi, Bakari M. M. Mwinyiwiwa

Tanzania Journal of Engineering and Technology (TJET)

Voltage disturbances are the most important power quality (PQ) complications that customers and power utilities face in this smart era. The growing adoption of sophisticated electronic equipment and integration of renewable energy sources (RES) into power grids has increased the susceptibility of power distribution networks (PDNs) to voltage sags, swells, interruptions, flicker, and voltage imbalance. These disturbances, mainly caused by upstream faults, switching operations, and RES integration, compromise voltage PQ and system reliability. Consequently, they accelerate equipment degradation, increase electronic waste (e-waste), raise reactive power demand and maintenance costs, increase power losses, and impose substantial economic losses on customers and …

Thermo-Entropic Analysis Of Unsteady Mhd Nanofluid Couette Flow: A Classical Spectral Approach With An Exploratory Quantum Linear-Solver Application, Serai Israel Mosala, Oluwole Daniel Makinde, Azwinndini Muronga

Mathematical Modelling and Numerical Simulation with Applications

This study presents a thermo-entropic analysis of unsteady MHD nanofluid Couette flow, combining a classical bivariate spectral quasilinearisation (BI-SQLM) and Crank--Nicolson solution with an exploratory application of quantum linear solvers. An incompressible, electrically conducting nanofluid flows between parallel plates under partial slip and convective heat exchange; the discretised linear systems are additionally solved via the Harrow--Hassidim--Lloyd (HHL) algorithm and the Variational Quantum Linear Solver (VQLS). Results are presented for Pure Water, Cu-Water, and Al₂O₃-Water across seven parameter variations. The Hartmann number dominates velocity suppression and entropy amplification, while velocity slip reduces upper-wall entropy generation more than …

Effects Of Co-Rotating Cylinders And Boundary Conditions On Heat And Mass Transfer In A Trapezoidal Enclosure, Olalekan Adebayo Olayemi, Olalekan Tajudeen Popoola, Leke Thaddeus Oladimeji, Isaac Kayode Adegun, Taofiq Omoniyi Amoloye, Oluwatayo Babatope Ojo, Benjamin Ohida Daniel, Michael Olabode Ibiwoye

Mathematical Modelling and Numerical Simulation with Applications

This study investigates the effects of rotational speed $(\Omega)$, Richardson number $\left(Ri\right)$, and Lewis number $\left(Le\right)$ on heat and mass transfer (HMT) around two co-rotating cylinders symmetrically situated in a trapezoidal enclosure. In the double-diffusive mechanism, the bottom wall of the trapezium is at a high temperature $\left(T_h\right)$ and concentration $\left(C_h\right)$, while the top wall is at a low temperature $\left(T_c\right)$ and concentration $\left(C_c\right)$. The side walls are perfectly insulated with zero mass flux, while the buoyancy ratio $\left(Br\right)$ is fixed at 1.0.  Two thermal-mass boundary scenarios are considered. Scenario 1: Cylinders are thermally insulated with zero mass flux, and …

Exact Solution Structure And Dynamic Preservation In A Nonstandard Finite-Difference Formulation Of The Sir Model, Ramsha Shafqat, Wutiphol Sintunavarat

Mathematical Modelling and Numerical Simulation with Applications

A new nonstandard finite difference method is developed for discretizing the classical Susceptible--Infected--Removed (SIR) epidemic model. In the proposed approach, the first derivatives in the continuous SIR system are replaced by two carefully constructed nonstandard denominator functions, yielding a dynamically consistent discrete model with improved numerical stability while preserving consistency with the original continuous system. It is proved that the proposed scheme preserves the fundamental dynamical properties and qualitative behavior of the continuous SIR model, including positivity and the long-term dynamics of the population classes. Moreover, the resulting discrete system admits explicit analytical solutions, providing further insight into the model …

Numerical Solution Of The Hjb Equation Using A Hybrid Bernstein–Fractional Chebyshev Spectral Collocation Method, Alvian Alif Hidayatullah, Subchan Subchan, Sena Safarina, Tahiyatul Asfihani, Dilshod Mashrabovich Akhmedov

Mathematical Modelling and Numerical Simulation with Applications

The Hamilton–Jacobi–Bellman (HJB) equation is fundamental in stochastic optimal control but is often difficult to solve accurately because of its nonlinear and multidimensional structure. This paper develops a hybrid tensor-product spectral collocation method based on shifted Bernstein polynomials in time and fractional-order Chebyshev polynomials in the state variables. The proposed shifted Bernstein–fractional Chebyshev (SBFC) framework allows the time and state directions to be approximated using different basis characteristics, while the fractional parameters provide adjustable resolution in the state domain. A weighted approximation framework is established, including density and convergence results for the mixed approximation space. The state collocation points are …

Neural Network-Based Analysis Of Heroin Epidemic Models With Modified Fractional Operators, M. A. El-Shorbagy, Sedat Pak, Mati Ur Rahman, Hossam A. Nabwey

Mathematical Modelling and Numerical Simulation with Applications

Heroin and synthetic narcotic abuse have become a major global concern, posing challenges to individuals, families, and communities. Their widespread availability and low cost have intensified the crisis. This study investigates a heroin transmission model using the modified Atangana--Baleanu--Caputo (mABC) fractional operator, with emphasis on non-zero solutions. Series solutions are derived by combining the Laplace transform with the Adomian decomposition method to address nonlinear components. Qualitative analysis is conducted through fixed-point theory, while stability is assessed using the T-Picard method. Numerical simulations explore the effects of different fractional orders and transmission parameters on the system. The study incorporates a deep …

Pattern Formation And Diffusion-Driven Instability Of A Modified Lengyel-Epstein Model, Müge Meyvacı, Hüseyin Merdan, Yafia Radouane

Mathematical Modelling and Numerical Simulation with Applications

We consider a modified Lengyel-Epstein model that is covered by a system of reaction-diffusion equations and investigate the sufficient conditions on parameters of the model at which the diffusion-driven instability occurs. We perform numerical simulations to support and extend the theoretical findings. Analytical results and numerical simulations show that the modified Lengyel-Epstein system with positive initial and non-flux boundary conditions presents Turing patterns under some conditions on its parameters.

When Three Points Aren't Enough: Parameter Estimation In Newton's Law Of Cooling, Alberto A. Condori, Cara D. Brooks, Madeline R. Goldberg

CODEE Journal

We begin with a paradoxical three-point problem where the standard parameter estimation formula fails because temperature data must satisfy a concavity condition reflecting Newton's Law of Cooling's exponential structure. Although a closed-form solution for A exists for equally-spaced measurements, high sensitivity to error motivates the use of overdetermined systems with many measurements. The "profiling over A" technique transforms this three-parameter nonlinear problem into a sequence of simple linear regressions, providing computational efficiency and conceptual transparency. This approach can be generalized to many parameter estimation problems in science and engineering, making it a valuable tool for undergraduates interested in applied mathematics. …

Electronic Structure Discretization And Compression Using Diagonal Basis Sets, Casey Lee Dowdle

Dartmouth College Ph.D Dissertations

Numerically solving the electronic structure problem is a fundamentally difficult problem due to the exponential growth in the dimension of the Hilbert space as the system size increases. In order to solve problems at a chemically relevant accuracy, both the choice of basis set and numerical method are important factors that are intrinsically connected.

In this thesis, we study the discretization and resulting compression of electronic Hamiltonians using diagonal basis sets. A diagonal basis set approximately diagonalizes the matrix and tensor representations of the one- and two-body potentials. This can reduce storage, simplify matrix-vector products, and lower the complexity of …

Human-Centered Modeling Of Traffic As A Complex System, Poorendra P. Ramlall

Discovery Day - Daytona Beach

Traffic systems are driven not only by motion, but by interaction: vehicles influence one another, drivers continuously adapt to surrounding behaviour, and cognitive processes shape decisions that can propagate through the flow of traffic. Understanding these layered interactions is essential for improving traffic safety and for designing the next generation of intelligent, connected, and automated transportation systems. This PhD research develops a multiscale, data-driven framework for identifying, modelling, and ultimately interpreting interaction structure in traffic systems. The work first established an information-theoretic basis for this problem, demonstrating how information flow can uncover directional relationships in traffic dynamics and help infer …

Analyzing Fungal Growth Dynamics Under Different Environmental Conditions Using A Lotka–Volterra Competition System, Maria Ordonez, Fabrio Araujo

Discovery Day - Daytona Beach

Fungi play a critical role in ecosystems as decomposers that recycle nutrients and maintain environmental balance. Their populations are influenced by multiple environmental factors such as temperature, humidity, nutrient availability, and interactions with other organisms. In this project, the Lotka–Volterra model is used to analyze how competing fungal species interact and how these interactions influence population dynamics over time. By modeling two fungal populations competing for the same limited resources, the equations illustrate how environmental conditions and competition coefficients determine whether one species dominates; both species coexist, or one species becomes extinct. The model provides insight into how changes in …

Modeling Stellar Structure: Comparing Numerical Solutions Of The Lane-Emden Equation, Jasman Jasmanjot, Bailey Dale, Ryan Dickey

Discovery Day - Daytona Beach

The Lane-Emden equation is a differential equation that is often used in astrophysics to describe the distribution of the density inside a star, and by extension, its pressure distribution. Analytic solutions of the Lane-Emden equation can only be found at polytropic indices n = 0,1,5, matching certain physical conditions. For all other polytropic values, a numerical solution is needed. In this work, a comparison between two numerical schemes for solving the Lane-Emden equation is presented, namely the Euler method and the classical 4th order Runge-Kutta method. The accuracy of these models is first compared to the cases with known analytical …

Numerical Modeling Of Thermo-Poroelasticity Using Finite Element Method, Maya Mckean

Discovery Day - Daytona Beach

We propose a numerical method for solving and modeling thermo-poroelasticity problems using a finite element formulation. Thermo-poroelasticity models describe the coupled interaction between mechanical deformation, fluid flow, and heat transfer in specific materials or environments over time. These models demonstrate the evolution of displacement, pressure, and temperature; we compute these fields in this work using backward Euler time discretization and Enriched Galerkin finite element spatial discretization. For these computations we used FreeFEM, a partial differential equation solver that uses the finite element method, which produced our numerical results. We then compared these values with the expected analytical solution. This was …

Numerical Methods For Nonlinear Problems Using The Finite Element Method, Logan Price

Discovery Day - Daytona Beach

Numerical Methods for Nonlinear Problems Using the Finite Element Method is a computational mathematics capstone that builds and tests finite element method (FEM) workflows for nonlinear partial differential equations in FreeFEM++, with ParaView used for visualization. Two nonlinear model problems are used to demonstrate the approach. The first is a semilinear reaction-diffusion equation with a cubic nonlinearity. A manufactured solution is used so accuracy can be checked at a fixed final time, and refinement studies in both time step and mesh size are run while nonlinear iteration counts are tracked to show solver effort. The second problem is the steady …

Classification Of Sequential Factors In Aviation Accident Cause Prediction, Sophia Nasca, Addyson Wolfe

Discovery Day - Daytona Beach

Uncovering the root causes of aviation accidents is a critical component of improving aviation safety. Traditional approaches are largely reactive, relying on post-incident analysis rather than proactively identifying risk factors. This project addresses the need for proactive safety by using a multi-source dataset that integrates aviation accident records, weather conditions, and maintenance data extracted from investigative reports. The objective of this work is to move beyond predicting broad probable causes and instead model the sequence of contributing factors that lead to aviation incidents. Using the Swiss Cheese Model, the study will capture layered failures across operational, environmental, and maintenance domains. …

Quantifying The Effect Of Metallicity On Stellar Properties And Evolutionary Timescales, Jacob Becker

Discovery Day - Daytona Beach

Stellar age estimations derived from asteroseismology depend on stellar models that are sensitive to metallicity (Z). Variations in this parameter could alter the agreement between gyrochronological and asteroseismic ages, as well as main sequence lifetimes and observational properties such as effective temperature and luminosity. We test how much typical metallicity differences (±0.2 dex) affect main-sequence models of solar-type stars. Using MESA, evolutionary tracks are created for 1 M☉ stars at three metallicities (Z = 0.009, 0.014, 0.022), and it is measured how these changes shift the positions of the zero-age main sequence and the corresponding main sequence lifetimes. This work …

Analytical And Numerical Solutions For The Hydrogen Atom, Kassidy Myers

Discovery Day - Daytona Beach

The Schrödinger equation is the foundational equation of non-relativistic quantum mechanics. The hydrogen atom is the simplest system for solving this equation, as it consists of only one proton and one electron. In this project, we work on the Schrödinger equation that models the spherically symmetric states of the hydrogen atom that depend only on the radial coordinate. We simplified and nondimensionalized the radial equation and solved the resulting equation using a power series (Frobenius) method. This approach revealed the physically meaningful solutions and led to quantized energy levels. In addition to finding the analytical solution, we numerically solve the …

Optimization Of Engine, Jordan Reed, Dev Shah

Discovery Day - Daytona Beach

A matrix-based framework for modeling and optimizing fluid and gas in feed systems to pressurize for propulsion applications using advanced linear algebra techniques will be used in this project. The governing equations are derived from conservation of mass, momentum, and energy and are formulated in state space form. This enables the system to be expressed as a set of coupled linear differential equations. These equations are assembled into structured system matrices that show the interactions between pressure, flow rate, and component dynamics. This representation allows for numerical implementation and scalability to complex systems with multiple components. System behavior is analyzed …

Bridging Discrete And Continuous Systems: Fibonacci Sequences And Exponential Growth From Odes, Martyna Wojcik

Discovery Day - Daytona Beach

Bridging Discrete and Continuous Systems: Fibonacci Sequences and Exponential Growth from Ordinary Differential Equations       It has been observed that nature often exhibits specific patterns of growth and structure in biological systems and spiral formations. The Fibonacci sequence, defined as a discrete recursive sequence where each term is generated as the sum of the two preceding terms, “has been applied extensively to understand some natural phenomena” (Pakdemirli, 2023). In contrast, exponential growth describes a continuous process in which the rate of change of a quantity is proportional to its current value. Such behavior is modeled using differential equations that “produce solutions …

Numerical Modeling Of Badminton Shuttlecock Trajectories, Lola G. Torres, Cassandra Pumphrey, Jadyn Peterson, Domenic Barsotti

Discovery Day - Daytona Beach

The Trajectory of a badminton Shuttlecock can vary significantly when compared to a classic projectile motion, primarily due to aerodynamic drag. This project aims to model the flight of the shuttlecock using Newton's second law for gravitational and drag related forces, resulting in a nonlinear system of a first order differential equation. The given parameters include the shuttlecock mass, cross-sectional area, air density, as well as the drag coefficient, determining the overall magnitude of the drag force. The resulting initial value problem is solved numerically using a multitude of Runge_Kutta methods to compare the accuracy and stability across different computational …

Numerical Analysis Of The Sir Model For Predicting Disease Spread, Victoria Gaibor, Isabel Tejada, Kate Moore

Discovery Day - Daytona Beach

This project, Numerical Solutions of the SIR Model for Predicting Disease Spread, investigates the application of numerical methods to analyze the dynamics of infectious diseases using the classical Susceptible–Infected–Recovered (SIR) model. The SIR model, a system of nonlinear ordinary differential equations, is widely used to describe how diseases such as COVID-19 propagate through a population. The primary objective of this study is to solve the SIR initial value problem using multiple numerical techniques, including Euler’s method, Runge–Kutta methods, and multistep methods, and to compare their accuracy and efficiency. The model is implemented using given initial conditions and parameters, and additional …

Numerical Investigation Of The Nonlinear Simple Pendulum And The Dependence Of Oscillation Period On Initial Angle, Kelly Wold, Aidan Hart, Patrick Gilliam

Discovery Day - Daytona Beach

Numerical Investigation of the Nonlinear Simple Pendulum and the Dependence of Oscillation Period on Initial Angle examines how the oscillation period of a simple pendulum varies with initial angular displacement and evaluates the accuracy of numerical methods in capturing this behavior. In classical treatments, the small-angle approximation simplifies the governing differential equation and predicts a constant period independent of amplitude; however, this assumption breaks down for larger angles, where the system exhibits nonlinear dynamics. The objective of this project is to model the full nonlinear equation of motion and quantify how the period depends on initial conditions. To achieve this, …

On A Classification Of Goy/Sabra-Type Shell Models Of Turbulence, Austin Ho

Theses and Dissertations

Since their inception in the 1970s, shell models have proved to be a fruitful testbed for studying theories of turbulence. This thesis considers variations of two classical shell models: the Gledzer–Ohkitani–Yamada (GOY) model and the Sabra model. Up to conjugacy classes, we classify all possible GOY-type and Sabra-type models. We numerically approximate the scaling exponents for the newly identified models and develop their mathematical well-posedness.

From Homogeneous To Heterogeneous Adaptive Bounded-Confidence Opinion Dynamics On Networks, Sally Hafez, Fatma R. Farag, Amira S. N. Tawadros

Northeast Journal of Complex Systems (NEJCS)

Adaptive bounded-confidence models (ABCMs) elucidate the coevolution of agent states and network structure via local interactions and rewiring mechanisms. Traditional formulations assume uniform interaction parameters, leading to distinct regime shifts encompassing fragmentation, polarization, and consensus. A symmetric heterogeneous extension of the adaptive bounded-confidence model is introduced, in which interaction parameters vary according to whether agents belong to the same or different groups. The model retains the original update and rewiring protocols but integrates within-group and between-group confidence bounds alongside tolerance thresholds. Initially, the classic homogeneous model is replicated to establish a reference point. Subsequently, the heterogeneous extension is assessed under …

Criticality In A Heterogeneous Neutron Transport Rod Model, Samuel Kaleb Crowford

All Graduate Reports and Creative Projects, Fall 2023 to Present

This work studies the stochastic behavior of neutron populations in a one-dimensional rod model using Monte Carlo simulation. The first part of this project reproduces the computational results of Dumonteil, Horton, Kyprianou, and Zoia (2025) by independently implementing the Monte Carlo algorithm described in their article, with the asymptotic behavior of the first moment analyzed in relation to the dominant eigenvalue and adjoint eigenfunction of the neutron transport operator. The model is then extended to a heterogeneous setting by introducing a central region where fission is suppressed. A global expectation over initial positions and directions is used to estimate the …

Sustainable Consumer Behavior Modeling: A Complex-Systems Approach To Neuromarketing, Preethi Nanjundan, Nupoor Sanjay Bhute, Ragini Topre, Lijo Thomas

Northeast Journal of Complex Systems (NEJCS)

This study examines the application of complex-systems modeling to neuromarketing data for gaining deeper insights into the mechanisms underlying sustainable consumer behavior. It investigates how neural and biometric responses, interpreted through a systems-based perspective, can uncover dynamic interactions, feedback mechanisms, and emergent behavioral patterns influencing sustainable purchase decisions. The research explores the impact of sustainability-oriented marketing stimuli on long-term behavioral intentions by emphasizing the interconnected roles of cognitive processing, emotional engagement, implicit associations, and collective consumer dynamics. Through simulation-based modeling and structural analysis, the study demonstrates how subconscious neural responses and affective mechanisms mediate the relationship between marketing interventions, consumer …

Provocable Forgiveness In Noisy Brand--Consumer Systems:An Agent-Based Study Of Repeated Interaction, Tahere Ahmadiyan, Hamidreza Navidi, Behbod Keshavarzi

Northeast Journal of Complex Systems (NEJCS)

Repeated brand--consumer exchange is often treated as a managerial problem of loyalty, recovery, and trust. It can also be read as a small complex system: many local decisions about cooperation, retaliation, and forgiveness accumulate into market-level selection. This study uses that perspective to examine which relational rules survive when communication is imperfect. Eight canonical Iterated Prisoner's Dilemma strategies are translated into marketing archetypes and evaluated through round-robin tournaments, a six-level noise sweep, proportional-fitness ecological dynamics, and finite-population Moran invasion tests. The tournament leaderboard is calculated without same-strategy self-play, so that reported payoffs reflect inter-archetype competition rather than homogeneous self-coordination. At …

Pseudodifferential Absorbing Boundary Conditions For Waves, Lauren Taylor

Mechanical Engineering Theses

Absorbing boundary conditions (ABCs) are required to truncate the computational domain when performing Finite Element Analyses of exterior acoustic problems where physical domain is unbounded. ABC is applied on a fictitious boundary containing the scatterer and ideally allows outgoing waves to leave without non-physical reflections. Preventing artificial reflections is essential to benefit from the accuracy of the numerical method used. Otherwise, ABC acts as a reflective surface, and artificial reflections distort the solution in the entire domain which is not recoverable by any type of refinement. Pseudodifferential ABCs were used to describe the Dirichlet-to-Neumann map which maps known boundary values …

Cdt-1d Cnn Integration With Simpson-Sobolev Regularization For High-Frequency Options Trading: With Fem-Based Heston Option Pricing, Daniel M. Margolis, Johannes Tausch, Arthur K. Selender

Mathematics Theses and Dissertations

This dissertation presents a computational framework for high-frequency options trading that combines Cross-Data-Type 1-D Convolutional Neural Networks (CDT-1D CNN) with Simpson-Sobolev regularization for directional prediction, and finite element methods (FEM) for realistic option pricing during backtesting. The core innovation lies in developing a mathematically rigorous regularization approach that maintains the adaptability of modern deep learning while enabling accurate evaluation through stochastic volatility models. The primary contribution is the Simpson-Sobolev regularization scheme, which extends traditional Sobolev regularization by incorporating Simpson’s rule for numerical integration. This approach achieves higher-order accuracy in approximating the Sobolev norms that control function smoothness. Simpson’s rule attains …