Open Access. Powered by Scholars. Published by Universities.®

Applied Mathematics Commons

Open Access. Powered by Scholars. Published by Universities.®

Faculty Publications

Discipline
Institution
Keyword
Publication Year

Articles 1 - 30 of 63

Full-Text Articles in Applied Mathematics

Analysis And Machine Learning Adaptation Of A Cognitive Model For Human Memory, Trevor Cross, Aihua W. Wood Jun 2026

Analysis And Machine Learning Adaptation Of A Cognitive Model For Human Memory, Trevor Cross, Aihua W. Wood

Faculty Publications

In this paper, we use the Duolingo SLAM dataset to analyze several cognitive models of second language acquisition and develop new approaches for enhanced performance. In particular, we consider the Predictive Performance Equation and some of its underlying power laws. Leveraging insights from machine learning, we develop simple one-feature models as building blocks for combined models that match or in certain cases outperform the existing models at much reduced computational cost. In addition, a neural network with one fully connected hidden layer is constructed that outperforms all other models on sufficiently large datasets.


Production Box Cost Estimating Relationships For Dod Avionics, Carla J. Cisneros, Edward D. White, Brandon M. Lucas, Jonathan D. Ritschel, Robert D. Fass, Shawn M. Valentine Oct 2024

Production Box Cost Estimating Relationships For Dod Avionics, Carla J. Cisneros, Edward D. White, Brandon M. Lucas, Jonathan D. Ritschel, Robert D. Fass, Shawn M. Valentine

Faculty Publications

The authors use historical information obtained from the Cost Assessment Data Enterprise to estimate recurring production unit cost for DoD avionics via cost estimating relationships (CERs). The specific modeled responses include mean unit cost, median unit cost, and the 100th production unit cost (T100) utilizing learning curve theory. For T100, the authors adopt both a multiplicative and an additive error for CER comparison. Recommended CERs consist of the mean unit cost and the T100 utilizing a multiplicative error. Moreover, results reveal that weight has a significant effect on cost as well as a potential underaccounting of real price change or …


Malware Classification Through Abstract Syntax Trees And L-Moments, Anthony J. Rose, Christine M. Schubert Kabban, Scott R. Graham, Wayne C. Henry, Christopher M. Rondeau Sep 2024

Malware Classification Through Abstract Syntax Trees And L-Moments, Anthony J. Rose, Christine M. Schubert Kabban, Scott R. Graham, Wayne C. Henry, Christopher M. Rondeau

Faculty Publications

The ongoing evolution of malware presents a formidable challenge to cybersecurity: identifying unknown threats. Traditional detection methods, such as signatures and various forms of static analysis, inherently lag behind these evolving threats. This research introduces a novel approach to malware detection by leveraging the robust statistical capabilities of L-moments and the structural insights provided by Abstract Syntax Trees (ASTs) and applying them to PowerShell. L-moments, recognized for their resilience to outliers and adaptability to diverse distributional shapes, are extracted from network analysis measures like degree centrality, betweenness centrality, and closeness centrality of ASTs. These measures provide a detailed structural representation …


Deterministic Global 3d Fractal Cloud Model For Synthetic Scene Generation, Aaron M. Schinder, Shannon R. Young, Bryan J. Steward, Michael L. Dexter, Andrew Kondrath, Stephen Hinton, Ricardo Davila May 2024

Deterministic Global 3d Fractal Cloud Model For Synthetic Scene Generation, Aaron M. Schinder, Shannon R. Young, Bryan J. Steward, Michael L. Dexter, Andrew Kondrath, Stephen Hinton, Ricardo Davila

Faculty Publications

This paper describes the creation of a fast, deterministic, 3D fractal cloud renderer for the AFIT Sensor and Scene Emulation Tool (ASSET). The renderer generates 3D clouds by ray marching through a volume and sampling the level-set of a fractal function. The fractal function is distorted by a displacement map, which is generated using horizontal wind data from a Global Forecast System (GFS) weather file. The vertical windspeed and relative humidity are used to mask the creation of clouds to match realistic large-scale weather patterns over the Earth. Small-scale detail is provided by the fractal functions which are tuned to …


Association Between Short Birth Spacing And Child Malnutrition In Bangladesh: A Propensity Score Matching Approach, Foyez Ahmmed, Nahid Hasan, Faruk Hossain, Tareq Ferdous Khan Mar 2024

Association Between Short Birth Spacing And Child Malnutrition In Bangladesh: A Propensity Score Matching Approach, Foyez Ahmmed, Nahid Hasan, Faruk Hossain, Tareq Ferdous Khan

Faculty Publications

Objectives: This study aimed to explore the effects of short birth spacing (SBS), which is defined as a period of less than 33 months between two successive births, on multiple concurrent forms of child malnutrition (MCFCM) and at least one form of child malnutrition (ALOFCM) using propensity score matching (PSM). Methods: This study used data extracted from the 2017- 18 Bangladesh Demographic and Health Survey. PSM with four different distance functions, including logistic regression, classification and regression tree, single hidden layer neural network and random forest, were performed to evaluate the effects of SBS on MCFCM and ALOFCM. We also …


Training Data Augmentation With Data Distilled By Principal Component Analysis, Nikolay Metodiev Sirakov, Tahsin Shahnewaz, Arie Nakhmmani Jan 2024

Training Data Augmentation With Data Distilled By Principal Component Analysis, Nikolay Metodiev Sirakov, Tahsin Shahnewaz, Arie Nakhmmani

Faculty Publications

This work develops a new method for vector data augmentation. The proposed method applies principal component analysis (PCA), determines the eigenvectors of a set of training vectors for a machine learning (ML) method and uses them to generate the distilled vectors. The training and PCA-distilled vectors have the same dimension. The user chooses the number of vectors to be distilled and augmented to the set of training vectors. A statistical approach determines the lowest number of vectors to be distilled such that when augmented to the original vectors, the extended set trains an ML classifier to achieve a required accuracy. …


Understanding The Implications Of Underreporting, Vaccine Efficiency And Social Behavior On The Post-Pandemic Spread Using Physics Informed Neural Networks: A Case Study Of China, Samiran Ghosh, Alonso Ogueda-Oliva, Aditi Ghosh, Malay Banerjee Nov 2023

Understanding The Implications Of Underreporting, Vaccine Efficiency And Social Behavior On The Post-Pandemic Spread Using Physics Informed Neural Networks: A Case Study Of China, Samiran Ghosh, Alonso Ogueda-Oliva, Aditi Ghosh, Malay Banerjee

Faculty Publications

In late 2019, the emergence of COVID-19 in Wuhan, China, led to the implementation of stringent measures forming the zero-COVID policy aimed at eliminating transmission. Zero COVID policy basically aimed at completely eliminating the transmission of COVID-19. However, the relaxation of this policy in late 2022 reportedly resulted in a rapid surge of COVID-19 cases. The aim of this work is to investigate the factors contributing to this outbreak using a new SEIR-type epidemic model with time-dependent level of immunity. Our model incorporates a time-dependent level of immunity considering vaccine doses administered and time-post-vaccination dependent vaccine efficacy. We find that …


Accurate Covariance Estimation For Pose Data From Iterative Closest Point Algorithm, Rick H. Yuan, Clark N. Taylor, Scott L. Nykl Jul 2023

Accurate Covariance Estimation For Pose Data From Iterative Closest Point Algorithm, Rick H. Yuan, Clark N. Taylor, Scott L. Nykl

Faculty Publications

One of the fundamental problems of robotics and navigation is the estimation of the relative pose of an external object with respect to the observer. A common method for computing the relative pose is the iterative closest point (ICP) algorithm, where a reference point cloud of a known object is registered against a sensed point cloud to determine relative pose. To use this computed pose information in downstream processing algorithms, it is necessary to estimate the uncertainty of the ICP output, typically represented as a covariance matrix. In this paper, a novel method for estimating uncertainty from sensed data is …


A Predator-Prey Biological Model With Combined Birth Rates, Self-Limitation And Competition Terms, Joon Hyuk Kang, Lucinda Ford Apr 2023

A Predator-Prey Biological Model With Combined Birth Rates, Self-Limitation And Competition Terms, Joon Hyuk Kang, Lucinda Ford

Faculty Publications

The purpose of this paper is to give sufficient conditions for the existence and uniqueness of positive solutions to a rather general type of elliptic system of the Dirichlet problem on a bounded domain Ω in Rn. Also considered are the effects of perturbations on the coexistence state and uniqueness. The techniques used in this paper are super-sub solutions method, eigenvalues of operators, maximum principles, spectrum estimates, inverse function theory, and general elliptic theory. The arguments also rely on some detailed properties for the solution of logistic equations. These results yield an algebraically computable criterion for the positive …


Evolution Of Coronal Magnetic Field Parameters During X5.4 Solar Flare, Seth H. Garland, Benjamin F. Akers, Vasyl B. Yurchyshyn, Robert D. Loper, Daniel J. Emmons Mar 2023

Evolution Of Coronal Magnetic Field Parameters During X5.4 Solar Flare, Seth H. Garland, Benjamin F. Akers, Vasyl B. Yurchyshyn, Robert D. Loper, Daniel J. Emmons

Faculty Publications

The coronal magnetic field over NOAA Active Region 11,429 during a X5.4 solar flare on 7 March 2012 is modeled using optimization based Non-Linear Force-Free Field extrapolation. Specifically, 3D magnetic fields were modeled for 11 timesteps using the 12-min cadence Solar Dynamics Observatory (SDO) Helioseismic and Magnetic Imager photospheric vector magnetic field data, spanning a time period of 1 hour before through 1 hour after the start of the flare. Using the modeled coronal magnetic field data, seven different magnetic field parameters were calculated for 3 separate regions: areas with surface |Bz| ≥ 300 G, areas of flare brightening seen …


A Mathematical Modeling Technique To Understand The Role Of Decoy Receptors In Ligand‑Receptor Interaction, Subrata Dey, Aditi Ghosh, Malay Banerjee Jan 2023

A Mathematical Modeling Technique To Understand The Role Of Decoy Receptors In Ligand‑Receptor Interaction, Subrata Dey, Aditi Ghosh, Malay Banerjee

Faculty Publications

The ligand‑receptor interaction is fundamental to many cellular processes in eukaryotic cells such as cell migration, proliferation, adhesion, signaling and so on. Cell migration is a central process in the development of organisms. Receptor induced chemo‑tactic sensitivity plays an important role in cell migration. However, recently some receptors identified as decoy receptors, obstruct some mechanisms of certain regular receptors. DcR3 is one such important decoy receptor, generally found in glioma cell, RCC cell and many various malignant cells which obstruct some mechanism including apoptosis cell‑signaling, cell inflammation, cell migration. The goal of our work is to mathematically formulate ligand‑receptor interaction …


Uniqueness Of Steady State Positive Solutions To A General Elliptic System With Dirichlet Boundary Conditions, Joon Hyuk Kang Dec 2022

Uniqueness Of Steady State Positive Solutions To A General Elliptic System With Dirichlet Boundary Conditions, Joon Hyuk Kang

Faculty Publications

The purpose of this paper is to give conditions for the uniqueness of positive solution to a rather general type of elliptic system of the Dirichlet problem on a bounded domain Ω in Rn. Also considered are the effects of perturbations on the coexistence state and uniqueness.


Training Rbf Neural Networks For The Solution Of Elliptic Boundary Value Problems, Andreas Karageorghis, C.S. Chen Nov 2022

Training Rbf Neural Networks For The Solution Of Elliptic Boundary Value Problems, Andreas Karageorghis, C.S. Chen

Faculty Publications

We propose a radial basis function (RBF) neural network method for solving two- and three–dimensional second and fourth order elliptic boundary value problems (BVPs). The neural network in question is trained by minimizing a nonlinear least squares functional, thus determining the optimal values of the various RBF parameters involved. The functional minimization is carried out using standard MATLAB® software efficiently. Several numerical experiments are presented to demonstrate the efficacy of the proposed method.


Obstructions To Shake Sliceness For Links, Anthony Bosman Apr 2022

Obstructions To Shake Sliceness For Links, Anthony Bosman

Faculty Publications

Shake slice generalizes the notion of a slice link, naturally extending the notion of shake slice knots to links. There is also a relative version, shake concordance, that generalizes link concordance. We show that if two links are shake concordant, then their zero surgery manifolds are homology cobordant. Then we give several obstructions to a link being shake slice; for instance, the Arf invariants vanish for both the link and each component. Finally we show that a shake slice link bounds disjoint disks in a homology 4-ball and hence each component is algebraically slice.


Image Classification Based On Sparse Representation In The Quaternion Wavelet Domain, Long H. Ngo, Nikolay Metodiev Sirakov, Marie Luong, Emmanuel Viennet Mar 2022

Image Classification Based On Sparse Representation In The Quaternion Wavelet Domain, Long H. Ngo, Nikolay Metodiev Sirakov, Marie Luong, Emmanuel Viennet

Faculty Publications

In this study, we propose a novel sparse representation learning method in the Quaternion Wavelet (QW) domain for multi-class image classification. The proposed method takes advantages from: i) the QW decomposition, which promotes sparsity and provides structural information about the image data while allowing approximate shift-invariance, to extract meaningful features from low-frequency QW subbands, ii) the dimensionality reduction method using Principal Component Analysis (PCA) for reducing the complexity of the problem, and iii) the sparse representation of the generated QW features to efficiently learn and capture the meaningful and compact information of this data. After the QW decomposition, the features …


Robust Error Estimation Based On Factor-Graph Models For Non-Line-Of-Sight Localization, O. Arda Vanli, Clark N. Taylor Jan 2022

Robust Error Estimation Based On Factor-Graph Models For Non-Line-Of-Sight Localization, O. Arda Vanli, Clark N. Taylor

Faculty Publications

This paper presents a method to estimate the covariances of the inputs in a factor-graph formulation for localization under non-line-of-sight conditions. A general solution based on covariance estimation and M-estimators in linear regression problems, is presented that is shown to give unbiased estimators of multiple variances and are robust against outliers. An iteratively re-weighted least squares algorithm is proposed to jointly compute the proposed variance estimators and the state estimates for the nonlinear factor graph optimization. The efficacy of the method is illustrated in a simulation study using a robot localization problem under various process and measurement models and measurement …


Survivals Of Two Cooperating Species Of Animals, Joon Hyuk Kang Dec 2021

Survivals Of Two Cooperating Species Of Animals, Joon Hyuk Kang

Faculty Publications

The purpose of this paper is to give conditions for the existence and uniqueness of positive solution to a rather general type of elliptic system of the Dirichlet problem on a bounded domain Ω" role="presentation" style="box-sizing: border-box; margin: 0px; padding: 0px; display: inline-block; line-height: normal; font-size: 16.2px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; border: 0px; position: relative;">Ω in Rn" role="presentation" style="box-sizing: border-box; margin: 0px; padding: 0px; display: inline-block; line-height: normal; font-size: 16.2px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: …


Extending Critical Infrastructure Element Longevity Using Constellation-Based Id Verification, Christopher M. Rondeau, Michael A. Temple, J. Addison Betances, Christine M. Schubert Kabban Jan 2021

Extending Critical Infrastructure Element Longevity Using Constellation-Based Id Verification, Christopher M. Rondeau, Michael A. Temple, J. Addison Betances, Christine M. Schubert Kabban

Faculty Publications

This work supports a technical cradle-to-grave protection strategy aimed at extending the useful lifespan of Critical Infrastructure (CI) elements. This is done by improving mid-life operational protection measures through integration of reliable physical (PHY) layer security mechanisms. The goal is to improve existing protection that is heavily reliant on higher-layer mechanisms that are commonly targeted by cyberattack. Relative to prior device ID discrimination works, results herein reinforce the exploitability of constellation-based PHY layer features and the ability for those features to be practically implemented to enhance CI security. Prior work is extended by formalizing a device ID verification process that …


A Learning Curve Model Accounting For The Flattening Effect In Production Cycles, Evan R. Boone, John J. Elshaw, Clay M. Koschnick, Jonathan D. Ritschel, Adedeji B. Badiru Jan 2021

A Learning Curve Model Accounting For The Flattening Effect In Production Cycles, Evan R. Boone, John J. Elshaw, Clay M. Koschnick, Jonathan D. Ritschel, Adedeji B. Badiru

Faculty Publications

We investigate production cost estimates to identify and model modifications to a prescribed learning curve. Our new model examines the learning rate as a decreasing function over time as opposed to a constant rate that is frequently used. The purpose of this research is to determine whether a new learning curve model could be implemented to reduce the error in cost estimates for production processes. A new model was created that mathematically allows for a “flattening effect,” which typically occurs later in the production process. This model was then compared to Wright’s learning curve, which is a popular method used …


Modeling And Simulation Techniques Used In High Strain Rate Projectile Impact, Derek G. Spear, Anthony N. Palazotto, Ryan A. Kemnitz Jan 2021

Modeling And Simulation Techniques Used In High Strain Rate Projectile Impact, Derek G. Spear, Anthony N. Palazotto, Ryan A. Kemnitz

Faculty Publications

A series of computational models and simulations were conducted for determining the dynamic responses of a solid metal projectile impacting a target under a prescribed high strain rate loading scenario in three-dimensional space. The focus of this study was placed on two different modeling techniques within finite element analysis available in the Abaqus software suite. The first analysis technique relied heavily on more traditional Lagrangian analysis methods utilizing a fixed mesh, while still taking advantage of the finite difference integration present under the explicit analysis approach. A symmetry reduced model using the Lagrangian coordinate system was also developed for comparison …


On Subdiagonal Rational Pade Approximations And The Brenner-Thomee Approximation Theorem For Operator Semigroups, Frank Neubrander, Koray Ozer, Lee Windsperger Dec 2020

On Subdiagonal Rational Pade Approximations And The Brenner-Thomee Approximation Theorem For Operator Semigroups, Frank Neubrander, Koray Ozer, Lee Windsperger

Faculty Publications

The computational powers of Mathematica are used to prove polynomial identities that are essential to obtain growth estimates for subdiagonal rational Pade approximations of the exponential function and to obtain new estimates of the constants of the Brenner-Thomee Approximation Theorem of Semigroup Theory.


Residual-Based Multi-Filter Methodology For All-Source Fault Detection, Exclusion, And Performance Monitoring, Juan Jurado, John F. Raquet, Christine M. Schubert Kabban Air Force Institute Of Technology, Jonathon S. Gipson Aug 2020

Residual-Based Multi-Filter Methodology For All-Source Fault Detection, Exclusion, And Performance Monitoring, Juan Jurado, John F. Raquet, Christine M. Schubert Kabban Air Force Institute Of Technology, Jonathon S. Gipson

Faculty Publications

All-source navigation has become increasingly relevant over the past decade with the development of viable alternative sensor technologies. However, as the number and type of sensors informing a system increases, so does the probability of corrupting the system with sensor modeling errors, signal interference, and undetected faults. Though the latter of these has been extensively researched, the majority of existing approaches have constrained faults to biases and designed algorithms centered around the assumption of simultaneously redundant, synchronous sensors with valid measurement models, none of which are guaranteed for all-source systems. As part of an overall all-source assured or resilient navigation …


Through-The-Wall Radar Detection Using Machine Learning, Aihua W. Wood, Ryan Wood, Matthew Charnley Aug 2020

Through-The-Wall Radar Detection Using Machine Learning, Aihua W. Wood, Ryan Wood, Matthew Charnley

Faculty Publications

This paper explores the through-the-wall inverse scattering problem via machine learning. The reconstruction method seeks to discover the shape, location, and type of hidden objects behind walls, as well as identifying certain material properties of the targets. We simulate RF sources and receivers placed outside the room to generate observation data with objects randomly placed inside the room. We experiment with two types of neural networks and use an 80-20 train-test split for reconstruction and classification.


An Ultra-Sparse Approximation Of Kinetic Solutions To Spatially Homogeneous Flows Of Non-Continuum Gas, Alexander Alekseenko, Amy Grandilli, Aihua W. Wood Feb 2020

An Ultra-Sparse Approximation Of Kinetic Solutions To Spatially Homogeneous Flows Of Non-Continuum Gas, Alexander Alekseenko, Amy Grandilli, Aihua W. Wood

Faculty Publications

We consider a compact approximation of the kinetic velocity distribution function by a sum of isotropic Gaussian densities in the problem of spatially homogeneous relaxation. Derivatives of the macroscopic parameters of the approximating Gaussians are obtained as solutions to a linear least squares problem derived from the Boltzmann equation with full collision integral. Our model performs well for flows obtained by mixing upstream and downstream conditions of normal shock wave with Mach number 3. The model was applied to explore the process of approaching equilibrium in a spatially homogeneous flow of gas. Convergence of solutions with respect to the model …


A Sequential Partial Information Bomber‐Defender Shooting Problem, Krishna Kalyanam, David W. Casbeer, Meir Pachter Feb 2020

A Sequential Partial Information Bomber‐Defender Shooting Problem, Krishna Kalyanam, David W. Casbeer, Meir Pachter

Faculty Publications

No abstract provided.


A Comparison Of The Trojan Y Chromosome Strategy To Harvesting Models For Eradication Of Non-Native Species, Jingjing Lyu, Pamela J. Schofield, Kristen M. Reaver, Matthew Beauregard, Rana D. Parshad Aug 2019

A Comparison Of The Trojan Y Chromosome Strategy To Harvesting Models For Eradication Of Non-Native Species, Jingjing Lyu, Pamela J. Schofield, Kristen M. Reaver, Matthew Beauregard, Rana D. Parshad

Faculty Publications

The Trojan Y Chromosome Strategy (TYC) is a promising eradication method for biological control of non-native species. The strategy works by manipulating the sex ratio of a population through the introduction of supermales that guarantee male offspring. In the current manuscript, we compare the TYC method with a pure harvesting strategy. We also analyze a hybrid harvesting model that mirrors the TYC strategy. The dynamic analysis leads to results on stability of solutions and bifurcations of the model. Several conclusions about the different strategies are established via optimal control methods. In particular, the results affirm that either a pure harvesting …


Large And Small Data Blow-Up Solutions In The Trojan Y Chromosome Model, Rana D. Parshad, Matthew Beauregard, Eric M. Takyi, Thomas Griffin, Landrey Bobo Jul 2019

Large And Small Data Blow-Up Solutions In The Trojan Y Chromosome Model, Rana D. Parshad, Matthew Beauregard, Eric M. Takyi, Thomas Griffin, Landrey Bobo

Faculty Publications

The Trojan Y Chromosome Strategy (TYC) is an extremely well investigated biological control method for controlling invasive populations with an XX-XY sex determinism. In [35, 36] various dynamical properties of the system are analyzed, including well posedness, boundedness of solutions, and conditions for extinction or recovery. These results are derived under the assumption of positive solutions. In the current manuscript, we show that if the introduction rate of trojan fish is zero, under certain large data assumptions, negative solutions are possible for the male population, which in turn can lead to finite time blow-up in the female and male populations. …


Quenching Estimates For A Non-Newtonian Filtration Equation With Singular Boundary Conditions, Matthew Beauregard, Burhan Selcuk Jul 2019

Quenching Estimates For A Non-Newtonian Filtration Equation With Singular Boundary Conditions, Matthew Beauregard, Burhan Selcuk

Faculty Publications

In this paper, the quenching behavior of the non-Newtonian filtration equation (φ(u))t = (|ux| r−2 ux)x with singular boundary conditions, ux (0, t) = u −p (0, t), ux (a, t) = (1 − u(a, t))−q is investigated. Various conditions on the initial condition are shown to guarantee quenching at either the left or right boundary. Theoretical quenching rates and lower bounds to the quenching time are determined when φ(u) = u and r = 2. Numerical experiments are provided to illustrate and provide additional validation of the theoretical estimates to the quenching rates and times.


Ergodicity For The 3d Stochastic Navier-Stokes Equations Perturbed By Lévy Noise, Manil T. Mohan, K. Sakthivel, Sivaguru S. Sritharan May 2019

Ergodicity For The 3d Stochastic Navier-Stokes Equations Perturbed By Lévy Noise, Manil T. Mohan, K. Sakthivel, Sivaguru S. Sritharan

Faculty Publications

In this work we construct a Markov family of martingale solutions for 3D stochastic Navier–Stokes equations (SNSE) perturbed by Lévy noise with periodic boundary conditions. Using the Kolmogorov equations of integrodifferential type associated with the SNSE perturbed by Lévy noise, we construct a transition semigroup and establish the existence of a unique invariant measure. We also show that it is ergodic and strongly mixing.
/="/">Abstract © Wiley.


A Variable Nonlinear Splitting Algorithm For Reaction Diffusion Systems With Self- And Cross-Diffusion, Matthew Beauregard, Joshua L. Padgett Jan 2019

A Variable Nonlinear Splitting Algorithm For Reaction Diffusion Systems With Self- And Cross-Diffusion, Matthew Beauregard, Joshua L. Padgett

Faculty Publications

Self- and cross-diffusion are important nonlinear spatial derivative terms that are included into biological models of predator-prey interactions. Self-diffusion models overcrowding effects, while cross-diffusion incorporates the response of one species in light of the concentration of another. In this paper, a novel nonlinear operator splitting method is presented that directly incorporates both self- and cross-diffusion into a computational efficient design. The numerical analysis guarantees the accuracy and demonstrates appropriate criteria for stability. Numerical experiments display its efficiency and accuracy