pith. sign in

arxiv: 1907.03796 · v1 · pith:G6ARKLBFnew · submitted 2019-07-08 · 🧮 math.AP

Quenching estimates for a non-Newtonian filtration equation with singular boundary conditions

Matthew A. Beauregard , Burhan Selcuk This is my paper

Pith reviewed 2026-05-25 00:44 UTC · model grok-4.3

classification 🧮 math.AP
keywords quenchingnon-Newtonian filtrationsingular boundary conditionsquenching ratesinitial conditionsparabolic equationsfiltration equation
0
0 comments X

The pith

Conditions on initial data determine whether solutions quench at the left or right boundary for the non-Newtonian filtration equation with singular boundaries.

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

The paper demonstrates that suitable restrictions on the initial condition for the equation (ϕ(u))_t = (|u_x|^{r-2} u_x)_x force the solution to quench either at x=0 or at x=a under the given singular boundary conditions. When ϕ(u)=u and r=2, explicit quenching rates and lower bounds on the quenching time are derived. These results matter for predicting where and how rapidly solutions lose regularity in models of nonlinear diffusion with flux singularities at the ends. The work includes numerical checks that match the predicted rates and times.

Core claim

Various conditions on the initial condition 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.

What carries the argument

The singular boundary conditions u_x(0,t)=u^{-p}(0,t) and u_x(a,t)=(1-u(a,t))^{-q} that drive the solution to quench by making the spatial derivative unbounded as u approaches its singular values.

If this is right

  • Quenching occurs at the left boundary under initial data small enough near x=0.
  • Quenching occurs at the right boundary under initial data small enough near x=a.
  • Explicit rates describe how fast u approaches its singular value when ϕ(u)=u and r=2.
  • Lower bounds on the quenching time T are obtained for the linear case.

Where Pith is reading between the lines

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

  • The same initial-data comparison technique might locate quenching in related equations with different nonlinearities.
  • Physical filtration models with these boundary fluxes could use the rate formulas to estimate failure times.
  • Numerical schemes for the special case could be tuned to the predicted rates for accuracy checks.

Load-bearing premise

Solutions exist globally until the quenching time and the singular boundary conditions can be handled without extra regularity issues.

What would settle it

A numerical solution started from initial data satisfying one of the paper's conditions that quenches at the opposite boundary from the one predicted would show the location claim is incorrect.

Figures

Figures reproduced from arXiv: 1907.03796 by Burhan Selcuk, Matthew A. Beauregard.

Figure 1
Figure 1. Figure 1: (a) A graph of u ′ 0 (x) (RED) and x a (1−u0(x))q (BLUE) for u0(x) = 1 4 −4x−4x 2 . It is clear that u ′ 0 (x) ≥ x a (1 − u0(x))−q is satisfied throughout the domain 0 ≤ x ≤ 1/8. (b) A graph of u ′ 0 (x) (RED) and a−x a (u0(x))−p (BLUE) for u0(x) = 1 4 + 4x − 2x 2 . It is clear that u ′ 0 (x) ≥ a−x a (u0(x))−p is satisfied throughout the domain 0 ≤ x ≤ 1/8. 4. Numerical Approximation and Experiments Let xj… view at source ↗
Figure 2
Figure 2. Figure 2: Loglog plots of the numerical observed (a) 1 [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
read the original abstract

In this paper, the quenching behavior of the non-Newtonian filtration equation $(\phi (u))_{t}=(\left \vert u_{x}\right \vert ^{r-2}u_{x})_{x}$ with singular boundary conditions, $u_{x}\left( 0,t\right) =u^{-p}(0,t)$, $u_{x}\left( a,t\right) =(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 $\phi(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.

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.

Referee Report

1 major / 0 minor

Summary. The paper investigates quenching behavior for the non-Newtonian filtration equation (ϕ(u))_t = (|u_x|^{r-2} u_x)_x subject to the singular boundary conditions u_x(0,t) = u^{-p}(0,t) and u_x(a,t) = (1-u(a,t))^{-q}. It establishes conditions on the initial data that force quenching to occur at either the left or right boundary. For the linear case ϕ(u)=u and r=2, explicit quenching rates and lower bounds on the quenching time are derived. Numerical experiments are presented to support the theoretical estimates.

Significance. If the foundational assumptions hold, the work supplies new location criteria for boundary quenching and explicit rate formulas in a tractable special case, together with numerical corroboration. These results would add to the literature on singular-boundary nonlinear diffusion problems by providing concrete, testable predictions for the linear regime.

major comments (1)
  1. [analysis of the linear case (ϕ(u)=u, r=2)] The central claims on quenching location (under various initial-condition hypotheses) and the explicit rates/lower bounds for the case ϕ(u)=u, r=2 rest on the unstated assumption that classical or weak solutions exist globally on [0,T) up to the quenching time T and retain sufficient regularity for comparison principles and energy identities to apply. No existence theorem, continuation argument, or regularity result is cited or proved that closes this loop for the given singular fluxes; failure of this assumption would invalidate both the location results and the rate formulas.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and for highlighting this foundational issue. We address the concern directly below and will incorporate the necessary clarifications in a revised version.

read point-by-point responses

  1. Referee: [analysis of the linear case (ϕ(u)=u, r=2)] The central claims on quenching location (under various initial-condition hypotheses) and the explicit rates/lower bounds for the case ϕ(u)=u, r=2 rest on the unstated assumption that classical or weak solutions exist globally on [0,T) up to the quenching time T and retain sufficient regularity for comparison principles and energy identities to apply. No existence theorem, continuation argument, or regularity result is cited or proved that closes this loop for the given singular fluxes; failure of this assumption would invalidate both the location results and the rate formulas.

    Authors: We agree that the manuscript does not explicitly address the existence and regularity of solutions up to the quenching time T. The claims for the linear case (ϕ(u)=u, r=2) rely on the applicability of comparison principles and energy methods, which presuppose sufficient regularity. In the revised manuscript we will add a dedicated remark (or short subsection) in the preliminaries section. This remark will (i) recall standard local existence results for the linear heat equation with nonlinear boundary flux via approximation or fixed-point methods, (ii) note that the singular boundary conditions remain integrable up to but not including T, and (iii) cite relevant literature on degenerate parabolic equations with singular data to cover the general (ϕ,r) setting. These additions will explicitly close the logical loop without altering the main quenching-location and rate arguments. revision: yes

Circularity Check

0 steps flagged

No circularity: quenching rates derived from PDE analysis, not by construction

full rationale

The paper states that theoretical quenching rates and lower bounds are determined when ϕ(u)=u and r=2, using conditions on initial data to locate quenching at boundaries. No quoted steps reduce a claimed prediction to a fitted parameter, self-citation chain, or definitional equivalence. Numerical experiments are presented separately for validation. The derivation chain remains self-contained against the stated PDE and singular BCs; global existence until quenching is an assumption but does not create circularity in the rate formulas themselves.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract, the paper relies on standard existence and regularity assumptions from parabolic PDE theory rather than introducing new free parameters or invented entities.

axioms (1)
  • domain assumption Solutions to the PDE exist up to the quenching time and satisfy the singular boundary conditions in an appropriate weak sense.
    Required for the quenching analysis and rate derivations to be meaningful.

pith-pipeline@v0.9.0 · 5661 in / 1284 out tokens · 29331 ms · 2026-05-25T00:44:04.056674+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages

  1. [1]

    A., Numerical solutions to singular reaction diffusion equation over elliptical domains , Appl

    Beauregard, M. A., Numerical solutions to singular reaction diffusion equation over elliptical domains , Appl. Math. Comput., Vol. 254, 75-91, 2015

    work page 2015

  2. [2]

    A., Sheng, Q., An adaptive splitting approach for the quenching solution of reaction-diffusion equations ove r nonuniform grids, J

    Beauregard, M. A., Sheng, Q., An adaptive splitting approach for the quenching solution of reaction-diffusion equations ove r nonuniform grids, J. Comput. Appl. Math., Vol. 241, 30-44, 2013

    work page 2013

  3. [3]

    Beauregard, M.A., Sheng, Q., Explorations and expectations of equidis- tribution adaptations for nonlinear quenching problems , Adv. Appl. Math. Mech., Vol. 5, 407-422, 2013

    work page 2013

  4. [4]

    A., Padgett, J., Parshad, R

    Beauregard, M. A., Padgett, J., Parshad, R. D., A nonlinear splitting algorithm for systems of partial differential equations wit h self-diffusion , J. Comput. Appl. Math., Vol. 31, 8-25, 2017. 16

    work page 2017

  5. [5]

    Chan, C.Y., Yuen, S.I., Parabolic problems with nonlinear absorptions and releases at the boundaries , Appl. Math. Comput., 121, 203-209, 2001

    work page 2001

  6. [6]

    C. Y. Chan and L. Ke, Parabolic quenching for nonsmooth convex do- mains, J. Math. Anal. Appl. 186 (1994), pp. 52-65

    work page 1994

  7. [7]

    C. Y. Chan, A quenching criterion for a multi-dimensional parabolic problem due to a concentrated nonlinear source, J. Comput. Appl. Math., 235 (2011), pp. 3724-3727

    work page 2011

  8. [8]

    Cheng, P

    H. Cheng, P. Lin, Q. Sheng and R.C.E. Tan, Solving degenerate reaction- diffusion equations via variable step Peaceman-Rachford sp litting, SIAM J. Sci. Comput. 25 (2003), pp. 1273-1292

    work page 2003

  9. [9]

    Comput., Vol

    Hundsdorfer, W., Unconditional convergence of some Crank-Nicolson LOD method for initial-boundary value problems , Math. Comput., Vol. 58, 35-53, 1992

    work page 1992

  10. [10]

    Kawarada, On solutions of initial-boundary problem for ut = uxx + 1 1− u , Publ

    H. Kawarada, On solutions of initial-boundary problem for ut = uxx + 1 1− u , Publ. Res. Inst. Math. Sci., 10 (1975), pp. 729-736

    work page 1975

  11. [11]

    Z. Li, C. Mu, Critical exponents for a fast diffusive polytropic filtratio n equation with nonlinear boundary flux , J. Math. Anal. Appl., Vol. 346, Iss. 1, pp. 55-64, 2008

    work page 2008

  12. [12]

    X. Li, C. Mu, Q. Zhang, and S. Zhou, Quenching for a Non-Newtonian Filtration Equation with a Singular Boundary Condition , Abstr. Appl. Anal., Vol. 2012, Article ID 539161, doi:10.1155/2012/539161

    work page doi:10.1155/2012/539161 2012

  13. [13]

    Y. Mi, X. Wang and C. Mu, Blow-up set for the non-Newtonian poly- tropic filtration equation subjected to nonlinear Neumann b oundary con- dition, Applicable Analysis (2013) Vol. 92, No. 6, 1332–1344

    work page 2013

  14. [14]

    Selcuk, B., Ozalp, N., The quenching behavior of a semilinear heat equa- tion with a singular boundary outflux , Quarterly of Applied Mathemat- ics, Vol. 72, Iss. 4, 747-752, 2014

    work page 2014

  15. [15]

    Selcuk, B., Ozalp, N., Quenching behavior of semilinear heat equations with singular boundary conditions , Electron. J. Diff. Equ., Vol. 2015, No. 311, 1-13, 2015. 17

    work page 2015

  16. [16]

    Sheng and A

    Q. Sheng and A. Q. M. Khaliq, A revisit of the semi-adaptive method for singular degenerate reaction-diffusion equations , East Asia J. Appl. Math., 2 (2012), pp. 185-203

    work page 2012

  17. [17]

    Co., Belmont, CA, ISBN: 0-534-09984, 112- 134, 1989

    Strikwerda, John C., Finite difference schemes and partial differential equations, Wadsworth Publ. Co., Belmont, CA, ISBN: 0-534-09984, 112- 134, 1989

    work page 1989

  18. [18]

    Ozalp N, Selcuk B, The quenching behavior of a nonlinear parabolic equation with a singular parabolic with a singular boundary condition, Hacettepe Journal of Mathematics and Statistics,44: 615–621, 2 015

    work page

  19. [19]

    Z. Wang, J. Yin, and C. Wang, Critical exponents of the non-Newtonian polytropic filtration equation with nonlinear boundary con dition, Applied Mathematics Letters, vol. 20, no. 2, pp. 142–147, 2007

    work page 2007

  20. [20]

    Y. Zhi, C. Mu, The quenching behavior of a nonlinear parabolic equation with a nonlinear boundary outflux , Appl. Math. Comput., 184, 624-630, 2007. 18

    work page 2007