arxiv: 2605.11925 · v1 · submitted 2026-05-12 · 🧮 math.AP

Recognition: 2 theorem links

· Lean Theorem

A degenerate reaction-diffusion SIR model in interconnected regions

Abderrahim Zafrar, Jawad Salhi, Omar Elamraoui

Pith reviewed 2026-05-13 05:19 UTC · model grok-4.3

classification 🧮 math.AP MSC 35K5792D30

keywords SIR modelreaction-diffusiondegenerate diffusionwell-posednessFaedo-Galerkininterconnected regionsboundary conditionsepidemic modeling

0 comments

The pith

A degenerate reaction-diffusion SIR model with switching boundaries admits unique positive weak solutions in interconnected regions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper builds a coupled SIR system for two regions linked by migration, where diffusion degenerates to zero when local population density vanishes and boundary conditions switch dynamically from Robin to Neumann type to capture policy interventions such as lockdowns. The Faedo-Galerkin method is used to prove existence, uniqueness, and positivity of weak solutions for this system. This matters because it supplies a mathematically consistent way to simulate how mobility restrictions alter epidemic curves near and across borders. Finite-volume numerics then illustrate the practical influence of the degeneracy and the switches.

Core claim

The authors claim that the time-dependent degenerate reaction-diffusion SIR system posed in two adjacent domains with switching boundary conditions possesses a unique positive weak solution. The proof proceeds by constructing Faedo-Galerkin approximations that respect the degeneracy induced by vanishing density and the Robin-to-Neumann transition, then passing to the limit while preserving non-negativity.

What carries the argument

Degenerate diffusion coefficient vanishing with population density, paired with time-dependent switching of boundary conditions from Robin to Neumann.

If this is right

The model supplies a rigorous basis for comparing epidemic outcomes with and without mobility restrictions between regions.
Positivity of the weak solutions guarantees that simulated population compartments remain biologically meaningful.
Existence and uniqueness justify applying standard numerical schemes such as finite-volume discretizations without fear of spurious negative values.
The dynamic boundary switch can be tuned to study the timing and severity of lockdown effects on cross-border transmission.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same degeneracy mechanism might naturally limit mixing at low-density borders and could be tested against observed quarantine data.
Extending the two-region setup to a network of regions would allow systematic study of optimal intervention timing.
Stability or long-time behavior analysis of the proven weak solutions could reveal whether the model predicts endemic equilibria under repeated policy switches.

Load-bearing premise

Population density is allowed to reach zero locally, which forces the diffusion coefficient to degenerate, and boundary conditions are permitted to switch in time according to external policy decisions.

What would settle it

A sequence of Faedo-Galerkin approximations that fails to converge to a non-negative limit when density reaches zero near the interface would falsify the existence and positivity claim.

Figures

Figures reproduced from arXiv: 2605.11925 by Abderrahim Zafrar, Jawad Salhi, Omar Elamraoui.

**Figure 2.** Figure 2: Diffusion coefficient function σ(y, t) for λ = 0.01, a = 0.01 and ta = 50. This form ensures that the diffusion coefficient decreases exponentially over time, as shown in the figure 2 starting from a peak at t = ta days. The term e −a·(t−ta) captures this temporal decay. The spatial variation is governed by the term (2−y)y, which is a quadratic function of y. This quadratic term ensures that the diffusion … view at source ↗

**Figure 3.** Figure 3: SIR model in Ω1 and Ω2 The table 2 demonstrates how varying the parameter λ influences key epidemic metrics across different regions. As λ increases, the peak number of infected individuals decreases significantly, from 3673 at λ = 10−5 to 2012 at λ = 1. Similarly, both the total recovered individuals and total population size decrease with higher λ values, indicating a more dispersed disease impact. Impo… view at source ↗

**Figure 4.** Figure 4: Heatmaps showing the spatial and temporal evolution of the SIR model, [PITH_FULL_IMAGE:figures/full_fig_p024_4.png] view at source ↗

**Figure 5.** Figure 5: impact of varying λi , for i = 1, 2 Conclusion The numerical simulations presented in this study highlight the rich spatio-temporal dynamics of the proposed degenerate SIR model. By simulating a one-dimensional domain divided into two interacting regions, we demonstrated how population movement, infection thresholds, and policy-driven boundary dynamics jointly influence the epidemic trajectory. The inclus… view at source ↗

read the original abstract

This paper presents a novel time-space SIR (Susceptible-Infected-Recovered) model for simulating infectious disease dynamics in two interconnected regions. The model is formulated as a coupled reaction-diffusion system with boundary conditions that dynamically switch from Robin to Neumann types, effectively modelling policy-driven interventions such as lockdowns. A key innovation lies in the incorporation of degenerate diffusion, arising from vanishing population density, which significantly influences transmission behaviour near regional borders. The wellposedness of the model is rigorously established using the Faedo-Galerkin method, ensuring the existence, uniqueness, and positivity of weak solutions. Numerical simulations, performed using the Finite Volume Method, validate the theoretical findings and demonstrate the impact of migration and mobility restrictions on epidemic progression. This framework offers valuable insights for understanding and controlling disease spread in spatially heterogeneous and interconnected settings.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper sets up a SIR model with degenerate diffusion and policy-driven boundary switching in two regions, but the Faedo-Galerkin wellposedness argument is likely to fail on uniform bounds once density vanishes.

read the letter

The paper sets up a coupled reaction-diffusion SIR system across two interconnected regions. Diffusion degenerates when local density drops to zero, and the boundary operator switches dynamically from Robin to Neumann to capture interventions such as lockdowns. They assert existence, uniqueness, and positivity of weak solutions by Faedo-Galerkin and illustrate the behavior with finite-volume numerics that track migration and restriction effects on the epidemic curves. The specific combination of vanishing-density degeneracy with time-dependent boundary changes in a multi-region setting is the clearest novelty; prior SIR diffusion models usually keep positive diffusion or fixed boundary types. The numerics are standard but give a concrete picture of how mobility restrictions alter spread near the interface. The main weakness is the wellposedness step. Standard Galerkin needs coercive estimates and compactness that survive the limit, yet degeneracy removes the lower bound on the diffusion coefficient precisely where density approaches zero, especially near borders. The switching boundary term adds further trouble for trace control and monotonicity. The abstract gives no indication that the a priori bounds stay uniform in the degeneracy parameter or that the limit identification works after the boundary change. Without those details the existence claim does not yet hold. This is for readers working on spatially structured epidemic models or degenerate parabolic systems who want an applied example. A PDE analyst focused on wellposedness would find the setup worth checking, while a modeler might use the formulation for further simulation. The work is coherent enough on its own terms to merit referee time, though the analysis section will probably need substantial revision. I would send it to referees who handle degenerate equations rather than desk-reject it.

Referee Report

3 major / 2 minor

Summary. The manuscript introduces a coupled reaction-diffusion SIR model on two interconnected regions that incorporates degenerate diffusion (vanishing with local population density) and time-dependent switching of boundary conditions from Robin to Neumann type to model policy interventions such as lockdowns. The central claims are that the wellposedness of weak solutions (existence, uniqueness, positivity) is rigorously established via the Faedo-Galerkin method and that Finite Volume Method simulations confirm the impact of migration and mobility restrictions on epidemic dynamics.

Significance. If the wellposedness result is correct, the work supplies a mathematically consistent framework for analyzing spatially heterogeneous epidemic spread under dynamic interventions, with the degeneracy capturing realistic density-dependent mobility near borders. The numerical component provides concrete illustrations of how switching boundary conditions affect transmission, which could be useful for policy modeling in interconnected populations.

major comments (3)

[Wellposedness section (Faedo-Galerkin estimates)] The Faedo-Galerkin construction for existence requires uniform-in-degeneracy a priori estimates to pass to the limit. When the total density S+I+R approaches zero near the interface, the diffusion coefficient vanishes and coercivity of the bilinear form is lost; the manuscript must exhibit explicit bounds (e.g., on the L^2(0,T;H^1) norm of the approximations) that remain independent of the degeneracy parameter and of the time-dependent boundary switch. Without these, the compactness argument and identification of the limit solution are not justified.
[Uniqueness and positivity argument] Uniqueness and positivity for the weak solution rely on monotonicity or Lipschitz properties of the reaction terms and on the trace operators under the switching boundary condition. The dynamic change from Robin to Neumann perturbs the boundary integrals; the proof must verify that the difference of two solutions still satisfies a Gronwall-type inequality that is uniform with respect to the switching times.
[Numerical section] The numerical Finite Volume scheme is presented as validation, yet no convergence analysis or error estimates are given that account for the degeneracy and the time-dependent boundary switch. It is therefore unclear whether the discrete solutions remain consistent with the weak formulation when density vanishes locally.

minor comments (2)

[Model formulation] Notation for the switching times and the precise definition of the time-dependent boundary operator should be stated explicitly at the beginning of the model section to avoid ambiguity when reading the weak formulation.
[Introduction] The abstract states that the model is 'rigorously established'; the introduction should clarify which theorem contains the full wellposedness statement (existence + uniqueness + positivity) rather than leaving it implicit.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough and helpful report. We agree with the points raised and will revise the manuscript to address them, as detailed in the point-by-point responses below.

read point-by-point responses

Referee: The Faedo-Galerkin construction for existence requires uniform-in-degeneracy a priori estimates to pass to the limit. When the total density S+I+R approaches zero near the interface, the diffusion coefficient vanishes and coercivity of the bilinear form is lost; the manuscript must exhibit explicit bounds (e.g., on the L^2(0,T;H^1) norm of the approximations) that remain independent of the degeneracy parameter and of the time-dependent boundary switch. Without these, the compactness argument and identification of the limit solution are not justified.

Authors: We thank the referee for this comment. We will revise the well-posedness section to provide more detailed a priori estimates that explicitly demonstrate uniformity with respect to the degeneracy parameter and the boundary condition switches. This will include bounding the L^2(0,T;H^1) norm independently of these aspects by using the structure of the model and additional regularization arguments if necessary. revision: yes
Referee: Uniqueness and positivity for the weak solution rely on monotonicity or Lipschitz properties of the reaction terms and on the trace operators under the switching boundary condition. The dynamic change from Robin to Neumann perturbs the boundary integrals; the proof must verify that the difference of two solutions still satisfies a Gronwall-type inequality that is uniform with respect to the switching times.

Authors: We will revise the uniqueness proof to include a verification that the Gronwall inequality holds uniformly with respect to the switching times, by carefully estimating the boundary terms arising from the dynamic conditions and confirming that the Lipschitz properties of the reaction terms yield a constant independent of the switching sequence. revision: yes
Referee: The numerical Finite Volume scheme is presented as validation, yet no convergence analysis or error estimates are given that account for the degeneracy and the time-dependent boundary switch. It is therefore unclear whether the discrete solutions remain consistent with the weak formulation when density vanishes locally.

Authors: We agree that providing a convergence analysis would be beneficial. The Finite Volume method used is monotone and consistent with the weak form by using upwind fluxes for the advection (migration) terms and centered differences for diffusion, with the diffusion coefficient set to zero when the local density is below a small threshold to handle degeneracy. The time-dependent boundary conditions are implemented by adjusting the boundary fluxes at each time step. While a full error estimate is beyond the scope of the current work, which focuses on the continuous model, we will add a subsection discussing the consistency of the scheme with the weak solutions in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No circularity: wellposedness established via standard Faedo-Galerkin on explicitly defined degenerate system

full rationale

The paper introduces an explicit coupled reaction-diffusion SIR system with density-dependent degeneracy and time-dependent Robin-to-Neumann switching, then invokes the classical Faedo-Galerkin procedure to obtain existence, uniqueness and positivity of weak solutions. No load-bearing step reduces to a self-citation, a fitted parameter renamed as prediction, or an ansatz imported from prior author work; the estimates and compactness arguments are derived directly from the stated equations and standard functional-analytic tools. The derivation is therefore self-contained against external benchmarks and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on mathematical analysis of a new PDE system; no free parameters or invented entities mentioned in abstract. Relies on domain assumptions for SIR models and degeneracy conditions.

axioms (1)

standard math Standard assumptions for existence of weak solutions in degenerate parabolic systems
Invoked in the wellposedness proof using Faedo-Galerkin.

pith-pipeline@v0.9.0 · 5436 in / 1335 out tokens · 74711 ms · 2026-05-13T05:19:23.342130+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear
well-posedness ... using the Faedo-Galerkin method, ensuring the existence, uniqueness, and positivity of weak solutions
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
degenerate diffusion arising from vanishing population density

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages

[1]

A contribution to the mathematical theory of epidemics.Proceedings of the royal society of lon- don

William Ogilvy Kermack and Anderson G McKendrick. A contribution to the mathematical theory of epidemics.Proceedings of the royal society of lon- don. Series A, Containing papers of a mathematical and physical character, 115(772):700–721, 1927

work page 1927
[2]

Springer, 2012

Fred Brauer, Carlos Castillo-Chavez, and Carlos Castillo-Chavez.Mathematical models in population biology and epidemiology, volume 2. Springer, 2012

work page 2012
[3]

An Introduction

James Dickson Murray.MathematicalBiology I. An Introduction. Springer, 2002

work page 2002
[4]

Fractional time-delay mathematical modeling of oncolytic virotherapy.Chaos, Solitons & Fractals, 150:111123, 2021

Pushpendra Kumar, Vedat Suat Erturk, Abdullahi Yusuf, and Sunil Kumar. Fractional time-delay mathematical modeling of oncolytic virotherapy.Chaos, Solitons & Fractals, 150:111123, 2021

work page 2021
[5]

Optimal control for variable order fractional hiv/aids and malaria mathematical models with multi-time delay.Alexandria Engineering Journal, 59(5):3149–3162, 2020

NH Sweilam, SM Al-Mekhlafi, ZN Mohammed, and D Baleanu. Optimal control for variable order fractional hiv/aids and malaria mathematical models with multi-time delay.Alexandria Engineering Journal, 59(5):3149–3162, 2020

work page 2020
[6]

Stochastic epidemic modeling

Priscilla E Greenwood and Luis F Gordillo. Stochastic epidemic modeling. InMathematical and statistical estimation approaches in epidemiology, pages 31–52. Springer, 2009

work page 2009
[7]

JAJ Metz. The epidemic in a closed population with all susceptibles equally vulnerable; some results for large susceptible populations and small initial in- fections.Acta biotheoretica, 27(1):75–123, 1978

work page 1978
[8]

Stochastic and deterministic models in epidemiology

Carlos Castillo-Chavez, Samuel Fridman, Xiaolong Luo, et al. Stochastic and deterministic models in epidemiology. 1993

work page 1993
[9]

Stochastic mathematical models in epidemiology.Information, 20:6185–6196, 2017

Meksianis Z Ndii and Asep K Supriatna. Stochastic mathematical models in epidemiology.Information, 20:6185–6196, 2017

work page 2017
[10]

Boundary optimal control of time–space sir model with nonlinear robin boundary condition.International Journal of Dynamics and Control, 10(4):1279–1290, 2022

E-H Essoufi and A Zafrar. Boundary optimal control of time–space sir model with nonlinear robin boundary condition.International Journal of Dynamics and Control, 10(4):1279–1290, 2022. 26

work page 2022
[11]

Spatio-temporal sir model with robin boundary condition and automatic lockdown policy.International Journal of Applied and Computational Mathematics, 9(1):1–18, 2023

O Elamraoui, E-H Essoufi, and A Zafrar. Spatio-temporal sir model with robin boundary condition and automatic lockdown policy.International Journal of Applied and Computational Mathematics, 9(1):1–18, 2023

work page 2023
[12]

Asymptotic profiles of the steady states for an sis epidemic reaction-diffusion model.Dis- crete & Continuous Dynamical Systems, 21(1):1, 2008

Linda JS Allen, Ben M Bolker, Yuan Lou, and Andrew L Nevai. Asymptotic profiles of the steady states for an sis epidemic reaction-diffusion model.Dis- crete & Continuous Dynamical Systems, 21(1):1, 2008

work page 2008
[13]

On spatially explicit models of cholera epidemics.Journal of the Royal Society Interface, 7(43):321–333, 2010

E Bertuzzo, Renato Casagrandi, Marino Gatto, I Rodriguez-Iturbe, and A Ri- naldo. On spatially explicit models of cholera epidemics.Journal of the Royal Society Interface, 7(43):321–333, 2010

work page 2010
[14]

John Wiley & Sons, 2004

Robert Stephen Cantrell and Chris Cosner.Spatial ecology via reaction-diffusion equations. John Wiley & Sons, 2004

work page 2004
[15]

International migration 2020 highlights

Clare Menozzi. International migration 2020 highlights. 2021

work page 2020
[16]

Models for transmission of dis- ease with immigration of infectives.Mathematical Biosciences, 171(2):143–154, 2001

Fred Brauer and Pauline van den Driessche. Models for transmission of dis- ease with immigration of infectives.Mathematical Biosciences, 171(2):143–154, 2001

work page 2001
[17]

An epidemic patchy model with entry–exit screening.Bulletin of mathematical biology, 77(7):1237–1255, 2015

Xinxin Wang, Shengqiang Liu, Lin Wang, and Weiwei Zhang. An epidemic patchy model with entry–exit screening.Bulletin of mathematical biology, 77(7):1237–1255, 2015

work page 2015
[18]

Global analysis of two tubercu- losis models.Journal of Dynamics and Differential Equations, 16(1):139–166, 2004

C Connell McCluskey and P van den Driessche. Global analysis of two tubercu- losis models.Journal of Dynamics and Differential Equations, 16(1):139–166, 2004

work page 2004
[19]

Global stability for an sei model of infectious disease with immigration.Applied Mathematics and Computation, 243:684–689, 2014

Ram P Sigdel and C Connell McCluskey. Global stability for an sei model of infectious disease with immigration.Applied Mathematics and Computation, 243:684–689, 2014

work page 2014
[20]

Disease dynamics for the hometown of migrant workers.Mathematical Biosciences & Engineering, 11(5):1175, 2014

Ram P Sigdel and C Connell McCluskey. Disease dynamics for the hometown of migrant workers.Mathematical Biosciences & Engineering, 11(5):1175, 2014

work page 2014
[21]

Reaction diffusion equations with boundary degeneracy

Huashui Zhan. Reaction diffusion equations with boundary degeneracy. 2016

work page 2016
[22]

Blow-up at the boundary for degenerate semilinear parabolic equa- tions.Archive for Rational Mechanics and Analysis, 114(1):57–77, 1991

MS Floater. Blow-up at the boundary for degenerate semilinear parabolic equa- tions.Archive for Rational Mechanics and Analysis, 114(1):57–77, 1991

work page 1991
[23]

A compact adaptive approach for degenerate singular reaction-diffusion equations.Numerical Methods for Partial Differential Equations, 34(4):1166–1187, 2018

Yongbin Ge, Zhiquan Cai, and Qin Sheng. A compact adaptive approach for degenerate singular reaction-diffusion equations.Numerical Methods for Partial Differential Equations, 34(4):1166–1187, 2018. 27

work page 2018
[24]

Indirect diffusion effect in de- generate reaction-diffusion systems.SIAM Journal on Mathematical Analysis, 52(5):4314–4361, 2020

Amit Einav, Jeffrey J Morgan, and Bao Q Tang. Indirect diffusion effect in de- generate reaction-diffusion systems.SIAM Journal on Mathematical Analysis, 52(5):4314–4361, 2020

work page 2020
[25]

A time-space sir model for disease spread across two regions: Analysis and numerical simulations.In- ternational Journal of Biomathematics, 18(05):2450006, 2025

O Elamraoui, EH Essoufi, H Rahnaoui, and A Zafrar. A time-space sir model for disease spread across two regions: Analysis and numerical simulations.In- ternational Journal of Biomathematics, 18(05):2450006, 2025

work page 2025
[26]

Springer Singapore, Singapore, 2013

Nicolas Privault.Spatial Poisson Processes, pages 225–239. Springer Singapore, Singapore, 2013

work page 2013
[27]

C. Wang. Approximate controllability of a class of semilinear systems with boundary degeneracy.Journal of Evolution Equations, 10(1):163–193, 2010

work page 2010
[28]

F. D. Araruna, B. S. V. Ara´ ujo, and E. Ferna´ andez-Cara. Carleman estimates for some two-dimensional degenerate parabolic pdes and applications.SIAM Journal on Control and Optimization, 57(6):3985–4010, 2019

work page 2019
[29]

W. Wu, Y. Hu, Y. Liu, and D. Yang. Carleman estimates for degenerate parabolic equations with single interior degeneracy point and applications. Mathematical Control and Related Fields, 2025

work page 2025
[30]

Roger Temam.Navier–Stokes equations: theory and numerical analysis, volume

work page
[31]

American Mathematical Society, 2024. 28

work page 2024