REVIEW 2 major objections 3 minor 17 references
Nonlinear evolution of weak discontinuity waves in Darcy-type porous media
T0 review · 2 major / 3 minor · reviewed 2026-05-24 · grok-4.3
Pith's one-line read In Darcy-type porous media, compressive weak discontinuities above a critical initial amplitude grow into shocks while smaller ones decay.
desk verdict The paper's characteristic derivation for weak discontinuities does not hold because Darcy flow produces a parabolic system without real characteristics. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Evolution equation for weak discontinuities in the characteristic plane, obtained by using the characteristics of the quasilinear system as the reference coordinate system.
What would settle it
A direct numerical integration of the original governing equations showing that a compressive disturbance whose initial amplitude exceeds the predicted critical value fails to form a shock.
Extended reading notes
Core claim
Any compressive disturbance with initial amplitude greater than the critical one always grows into a shock wave, while the initial amplitude less than the critical one always decays. The breaking point at the wave front is determined, and porosity affects the steepening and flattening processes for waves with planar, cylindrical, and spherical symmetry.
Load-bearing premise
Weak discontinuities propagate along the characteristic paths of the governing quasilinear system.
Editorial extensions
If this is right
- The location where an acceleration wave breaks can be calculated explicitly from initial data and porosity.
- Increasing porosity delays or prevents shock formation for a given initial amplitude.
- The same critical-amplitude threshold applies across planar, cylindrical, and spherical geometries.
- Disturbances below the critical amplitude decay without forming discontinuities.
Reading between the lines
- The model could be tested by measuring shock onset times in laboratory flows through packed beds with controlled porosity.
- The derived evolution equation may serve as a reduced-order description for wave steepening in other quasilinear hyperbolic systems with damping.
- Extensions to radially symmetric flows in higher dimensions would follow the same characteristic reduction if the governing equations remain quasilinear.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes propagation of nonlinear weak discontinuity waves in one-dimensional unsteady compressible flow through Darcy-type porous media. It assumes weak discontinuities propagate along characteristics of the governing quasilinear system (used as reference coordinates), derives an evolution equation in the characteristic plane, determines the breaking point at the wave front, and examines how porosity affects steepening/flattening of acceleration waves for planar, cylindrical, and spherical symmetry. A critical initial amplitude is identified such that compressive disturbances above it form shocks while those below decay.
Significance. If the system is shown to be strictly hyperbolic, the explicit critical-amplitude thresholds and porosity dependence for different geometries would constitute a concrete, testable prediction for wave breaking in porous-media flows. The approach follows standard techniques for hyperbolic systems and could inform applications in geophysics or filtration if the hyperbolicity assumption holds.
major comments (2)
- [Abstract / foundational modeling choice] Abstract (modeling assumption): The central step assumes the governing equations form a quasilinear hyperbolic system whose characteristics can serve as the reference coordinate system for the evolution equation of weak-discontinuity amplitude. Standard Darcy drag replaces momentum balance with an algebraic relation, producing a parabolic pressure equation with no real characteristics or finite propagation speed. The manuscript must exhibit the full governing equations (continuity + momentum) and prove they remain strictly hyperbolic; absent this verification the derived evolution equation and critical-amplitude result rest on an unsupported premise.
- [Abstract] Abstract (critical-amplitude claim): The statement that “any compressive disturbance with initial amplitude greater than the critical one always grows into a shock wave, while the initial amplitude less than the critical one always decays” is obtained only after the characteristic evolution equation is derived. Because that equation depends on the unverified hyperbolicity, the claim cannot be regarded as established; a direct check against the original system (or a numerical solution of the full equations) is required before the result can be accepted.
minor comments (3)
- [Abstract] Abstract: “propagate long the characteristic path” should read “propagate along the characteristic path.”
- [Abstract] Abstract: “flattering of acceleration waves” should read “flattening of acceleration waves.”
- [Abstract] Abstract: The phrase “one dimensional space, unsteady and compressible flow” is ambiguous; clarify whether the analysis is strictly one-dimensional or employs symmetry reductions for cylindrical/spherical cases.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive comments. We address each major comment below, clarifying the modeling assumptions and proposing targeted revisions.
read point-by-point responses
-
Referee: [Abstract / foundational modeling choice] Abstract (modeling assumption): The central step assumes the governing equations form a quasilinear hyperbolic system whose characteristics can serve as the reference coordinate system for the evolution equation of weak-discontinuity amplitude. Standard Darcy drag replaces momentum balance with an algebraic relation, producing a parabolic pressure equation with no real characteristics or finite propagation speed. The manuscript must exhibit the full governing equations (continuity + momentum) and prove they remain strictly hyperbolic; absent this verification the derived evolution equation and critical-amplitude result rest on an unsupported premise.
Authors: We agree that explicit display of the equations and verification of hyperbolicity will strengthen the manuscript. In the revised version we will present the full continuity and momentum equations in Section 2. The momentum equation retains the unsteady and convective inertial terms together with the Darcy drag (treated as a source term). The quasilinear principal part is identical to that of one-dimensional gas dynamics; the eigenvalues are u ± a (a the local sound speed), which are real and distinct, establishing strict hyperbolicity. The corresponding right and left eigenvectors will also be supplied. revision: yes
-
Referee: [Abstract] Abstract (critical-amplitude claim): The statement that “any compressive disturbance with initial amplitude greater than the critical one always grows into a shock wave, while the initial amplitude less than the critical one always decays” is obtained only after the characteristic evolution equation is derived. Because that equation depends on the unverified hyperbolicity, the claim cannot be regarded as established; a direct check against the original system (or a numerical solution of the full equations) is required before the result can be accepted.
Authors: Once strict hyperbolicity is verified (as outlined above), the evolution equation for the discontinuity amplitude follows from the standard theory of weak discontinuities propagating along characteristics. The critical initial amplitude is obtained directly from the sign of the quadratic coefficient in that Riccati equation and is therefore established within the hyperbolic framework. We do not consider a separate numerical integration of the original system necessary to validate the analytical result, but we will add an explicit statement in the revised abstract and conclusions that all claims presuppose the verified hyperbolicity of the governing system. revision: partial
Circularity Check
No circularity: standard characteristic analysis applied to quasilinear system without reduction to inputs or self-citations.
full rationale
The derivation begins from the governing quasilinear system for compressible flow in Darcy-type porous media, assumes weak discontinuities propagate along characteristics (a conventional modeling step for hyperbolic systems), and derives an evolution equation in the characteristic plane to determine breaking points and critical amplitudes. No equations or claims reduce by construction to fitted parameters, self-defined quantities, or load-bearing self-citations; the critical-amplitude result follows directly from solving the derived transport equation along characteristics. The approach is self-contained against external benchmarks for hyperbolic PDE methods and does not rename known results or smuggle ansatzes via citation.
Assumptions & free parameters
assumptions (2)
- domain assumption The flow is governed by a quasilinear hyperbolic system describing unsteady compressible flow in Darcy-type porous media.
- domain assumption Weak discontinuities propagate along the characteristic paths of the system.
Cite this review
Pith. "Pith review of Nonlinear evolution of weak discontinuity waves in Darcy-type porous media." pith.science (2026). https://pith.science/paper/YANNYCYN
@misc{pith2026190709718,
author = {Pith},
title = {Pith review of: Nonlinear evolution of weak discontinuity waves in Darcy-type porous media},
year = {2026},
howpublished = {\url{https://pith.science/paper/YANNYCYN}},
note = {Machine review of arXiv:1907.09718}
}
read the original abstract
The propagation of nonlinear waves in one dimensional space, unsteady and compressible flow in Darcy-type porous media is analyzed. It is assumed that the weak discontinuity propagate long the characteristic path using the characteristics of the governing quasilinear system as the reference coordinate system. Evolution equation in the characteristic plane is derived. As an application of the theory the breaking point at the wave front is determined. It is assessed as to how the porosity of the medium affects the process of steepening and flattering of acceleration waves with planar, cylindrical, and spherical symmetry. The critical amplitude of the initial disturbance has been determined such that any compressive disturbance with initial amplitude greater than the critical one always grows into a shock wave, while the initial amplitude less than the critical one always decays.
Figures
Reference graph
Works this paper leans on
-
[1]
Morse, P. M. and Ingard, K.U. (1968) Theoretical Acoustic, New York: McGraw- Hill
work page 1968
-
[2]
Pascal, H., (1986) Pressure wave propagation in a fluid flowing through a porous medium and problem related to interpretation of Stoneley’s w ave attenuation in acoustical and well logging, I. J. Eng. Sc., 24, pp. 1553-1570
work page 1986
-
[3]
Nield, D. A. and Bejan, A., (1999) Convection in Porous Media, Berlin: Springer
work page 1999
-
[4]
Jordan, P.M., (2005) Growth and decay of acoustic acceleration waves in Darcy - type porous media, Proceed. Roy. Soc. A, 461, pp. 2749-2766
work page 2005
-
[5]
Ville, A.D., (1996) On the properties of compressible gas flow in a porous m edia, Transport in Porous Media, 22, pp. 287-306
work page 1996
-
[6]
Hsiao, L., (1999) Initial boundary value problem for the system of compressible adiabatic flow through porous media, J. Diff. Eq., 159, pp. 280-305
work page 1999
-
[7]
Pan, R., (2006) Darcy’s law as long -time limit of adiabatic porous media flow, J. Diff. Eq., 220, pp. 121-146
work page 2006
-
[8]
Jeffrey, A., (1976) Quasilinear hyperbolic system and wave, London
work page 1976
Show all 17 references
-
[9]
and Ruggeri, T., (1965) On the evolution law of weak discontinuity for hyperbolic quasilinear systems, Wave Motion, 1, pp.149-152
Boillat, G. and Ruggeri, T., (1965) On the evolution law of weak discontinuity for hyperbolic quasilinear systems, Wave Motion, 1, pp.149-152
1965
-
[10]
and Ferraioli, F., (1982) Evolution of weak discontinui ty waves in self- similar flows and form ation of s econdary sh ocks
Virgopia, N. and Ferraioli, F., (1982) Evolution of weak discontinui ty waves in self- similar flows and form ation of s econdary sh ocks. The point explosion model, J. Appl. Math. Phys. 33, pp. 63-80
1982
-
[11]
Ram, R ., (1978) Effect of radiative heat transfer on the growth and decay of acceleration waves, Appl. Sc. Res., 34, pp. 93-104
1978
-
[12]
D., (2007) Interaction of a characteristic shock with a weak discontinuity in a non-ideal gas, Wave Motion, 44, pp
Pandey, Manoj and Sharma, V. D., (2007) Interaction of a characteristic shock with a weak discontinuity in a non-ideal gas, Wave Motion, 44, pp. 346-354
2007
-
[13]
Sugiyama, M
Mentrelli, A., Ruggeri, T. Sugiyama, M. Zhao, N., (2008) Interaction between a shock and acceleration waves in a perfect gas for increasing shock strength, Wave Motion, 45, pp. 498–517
2008
-
[14]
and Husain, A., (2010) Propagation of nonlinear travelling waves in Darc y-type porous media, Acta Astronautica , 67, pp
Singh, M., Singh, L.P. and Husain, A., (2010) Propagation of nonlinear travelling waves in Darc y-type porous media, Acta Astronautica , 67, pp. 1053- 1058
2010
-
[15]
Safargulova, S.I., (1991) Propagation of weak disturbances in the combustion of compressible porous fuels, Combustion, Explosion and Shock Wave, 27 (2), pp
Smirnov, N.N. Safargulova, S.I., (1991) Propagation of weak disturbances in the combustion of compressible porous fuels, Combustion, Explosion and Shock Wave, 27 (2), pp. 26-34
1991
-
[16]
Courant, R., and Friedrichs, K.O., (1948) Supersonic Flow and Shock Waves, New York: Wiley Interscience
1948
-
[17]
0.0 0.5 1.0 1.5 2.0 4 2 0 2 4 6 Fig.1 Evolution of the amplitude of acceleration waves influenced by Darcy -type and non-Darcy type porous media for plane ( 0m = ) flows
Whitham, G.B., (1974) Linear and Nonlinear Waves, New York: Wiley. 0.0 0.5 1.0 1.5 2.0 4 2 0 2 4 6 Fig.1 Evolution of the amplitude of acceleration waves influenced by Darcy -type and non-Darcy type porous media for plane ( 0m = ) flows. 0.0 0.5 1.0 1.5 2.0 2.5 3.0 4 2 0 2 4 F...
1974
Reviewed May 24, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.