On the Numerical Solution of the Multiscale Models Characterising Noroviral Infectious Disease Systems Using the Multistage Spectral the Relaxation Method and the Nonstandard Finite Difference Method

Multiscale models for infectious disease systems are highly nonlinear. This presents substantial difficulties for both analysis and computation. Consequently, there is a continuous demand for the development of efficient numerical methods that provide reliable solutions. Over time, various numerical techniques have been developed for single-scale models. However, multiscale models have increasingly relied on built-in solvers and, more recently, the nonstandard finite difference method. The primary aim of this study is to offer an in-depth examination of the multistage spectral relaxation method (MSRM) applied to the multiscale model characterising norovirus infection. This study also provides a comprehensive comparison with the nonstandard finite difference method (NSFDM) in terms of convergence and accuracy, as well as their handling of highly nonlinear systems of ordinary differential equations. Both schemes were validated against an adaptive ODE45 scheme in MATLAB2025a, which served as the reference solver. The numerical outcomes indicate that the MSRM achieves a better accuracy level, with relative L2 errors ranging from 10−9 to 10−7 for Δτ≤0.10 and Nch=8. In contrast, the NSFDM provides solid first-order accuracy, yielding L2 errors around 10−4, while ensuring strict unconditional positivity and stability across all tested step sizes from h=0.01 to 2.0. This study features a parameter sensitivity analysis that varies βH,δH,αh,μV,μH, as well as semi-log error graphs, convergence rate visuals, and CPU performance benchmarks. The MSRM costs about 24 times more per run, but it offers better accuracy for between-host variables. This makes it the preferred method for high-fidelity applications. In contrast, the NSFDM is the best option for quick parameter sweeps and long simulations that require guaranteed positivity.

Authors

Institutions

Publication Details

Journal
Mathematical and Computational Applications
Published
2026-08-03
DOI
https://doi.org/10.3390/mca31040154
Primary Topic
Advanced Mathematical Modeling in Engineering
Type
article
Field-Weighted Citation Impact
0.00

Funders

Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

On the Numerical Solution of the Multiscale Models Characterising Noroviral Infectious Disease Systems Using the Multistage Spectral the Relaxation Method and the Nonstandard Finite Difference Method

Kizito Muzhinji
Mathematical and Computational Applications
Advanced Mathematical Modeling in Engineering
article

On the Numerical Solution of the Multiscale Models Characterising Noroviral Infectious Disease Systems Using the Multistage Spectral the Relaxation Method and the Nonstandard Finite Difference Method

Kizito Muzhinji
article en

Abstract

Multiscale models for infectious disease systems are highly nonlinear. This presents substantial difficulties for both analysis and computation. Consequently, there is a continuous demand for the development of efficient numerical methods that provide reliable solutions. Over time, various numerical techniques have been developed for single-scale models. However, multiscale models have increasingly relied on built-in solvers and, more recently, the nonstandard finite difference method. The primary aim of this study is to offer an in-depth examination of the multistage spectral relaxation method (MSRM) applied to the multiscale model characterising norovirus infection. This study also provides a comprehensive comparison with the nonstandard finite difference method (NSFDM) in terms of convergence and accuracy, as well as their handling of highly nonlinear systems of ordinary differential equations. Both schemes were validated against an adaptive ODE45 scheme in MATLAB2025a, which served as the reference solver. The numerical outcomes indicate that the MSRM achieves a better accuracy level, with relative L2 errors ranging from 10−9 to 10−7 for Δτ≤0.10 and Nch=8. In contrast, the NSFDM provides solid first-order accuracy, yielding L2 errors around 10−4, while ensuring strict unconditional positivity and stability across all tested step sizes from h=0.01 to 2.0. This study features a parameter sensitivity analysis that varies βH,δH,αh,μV,μH, as well as semi-log error graphs, convergence rate visuals, and CPU performance benchmarks. The MSRM costs about 24 times more per run, but it offers better accuracy for between-host variables. This makes it the preferred method for high-fidelity applications. In contrast, the NSFDM is the best option for quick parameter sweeps and long simulations that require guaranteed positivity.

Mathematical and Computational ApplicationsVol. 31(4)
University of Venda (ZA)
University of Venda
Openalex Percentile: Top 7%
Advanced Mathematical Modeling in Engineering
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.