pith. sign in

arxiv: 2506.22225 · v2 · submitted 2025-06-27 · 🧮 math.AP

Existence and uniqueness of solutions of unsteady Darcy-Brinkman problem for modelling miscible reactive flows in porous media

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

classification 🧮 math.AP
keywords Darcy-Brinkman equationporous mediaweak solutionsfinite-time blow-upexistence and uniquenessreactive flowsKorteweg stressesadvection-reaction-diffusion
0
0 comments X

The pith

Weak solutions to a Darcy-Brinkman model of reactive porous media flows exist globally only if the initial concentration is at most 1.

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

The paper proves existence of weak solutions for a system that couples an advection-reaction-diffusion equation for solute concentration to an unsteady Darcy-Brinkman equation including Korteweg stresses, under Neumann boundary conditions for concentration and no-flow for velocity. When a second-order nonlinear reaction is present, the long-time behavior hinges on the initial concentration value: solutions remain global in time for initial concentrations between 0 and 1 but develop singularities in finite time when the initial concentration exceeds 1. Uniqueness of solutions is established specifically in two-dimensional domains. Numerical experiments using finite elements confirm the non-negativity of concentration, its decay over long times, and the occurrence of blow-up as predicted by the analysis.

Core claim

For the coupled system of advection-reaction-diffusion for concentration and unsteady Darcy-Brinkman for velocity with Korteweg stresses, weak solutions exist for broad initial data. With the nonlinear reaction, global existence holds when 0 ≤ C0 ≤ 1 while finite-time blow-up occurs for C0 > 1, and uniqueness is shown in 2D.

What carries the argument

The weak formulation of the coupled PDE system, analyzed via energy methods and compactness to establish existence, together with comparison principles for the concentration to determine global existence versus blow-up depending on initial value C0.

If this is right

  • If 0 ≤ C0 ≤ 1, the concentration stays non-negative and bounded for all future times.
  • For C0 > 1, a singularity forms in finite time, limiting the model's applicability to short-time predictions.
  • Solutions are unique in two dimensions, enabling precise tracking of flow and concentration evolution in planar settings.
  • The model applies to heterogeneous porous media such as fractured karst reservoirs where permeability varies spatially.

Where Pith is reading between the lines

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

  • The blow-up phenomenon for high initial concentrations may physically correspond to rapid solute accumulation or phase separation in real porous media applications.
  • Extending the uniqueness result to three dimensions would likely require stronger regularity on the permeability or additional dissipation terms.
  • Numerical schemes could be designed to adaptively detect the blow-up time for C0 slightly above 1 to validate the threshold.

Load-bearing premise

The permeability must be sufficiently smooth and the Korteweg stress must appear in the exact form included in the unsteady Darcy-Brinkman equation.

What would settle it

A computation starting with initial concentration C0 = 1.5 that remains bounded and smooth for arbitrarily long times would falsify the finite-time blow-up result.

read the original abstract

In this work, we investigate a model describing flow through porous media with permeability heterogeneity, combining an advection-reaction-diffusion equation for solute concentration with an unsteady Darcy-Brinkman equation with Korteweg stresses in the presence of external body forces for the flow field. Such models are appropriate in describing flows in fractured karst reservoirs, mineral wool, industrial foam, coastal mud, etc. These equations are coupled with Neumann boundary conditions for the solute concentration and no-flow conditions for the fluid velocity. For a broad class of initial data, we proved the existence of weak solutions. In the presence of a second-order nonlinear reaction, we show that the long-time behaviour of the solution depends on the initial concentration \(C_0\). More precisely, the solution exists for all time if \(0\leq C_0\leq 1\), and blows up at finite time if $C_0>1$. Furthermore, the uniqueness of the solution is proved for a two-dimensional domain. Finally, numerical simulations based on the finite element method have been presented that illustrate non-negativity of the concentration, long-time decay, and finite-time blow-up in agreement with theoretical estimates.

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

2 major / 2 minor

Summary. The manuscript proves existence of weak solutions to a coupled system of an advection-reaction-diffusion equation for solute concentration C and an unsteady Darcy-Brinkman equation with Korteweg stresses for the velocity field in heterogeneous porous media, subject to Neumann boundary conditions on C and no-flow conditions on velocity. For a quadratic nonlinear reaction, it establishes global existence when 0 ≤ C0 ≤ 1 and finite-time blow-up when C0 > 1; uniqueness is shown in two space dimensions; finite-element numerical simulations are presented to illustrate non-negativity, long-time decay, and blow-up.

Significance. If the estimates close rigorously, the results supply a mathematical basis for predicting global versus explosive long-time behavior in models of reactive flows through fractured or heterogeneous porous media, with direct relevance to applications such as karst reservoirs. The explicit dependence on initial concentration provides a clear, testable criterion, and the combination of analysis with numerical validation is a positive feature.

major comments (2)
  1. [Blow-up analysis (Theorem on finite-time blow-up)] Blow-up section (proof that solutions blow up for C0 > 1): the differential inequality for finite-time explosion is obtained by integrating the concentration equation against the constant test function 1. However, the existence theorem only yields weak solutions satisfying C ∈ L²(0,T;H¹(Ω)) ∩ L^∞(0,T;L²(Ω)) (without an a-priori L^∞ bound independent of T). It is therefore unclear whether the quadratic term C² remains integrable enough to justify the multiplication step or to pass to the limit inside the integrated form; an additional truncation or approximation argument appears necessary to close the estimate inside the exact function space delivered by the existence result.
  2. [Existence of weak solutions (Theorem 3.1)] Existence proof (Galerkin or compactness argument): after passage to the limit, the reaction term must lie in the dual of the test-function space for the weak formulation to be satisfied. The manuscript should explicitly verify that the quadratic nonlinearity maps the constructed solution class into this dual space uniformly in the approximation parameter, especially when the velocity field is only weakly regular from the Darcy-Brinkman equation.
minor comments (2)
  1. [Section 2 (model equations)] Notation for the Korteweg stress tensor and the precise form of the body-force term should be restated once in the weak formulation to avoid cross-referencing between the strong and weak statements.
  2. [Numerical section] Figure captions for the numerical blow-up examples should include the precise value of C0 used and the observed blow-up time to allow direct comparison with the theoretical criterion.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments. We address each major comment below and indicate the changes we will make to strengthen the rigor of the proofs.

read point-by-point responses
  1. Referee: [Blow-up analysis (Theorem on finite-time blow-up)] Blow-up section (proof that solutions blow up for C0 > 1): the differential inequality for finite-time explosion is obtained by integrating the concentration equation against the constant test function 1. However, the existence theorem only yields weak solutions satisfying C ∈ L²(0,T;H¹(Ω)) ∩ L^∞(0,T;L²(Ω)) (without an a-priori L^∞ bound independent of T). It is therefore unclear whether the quadratic term C² remains integrable enough to justify the multiplication step or to pass to the limit inside the integrated form; an additional truncation or approximation argument appears necessary to close the estimate inside the exact function space delivered by the existence result.

    Authors: We agree that the blow-up argument, as currently written, requires additional justification to rigorously handle the quadratic term within the weak-solution regularity class. We will revise the blow-up section by inserting a truncation (or mollification) procedure: we first derive the integrated inequality for a truncated version of the concentration, pass to the limit using the available integrability, and then obtain the differential inequality for the total mass. This closes the argument inside the precise function space delivered by the existence theorem and will be added to the revised manuscript. revision: yes

  2. Referee: [Existence of weak solutions (Theorem 3.1)] Existence proof (Galerkin or compactness argument): after passage to the limit, the reaction term must lie in the dual of the test-function space for the weak formulation to be satisfied. The manuscript should explicitly verify that the quadratic nonlinearity maps the constructed solution class into this dual space uniformly in the approximation parameter, especially when the velocity field is only weakly regular from the Darcy-Brinkman equation.

    Authors: We thank the referee for this observation. In the existence proof we use a Galerkin scheme and compactness to pass to the limit, but we acknowledge that an explicit verification of the reaction term’s membership in the dual space (uniformly with respect to the approximation parameter) is not stated in sufficient detail. We will add a short lemma after the a-priori estimates that confirms the quadratic nonlinearity maps the solution class into the dual of the test space, exploiting the Sobolev embeddings available from the Darcy-Brinkman regularity and the uniform bounds on the Galerkin approximations. This verification will be included in the revised version. revision: yes

Circularity Check

0 steps flagged

No circularity: direct PDE analysis with standard existence/uniqueness techniques

full rationale

The paper establishes existence of weak solutions for the coupled advection-reaction-diffusion and unsteady Darcy-Brinkman system via Galerkin approximation, a priori estimates, and compactness arguments applied to the given PDEs, boundary conditions, and initial data. The long-time behavior (global existence for 0 ≤ C0 ≤ 1 versus finite-time blow-up for C0 > 1 under quadratic reaction) follows from energy inequalities and comparison principles derived directly from the equations without parameter fitting or redefinition. Uniqueness in 2D is obtained from standard difference estimates. No load-bearing self-citations, ansatzes smuggled via prior work, or renamings of known results appear; the derivation chain remains self-contained against the stated model assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard functional-analytic assumptions for weak solutions of parabolic and Navier-Stokes-type systems; no free parameters or new postulated entities are introduced.

axioms (2)
  • standard math Existence of weak solutions is established via standard Galerkin approximation and compactness arguments in appropriate Sobolev spaces.
    Typical background framework for PDE existence proofs; invoked implicitly for the coupled system.
  • domain assumption The permeability heterogeneity is sufficiently regular to allow the Darcy-Brinkman operator to be well-defined.
    Required for the flow equation to make sense in the weak formulation.

pith-pipeline@v0.9.0 · 5740 in / 1610 out tokens · 37444 ms · 2026-05-19T07:59:09.709342+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

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

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages

  1. [1]

    (u, C) ∈ L2 (0, T; V1) , L2 (0, T; V2) ,

  2. [2]

    u (·, 0) = u0 and C (·, 0) = C0 a.e. in Ω,

  3. [3]

    on (0, T) , (10) ∂C (t) ∂t , B + d (∇C (t) , ∇B) + (u (t) · ∇C (t) , B) + κ (C (t) (1 − C (t)) , B) = 0 a.e

    For everyv ∈ V1 and for everyB ∈ V2, ∂u(t) ∂t , v + µe (∇u (t) , ∇v) + (F (C (t)) u (t) , v) − ⟨∇ · T (C (t)) , v⟩ + (f (t) , v) = 0 a.e. on (0, T) , (10) ∂C (t) ∂t , B + d (∇C (t) , ∇B) + (u (t) · ∇C (t) , B) + κ (C (t) (1 − C (t)) , B) = 0 a.e. on (0, T) . (11) B. Assumptions A1. The function F is assumed to be a mappingF : L2(0, T; V2) → L2(0, T; H1(Ω)...

  4. [4]

    We first construct approximate solutions using the semi-discrete Galerkin method

    A Priori Estimates The proof of Theorem III.1 is based on the well-known Faedo-Galerkin method, as detailed in the works of Lions [41], Temam [42, 43], and Ladyzhenskaya [44]. We first construct approximate solutions using the semi-discrete Galerkin method. We consider {wk}k∈N and {zk}k∈N to be basis ofV1 and V2, respectively, wherezk’s are eigen vectors ...

  5. [5]

    Consequently, we establish the uniform boundedness of∂Cn ∂t in L2 0, T; L2(Ω) . Lemma 4. ∂un ∂t is uniformly bounded inL2 (0, T; V ∗ 1 ). 8 Proof. Let w ∈ V1 be any non-zero vector. Then w can be written as w = w1 + w2 where w1 ∈ (V1)n , w2 ∈ (V1)⊥ n . It can be shown that ∂un (t) ∂t , w2 = 0. Then from equation (14) we get, ∂un (t) ∂t , w = ∂un (t) ∂t , ...

  6. [6]

    ApplyingAubin-Lionslemma( [46, Corollary 4])wegetstronglyconvergentsubsequencesof (un, Cn)in L2 (0, T; S1)× L2 (0, T; S2)

    Existence Now from the results we obtained in Lemma 1–4 approaching similarly as [20, 37, 38] we obtain uniformly convergent subsequence of un, Cn (again using the same notation), Cn ⇀ C in L 2 (0, T; V2) , Cn ∗ ⇀ C in L ∞ (0, T; S2) , ∂Cn ∂t ⇀ ∂C ∂t in L 2 (0, T; V ∗ 2 ) , un ⇀ u in L 2 (0, T; V1) , un ∗ ⇀ u in L ∞ (0, T; S1) , ∂un ∂t ⇀ ∂u ∂t in L 2 (0, ...

  7. [7]

    We define,u = u1 − u2, C = C1 − C2

    Uniqueness Let (u1, C1) and (u2, C2) be two solutions. We define,u = u1 − u2, C = C1 − C2. Then we have the following system of equations, ∂C (t) ∂t , B + d (∇C (t) , ∇B) + (u1 (t) · ∇C1 (t) − u2 (t) · ∇C2 (t) , B) +κ (C (t) , B) − κ ((C1 (t) + C2 (t)) C (t) , B) = 0 ∀B ∈ V2 a.e. on (0, T) , (42) ∂u (t) ∂t , v + µe (∇u (t) , ∇v) + (F (C1 (t)) u1 (t) − F (...

  8. [8]

    Positiveness of Solute Concentration We now show that if the initial concentrationC0 is positive almost everywhere onΩ then the solute concentration C remains non negative almost everywhere onΩ × (0, T). For that we takeC − = max {0, −C} as test function in equation (11) and obtain the following equation: 1 2 d dt ∥C −(t)∥2 L2 + d∥∇C −∥2 L2 = −κ∥C −(t)∥2 ...

  9. [9]

    Darcy, Les fontaines publiques de la ville de dijon, (1856)

    H. Darcy, Les fontaines publiques de la ville de dijon, (1856)

  10. [10]

    H. E. Huppert and J. A. Neufeld, The fluid mechanics of carbon dioxide sequestration, Annual Review of Fluid Mechanics 46, 255 (2014)

  11. [11]

    C. E. Neuzil, Groundwater flow in low-permeability environments, Water Resources Research22, 1163 (1986)

  12. [12]

    G. M. Homsy, Viscous fingering in porous media, Annual Review of Fluid Mechanics19, 271 (1987)

  13. [13]

    Zheng, K

    X. Zheng, K. Zhao, T. Jackson, and J. Lowengrub, Tumor growth towards lower extracellular matrix conductivity regions under Darcy’s law and steady morphology, Journal of mathematical biology85, 5 (2022)

  14. [14]

    F. J. Molz, M. Widdowson, and L. Benefield, Simulation of microbial growth dynamics coupled to nutrient and oxygen transport in porous media, Water Resources Research22, 1207 (1986)

  15. [15]

    Khalifa, I

    M. Khalifa, I. Wahyudi, and P. Thomas, New extension of Darcy’s law to unsteady flows, Soils and Foundations42, 53 (2002)

  16. [16]

    McCurdy, N

    M. McCurdy, N. Moore, and X. Wang, Convection in a coupled free flow-porous media system, SIAM Journal on Applied Mathematics 79, 2313 (2019)

  17. [17]

    H. C. Brinkman, A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles, Flow, Turbulence and Combustion 1, 27 (1949)

  18. [18]

    Vafai,Handbook of Porous Media(CRC Press, 2015) p

    K. Vafai,Handbook of Porous Media(CRC Press, 2015) p. 110

  19. [19]

    F. A. Morales and R. E. Showalter, A Darcy–Brinkman model of fractures in porous media, Journal of Mathematical Analysis and Applications452, 1332 (2017)

  20. [20]

    Ebenbeck and H

    M. Ebenbeck and H. Garcke, On a Cahn–Hilliard–Brinkman model for tumor growth and its singular limits, SIAM Journal on Mathematical Analysis51, 1868 (2019)

  21. [21]

    W. J. Layton, F. Schieweck, and I. Yotov, Coupling fluid flow with porous media flow, SIAM Journal on Numerical Analysis 40, 2195 (2002)

  22. [22]

    K. G. Link, M. G. Sorrells, N. A. Danes, K. B. Neeves, K. Leiderman, and A. L. Fogelson, A mathematical model of platelet aggregation in an extravascular injury under flow, Multiscale Modeling & Simulation18, 1489 (2020)

  23. [23]

    Babuška and G

    I. Babuška and G. N. Gatica, A residual-based a posteriori error estimator for the Stokes–Darcy coupled problem, SIAM Journal on Numerical Analysis48, 498 (2010)

  24. [24]

    Cimolin and M

    F. Cimolin and M. Discacciati, Navier–Stokes/Forchheimer models for filtration through porous media, Applied Numerical Mathematics 72, 205 (2013)

  25. [25]

    D. L. Koch and R. J. Hill, Inertial effects in suspension and porous-media flows, Annual Review of Fluid Mechanics33, 619 (2001)

  26. [26]

    S. B. Naqvi, D. Śnieżek, D. Strzelczyk, M. Mądrala, and M. Matyka, Inertial effects on fluid flow through low-porosity media, arXiv preprint arXiv:2505.10418 (2025). 16

  27. [27]

    Strzelczyk, G

    D. Strzelczyk, G. Kosec, and M. Matyka, The role of porosity in the transition to inertial regime in porous media flows, arXiv preprint arXiv:2505.10087 (2025)

  28. [28]

    Kundu, S

    S. Kundu, S. N. Maharana, and M. Mishra, Existence and uniqueness of solution to unsteady Darcy-Brinkman problem with Korteweg stress for modelling miscible porous media flow, Journal of Mathematical Analysis and Applications539, 128532 (2024)

  29. [29]

    J. Hou, M. Qiu, X. He, C. Guo, M. Wei, and B. Bai, A dual-porosity-Stokes model and finite element method for coupling dual-porosity flow and free flow, SIAM Journal on Scientific Computing 38, B710 (2016), https://doi.org/10.1137/15M1044072

  30. [30]

    Evje and H

    S. Evje and H. Wen, A Stokes two-fluid model for cell migration that can account for physical cues in the microenvironment, SIAM Journal on Mathematical Analysis50, 86 (2018), https://doi.org/10.1137/16M1078185

  31. [31]

    Carbonell and S

    R. Carbonell and S. Whitaker, Heat and mass transfer in porous media (1984) pp. 121–198

  32. [32]

    J. M. Delgado,Heat and mass transfer in porous media, Vol. 13 (Springer Science & Business Media, 2011)

  33. [33]

    Pramanik and M

    S. Pramanik and M. Mishra, Linear stability analysis of Korteweg stresses effect on miscible viscous fingering in porous media, Physics of Fluids25 (2013)

  34. [34]

    Dehghan and Z

    M. Dehghan and Z. Gharibi, A fully mixed virtual element method for Darcy–Forchheimer miscible displacement of incompressible fluids appearing in porous media, IMA Journal of Numerical Analysis44, 797 (2024)

  35. [35]

    Caucao, G

    S. Caucao, G. N. Gatica, R. Oyarzúa, and N. Sánchez, A fully-mixed formulation for the steady double-diffusive convection system based upon Brinkman–Forchheimer equations, Journal of Scientific Computing85, 44 (2020)

  36. [36]

    Caucao, G

    S. Caucao, G. N. Gatica, R. Oyarzúa, and P. Zúñiga, A posteriori error analysis of a mixed finite element method for the coupled Brinkman–Forchheimer and double-diffusion equations, Journal of Scientific Computing93, 50 (2022)

  37. [37]

    D. Korteweg, On the form which the equations of motion of fluids take if we take into account the capillary forces caused by considerable but known variations of density and on the theory of capillarity in the hypothesis of a continuous variation of the density, Dutch Arch. Exact Nat. Sci6, 1 (1901)

  38. [38]

    C.-Y. Chen, L. Wang, and E. Meiburg, Miscible droplets in a porous medium and the effects of Korteweg stresses, Physics of Fluids 13, 2447 (2001)

  39. [39]

    Swernath, B

    S. Swernath, B. Malengier, and S. Pushpavanam, Effect of Korteweg stress on viscous fingering of solute plugs in a porous medium, Chemical Engineering Science65, 2284 (2010)

  40. [40]

    De Wit, Chemo-hydrodynamic patterns in porous media, Philosophical Transactions of the Royal Society A: Mathe- matical, Physical and Engineering Sciences374, 20150419 (2016)

    A. De Wit, Chemo-hydrodynamic patterns in porous media, Philosophical Transactions of the Royal Society A: Mathe- matical, Physical and Engineering Sciences374, 20150419 (2016)

  41. [41]

    C. Rana, S. Pramanik, M. Martin, A. De Wit, and M. Mishra, Influence of Langmuir adsorption and viscous fingering on transport of finite size samples in porous media, Physical Review Fluids4, 104001 (2019)

  42. [42]

    De Wit, Fingering of chemical fronts in porous media, Phys

    A. De Wit, Fingering of chemical fronts in porous media, Phys. Rev. Lett.87, 054502 (2001)

  43. [43]

    Nagatsu, Y

    Y. Nagatsu, Y. Ishii, Y. Tada, and A. De Wit, Hydrodynamic fingering instability induced by a precipitation reaction, Physical Review Letters113, 024502 (2014)

  44. [44]

    Migórski and P

    S. Migórski and P. Szafraniec, Nonmonotone slip problem for miscible liquids, Journal of Mathematical Analysis and Applications 471, 342 (2019)

  45. [45]

    Allali, Existence and uniqueness of solutions to a model describing miscible liquids, Comptes Rendus

    K. Allali, Existence and uniqueness of solutions to a model describing miscible liquids, Comptes Rendus. Mathématique 355, 1148 (2017)

  46. [46]

    Kostin, M

    I. Kostin, M. Marion, R. Texier-Picard, and V. A. Volpert, Modelling of miscible liquids with the Korteweg stress, ESAIM: Mathematical Modelling and Numerical Analysis37, 741 (2003)

  47. [47]

    D. D. Joseph, Fluid dynamics of two miscible liquids with diffusion and gradient stresses, European Journal of Mechanics - B/Fluids 9, 565 (1990)

  48. [48]

    L. C. Evans,Partial Differential Equations, Graduate studies in mathematics (American Mathematical Society, 2010)

  49. [49]

    J. L. Lions,Quelques méthodes de résolution des problèmes aux limites non linéaires(Gauthier-Villars, Paris, 1969)

  50. [50]

    Temam,Navier-Stokes Equations

    R. Temam,Navier-Stokes Equations. Theory and Numerical Analysis, Studies in Mathematics and its Applications, Vol. 2 (North-Holland Publishing Co., Amsterdam–New York, 1979)

  51. [51]

    Temam,Navier-Stokes Equations and Nonlinear Functional Analysis(SIAM, 1983)

    R. Temam,Navier-Stokes Equations and Nonlinear Functional Analysis(SIAM, 1983)

  52. [52]

    O. A. Ladyzhenskaya,Mathematical Theory of Viscous Incompressible Flow(Gordon and Breach, 1963)

  53. [53]

    Kesavan,Topics in Functional Analysis and Applications(New Age International Pvt Limited Publishers, 2019)

    S. Kesavan,Topics in Functional Analysis and Applications(New Age International Pvt Limited Publishers, 2019)

  54. [54]

    Simon, Compact sets in the spaceLp(0, T ; B), Annali di Matematica Pura ed Applicata146, 65 (1986)

    J. Simon, Compact sets in the spaceLp(0, T ; B), Annali di Matematica Pura ed Applicata146, 65 (1986)

  55. [55]

    Boyer and P

    F. Boyer and P. Fabrie,Mathematical Tools for the Study of the Incompressible Navier-Stokes Equations and Related Models (Springer Science & Business Media, 2012)

  56. [56]

    L. C. Morrow, C. Cuttle, and C. W. MacMinn, Gas compression systematically delays the onset of viscous fingering, Physical Review Letters131, 224002 (2023)