REVIEW 2 major objections 3 minor 24 references
Two-dimensional equilibrium configurations in Korteweg fluids
T0 review · 2 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read For a two-dimensional Korteweg fluid at constant temperature, the paper shows that all mechanical equilibrium configurations are governed by a single nonlinear elliptic equation that identically satisfies the overdetermined equilibrium…
desk verdict Honest, compact reduction of 2D Korteweg equilibrium to one scalar PDE; the result is probably right but the key algebra is asserted and the numerics contain a boundary mismatch. 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
The load-bearing object is the compatibility condition (3.3), $\alpha_3^2-\alpha_1\,\partial\alpha_3/\partial\rho+2\alpha_2\alpha_3=0$, combined with the separated equilibrium entropy (3.4), $s_0(\rho,\varepsilon)=s_{01}(\rho)+s_{02}(\varepsilon)$. When both hold, the two equilibrium equations in (3.2) are no longer independent: they become derivatives of a single scalar expression, so the system reduces to equation (3.6). That single nonlinear elliptic equation is the mechanism carrying the argument; its dimensionless version (3.9) is the object actually solved, and the power-law specialization (3.7) turns it into a boundary-value problem amenable to finite-difference computation.
What would settle it
Numerically solve the dimensionless equilibrium equation (3.9) for a fixed parameter set on a grid, substitute the result into both components of the original equilibrium system (3.2), and evaluate the residuals; the claim that (3.6) identically satisfies the system stands or falls on whether those residuals vanish at the level of the numerical tolerance.
Extended reading notes
Core claim
The central discovery is a reduction: for a third-grade Korteweg fluid at constant temperature $\theta_0$ in two dimensions, the overdetermined equilibrium system (3.2) collapses to the single scalar equation $2\rho s_1(\rho)(\rho_{xx}+\rho_{yy}) + \frac{d(\rho s_1)}{d\rho}(\rho_x^2+\rho_y^2) - \frac{d(\rho s_{01})}{d\rho} + \frac{g}{\theta_0}y - \kappa = 0$, where $s_1(\rho)\le 0$ is the coefficient of $|\nabla\rho|^2$ in the entropy and $s_{01}(\rho)$ is the density-dependent part of the equilibrium entropy. Any sufficiently regular solution of this equation automatically satisfies both components of (3.2). The reduction relies on the compatibility condition (3.3), which the thermodynamically derived constitutive relations satisfy precisely under the separated-entropy hypothesis (3.4). With the power-law choices $s_{01}=\kappa_1\rho^m$ and $s_1=-\kappa_2\rho^n$, the dimensionless form (3.9) is obtained; for $m=1,n=-1$ it becomes a Poisson equation, and for $m=-1,n=-1$ a Laplace equation, both admitting separable analytical solutions.
Load-bearing premise
The derivation requires that the equilibrium entropy split as $s_0(\rho,\varepsilon)=s_{01}(\rho)+s_{02}(\varepsilon)$ and that the equilibrium temperature be constant; if a real fluid has a non-separable equilibrium entropy or a non-uniform temperature field, equation (3.6) does not describe its equilibrium density.
Editorial extensions
If this is right
- Equilibrium density fields in two dimensions can be obtained by solving one scalar elliptic equation with boundary data instead of an overdetermined pair of PDEs.
- For the parameter choices $m=1$, $n=-1$ and $m=-1$, $n=-1$, the reduced equation is linear and becomes respectively a Poisson equation and a Laplace equation, so analytical separable solutions exist.
- Because the compatibility condition (3.3) is met, the classical restriction of equilibrium phase boundaries to spherical, cylindrical, or planar shapes is lifted, making more general interface geometries admissible.
- The preliminary finite-difference solutions of the Dirichlet problem show that the reduced equation is numerically tractable in both linear and nonlinear regimes.
- The single-equation formulation provides a starting point for three-dimensional equilibria and for comparisons with laboratory experiments, both listed by the authors as future work.
Reading between the lines
- The reduction is demonstrated only for power-law choices of $s_{01}$ and $s_1$; if the compatibility condition is the only essential ingredient, the same collapse to a single equation should occur for other functional forms, which could be tested by solving (3.6) with non-power-law data and checking the residual of (3.2).
- The separability assumption likely marks the boundary of the single-equation description: for non-separable equilibrium entropy, the overdetermined system should generically admit only the simple phase-boundary geometries, so a physical test would be to look for non-spherical equilibrium interfaces in a Korteweg fluid with a known non-separable entropy.
- The constant-temperature hypothesis fixes $\theta_0$; allowing $\theta=\theta(x,y)$ gives two unknowns and two equations, so the authors' stated future problem of non-isothermal stationary solutions is a different regime in which the single-equation shortcut may not hold.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies mechanical equilibrium of a third-grade Korteweg fluid in two dimensions at constant temperature. After reviewing constitutive relations from a prior extended-Liu-procedure paper [20], it assumes that the equilibrium entropy separates as s0=s01(ρ)+s02(ε) and asserts that the overdetermined equilibrium system (3.2) reduces to the single scalar PDE (3.6) (or its dimensionless form (3.9)). This claim is the paper's central result. The paper then analyzes linear special cases, proposes a Dirichlet boundary value problem, and reports preliminary numerical solutions obtained with finite differences and a Matlab solver.
Significance. If the reduction is correct, the paper provides a practical tool: equilibrium density fields of Korteweg fluids satisfying Serrin's geometric condition can be computed by solving one elliptic PDE instead of an overdetermined system. Strengths include the explicit use of constitutive inputs from [20], a consistent nondimensionalization, and correct linear-case reductions. The authors are also appropriately cautious, labeling the numerics as preliminary and deferring three-dimensional and non-uniform-temperature problems. The main weaknesses are the missing derivation of the central reduction and an inconsistent numerical boundary condition.
major comments (2)
- [Section 3, Eq. (3.6)] The central reduction from the overdetermined equilibrium system (3.2) to the single scalar equation (3.6) is stated without derivation; the text between (3.5) and (3.6) says only 'after simple algebraic manipulations' and introduces an arbitrary constant κ. Since the paper's main claim is precisely that every solution of (3.6) identically satisfies (3.2), the authors should provide the complete derivation, including the role of the Serrin condition (3.3), the handling of the arbitrary function of y that arises upon integrating the x-component of (3.2), and any regularity or sign assumptions (ρ>0, s1≠0) needed for divisions. Without this, the central equivalence is an assertion rather than a demonstrated result.
- [Section 3.1, Eq. (3.11) and Figures 1-2] The Dirichlet data in (3.11) are discontinuous at the upper corners: at (0,d) and (1,d), the top boundary condition u(x,d)=ρ0−x^2(1−x)^2 gives u=ρ0, while the side conditions u(0,y)=u(1,y)=(ρ1−ρ0)y/d+ρ0 give u=ρ1 at y=d. Since ρ0≠ρ1, no continuous solution to the stated boundary value problem exists, so the numerical solutions in Figures 1 and 2 cannot be solutions of the stated problem. The boundary data should be made compatible (or the corner singularities should be discussed and treated) and the computations repeated.
minor comments (3)
- [Section 3.1, before Eq. (3.11)] The boundary data are written in terms of u(x,0), u(x,d), u(0,y), u(1,y), but the unknown in (3.11) is ρ; the notation should be made consistent.
- [Section 3.1, text after Eq. (3.10)] The equation ρxx+ρyy+2αρ=0 is a Helmholtz equation, not a Poisson equation, since it contains the zero-order term 2αρ.
- [Section 3.1, text before Eq. (3.12)] The text says first and second derivatives are approximated by second-order and fourth-order finite differences, respectively, but the formulas displayed in (3.12) are both standard second-order central differences.
Circularity Check
No significant circularity: the scalar equilibrium equation is a genuine reduction, not a disguised restatement of the paper's inputs.
full rationale
The central claim—that every solution of the scalar PDE (3.6) identically satisfies the overdetermined equilibrium system (3.2)—is a mathematical reduction, not a definition or a fit. The paper combines the equilibrium condition (3.1)-(3.2) with the constitutive relations (2.5) inherited from [20] and the separability condition (3.4), and it introduces the constant κ as an integration constant. No equation is defined in terms of the target result, and no parameter is fitted to data and then renamed as a prediction: the forms s01=κ1ρ^m and s1=-κ2ρ^n in (3.7) are free choices, and the numerical solutions are explicitly labeled preliminary. The self-citation to [20] is essential input, but [20] derives the constitutive functions from an extended Liu procedure that does not include the equilibrium PDE (3.6); under the review rules this is independent support, not circularity. The Serrin/Pucci theorems are external, not self-citations. The paper's acknowledged limitations—constant temperature θ0, separable s0, and deferral of non-uniform temperature to future work—are scope restrictions, not circular reasoning. Two non-circular weaknesses are flagged: the phrase 'after simple algebraic manipulations' omits the details of the reduction to (3.6), and the Dirichlet data in (3.11) are inconsistent at the corners when ρ0≠ρ1; both are correctness/transparency concerns, not circularity.
Assumptions & free parameters
free parameters (9)
- m =
1 or -1 in linear cases; 1 in nonlinear examples
- n =
-1 in linear cases; 0, 1, -2, -3 in nonlinear examples
- κ1 =
not numerically specified
- κ2 =
not numerically specified
- α =
1.0 in all numerical plots
- β =
-1.2 in all numerical plots
- γ =
±1 in numerical plots
- ρ0, ρ1 =
ρ0=1.4, ρ1=1.3 in numerical plots
- d =
3 in numerical plots
assumptions (4)
- domain assumption Constitutive relations (2.5) from [20] are accepted as correct.
- standard math Serrin's compatibility condition (3.3) is needed for non-simple phase boundaries.
- ad hoc to paper The equilibrium entropy separates as s0=s01(ρ)+s02(ε).
- domain assumption Temperature is constant at equilibrium (θ0).
Cite this review
Pith. "Pith review of Two-dimensional equilibrium configurations in Korteweg fluids." pith.science (2026). https://pith.science/paper/NALLJ4NZ
@misc{pith2026250522026,
author = {Pith},
title = {Pith review of: Two-dimensional equilibrium configurations in Korteweg fluids},
year = {2026},
howpublished = {\url{https://pith.science/paper/NALLJ4NZ}},
note = {Machine review of arXiv:2505.22026}
}
read the original abstract
In this paper, after reviewing the form of the constitutive equations for a third grade Korteweg fluid, recently derived by means of an extended Liu procedure, an equilibrium problem is investigated. By considering a two--dimensional setting, it is derived a single nonlinear elliptic equation such that the equilibrium conditions are identically satisfied. Such an equation is discussed both analytically and numerically. Moreover, by considering a particular boundary value problem of Dirichlet type, some preliminary numerical solutions are presented.
Figures
Reference graph
Works this paper leans on
-
[20]
M. Gorgone, P. Rogolino, On the characterization of constitutive equa- tions for third-grade viscous Korteweg fluids , Phys. Fluids 33 (2021), 043107
work page 2021
-
[1]
D. J. Korteweg, Sur la forme qui prennent les ´ equations du mouvement des fluids si l’on tient compte des forces capillaires par des variations de densit´ e consid´ erables mais continues et sur la th´ eorie de la capil- larit´ e dans l’hypoth` ese d’une variation continue de la densit´ e, Archives N´ eerlandaises des sciences exactes et naturelles 6 (190...
work page 1901
- [2]
-
[3]
J. E. Dunn, K. R. Rajagopal, Fluids of differential type: critical review and thermodynamic analysis , Int. J. Eng. Sci. 33 (1995), 689–729
work page 1995
-
[4]
C. Truesdell, K. R. Rajagopal, An Introduction to the Mechanics of Fluids, Birkh¨ auser, Boston–Basel–Berlin, 2000
work page 2000
-
[5]
J. E. Dunn, J. Serrin, On the thermomechanics of the interstitial work- ing, Arch. Rat. Mech. Anal. 88 (1985), 95–133
work page 1985
-
[6]
J. E. Dunn, Interstitial working and a nonclassical continuum thermo- dynamics, in: J. Serrin (ed.), New Perspectives in Thermodynamics , Springer Verlag, Berlin, 1986, 11, 187–222
work page 1986
-
[7]
Truesdell, Rational Thermodynamics
C. Truesdell, Rational Thermodynamics. 2nd ed. , Springer Verlag, Berlin, 1984
work page 1984
Show all 24 references
-
[8]
M¨ uller,On the Entropy Inequality , Arch
I. M¨ uller,On the Entropy Inequality , Arch. Rat. Mech. Anal. 26 (1967), 118–141
1967
-
[9]
V. A. Cimmelli, A. Sellitto, V. Triani, A new thermodynamic framework for second–grade Korteweg–type viscous fluids, J. Math. Phys. 50 (2009), 053101
2009
-
[10]
V. A. Cimmelli, A. Sellitto, V. Triani, A new perspective on the form of the first and second laws in rational thermodynamics: Kortew eg fluids as an example , J. Non–Equil. Thermodyn. 35 (2010), 251–265. 13
2010
-
[11]
V. A. Cimmelli, F. Oliveri, A. R. Pace, On the thermodynamics of Ko- rteweg fluids with heat conduction and viscosity , J. Elast. 104 (2011), 115–131
2011
-
[12]
V. A. Cimmelli, F. Oliveri, A. R. Pace, Thermodynamical setting for gradient continuum theories with vectorial internal varia bles: application to granular materials , Int. J. Non–Linear Mech. 49 (2013), 72–76
2013
-
[13]
V. A. Cimmelli, F. Oliveri, A. R. Pace, A nonlocal phase–field model of Ginzburg–Landau–Korteweg fluids , Contin. Mech. Thermodyn. 27 (2015), 367–378
2015
-
[14]
Oliveri, A
F. Oliveri, A. Palumbo, P. Rogolino, On a model of mixtures with inter- nal variables: extended Liu procedure for the exploitation of the entropy principle, Atti Accad. Pelorit. Pericol. Cl. Sci. Fis. Mat. Nat. 94 (2016), A2–1–17
2016
-
[15]
V. A. Cimmelli, M. Gorgone, F. Oliveri, A. R. Pace, Weakly nonlocal thermodynamics of binary mixtures of Korteweg fluids with tw o velocities and two temperatures , Eur. J. Mech. B/Fluids 83 (2020), 58–65
2020
-
[16]
Gorgone, F
M. Gorgone, F. Oliveri, P. Rogolino, Continua with non–local constitu- tive laws: Exploitation of entropy inequality , Int. J. Non–Linear Mech. 126 (2020), 103573
2020
-
[17]
Gorgone, F
M. Gorgone, F. Oliveri, P. Rogolino, Thermodynamical analysis and constitutive equations for a mixture of viscous Korteweg flu ids, Phys. Fluids 33 (2021), 093102
2021
-
[18]
Gorgone, P
M. Gorgone, P. Rogolino, A thermodynamical description of third grade fluid mixtures , J. Non–Equil. Thermodyn. 47 (2022), 133–142
2022
-
[19]
I–Shih Liu, Method of Lagrange multipliers for exploitation of the en- tropy principle , Arch. Rat. Mech. Anal. 46 (1972), 131–148
1972
-
[21]
V. A. Cimmelli, An extension of Liu procedure in weakly nonlocal ther- modynamics, J. Math. Phys. 48 (2007), 113510. 14
2007
-
[22]
Serrin, The form of interfacial surfaces in Korteweg’s theory of pha se equilibria, Quart
J. Serrin, The form of interfacial surfaces in Korteweg’s theory of pha se equilibria, Quart. Appl. Math. 41 (1983), 357–364
1983
-
[23]
Pucci, An overdetermined system , Quart
P. Pucci, An overdetermined system , Quart. Appl. Math. 41 (1983), 365–367
1983
-
[24]
A. R. Mitchell, D. F. Griffiths, The Finite Difference Method in Partial Differential Equations , John Wiley & Sons, New York, 1980. 15
1980
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.