Titles and Abstracts
10th Annual UMKC Math and Stat Research Day
Division of Computing, Analytics & Mathematics,
School of Science and Engineering
University of Missouri - Kansas City
A Decoupled Scheme for a Single-Phase Model in Ferrohydrodynamics
Lucas Delibas
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
In this research, we study a numerical approximation of a simplified PDE model for single phase ferrofluid flow, which consists of the Navier–Stokes equations coupled with the Poisson equation. Using existing spatial and temporal discretization techniques for the weak formulation, we construct a fully discretized coupled and decoupled finite element schemes with Newton linearization to solve the nonlinear and coupled multi-physics PDE system in Ferrohydrodynamics (FHD). Using MATLAB code packages, we perform numerical simulations for both coupled and decoupled schemes to verify convergence of the numerical solutions. Ferrofluids are colloidal solutions made of ferromagnetic nanoparticles suspended in a liquid, also known magnetic liquid. Such fluid can be controlled directly by the application of a magnetic field, its magnetization properties vanish entirely when the applied magnetic field is removed. This magnetizing property of ferrofluids widely utilized in many scientific and engineering fields requiring precise control, including drug delivery, electronic, biomedical engineering, nanotechnologies, assembly of nanoparticles, fluid transport and control, in optics, electronic and devices, etc.
Analyzing Biosecurity Adherence and Antibiotic Resistance Dynamics in Animals and Humans through Mathematical Modeling
Arash Arjmand1,
Joint work with Majid Bani-Yaghoub1, Kiel Corkran1, Pranav S. Pandit2, and Sharif Aly2
1Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
2Department of Population Health and Reproduction, University of California, Davis School of Veterinary Medicine, Davis, California, United States of America
In addressing the growing challenge of antimicrobial resistance, our study investigates the impact of biosecurity measures on Salmonella transmission among dairy cattle. Through a framework of ODE, we examine how different levels of biosecurity compliance (Excellent, Good, Marginal and Low) affect the spread of pathogens within animal populations. Our findings highlight that while individual biosecurity actions are important, no single method completely stops outbreaks. Instead, combining measures, especially those reducing direct animal-to-animal transmission, significantly lowers infection rates. Our study underscores the need for comprehensive biosecurity to safeguard livestock and community health.
Approximating the periodic solutions of an Hes1-mRNA model
Mohammed Alanazi
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
Hes1 (Hairy and Enhancer of Split 1) is a gene that plays a crucial role in embryonic development and cellular differentiation. mRNA plays a crucial role in the process of gene expression. The importance of Hes1 mRNA interaction lies in its role as a regulator of cell fate determination and differentiation during development. This research outlines the preliminary steps towards approximating periodic solutions for a model governing the dynamics of Hes1 mRNA interactions. By linearizing this model, we have identified a critical value (i.e., a bifurcation point) resulting in loss of linear stability. This bifurcation gives rise to periodic solutions as verified via numerical simulations. The subsequent phase of this research will apply the Poincare method to drive an approximation of the periodic solutions
An Extended Model to Study the Effects of Healthy Cell Growth Behaviors on the Dynamics of Vascular Cancer
Priscilla Owusu Sekyere
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
This study investigates the effectiveness of chemotherapy and anti-angiogenic therapy in controlling the growth of cancer cells within the context of an extended cancer model. We extend a cancer model consisting of a system of five ordinary differential equations that simulates the interactions between normal cells, cancer cells, endothelial cells, chemotherapy agent and anti-angiogenic agent in tumor growth. Numerical simulations with varying initial conditions and combinations are used to assess the impact of healthy cell growth behaviors on the dynamics of vascular cancer. The numerical method used is the Runge-Kutta method. It was used through Matlab ODE 45. These findings help guide the development of targeted treatments for vascular cancer.
Comparative Analysis of FTCS and FECD for Numerical Simulation of SIR RD Model
Adriana Martínez Cappello
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
This study aims to compare the computational speed of two explicit numerical methods, Forward Time Centered Scheme (FTCS) and Forward Euler Centered Difference (FECD), for simulating the Susceptible-Infectious-Removed (SIR) Reaction-Diffusion (RD) model. The SIR model is a classical tool in epidemiology used to understand disease spread within a population. Both methods utilize finite difference schemes to approximate spatial derivatives by discretizing space and time. We evaluate the efficiency of FTCS and FECD by comparing their computational speed and scalability performance. Results demonstrate that FECD exhibits greater computational efficiency than FTCS using the above metrics. Further, increasing the disease transmission rate (β) leads to a rise in the infected population concentration, which spreads to the susceptible population, regardless of the chosen method.
Probabilistic Learning on Manifold for Unsteady Fluid Dynamics Simulations under Uncertainties
Long Dang
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
In the field of computational fluid dynamics (CFD), simulating unsteady fluid flows, particularly those governed by the Navier-Stokes equations, presents significant challenges due to the complexity of nonlinear convection, diffusion, and pressure velocity coupling.
We solve the unsteady 2D Navier Stokes equations for incompressible fluid started with fixing Reynold's number using Chorin’s Splitting Method. This method separates the pressure-velocity coupling problem into distinct, more manageable sub-problems. The stability, consistency, and convergence conditions will be used to analyze this method.
This research aims to employ Probabilistic Learning on Manifold(PLoM) to explore the fluid dynamics under uncertainty, where the Reynold number is treated as the random control parameter. The Karhunen-Loeve (KL) expansion is implemented to reduce the dimensionality of the obtained velocity to capture essential features. Then PLoM is used to generate additional realization (sample) to robust model predictions and uncertainty quantification
Investigation of a numerical approach replacement for domain discretizing in wall bounded laminar flows
Milad Mohammadi
Division of Energy, Matter and Systems, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
Laminar flows carry on simple physical features compared to turbulent flows. One of the major applications of these flows is in medical fields (Cerebral spinal fluid and vessel diseases). Since the
governing equations of the laminar flow is much simpler than turbulent flows, maybe the reasonable approach to discretize the equations would be a method with computational costs and numerical challenges less than original finite volume method (FVM). In this study, we are investigating a possible replacement of FVM with FDM to get the desirable results with less computational cost and less
numerical challenges. We simulated 18 cases of channel flow with different step height. For grid generation, the number of grid points for both approaches considered to be the same in order to make a reasonable comparison. Simulations carried out using Ansys Fluent software and MATLAB software for the FVM and FDM methods, respectively. The velocity fields from the FDM and FVM methods for
18 cases were compared. The results show that with FDM approach, the computational costs and simulation challenges (such as using different software for geometry and grid generation) reduces,
compared to FVM method. Although, there is still the weakness of FDM method regarding complex geometries to deal with, which diminishes its range of applicability.
Evaluating PDE-Refiner: Comparing a Neural PDE Solver Against Traditional Numeric Methods for
Long Rollout of the 1D Kuramoto-Sivashinsky Equation
David Wagner
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
In this paper we aim to reproduce and further quantify the results comparing the PDE-Refiner neural PDE solver to the traditional finite difference method. First a comparison with the results of the original paper will be made to quantify the reproduction of the experiment. The rollouts between the two methods will be compared visually with their outputs alongside a diff of the outputs. Then the compute tradeoffs between the methods, the stability and consistency of
rollouts, as well as long term accuracy of rollouts will be assessed. In evaluating the strengths and weaknesses of the neural PDE solver we will identify areas for future research in improving performance of neural PDE solvers for long rollouts.
Social and Economic Contributions to COVID-19 Prevalence in Nursing Homes
Julia Pluta
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
The COVID-19 pandemic struck long-term care facilities hard in the first year of the pandemic, prior to the widespread availability of vaccines. In the general population, COVID-19 generally hit poorer communities and communities of color harder than affluent communities or communities with a higher percentage of white residents. The Regional Planning Association and American Planning Association use a concept of “Megaregions” to describe large, multicity areas connected through common transportation and commercial hubs. The placement of Kansas City in those Megaregions is inconsistent.
In prior work, clusters of COVID risk roughly followed Megaregion patterns in the early phases of the pandemic. However, national-scale analysis proved unwieldy. This study was meant to serve as a proof-of-concept for analyzing the pandemic based on Megaregion grouping, while trying to use Kansas City to determine which Megaregion model worked best. At the same time, this was also serving as a small-scale test to determine whether Social Determinants of Health could be combined with factors unique to a given nursing home to determine whether they could be a suitable basis for predicting pandemic risk.
Numerical Explorations of the Competitive Reaction-Diffusion Lotka-Volterra Model with Seasonal Growth Rates
Most Shewly Aktar
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
This research explores the spatio-temporal dynamics of competing species within the framework of the Lotka-Volterra model. In the absence of diffusion, the parameter space for the classical model is divided into four regions of coexistence, competitive exclusion or founder control. The main objective of this research is to determine whether the presence of diffusion can alter any of these regions. We use Matlab pdepe to numerically explore the solutions of the PDE model with Neuman boundary conditions. The numerical simulations suggest that diffusion has stabilizing effects and the outcomes of the PDE model are consistent with the corresponding ODE model. Also, the seasonal growth rates did not have any effect on the convergence of the PDE solution to the equilibria
Comparative analysis of ODE solutions by using Fast and well-conditioned Spectral and Numerical methods
Farzana Sultana Rafi
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
Spectral methods are numerical techniques utilized for solving differential equations with high accuracy and efficiency. Fast and well-conditioned spectral methods are particularly notable for their ability to swiftly converge to accurate solutions.
In this research, I apply fast and well-conditioned spectral methods to solve ordinary differential equation (ODE) problems. Specifically, here I employ spectral methods with increasing Chebyshev nodes to approximate the solutions. The accuracy of these approximations is compared with both exact solutions and solutions obtained using the ODE45.
My research observes how spectral methods perform when increasing the number of Chebyshev nodes. For 1st example as the number of Chebyshev nodes increases from N=1 to N=3, the solutions provided by the spectral method gradually converge towards the exact solution. However, at N=4 and N=5, differences arise between the spectral method's solution and the exact one, indicating potential limitations of the approach for higher node counts. Nevertheless, for the 2nd ODE problem, it observes that increasing the number of Chebyshev nodes to N=10, N=20, and N=30 leads to significantly improved solutions compared to lower node counts.
For one ODE, the results vary at N=4 and 5, while for another, the spectral method improves as Chebyshev node count increases. Further research is needed to confirm if this trend extends to other ODEs.
Converging of PDE Solutions to Travelling Wave Trains of Competitive Lotka-Volterra Model
Barsha Saha
Division of Computing, Analytics & Mathematics, School of Science and Engineering, University of Missouri - Kansas City, MO, USA
Lotka-Volterra (LV) Models have been used to study population dynamics of interacting species. The ODE version of Competitive Lotka-Volterra (CLV) models exhibits several periodic solutions. Three stable periodic limit cycles have been identified within Zeeman class of 29. The main focus of this study is to numerically explore how the solutions of the corresponding PDE CLV model behave with Neumann boundary conditions near the identified limit cycles. The numerical simulations suggest that the PDE solutions converge to the traveling wave trains of the model. This study confirms the previous theoretical results.