Non-stationary adiabatic filtration of gases in porous media

Anna Duyunova; Sergey Tychkov; Valentin Lychagin

REVIEW 3 major objections 5 minor 5 references

Non-stationary adiabatic filtration of gases in porous media

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper computes the full point symmetry algebra for non-stationary isentropic gas filtration in porous media and gives explicit invariant solutions for ideal and van der Waals gases.

desk verdict A solid, contained symmetry-classification paper whose printed Theorem 4 has a missing-factor typo that must be fixed, but the underlying reduction and solutions look recoverable and worth peer review. read the letter →

arxiv 1908.09316 v1 pith:Q53E5KUE submitted 2019-08-25 math-ph math.MP

Anna Duyunova , Valentin Lychagin , Sergey Tychkov This is my paper · ORCID

classification math-phmath.MP MSC 35B0676S0558J70

keywords gasfiltrationporousmediaisentropicflowLiepointsymmetriesinvariantsolutionsvanderWaalsphasetransitionsMassieu-Planckpotential

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper treats non-stationary filtration of gases through porous media as a PDE system combining Darcy's law, mass conservation, entropy transport, and thermodynamic equations of state. It establishes that for any porous medium the point symmetry algebra is generated by eight vector fields—translations in space and time, rotations, and one scaling—and that for ideal gases one or two additional symmetries appear exactly when the transport coefficient $\mu(v,T)$ has one of five listed forms. The authors then reduce the system by the rotation-scaling subalgebra $g_4$, obtain an explicit $g_4$-invariant solution for ideal gases, build first-order van der Waals corrections to it, and locate the phase-transition curves for a worked methane example. A sympathetic reader would care because explicit non-stationary solutions of filtration equations are rare, and the symmetry classification gives a principled way to choose media for which exact or asymptotic solutions exist.

What carries the argument

The load-bearing object is the PDE system $E$ consisting of Darcy's law $u=-\mu(v,T)\,\mathrm{grad}\,p$, the mass balance $q v_t+u\cdot\mathrm{grad}\,v=v\,\mathrm{div}\,u$, entropy transport $s_t+u\cdot\mathrm{grad}\,s=0$, and the Massieu-Planck equations of state $p=RT\varphi_v$, $\epsilon=RT^2\varphi_T$, $s=R(\varphi+T\varphi_T)$. Point symmetries are computed by requiring the second prolongation of a vector field to be tangent to $E^{(2)}$; the resulting Lie algebra is the classification device. For invariant solutions the paper uses the subalgebra $g_4=\mathrm{so}(3)\oplus\langle X_8\rangle$, whose orbits are three-dimensional and whose invariant is $r^2=(x^2+y^2+z^2)/t$; reduction by $g_4$ turns the PDE system into an ODE system that can be integrated on the constant-entropy factor. The phase-transition analysis uses the Massieu-Planck coexistence equations $\varphi_v(v_1,T)=\varphi_v(v_2,T)$ and $\varphi(v_2,T)-\varphi(v_1,T)+v_1\varphi_v(v_1,T)-v_2\varphi_v(v_2,T)=0$.

What would settle it

Substitute the Theorem 4 solution for a non-constant entropy profile $s(t,x)$ into the original system $E$; if the equation $s_t+u\cdot\mathrm{grad}\,s=0$ holds for some nonconstant $s$, the constant-entropy restriction is not necessary, whereas if it forces $\mathrm{grad}\,s=0$, the solution class is exactly isentropic. Alternatively, run a symmetry computation for $\mu(v,T)=v^2+T^2$; if the algebra exceeds eight dimensions, the 'arbitrary porous medium' part of Theorem 2 fails.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is a complete classification of the point symmetries of the isentropic filtration system $E$. Theorem 2 states that for an arbitrary porous medium the symmetry Lie algebra is generated by $X_1=\partial_x$, $X_2=\partial_y$, $X_3=\partial_z$, $X_4=\partial_t$, the three rotation fields $X_5,X_6,X_7$, and the scaling field $X_8=2t\partial_t+x\partial_x+y\partial_y+z\partial_z$. Theorem 3 states that for ideal gases, with state potential $\varphi(v,T)=\frac{n}{2}\ln T+\ln v$, the algebra is enlarged by one or two extra symmetries precisely for the listed special forms of $\mu(v,T)$; for example $\mu(v,T)=f(v)T^\alpha$ admits $X_9=(1+\alpha)t\partial_t-T\partial_T$. Theorem 4 gives the $g_4$-invariant ideal-gas solution $v(r)=RC_1 r^{3/(1-q)}$, $p(r)=-\frac12\int^r \mu(v,T)\,dr$, $T(r)=p(r)v(r)/R$, where $r=\sqrt{x^2+y^2+z^2}/t$. The reduction notes that the second factor in the first reduced ODE corresponds to constant entropy, so the explicit solution lives on the isentropic branch. The paper also constructs first-order corrections for a van der Waals gas and studies where the resulting solution is physically admissible.

Load-bearing premise

The explicit invariant solutions are derived on the branch where specific entropy is constant, so the answer only covers globally isentropic flows; if entropy varies, the formulas need not satisfy the full adiabatic system.

Editorial extensions

If this is right

For an arbitrary porous medium, any point symmetry of the isentropic filtration system is a combination of space-time translations, rotations, and the scaling $X_8$; no other point symmetries exist unless the medium satisfies one of the special conditions.
For ideal gases, the admissible extra symmetries are completely classified by the form of $\mu(v,T)$, so the symmetry algebra can be read off directly from the transport coefficient.
The explicit $g_4$-invariant solution provides closed-form pressure, temperature, and volume profiles for media with $\mu=\alpha(T/v)^\beta$, $\mu=\alpha v/T$, $\mu=\alpha v^\beta T^\gamma$, and $\mu=\alpha v^\beta/T$.
First-order van der Waals corrections to the ideal-gas solution are obtained by solving a linear system of ODEs, giving concrete temperature corrections $T_1(r)$ and $T_2(r)$ for the worked example.
Phase-transition curves for the constructed solution can be plotted on the distance-time plane, and in the methane example they lie outside the region where the solution is physically admissible.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The same reduction-by-symmetry strategy could be applied to the non-isentropic case by replacing $s_t+u\cdot \mathrm{grad}\,s=0$ with a full energy equation; the extra symmetries would likely depend on the form of the heat-flux law, giving a parallel classification.
Because the explicit solution is built on the constant-entropy branch, the title's 'adiabatic' framing is broader than the proved result; a testable next step is to integrate the other factor of the reduced ODE system to see whether nonconstant-entropy invariant solutions exist and how they differ.
The list of special $\mu(v,T)$ forms may serve as a diagnostic: if experimental pressure profiles exhibit one of the extra scaling symmetries, that constrains the functional form of permeability and viscosity, effectively using symmetry as a measurement tool.
The virial-asymptotic construction for van der Waals gases suggests a recursion in the virial coefficients $A_k(T)$; truncating at higher order would produce corrections for more realistic equations of state, not just the $a,b$ model.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Bottom line: this is a legitimate extension of the steady-state analysis in [4] to non-stationary isentropic filtration, and it gives explicit invariant solutions that should be useful as benchmark cases for numerical simulators. It deserves a real referee, but only after the printed Theorem 4 is corrected.

What is genuinely new: the 8-dimensional universal symmetry algebra for arbitrary porous media, the classification of extra symmetries for special ideal-gas media, and the g4-invariant reduction to a solvable ODE system leading to closed-form pressure profiles for several viscosity models. The van der Waals first-order corrections and the phase-transition analysis for the time-dependent solution are also not in the earlier steady-state paper. I want to give credit where it is earned: the reduction idea is clean, and the explicit profiles for the four special cases look like they came from a real computation, not from fitting.

The soft spots are real but localized. Most importantly, Theorem 4 as stated is internally inconsistent. From the first factor of the reduced ODE system you get p_r = -r/(2µ), which means p(r) = C - (1/2)∫ r µ dr. The theorem prints p(r) = -(1/2)∫ µ dr, and the Gaussian example µ = αv/T confirms the missing factor r: their p = C2 exp(-r²/(4αR)) satisfies p' = -r p/(2αR) = -r/(2µ), not p' = -µ/2. The four listed cases seem to match the corrected relation, so this is likely a typo in the theorem statement rather than a fatal flaw, but it has to be fixed.

The reader's claim that the solution sits on the constant-entropy branch is wrong. The constant-entropy branch is the second factor, 2T v_r + n v T_r = 0; the theorem's solution uses the first factor. So the concern should be redirected to the printed formula, not the branch choice.

Two smaller issues. The symmetry classification is outsourced to Maple files; that is acceptable practice, but the paper should make clear which files contain which computations, and the authors should confirm the files are actually available at the cited URL. Also, the abstract says classification of media "with respect to admissible symmetries," but the full classification is only for ideal gases; arbitrary media get only the universal algebra. That is a scope overstatement, easy to fix in both abstract and introduction.

Who is this for? Researchers working on gas filtration, porous media flows, or applications of symmetry reduction to PDE systems. It is not a landmark, but it is a competent technical contribution. My recommendation: engage with it, send it to peer review, and require the authors to fix Theorem 4 and clarify the classification scope. I would not cite it in my own next paper, but I would not be embarrassed to have it in the literature once corrected.

Referee Report

3 major / 5 minor

Summary. This paper studies the PDE system governing non-stationary adiabatic filtration of gases in a porous medium of constant porosity, comprising Darcy's law, mass conservation, entropy advection along the flow, and equations of state expressed through a Massieu-Planck potential. The authors compute the point-symmetry algebra for the general system (Theorem 2), classify ideal-gas media according to the additional admissible symmetries (Theorem 3 and the table in Section 3), and construct a rotation- and scaling-invariant (g4) family of self-similar solutions for ideal gases, with explicit pressure profiles for several viscosity/permeability models (Theorem 4 and Section 4.1). They then compute first-order van der Waals corrections to this solution and analyze phase coexistence curves on the (distance,time) plane. The computations are supported by Maple files referenced at d-omega.org.

Significance. If the symmetry classification is correct, the paper gives a useful and fairly complete Lie-symmetry analysis for a physically motivated nonlinear filtration model, including a media classification and explicit invariant solutions with estimates of their physical domain of applicability. A notable strength is that the reduction and classification are computational and reproducible through the cited Maple files, and no free parameters are fitted to force the advertised solutions. The main limitation is that the central invariant-solution formula, Theorem 4, is misstated as printed; the paper should not be evaluated on the basis of that formula until it is corrected.

major comments (3)

[§4.1, Theorem 4] The printed formula p(r) = -1/2 ∫^r μ(v,T) dr is not the integral of the reduced ODE. Setting the first factor of the first reduced equation to zero gives 2Rμ(v_rT - vT_r) - r v^2 = 0; since p = RT/v, this is equivalent to p_r = -r/(2μ), so the correct relation is p(r) = C - (1/2)∫^r rμ(v,T) dr. The missing factor r is not cosmetic: for μ = αv/T = αR/p (ideal gas), the printed relation would give p^2 = C - αRr, whereas the paper's own case 2 gives p = C2 exp(-r^2/(4αR)), which satisfies the corrected equation p_r = -r/(2μ). The four listed cases are consistent with the corrected formula; the theorem and the surrounding discussion must be amended.
[§3 and §4.1] The completeness of the symmetry algebra (Theorem 2), the classification table (Theorem 3), and the invariant reduction leading to the ODE system (Theorem 4) are asserted with the proof deferred entirely to a Maple file. Since these results are the central claims of the paper, the authors should include at least the determining equations and a reproducible derivation sketch in the text, or state the exact Maple commands and output that verify the classification, so a reader does not have to trust an external computation blindly.
[§4.1, entropy branch] The invariant solution of Theorem 4 is obtained from the first factor of the first reduced equation, not from the second factor. The sentence 'Note that the second factor in the first equation corresponds the case when the entropy s is constant' is correct, but the paper should explicitly state that the constructed solution lies on the first branch, where entropy varies along the flow. This distinction matters for what the solution represents relative to the title's 'adiabatic' wording and should not be left implicit.

minor comments (5)

[Abstract and Section 1] The terms 'isentropic' and 'adiabatic' are used interchangeably, but the equation s_t + u·grad s = 0 only expresses entropy conservation along particle trajectories; please align the terminology.
[Theorem 4] The constants C1 and C2 are introduced without explicit definitions, and 'C ∈ R' should read 'C1 ∈ R'; the status of R as the specific gas constant should be stated.
[Section 3, table] The classification table would be more complete if each row stated the domain of the function (e.g. β ≠ -1, positive α, q ≠ 1) and if a sentence explained why these cases exhaust the possibilities for the ideal-gas model.
[Section 4.2] The linear system for the first-order van der Waals corrections is not displayed in the text; please either display the system or identify the specific Maple file and computation so that the printed expressions for T1 and T2 can be checked.
[Throughout] There are several grammatical and typographical issues (e.g. 'This gives us understanding when the solution is applicable') that should be corrected in a final revision.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: symmetry computation and invariant solution are derived from the explicit PDE system without fitted inputs.

full rationale

The paper's claimed results — the symmetry algebra in Theorems 2 and 3, and the g4-invariant solution in Theorem 4 — are obtained by direct prolongation computation and by reduction of the explicit PDE system (1)-(3) with equations of state (6) to an ODE system, followed by explicit integration. No parameter is fitted to a data subset and then renamed a prediction; the constants C, C1, C2, C3, C4 are free integration constants. The author-overlap citations [3], [4], [5] are used only to supply the Massieu-Planck representation of thermodynamic states and the phase-coexistence equations; these are exogenous inputs to the model, not consequences of the symmetry or solution results, and no uniqueness theorem or prior ansatz is invoked to force the final answer. I therefore find no circular step. (The printed p-formula in Theorem 4 appears to omit a factor r and to that extent is internally inconsistent, but that is a correctness defect, not a circularity.)

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The central claim uses no curve-fitted constants: C1, C2, C3, C4 are integration constants, and q, α, a, b are material or gas parameters assigned in the example. The main external inputs are the thermodynamic formalism from the authors' prior papers and a black-box Maple computation, which are the axioms listed above.

assumptions (5)

domain assumption The Darcy law, mass conservation, and entropy advection equations (1)-(3) with constant porosity q constitute a valid model for non-stationary isentropic filtration.
Section 1 states these as the governing equations; the entire symmetry analysis is conditional on them.
domain assumption Thermodynamic states of real gases are Legendrian/Lagrangian manifolds represented by a Massieu-Planck potential φ(v,T) through equations (6).
Section 2 imports Theorem 1 from the authors' prior work [3] and uses the virial expansion (7) without proving it here.
domain assumption For ideal gas, φ=(n/2)lnT+lnv; for van der Waals gas, φ=(n/2)lnT+ln(v-b)+a/(RvT)+o2 approximates the potential.
Section 4 builds the explicit solutions on these potentials, restricting the results to those gas models.
ad hoc to paper The Maple DifferentialGeometry package computes the determining equations for the symmetry algebra correctly and completely.
The proofs of Theorems 2 and 3 are deferred to Maple files; no independent derivation appears in the text.
domain assumption The first-order asymptotic expansion in a and b (equation (10)) with the ideal-gas solution as the zero-order term is a valid approximation for the van der Waals system.
Section 4.2 assumes a and b are small and that first-order corrections capture the real-gas behavior; convergence is not analyzed.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Non-stationary adiabatic filtration of gases in porous media." pith.science (2026). https://pith.science/paper/Q53E5KUE

@misc{pith2026190809316,
  author       = {Pith},
  title        = {Pith review of: Non-stationary adiabatic filtration of gases in porous media},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Q53E5KUE}},
  note         = {Machine review of arXiv:1908.09316}
}

read the original abstract

A non-stationary isentropic filtration of gases in porous media is considered. Thermodynamics in terms of contact and symplectic geometries is briefly discussed. Algebra of symmetries for the PDE system is found, and the classification of media with respect to admissible symmetries is given. Solution for one class of media is found and the phase transitions for this solution are studied. Keywords: gas filtration, porous media, symmetry algebra, phase

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

5 extracted references · 4 canonical work pages

[4]

arXiv:1903.00276

Lychagin V., Roop M., Phase transitions in ﬁltration of real gases (2019). arXiv:1903.00276

arXiv 2019
[1]

The physics of ﬂow through porous media

Scheidegger A. The physics of ﬂow through porous media. Revised edition, The Macmillan Co., New York, 1960. 10 Figure 4

work page 1960
[2]

and Torre, Charles G., The Dif- ferential Geometry Package (2016)

Anderson, Ian M. and Torre, Charles G., The Dif- ferential Geometry Package (2016). Downloads. Paper 4. http://digitalcommons.usu.edu/dg downloads/4

work page 2016
[3]

Proceedings of the Wisla 18 Summer School (2019) 354, Springer Nature, Switzerland

Lychagin V., Contact Geometry, Measurement and Thermody- namics, in: Nonlinear PDEs, Their Geometry and Applications. Proceedings of the Wisla 18 Summer School (2019) 354, Springer Nature, Switzerland

work page 2019
[5]

Diﬀerential invariants for plane ﬂows of viscid ﬂuids

Duyunova A., Lychagin V., Tychkov S. Diﬀerential invariants for plane ﬂows of viscid ﬂuids. Lobachevskii Journal of Mathematics, 2017, Vol. 38, No. 4, 644-652. 11

work page 2017

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [4]

arXiv:1903.00276

Lychagin V., Roop M., Phase transitions in ﬁltration of real gases (2019). arXiv:1903.00276

arXiv 2019

[2] [1]

The physics of ﬂow through porous media

Scheidegger A. The physics of ﬂow through porous media. Revised edition, The Macmillan Co., New York, 1960. 10 Figure 4

work page 1960

[3] [2]

and Torre, Charles G., The Dif- ferential Geometry Package (2016)

Anderson, Ian M. and Torre, Charles G., The Dif- ferential Geometry Package (2016). Downloads. Paper 4. http://digitalcommons.usu.edu/dg downloads/4

work page 2016

[4] [3]

Proceedings of the Wisla 18 Summer School (2019) 354, Springer Nature, Switzerland

Lychagin V., Contact Geometry, Measurement and Thermody- namics, in: Nonlinear PDEs, Their Geometry and Applications. Proceedings of the Wisla 18 Summer School (2019) 354, Springer Nature, Switzerland

work page 2019

[5] [5]

Diﬀerential invariants for plane ﬂows of viscid ﬂuids

Duyunova A., Lychagin V., Tychkov S. Diﬀerential invariants for plane ﬂows of viscid ﬂuids. Lobachevskii Journal of Mathematics, 2017, Vol. 38, No. 4, 644-652. 11

work page 2017