Pith. sign in

REVIEW 3 major objections 5 minor 24 references

Compressibility and volume variations due to composition in multicomponent fluids

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

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

Pith's one-line read For a multicomponent fluid, zero isothermal compressibility does not force a constant density: the paper defines $k=\beta(T,\mu_t)-\beta(T,y)$, built from derivatives of the Gibbs function, which measures composition-driven volume change…

desk verdict A clean derivation of a new composition-compressibility coefficient, with honest but model-dependent numerics that need stability checks. read the letter →

arxiv 2411.17742 v2 pith:ODVQUV42 submitted 2024-11-24 physics.flu-dyn

Pierre-Etienne Druet This is my paper · ORCID
classification physics.flu-dyn MSC 76T3076N9980A1735Q35
keywords multicomponentfluidcompressibilityGibbsfunctionincompressiblelimitconstantdensityapproximationvolumeofmixingchemicalpotentialaqueoussolutions
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

For a single-component fluid, zero isothermal compressibility forces a constant density. The paper shows that the same conclusion does not hold for mixtures: even when $\beta(T,y)=0$, the specific volume can still change because composition changes. To make that effect quantitative, the paper defines a coefficient $k=\beta(T,\mu_t)-\beta(T,y)$, expressed through derivatives of the Gibbs function, which treats volume change on mixing as a true compressibility. The ratio $k/\beta(T,y)$ then gives a criterion for when the constant-density approximation is justified. Three aqueous solutions, salt water, sucrose water, and ethanol-water, illustrate the sizes of this ratio in practice.

What carries the argument

The load-bearing object is the coefficient $k=\beta(T,\mu_t)-\beta(T,y)$, constructed from the Gibbs function $g(T,p,y)$ and its tangential compositional Hessian $\mathrm{D}^2_{yy}g$. Explicitly, $k=\partial_p g\,(\mathrm{D}^2_{yy}g)^{-1}(\partial_y v/v)\cdot(\partial_y v/v)$, with derivatives taken along the hypersurface $\sum_i y_i=1$. Because $\beta(T,\mu_t)$ is the pressure derivative of volume at fixed relative chemical potentials, $k$ has the same dimension and meaning as a compressibility; the inverse Hessian converts the composition response into this pressure-like measure. The derivation also uses the entropic-variable transformation of the full mixture equations to identify $\beta(T,\mu_t)$, and the ideal-mixture limit yields the closed form $k=\frac{1}{RT}\sum_i M_i\rho_i(1/\varrho-v_i)^2$, which makes the vanishing cases, a pure fluid or equal specific volumes, explicit.

What would settle it

Take a binary mixture such as sucrose-water, measure density and isothermal compressibility over a range of temperature, pressure, and composition, construct the Gibbs function from these data, then compare $\beta(T,\mu_t)-\beta(T,y)$ from eq. (55) with the measured difference between compressibility at fixed relative chemical potential and at fixed composition; a systematic mismatch would show the claimed identity does not hold.

Watch

Extended reading notes

Core claim

The central claim is that composition-induced volume changes in a multicomponent fluid can and should be quantified as a compressibility. If $v(T,p,y)$ is the specific volume and $\beta(T,y)=-(1/v)\partial_p v$ is the isothermal compressibility at fixed composition, the paper proves that the compressibility measured at fixed relative chemical potentials is $\beta(T,\mu_t)=\beta(T,y)+k$, where $k$ is defined from the Gibbs function by $k=\partial_p g\,(\mathrm{D}^2_{yy}g)^{-1}(\partial_y v/v)\cdot(\partial_y v/v)$. Thus $k$ is the compressibility of mixing. In the incompressible limit $\beta(T,y)\to 0$ the coefficient $k$ need not vanish, so an incompressible mixture can still have a density that varies with composition; the density is constant only when $k$ also vanishes. For ideal mixtures the paper derives the explicit formula $\beta(T,\mu_t)-\beta(T,y)=\frac{1}{RT}\sum_i M_i\rho_i(1/\varrho-v_i)^2$, which vanishes only for a pure fluid or when all species have the same specific volume. The ratio $k/\beta(T,y)$ is proposed as the quantitative criterion for the validity of the constant-density approximation.

Load-bearing premise

The derivation assumes the mixture is thermodynamically stable, so the curvature of the Gibbs function with respect to composition is positive and invertible; near a spinodal or unstable state this inverse fails, $k$ diverges, and the split into $\beta(T,y)$ and $k$ no longer describes a physical compressibility.

Editorial extensions

If this is right

  • In a multicomponent fluid with $\beta(T,y)=0$, the density can still vary with composition through $k$; only $k=0$ gives a truly constant density.
  • The ratio $k/\beta(T,y)$ is a quantitative test for the constant-density approximation: when it is large, composition effects dominate and the approximation fails.
  • For ideal mixtures, $k$ is nonnegative and vanishes only for a pure fluid or when all species have identical specific volumes, so mixing alone always contributes a positive compressibility-like term.
  • The incompressible limit splits into two distinct regimes: $\beta(T,y)\to 0$ with the constraint $\sum_i \bar v_i\rho_i=1$ (variable density), versus $\beta(T,\mu_t)\to 0$ with $\varrho=$ constant.
  • In the aqueous examples, the ratio exceeds about 6 for concentrated salt water and 10 for sucrose water, while for ethanol-water it stays near order one, suggesting the constant-density approximation is unsafe for the first two but may be acceptable for the third.

Reading between the lines

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

  • [Editorial inference] The same coefficient $k$ could be extracted from tabulated mixture-density data by finite differencing, even without a full Gibbs function, which would make the constant-density test routinely available for engineering mixtures.
  • [Editorial inference] Because the appendix shows the compressive part of the pressure gradient is controlled by $\beta(T,\mu_t)$ rather than $\beta(T,y)$, sound-speed measurements in a mixture could provide a direct experimental route to $k$.
  • [Editorial inference] In a ternary reactive model such as the ethanol-water one, $k$ absorbs the volume-of-mixing defect into an effective species, so $k$ could serve as a general measure of non-ideal volume behaviour in mixtures with chemical association.
  • [Editorial inference] Since $k$ involves the inverse compositional Hessian, it should diverge as a mixture approaches a spinodal, potentially making $k$ an early-warning indicator of phase instability; the paper does not pursue this.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper introduces a coefficient k defined as the difference between the isothermal compressibility at fixed chemical potentials, β(T,μ_t), and the usual isothermal, isocompositional compressibility β(T,y). The central identity is β(T,μ_t) = β(T,y) + k, where k is expressed via the Gibbs function and the inverse of its compositional Hessian (Eqs. (2), (54), (55)). The author proves this identity by transforming variables from (T,p,y) to (T,p,μ_t) and argues that in the incompressible limit β(T,y)→0, k can remain nonzero, so an incompressible multicomponent fluid need not have constant density. The ratio k/β(T,y) is proposed as a quantitative criterion for the validity of the constant-density approximation. Three aqueous solutions (seawater, sucrose, and ethanol) are used as illustrations, with the ethanol case treated by both a binary ideal model and a ternary ideal model with a chemical equilibrium constraint (Appendix C).

Significance. If the decomposition is valid and can be computed reliably, the paper provides a useful conceptual distinction between 'incompressible' and 'constant density' for multicomponent fluids, with a concrete criterion for when the constant-density approximation fails. The calculus derivation in §3.2 is transparent and internally consistent, and the explicit ideal-mixture formula (57) is a clean closed-form result. The paper is also honest in acknowledging that wider data sets and full Gibbs functions are needed (conclusion). However, the numerical support is not yet sufficient to establish the quantitative claims: the seawater ideal model has a 27.9% error in the very compressibility used in the ratio, the sucrose and ethanol ideal models have errors of 7.2% and 35% respectively, and the ternary water-ethanol model is fitted with only two parameters and its thermodynamic stability is never verified. Thus the theoretical decomposition is sound, but the practical claim that k/β(T,y) provides a quantitative criterion for real multicomponent fluids is not established.

major comments (3)
  1. [§3.2, Eqs. (53)–(55)] The decomposition β(T,μ_t) = β(T,y)+k requires the compositional Hessian D²_yy g to be positive definite (hence invertible) on the relevant domain. The paper states this stability condition in §2.1 but never verifies it for the fitted models used in §4. In the ternary water-ethanol model of Appendix C, g is an effective Gibbs function obtained after imposing chemical equilibrium (Eq. (79)); the effective reduced Hessian is not computed, and the parameters D(T), D′(T) are fitted to one reference density and one reference compressibility (Eqs. (81)–(82)). If an eigenvalue of D²_yy g approaches zero or the effective model is not convex on the plotted composition range, the inverse in Eqs. (54)–(55) fails or becomes indefinite, so k loses its interpretation as a physical compressibility and the ratios shown in Figs. 5–6 become artefacts. This issue is load-bearing because the proposed criterion k/β(T,y) presupposes that k is well defined.
  2. [§4.1, Eq. (65) and following] The ideal binary model for seawater has relative error Err(β(T,y)) = 0.279 in the isothermal compressibility, yet this model is used to compute k = β(q)−β(T,y) and to compare it with β(T,y). Figure 2 shows that for high salinity the ideal model yields a k about six times larger than β(T,y), while the Gibbs seawater library attains a significantly smaller maximum. With a 27.9% error in the very quantity that defines the ratio k/β(T,y), the quantitative estimates for seawater are not reliable. The paper should either compute k from the full Gibbs function provided by the gsw library (as it already does for comparison) or provide uncertainty bounds before drawing conclusions about the incompressible-limit criterion.
  3. [§4.2 and §4.3] The same reliability issue affects the sucrose and ethanol examples. For sucrose, Err(β(T,y)) = 0.072 (§4.2); for the binary ethanol model, Err(β(T,y)) = 0.35 and Err(ϱ) = 0.14 (§4.3). The ternary ethanol model improves the fits but is calibrated with only two parameters D(T) and D′(T) to one reference density and one reference compressibility (Eqs. (81)–(82)). Consequently, the reported maximum ratios (more than 10 for sucrose, slightly above 1 for ethanol) are model outcomes, not validated physical predictions. The paper's own caveats in §4.1 ('Further investigations will be necessary') and in the conclusion ('wider data sets and ... full Gibbs functions are available') are appropriate, but they are in tension with the abstract's statement that the phenomenon becomes quantitatively comparable through k.
minor comments (5)
  1. [Eq. (2) and §2.1] The notation D²_{y,y}g is used in Eq. (2) before the tangential derivative is defined in §2.1; please define the notation at first use.
  2. [Page 5, first paragraph] The phrase 'completely equivalent ot (8)' contains a typo: 'ot' should be 'to'.
  3. [§2.1, paragraph on entropy] The sentence 'We let ϱs be the bulk density of the mixture entropy possesses of the special form' is grammatically awkward and should be rephrased.
  4. [§4.1, Eq. (65)] The error metric Err(·) is defined as the maximum absolute deviation normalized by the data range; this makes the values Err(ϱ)=0.029 and Err(β)=0.279 difficult to interpret. Please also report a typical or root-mean-square error, which would give a clearer picture of the model's accuracy.
  5. [Appendix B] The wave-equation analysis is only linearized and heuristic; the statement that diffusion processes 'accelerate' compression waves (around Eq. (73)) should be qualified as a linearized, near-equilibrium result.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central identity is a self-contained calculus derivation, and the numerical illustrations are openly labelled fitted models rather than disguised predictions.

full rationale

The paper's central claim, Eq. (3) together with Eqs. (53)–(55), is derived in Section 3.2 as a chain-rule identity relating derivatives at fixed composition y to derivatives at fixed chemical-potential projection μt. The derivation uses only the Gibbs-function identities (20)–(22), the definition of β(T,y) in (25), and the positivity of the compositional Hessian as a stability assumption. No fitted coefficient enters the derivation, and no step in the paper reduces the target result to its own input by construction. The numerical sections fit ideal-model parameters to reference density or compressibility data, but the paper explicitly presents these as approximate ideal models and compares the resulting k values against external libraries (seawater, Gibbs library) and experimental data. Where disagreements occur, they are reported honestly, e.g. the relative error Err(β(T,y)) = 0.279 for seawater and the statement that the Gibbs-library value attains a significantly smaller maximum than the ideal model. The same honest limitation appears in the conclusion, which calls for wider data sets and full Gibbs functions. The unresolved question of whether D²_yy g remains positive definite for the fitted ternary water-ethanol model is a validity condition and a correctness risk, not circularity: the paper does not use a positive-definite check as evidence for its central identity. Self-citations to [4], [5], [9], and [10] supply framework and earlier mathematical lemmas, but the key decomposition is re-derived in this paper and does not depend on the target result being assumed. Therefore the derivation chain is self-contained, and the fitted examples do not amount to predictions forced by the fit.

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

The central theoretical result is derived from a standard Gibbs-function thermodynamic framework; no new entities are postulated. The empirical sections introduce four fitted model parameters (two rho_S values and two pressure-deviation functions) and rely on ideal-mixture and ternary-reaction model assumptions. The framework is therefore not parameter-free when applied to real solutions.

free parameters (4)
  • rho_S (seawater) = not reported (constant fitted at T0=24.3 C, p0=20.21 dbar, y0=0.03956)
    Solute density in the binary ideal model; fit so that model density (61) matches seawater density at reference conditions. Used in eqs. (62)-(64) to compute beta and k.
  • rho_S (sucrose) = not reported (fitted at T0=20 C, p0=10.1325 dbar, y0=0.8)
    Solute density in the binary ideal model for sugar water; fitted to match measured density at reference. Appears in k via (64).
  • D(T) = determined by eq. (81) from V0(T) at reference composition and pressure
    Pressure-linear deviation term in the ternary ideal model for water-ethanol; fitted to the measured molar volume of the mixture at reference state. Appears in (80).
  • D'(T) = determined by eq. (82) from reference isothermal compressibility
    Quadratic pressure deviation coefficient in the ternary model; fitted to match the measured compressibility at reference state. Appears in (82).
assumptions (5)
  • domain assumption The entropy function h is of Legendre type, strictly convex, so the conjugate h* is well-defined and D^2h* is positive definite.
    Invoked in Section 2.2 to justify the variable transformation and in Section 3.1 to guarantee positivity of denominators such as D^2h* xi dot xi.
  • domain assumption Thermodynamic stability: {d^2_{T,p}g} < 0 and {d^2_{y_i,y_j}g}^t > 0, so the compositional Hessian is invertible.
    Needed in Section 3.2 to define (D^2_{y,y}g)^(-1) in eqs. (54)-(55) and to ensure k is positive.
  • domain assumption The mixture is an ideal mixture, with chemical potentials mu_i = g_i(T,p) + (RT/M_i) ln x_i.
    Used in Section 3.3 to derive the explicit formula (57) and throughout Section 4 for the binary and ternary ideal models.
  • ad hoc to paper Seawater can be treated as a binary ideal mixture with an incompressible solute (beta_S^T = 0).
    Stated in Section 4.1; enables formulas (61)-(64) for seawater, though the paper notes no pure-salt data are available.
  • ad hoc to paper Water-ethanol can be modeled as a ternary ideal mixture with a single chemical reaction A1 + gamma2 A2 -> A3 in chemical equilibrium, with a quadratic-in-pressure ansatz for the reaction Gibbs energy.
    Introduced in Appendix C; the deviation g tilde(T,p) = (p-p0)(D(T)+D'(T)(p-p0)/2) is a modeling choice, not derived.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Compressibility and volume variations due to composition in multicomponent fluids." pith.science (2026). https://pith.science/paper/ODVQUV42

@misc{pith2026241117742,
  author       = {Pith},
  title        = {Pith review of: Compressibility and volume variations due to composition in multicomponent fluids},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ODVQUV42}},
  note         = {Machine review of arXiv:2411.17742}
}
read the original abstract

For single components fluids, vanishing isothermal compressibility implies that the mass density is constant, but the same conclusion is unknown for multicomponent fluids. Here the volume remains affected by changes of the composition. In the present paper we discuss an apparently natural way to conceptualise, based on derivatives of the Gibbs function g, this 'volume change due to composition' as a compression. In this way the phenomenon becomes quantitatively comparable to the usual coefficients of compressibility. This is a first step to investigate the range of validity of the constant density approximation of multicomponent fluids. As an illustration, three different aqueous solutions are discussed.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 23 canonical work pages

  1. [1]

    R. D. Barbosa. High pressure and temperature dependence of thermodynamic properties of model food so- lutions obtained from in situ ultrasonic measurements. PhD thesis, University of Florida, 2003

    work page 2003

  2. [2]

    G. C. Benson and O. Kiyohara. Thermodynamics of aqueous mixtures of nonelectrolytes. i. excess volumes of water-n-alcohol mixtures at several temperatures. Journal of Solution Chemistry , 9:791–804, 1980. doi: 10.1007/BF00646798

    work page doi:10.1007/bf00646798 1980

  3. [3]

    Bothe and W

    D. Bothe and W. Dreyer. Continuum thermodynamics of chemically reacting fluid mixtures. Acta Mech. , 226(6):1757–1805, 2015. doi:10.1007/s00707-014-1275-1

    work page doi:10.1007/s00707-014-1275-1 2015

  4. [4]

    Bothe, W

    D. Bothe, W. Dreyer, and P .-E. Druet. Multicomponent incompressible fluids – an asymptotic study.ZAMM, Z. Angew. Math. Mech., 103(7):48, 2023. Id/No e202100174. doi:10.1002/zamm.202100174

    work page doi:10.1002/zamm.202100174 2023

  5. [5]

    Bothe and P .-E

    D. Bothe and P .-E. Druet. Well-posedness analysis of multicomponent incompressible flow models. J. Evol. Equ., 21(4):4039–4093, 2021. doi:10.1007/s00028-021-00712-3

    work page doi:10.1007/s00028-021-00712-3 2021

  6. [6]

    Bothe and P .-´E

    D. Bothe and P .-´E. Druet. On the structure of continuum thermodynamical diffusion fluxes – a novel closure scheme and its relation to the Maxwell-Stefan and the Fick-Onsager approach.Int. J. Eng. Sci., 184:33, 2023. Id/No 103818. doi:10.1016/j.ijengsci.2023.103818

    work page arXiv 2023

  7. [7]

    Dreyer, C

    W. Dreyer, C. Guhlke, and R. M¨uller. Overcoming the shortcomings of the nernst-planck model. Phys. Chem. Chem. Phys., 15:7075–7086, 2013. doi:10.1039/C3CP44390F

    work page doi:10.1039/c3cp44390f 2013

  8. [8]

    P .-E. Druet. Global-in-time existence for liquid mixtures subject to a generalised incompressibility constraint. J. Math. Anal. Appl., 499(2):56, 2021. Id/No 125059. doi:10.1016/j.jmaa.2021.125059

    work page arXiv 2021

Show all 24 references
  1. [9]

    P .-E. Druet. Maximal mixed parabolic-hyperbolic regularity for the full equations of multicomponent fluid dy- namics. Nonlinearity, 35(7):3812–3882, 2022. doi:10.1088/1361-6544/ac5679

    2022 doi

  2. [10]

    P .-E. Druet. Incompressible limit for a fluid mixture. Nonlinear Anal., Real World Appl. , 72:68, 2023. Id/No 103859. doi:10.1016/j.nonrwa.2023.103859. 26 P .-E. DRUET

    2023

  3. [11]

    N. P . Fofonoff and R. Millard Jr. Algorithms for the computation of fundamental properties of seawater . UN- ESCO Technical Papers in Marine Sciences; 44. Unesco, Paris, 1983. URL: http://hdl.handle.net/ 11329/10910.25607/OBP-1450

    1983

  4. [12]

    Giovangigli

    V. Giovangigli. Multicomponent flow modeling. Sci. China, Math. , 55(2):285–308, 2012. doi:10.1007/ s11425-011-4346-y

    2012

  5. [13]

    Giovangigli and L

    V. Giovangigli and L. Matuszewski. Supercritical fluid thermodynamics from equations of state. Physica D, 241(6):649–670, 2012. doi:10.1016/j.physd.2011.12.002

    2012 doi

  6. [14]

    Gouin, A

    H. Gouin, A. Muracchini, and T. Ruggeri. On the m¨uller paradox for thermal-incompressible media.Continuum Mechanics and Thermodynamics, 24:505–513, 2012. doi:10.1007/s00161-011-0201-1

    2012 doi

  7. [15]

    W. E. Haynes. CRC Handbook of Chemistry and Physics (95th ed.) . CRC Press, 2014. doi:10.1201/ b17118

    2014

  8. [16]

    The international thermodynamic equation of seawater – 2010: Calculation and use of thermodynamic properties

    IOC, SCOR, and IAPSO. The international thermodynamic equation of seawater – 2010: Calculation and use of thermodynamic properties . Intergovernmental Oceanographic Commission, Manuals and Guides No

    2010

  9. [17]

    B. S. Krumgalz and F . J. Millero. Physico-chemical study of dead sea waters: Ii. density measurements and equation of state of dead sea waters at 1 atm. Marine Chemistry , 11(5):477–492, 1982. doi:10.1016/ 0304-4203(82)90012-3

    1982

  10. [18]

    M ¨uller

    I. M ¨uller. Thermodynamics. Interaction of mechanics and mathematics series. Pitman, 1985. URL: https: //books.google.de/books?id=-Jd6AAAAIAAJ

    1985

  11. [19]

    M ¨uller

    I. M ¨uller. Grundz¨uge der Thermodynamik: mit historischen Anmerkungen . Springer-Verlag, 2013. doi:10. 1007/978-3-642-56474-1

    2013

  12. [20]

    C. R. Randall. Ultrasonic measurements of the compressibility of solutions and of solid particles in suspension. Bureau of Standards Journal of Research , 8(1):79, 1932. URL: https://nvlpubs.nist.gov/nistpubs/ jres/8/jresv8n1p79_a2b.pdf

    1932

  13. [21]

    Roux and J

    A. Roux and J. Desnoyers. Association models for alcohol-water mixtures. In Proceedings of the Indian academy of sciences-chemical sciences , volume 98, pages 435–451. Springer, 1987. doi:10.1007/ BF02861539

    1987

  14. [22]

    G. A. Schmidt, G. R. Bigg, and E. J. Rohling. Global seawater oxygen-18 database - v1.22. 1999. URL: https://data.giss.nasa.gov/o18data

    1999

  15. [23]

    Tanaka, T

    Y . Tanaka, T. Y amamoto, Y . Satomi, H. Kubota, and T. Makita. Specific volume and viscosity of ethanol-water mixtures under high pressure. The Review of Physical Chemistry of Japan, 47(1):12–24, 1977. URL: https: //repository.kulib.kyoto-u.ac.jp/dspace/bitstream/2433/47042/1...

    1977

  16. [56]

    URL: https://www

    UNESCO: United Nations Educational, Scientific and Cultural Organisation, 2010. URL: https://www. teos-10.org

    2010

Pith tools

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