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
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.
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
- 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.
Referee Report
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)
- [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.
- [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)
- [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.
- [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
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
-
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
-
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
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
axioms (2)
- standard math Existence of weak solutions is established via standard Galerkin approximation and compactness arguments in appropriate Sobolev spaces.
- domain assumption The permeability heterogeneity is sufficiently regular to allow the Darcy-Brinkman operator to be well-defined.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We establish the existence and uniqueness of a weak solution of reactive miscible flows in porous media... viscosity and permeability... positive functions of the concentration... second-order chemical reaction.
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
-
[1]
(u, C) ∈ L2 (0, T; V1) , L2 (0, T; V2) ,
-
[2]
u (·, 0) = u0 and C (·, 0) = C0 a.e. in Ω,
-
[3]
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]
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]
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]
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]
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]
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 ...
work page 2021
-
[9]
Darcy, Les fontaines publiques de la ville de dijon, (1856)
H. Darcy, Les fontaines publiques de la ville de dijon, (1856)
-
[10]
H. E. Huppert and J. A. Neufeld, The fluid mechanics of carbon dioxide sequestration, Annual Review of Fluid Mechanics 46, 255 (2014)
work page 2014
-
[11]
C. E. Neuzil, Groundwater flow in low-permeability environments, Water Resources Research22, 1163 (1986)
work page 1986
-
[12]
G. M. Homsy, Viscous fingering in porous media, Annual Review of Fluid Mechanics19, 271 (1987)
work page 1987
- [13]
-
[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)
work page 1986
-
[15]
M. Khalifa, I. Wahyudi, and P. Thomas, New extension of Darcy’s law to unsteady flows, Soils and Foundations42, 53 (2002)
work page 2002
-
[16]
M. McCurdy, N. Moore, and X. Wang, Convection in a coupled free flow-porous media system, SIAM Journal on Applied Mathematics 79, 2313 (2019)
work page 2019
-
[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)
work page 1949
-
[18]
Vafai,Handbook of Porous Media(CRC Press, 2015) p
K. Vafai,Handbook of Porous Media(CRC Press, 2015) p. 110
work page 2015
-
[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)
work page 2017
-
[20]
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)
work page 2019
-
[21]
W. J. Layton, F. Schieweck, and I. Yotov, Coupling fluid flow with porous media flow, SIAM Journal on Numerical Analysis 40, 2195 (2002)
work page 2002
-
[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)
work page 2020
-
[23]
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)
work page 2010
-
[24]
F. Cimolin and M. Discacciati, Navier–Stokes/Forchheimer models for filtration through porous media, Applied Numerical Mathematics 72, 205 (2013)
work page 2013
-
[25]
D. L. Koch and R. J. Hill, Inertial effects in suspension and porous-media flows, Annual Review of Fluid Mechanics33, 619 (2001)
work page 2001
- [26]
-
[27]
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]
-
[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]
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]
R. Carbonell and S. Whitaker, Heat and mass transfer in porous media (1984) pp. 121–198
work page 1984
-
[32]
J. M. Delgado,Heat and mass transfer in porous media, Vol. 13 (Springer Science & Business Media, 2011)
work page 2011
-
[33]
S. Pramanik and M. Mishra, Linear stability analysis of Korteweg stresses effect on miscible viscous fingering in porous media, Physics of Fluids25 (2013)
work page 2013
-
[34]
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)
work page 2024
- [35]
- [36]
-
[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)
work page 1901
-
[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)
work page 2001
-
[39]
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)
work page 2010
-
[40]
A. De Wit, Chemo-hydrodynamic patterns in porous media, Philosophical Transactions of the Royal Society A: Mathe- matical, Physical and Engineering Sciences374, 20150419 (2016)
work page 2016
-
[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)
work page 2019
-
[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)
work page 2001
-
[43]
Y. Nagatsu, Y. Ishii, Y. Tada, and A. De Wit, Hydrodynamic fingering instability induced by a precipitation reaction, Physical Review Letters113, 024502 (2014)
work page 2014
-
[44]
S. Migórski and P. Szafraniec, Nonmonotone slip problem for miscible liquids, Journal of Mathematical Analysis and Applications 471, 342 (2019)
work page 2019
-
[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)
work page 2017
- [46]
-
[47]
D. D. Joseph, Fluid dynamics of two miscible liquids with diffusion and gradient stresses, European Journal of Mechanics - B/Fluids 9, 565 (1990)
work page 1990
-
[48]
L. C. Evans,Partial Differential Equations, Graduate studies in mathematics (American Mathematical Society, 2010)
work page 2010
-
[49]
J. L. Lions,Quelques méthodes de résolution des problèmes aux limites non linéaires(Gauthier-Villars, Paris, 1969)
work page 1969
-
[50]
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)
work page 1979
-
[51]
Temam,Navier-Stokes Equations and Nonlinear Functional Analysis(SIAM, 1983)
R. Temam,Navier-Stokes Equations and Nonlinear Functional Analysis(SIAM, 1983)
work page 1983
-
[52]
O. A. Ladyzhenskaya,Mathematical Theory of Viscous Incompressible Flow(Gordon and Breach, 1963)
work page 1963
-
[53]
S. Kesavan,Topics in Functional Analysis and Applications(New Age International Pvt Limited Publishers, 2019)
work page 2019
-
[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)
work page 1986
-
[55]
F. Boyer and P. Fabrie,Mathematical Tools for the Study of the Incompressible Navier-Stokes Equations and Related Models (Springer Science & Business Media, 2012)
work page 2012
-
[56]
L. C. Morrow, C. Cuttle, and C. W. MacMinn, Gas compression systematically delays the onset of viscous fingering, Physical Review Letters131, 224002 (2023)
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.