{"id":"b6b97ddb-5779-4f54-8a03-25014b845a22","arxiv_id":"1908.09316","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"A symmetry classification and explicit non-stationary similarity solutions are derived for isentropic gas filtration in porous media, with phase-transition curves for a van der Waals gas.","lead":"This paper uses symmetry methods to find exact time-dependent solutions for gas flow through porous rocks. It classifies when such flows have extra symmetries and gives explicit pressure and temperature profiles for special gas models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 4 as printed is internally inconsistent: the reduced ODE gives p_r = -r/(2μ), but the theorem states p = -½∫ μ dr; the Gaussian example confirms the missing factor r.","rationale":"The most load-bearing issue is the false formula in Theorem 4, since that theorem is the explicit invariant solution underpinning Section 4 and the phase-transition study. The paper's own examples show the intended corrected relation, so the flaw is localized and fixable, not a collapse of the symmetry classification. The reader correctly noticed that the Gaussian profile contradicts the printed integral relation, and this is enough to keep the verdict conditional. However, the reader's stated weakest assumption is wrong: the solution is not obtained on the constant-entropy branch. The second factor 2Tv_r + nvT_r = 0 is the constant-entropy branch, while the solution satisfies the first factor (non-constant entropy branch). The paper's note about the second factor does not imply that branch was chosen. Thus the reader's main criticism should be replaced by the concrete missing-r inconsistency. A single algebraic substitution into the reduced ODE settles the matter, so the verdict should remain CONDITIONAL pending correction of Theorem 4 and verification of the Maple computations behind Theorems 2 and 3.","tokens_in":5234,"tokens_out":27181,"duration_ms":265250,"concrete_test":"Re-derive the first reduced ODE from Eqs. (1)–(3) using the g4-invariant ansatz with r = √((x^2+y^2+z^2)/t). Substitute the theorem's printed relation p_r = −μ/2 and the proposed v = RC1 r^{3/(1−q)} into the first factor 2Rμ(v_rT − vT_r) − rv^2 = 0; the left-hand side fails to vanish unless r = 1, while the corrected relation p_r = −r/(2μ) makes it vanish identically. Equivalently, for μ = αv/T, compute p from the printed integral (p^2 = C − αR r) and compare with the paper's Gaussian Example (p = C2 exp(−r^2/(4αR))); the mismatch settles the missing factor.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"On the non-constant-entropy branch (first factor of the first reduced ODE), the invariant solution satisfies 2Rμ(v_rT − vT_r) − rv^2 = 0. Since p = RT/v, this is equivalent to p_r = −r/(2μ), hence p(r) = C − ½∫^r r μ(v,T) dr. Theorem 4 prints p(r) = −½∫^r μ dr, omitting the factor r. This is not cosmetic: for μ = αv/T (ideal gas, so v/T = R/p), the printed relation gives p′ = −αR/(2p), whose solution is p^2 = C − αR r; the paper's own Example gives p = C2 exp(−r^2/(4αR)), which satisfies p′ = −r p/(2αR) = −r/(2μ). Thus the theorem statement is false as written, while the four listed cases agree with the corrected relation. The reader's claim that the solution is built on the constant-entropy branch is also incorrect: the constant-entropy branch is the second factor 2Tv_r + nvT_r = 0, whereas Theorem 4's solution satisfies the first factor. The central solution construction is therefore recoverable, but the printed theorem needs correction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":5523,"tokens_out":7191,"duration_ms":74827,"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":[{"comment":"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.","section":"§4.1, Theorem 4"},{"comment":"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.","section":"§3 and §4.1"},{"comment":"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.","section":"§4.1, entropy branch"}],"minor_comments":[{"comment":"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.","section":"Abstract and Section 1"},{"comment":"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":"Theorem 4"},{"comment":"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":"Section 3, table"},{"comment":"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.","section":"Section 4.2"},{"comment":"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.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The factor-r error in Theorem 4 is a fixable typo, but it currently makes the theorem false as stated, so the paper needs revision rather than acceptance. The reliance on the authors' earlier thermodynamic formalism [3,4] is natural and not a cause for concern; the symmetry computation itself is independent of that formalism. The main editorial question is whether deferring all derivations to external Maple files is acceptable for this journal; if not, that should be communicated to the authors."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nWhat 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.\n\nThe 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.\n\nThe 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.\n\nTwo 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.\n\nWho 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.","headline":"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.","tokens_in":6018,"tokens_out":2014,"would_cite":false,"duration_ms":22822,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B06","76S05","58J70"],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["gas filtration","porous media","isentropic flow","Lie point symmetries","invariant solutions","van der Waals gas","phase transitions","Massieu-Planck potential"],"falsifier":"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.","tokens_in":5038,"feed_emoji":"💨","tokens_out":7371,"duration_ms":67147,"temperature":0.7,"pith_summary":"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.","feed_headline":"Gas filtration symmetries reduced to eight vector fields","feed_subtitle":"Same eight space-time motions govern every porous medium; special gas-medium pairs add one or two symmetries and yield exact solutions.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the standard Darcy-flow and conservation-law model that the PDE system $E$ is built from.","marker":"[1]"},{"why":"The computational differential-geometry package used to obtain the quoted determining equations and reductions.","marker":"[2]"},{"why":"Provides the contact-geometric thermodynamics and the Massieu-Planck phase-equilibrium equations used throughout.","marker":"[3]"},{"why":"The steady-filtration predecessor whose symmetry classification the paper extends to non-stationary processes.","marker":"[4]"}],"fun_headline_variants":["Eight symmetries govern every porous-medium gas flow","Gas filtration symmetries fully classified; ideal gases add two","Exact self-similar solution for special gas in porous media","Symmetry algebra of gas filtration: eight fields plus extras","Phase transitions in new exact gas filtration solution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Eight symmetries govern every porous-medium gas flow","Gas filtration symmetries fully classified; ideal gases add two","Exact self-similar solution for special gas in porous media","Symmetry algebra of gas filtration: eight fields plus extras","Phase transitions in new exact gas filtration solution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000621,"raw_usage":{"total_tokens":2870,"prompt_tokens":931,"completion_tokens":1939,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":1862}},"tokens_in":547,"tokens_out":1939,"duration_ms":15698,"temperature":1.0,"reasoning_tokens":1862,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:16:33.932203+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"The physics of flow through porous media","cited_arxiv_id":null,"evidence_quote":"Supplies the standard Darcy-flow and conservation-law model that the PDE system $E$ is built from."},{"cited_title":"and Torre, Charles G., The Dif- ferential Geometry Package (2016)","cited_arxiv_id":null,"evidence_quote":"The computational differential-geometry package used to obtain the quoted determining equations and reductions."},{"cited_title":"Proceedings of the Wisla 18 Summer School (2019) 354, Springer Nature, Switzerland","cited_arxiv_id":null,"evidence_quote":"Provides the contact-geometric thermodynamics and the Massieu-Planck phase-equilibrium equations used throughout."}],"review_version":1}