pith. sign in

arxiv: 2512.21528 · v2 · submitted 2025-12-25 · ⚛️ physics.flu-dyn · cond-mat.soft

Under pressure: poroelastic regulation of flow in espresso brewing

Pith reviewed 2026-05-16 20:01 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn cond-mat.soft
keywords espresso brewingporoelasticityporous media flowdissolutioncoffee extractionnonlinear flowfluid dynamics
0
0 comments X

The pith

Interplay of elasticity and porosity in the coffee puck governs long-time espresso flow rate and solubles concentration.

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

The paper establishes that during espresso brewing the coffee grounds form a porous elastic medium whose permeability and stiffness change as soluble material dissolves under pressure. This coupled poroelastic response produces a non-linear relation between applied pressure and flow rate that persists over the extraction time. Consequently the amount of solubles carried into the cup is set by the time history of that regulated flow. A minimal model is introduced and shown to reproduce the observed temporal evolution when tested on a standard espresso machine. Dissolution is identified as the dominant process driving the changes in flow.

Core claim

The central claim is that the interplay between elasticity and porosity governs the long-time flow rate during espresso extraction and, consequently, the concentration of solubles in the final beverage. Using controlled experiments on a café-grade machine the authors show that a minimal model capturing the non-linear pressure-flow relationship, driven by dissolution-induced changes in the puck, reproduces the time-dependent behaviour of the brewing process.

What carries the argument

A minimal poroelastic model of the coffee puck in which permeability and elastic moduli evolve through dissolution under pressure.

If this is right

  • The long-time flow rate is set by the evolving poroelastic state of the puck rather than by a fixed permeability value.
  • Solubles concentration in the beverage is determined by the time-integrated flow under this poroelastic regulation.
  • The model reproduces the full temporal evolution of flow observed during extraction.
  • Dissolution dynamics are the central driver of changes in flow rate over time.

Where Pith is reading between the lines

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

  • Varying pressure profiles during extraction could exploit the elastic response to reach targeted concentrations.
  • The same poroelastic regulation may operate in other pressure-driven flows through dissolving porous materials.
  • Direct measurements of puck deformation during brewing would provide an independent check on the elastic component.

Load-bearing premise

The coffee puck can be treated as a homogeneous poroelastic medium whose permeability and elastic moduli evolve primarily through dissolution in a way captured by a minimal model.

What would settle it

A measured flow-rate curve under constant pressure that deviates substantially from the minimal model's prediction, or direct imaging showing that local heterogeneity rather than bulk poroelastic evolution dominates the regulation.

Figures

Figures reproduced from arXiv: 2512.21528 by Franciszek Myck, {\L}ukasz Bia{\l}as, Maciej Lisicki, Maria Puciata-Mroczynska, Micha{\l} Dzikowski, Piotr Szymczak, Radost Waszkiewicz.

Figure 1
Figure 1. Figure 1: The espresso brewing process. First the weighed grounds are placed in the brewing chamber (called portafilter), then clumps are removed with a needle tool (panel (A), necessary for good reproducibility). Then the grounds are tamped level with a tamper (panel (B), manual tamper pictured for demonstration only) to ensure consistent preparation pressure, the espresso can then be brewed (panel (C)). (D) Simpli… view at source ↗
Figure 2
Figure 2. Figure 2: Brewing assembly pressure drop calibra￾tion. (A) Espresso machine with empty basket with second pressure sensor and control valve added below the portafilter. (B) Combined measurements from multiple series and fitted quadratic approximation. Two with a control valve sweep and two with pressure regulator sweep. Brewer pressure drop de￾pends only on the flow rate through it. Brewer pressure drop, p1 − p2, as… view at source ↗
Figure 4
Figure 4. Figure 4: Theoretical equilibrium flow rate as a func￾tion of pressure Full theoretical expression (15) for nor￾malised flow rate Qˆ as a function of normalised pressure Pˆ at several values of control parameter Φ together with limit Φ → 0. Within the range of possible values of Φ the depen￾dence is negligible. which we integrate along the coffee bed to get µu(h0 − z) + C = K(σ), K(σ) = Z σ 0 k(ϕ(σ ′ ))dσ ′ (5) wher… view at source ↗
Figure 5
Figure 5. Figure 5: Experimental data on espresso flow at controlled pressures. (A) Basket pressure as function of time. (B) Simultaneously measured mass in cup. (C) Derived mass flow rate as a function of time. At low pressures (dark blue) increase in pressure results in an increase in flow rate, this is contrasted with opposite trend for large pressures (light yellow). Each line is an average result of multiple experiments … view at source ↗
Figure 6
Figure 6. Figure 6: Experimental long-time flow rate as a func￾tion of pressure. Symbols indicate flow rate values 110 − 120s after brewer activation, averaged over multiple experi￾mental runs (error bars show the standard error of the mean), together with the fitted theoretical model of Eq. (16) that as￾sumes pressure-induced compactification of the coffee bed and predicts flow rate saturation at higher pressures, while repr… view at source ↗
Figure 7
Figure 7. Figure 7: Determination of dissolution kinetics (A) A single espresso split into 14 fractions of 5 seconds each. Flow rate during first 20 seconds is much smaller than in later intervals. (B) Total dissolved solids (TDS) derived from refractometric measurements together with fitted approximation of Eq. (19). (C) Total flow rate (blue) and solutes flow rate (orange) obtained by multiplying total flow rate by the smoo… view at source ↗
Figure 8
Figure 8. Figure 8: Dissolution induced changes in porosity. (A) Comparison of time dependent flow profile derived from the exponential dissolution model for the most typical (9 bar) brewing pressure. (B) Comparison of several time dependent trajectories at different pressures. 1 cm (A) Pre-brew (B) Post-brew [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: X-ray tomographic slice through an espresso puck. (A) Tamped puck before brewing. Glass microbeads (dark) are added for visualisation. (B) Coffee puck after brewing. Considerable swelling is seen as compared to the pre-brewing puck, as well as numerous horizontal delamina￾tions. The coffee puck as a whole was also lifted from the bottom filter mesh, leaving a layer of air in between. V. DISCUSSION A. Groun… view at source ↗
read the original abstract

The sensory richness of coffee is widely recognised and arises from the complex chemistry and immersion in cultural practices of coffee preparation. In contrast, the physical complexity of espresso has received less attention. The multiphase reactive flow through a dissolving, elastic porous medium remains challenging to describe. Using a controlled experimental setup based on a caf\'e-grade espresso machine, we demonstrate that the interplay between elasticity and porosity governs the long-time flow rate during espresso extraction and, consequently, the concentration of solubles in the final beverage. We introduce a minimal model that captures the resulting non-linear pressure-flow relationship and propose a methodology capable of reproducing the time-dependent behaviour of the espresso brewing process. Finally, we show that dissolution dynamics play a central role in determining the temporal evolution of flow during extraction.

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 reports experiments on a café-grade espresso machine showing that long-time flow rate during extraction is regulated by the interplay of elasticity and porosity in the coffee puck, with dissolution driving the evolution. A minimal poroelastic model is introduced that reproduces the observed non-linear pressure-flow relationship and the time-dependent flow behavior, attributing the dynamics primarily to dissolution effects on permeability and elastic properties.

Significance. If the central claim is substantiated, the work supplies a physically motivated minimal model linking poroelastic mechanics to espresso flow regulation and solute extraction, which could inform reproducible brewing protocols and beverage quality control. The emphasis on a simple, dissolution-coupled description is a positive feature, though its value hinges on demonstrating that the model is not merely descriptive of the fitted data.

major comments (2)
  1. [Abstract and model description] The minimal model parameters (effective elastic modulus and permeability evolution rate) are constrained exclusively by the same pressure-flow time series used to support the central claim. No independent measurements of k(t) or E (e.g., via permeability cells or uniaxial compression on extracted pucks) are reported, leaving open whether the non-linear relation arises from the posited elasticity-porosity coupling or from unmodeled effects such as particle rearrangement or fines migration.
  2. [Experimental setup and model] The assumption that the coffee puck behaves as a homogeneous poroelastic medium whose permeability and elastic moduli evolve primarily through dissolution is load-bearing for the claim that dissolution dynamics govern long-time flow. The manuscript provides no controls or auxiliary data to rule out competing mechanisms, and the model equations are not shown to be predictive rather than fitted post hoc to the observed curves.
minor comments (2)
  1. [Abstract] The abstract states that the model 'captures' the non-linear relationship but does not indicate whether the governing equations are solved analytically, numerically, or via fitting; adding a brief outline of the solution procedure would improve clarity.
  2. [Model description] Notation for the time-dependent permeability k(t) and elastic modulus E should be defined explicitly when first introduced, including any functional forms assumed for their dissolution-driven evolution.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important aspects of model validation and experimental controls that we address point by point below. We agree that further clarification on parameter determination and discussion of alternative mechanisms will improve the manuscript.

read point-by-point responses
  1. Referee: [Abstract and model description] The minimal model parameters (effective elastic modulus and permeability evolution rate) are constrained exclusively by the same pressure-flow time series used to support the central claim. No independent measurements of k(t) or E (e.g., via permeability cells or uniaxial compression on extracted pucks) are reported, leaving open whether the non-linear relation arises from the posited elasticity-porosity coupling or from unmodeled effects such as particle rearrangement or fines migration.

    Authors: We acknowledge that the effective elastic modulus and permeability evolution rate are obtained by fitting the minimal model to the measured pressure-flow time series. This is inherent to the construction of a minimal model intended to isolate the essential poroelastic-dissolution coupling with the smallest number of parameters. Independent measurements of instantaneous permeability or modulus during active extraction are technically difficult because the puck properties change continuously under flow and pressure. We will revise the manuscript to include an expanded section on the fitting procedure, sensitivity analysis, and explicit discussion of how competing mechanisms (particle rearrangement, fines migration) could produce similar phenomenology, thereby clarifying the scope of the claim. revision: partial

  2. Referee: [Experimental setup and model] The assumption that the coffee puck behaves as a homogeneous poroelastic medium whose permeability and elastic moduli evolve primarily through dissolution is load-bearing for the claim that dissolution dynamics govern long-time flow. The manuscript provides no controls or auxiliary data to rule out competing mechanisms, and the model equations are not shown to be predictive rather than fitted post hoc to the observed curves.

    Authors: The homogeneous poroelastic description is adopted as the simplest framework that couples elasticity, porosity, and dissolution to reproduce the observed non-linear pressure-flow relation and the temporal evolution of flow rate. While auxiliary controls that isolate particle rearrangement or fines migration are not reported, the model’s ability to capture both the steady-state non-linearity and the transient flow behavior with a single set of dissolution-driven parameters provides evidence that these effects are dominant on the timescales of interest. To address the post-hoc fitting concern, we will add cross-validation by withholding portions of the time series for forward prediction and will include a brief discussion of possible heterogeneities in the revised text. revision: partial

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper introduces a minimal poroelastic model to describe the interplay of elasticity and porosity in espresso flow, fitted to pressure-flow time series from a commercial machine. No load-bearing step reduces by construction to its inputs: the model equations are derived from Darcy flow coupled to poroelastic constitutive relations and a dissolution ODE, with parameters constrained by data but not shown to be tautological renamings or self-citations. The central claim that dissolution dynamics govern long-time flow rests on physical assumptions and experimental reproduction rather than self-definition or fitted-input predictions. The derivation chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is inferred from the stated minimal model; it likely rests on standard poroelastic constitutive relations plus an assumption that dissolution dominates permeability evolution.

free parameters (1)
  • effective elastic modulus and permeability evolution rate
    These quantities are expected to be adjusted to match measured flow curves in the minimal model.
axioms (1)
  • domain assumption coffee puck behaves as a homogeneous poroelastic medium whose properties change primarily via dissolution
    Invoked to justify the minimal model that links elasticity, porosity, and flow.

pith-pipeline@v0.9.0 · 5462 in / 1373 out tokens · 36061 ms · 2026-05-16T20:01:12.356818+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

48 extracted references · 48 canonical work pages

  1. [1]

    used time-resolved CT to track the infiltration front of the water into the dry coffee bed, and proposed a one- dimensional unsaturated porous medium flow model to explain the dynamics of initial wetting. Since the preparation of an espresso in a pressurised container requires a consolidated bed of finely ground coffee, an intuitive modelling choice to qu...

  2. [2]

    Farah, Coffee constituents, inCoffee(John Wiley & Sons, Ltd, 2012) Chap

    A. Farah, Coffee constituents, inCoffee(John Wiley & Sons, Ltd, 2012) Chap. 2, pp. 21–58

  3. [3]

    J.Morris,MakingItalianEspresso, MakingEspressoItal- ian, Food and History8, 155–183 (2010)

  4. [4]

    J. Morris, Why espresso? explaining changes in euro- pean coffee preferences from a production of culture per- spective,EuropeanReviewofHistory: Revueeuropéenne d’histoire20, 881–901 (2013)

  5. [5]

    A. N. Gloess, B. Schönbächler, B. Klopprogge, L. D‘Ambrosio, K. Chatelain, A. Bongartz, A. Strittmat- ter, M. Rast, and C. Yeretzian, Comparison of nine com- moncoffeeextractionmethods: instrumentalandsensory analysis, European Food Research and Technology236, 607–627 (2013)

  6. [6]

    Hoffmann,The World Atlas of Coffee(Octopus Pub- lishing Group, 2018)

    J. Hoffmann,The World Atlas of Coffee(Octopus Pub- lishing Group, 2018)

  7. [7]

    M. I. Cameron, D. Morisco, D. Hofstetter, E. Uman, J. Wilkinson, Z. C. Kennedy, S. A. Fontenot, W. T. Lee, C. H. Hendon, and J. M. Foster, Systematically improv- ing espresso: Insights from mathematical modeling and experiment, Matter2, 631 (2020)

  8. [8]

    W. T. Lee, A. Smith, and A. Arshad, Uneven ex- traction in coffee brewing, Physics of Fluids35, 10.1063/5.0138998 (2023)

  9. [9]

    Andueza, L

    S. Andueza, L. Maeztu, B. Dean, M. P. de Peña, J. Bello, and C. Cid, Influence of water pressure on the final qual- ity of arabica espresso coffee. application of multivariate analysis, Journal of Agricultural and Food Chemistry50, 7426–7431 (2002)

  10. [10]

    Navarini and D

    L. Navarini and D. Rivetti, Water quality for espresso coffee, Food Chemistry122, 424–428 (2010)

  11. [11]

    Andueza, L

    S. Andueza, L. Maeztu, L. Pascual, C. Ibáñez, M. P. de Peña, and C. Cid, Influence of extraction tempera- ture on the final quality of espresso coffee, Journal of the Science of Food and Agriculture83, 240–248 (2003)

  12. [12]

    C. A. Salamanca, N. Fiol, C. González, M. Saez, and I. Villaescusa, Extraction of espresso coffee by using gra- dient of temperature. effect on physicochemical and sen- sorial characteristics of espresso, Food Chemistry214, 622–630 (2017)

  13. [13]

    Moroney, W

    K. Moroney, W. Lee, F. Suijver, J. Marra,et al., Mod- elling of coffee extraction during brewing using multiscale methods: An experimentally validated model, Chemical Engineering Science137, 216 (2015)

  14. [14]

    Smrke, A

    S. Smrke, A. Eiermann, and C. Yeretzian, The role of fines in espresso extraction dynamics, Scientific Reports 14, 10.1038/s41598-024-55831-x (2024)

  15. [15]

    World Coffee Research, Sensory lexicon,https://worl dcoffeeresearch.org/resources/sensory-lexicon (2017)

  16. [16]

    Fasman, Defining the ever-changing espresso, 25 Mag- azine3(2018)

    D. Fasman, Defining the ever-changing espresso, 25 Mag- azine3(2018)

  17. [17]

    Spencer, E

    M. Spencer, E. Sage, M. Velez, and J. Guinard, Using single free sorting and multivariate exploratory methods to design a new coffee taster’s flavor wheel, Journal of Food Science81, 10.1111/1750-3841.13555 (2016). 12

  18. [18]

    Fasano, F

    A. Fasano, F. Talamucci, and M. Petracco, The espresso coffee problem, inComplex flows in industrial processes, edited by A. Fasano (Springer, 2000) pp. 241–280

  19. [19]

    A. F. V. Matias, D. F. Valente-Matias, N. R. Neng, J. M. F. Nogueira, J. S. Andrade, R. C. V. Coelho, and N. A. M. Araújo, Continuum model for extraction and retention in porous media, Physics of Fluids35, 10.1063/5.0176273 (2023)

  20. [20]

    Mahadevan, A

    A. Mahadevan, A. V. Orpe, A. Kudrolli, and L. Mahade- van, Flow-induced channelization in a porous medium, EPL (Europhysics Letters)98, 58003 (2012)

  21. [21]

    Giacomini, G

    J. Giacomini, G. Khamitova, P. Maponi, S. Vittori, and L. Fioretti, Water flow and transport in porous media for in-silico espresso coffee, International Journal of Mul- tiphase Flow126, 103252 (2020)

  22. [22]

    C. Mo, R. Johnston, L. Navarini, F. S. Liverani, and M. Ellero, Exploring the link between coffee matrix mi- crostructure and flow properties using combined x-ray microtomography and smoothed particle hydrodynamics simulations, Scientific Reports13, 16374 (2023)

  23. [23]

    C. Mo, R. Johnston, L. Navarini, and M. Ellero, Model- ing the effect of flow-induced mechanical erosion during coffee filtration, Physics of Fluids33, 10.1063/5.0059707 (2021)

  24. [24]

    A. F. Matias, R. C. Coelho, J. S. Andrade, and N. A. Araújo, Flow through time–evolving porous media: Swelling and erosion, Journal of Computational Science 53, 101360 (2021)

  25. [25]

    K. M. Moroney, K. O’Connell, P. Meikle-Janney, S. B. O’Brien, G. M. Walker, and W. T. Lee, Analysing ex- traction uniformity from porous coffee beds using math- ematicalmodellingandcomputationalfluiddynamicsap- proaches, PLoS One14, e0219906 (2019)

  26. [26]

    Foster, W

    J. Foster, W. Lee, K. Moroney, D. Prjamkov, M. Sala- mon, A.Smith, J.Petrassem-deSousa,andM.Vynnycky, Dynamics of liquid infiltration into an espresso bed us- ing time-resolved micro-computed tomography: Insights from experiment and modeling, Physics of Fluids37, 10.1063/5.0245167 (2025)

  27. [27]

    B. R. Corrochano, J. Melrose, A. Bentley, P. Fryer, and S. Bakalis, A new methodology to estimate the steady- state permeability of roast and ground coffee in packed beds, J. Food Eng.150, 106 (2015)

  28. [28]

    Gianino, Experimental analysis of the italian coffee pot “moka”, American Journal of Physics75, 43 (2007)

    C. Gianino, Experimental analysis of the italian coffee pot “moka”, American Journal of Physics75, 43 (2007)

  29. [29]

    W. D. King, The physics of a stove-top espresso machine, American Journal of Physics76, 558–565 (2008)

  30. [30]

    Navarini, E

    L. Navarini, E. Nobile, F. Pinto, A. Scheri, and F. Suggi- Liverani, Experimental investigation of steam pressure coffee extraction in a stove-top coffee maker, Applied Thermal Engineering29, 998–1004 (2009)

  31. [31]

    Hedrick, This changes everything

    L. Hedrick, This changes everything... (Wdt, Autocomb, Moonraker, NCD, etc),https://www.youtube.com/wa tch?v=OyUaCuoqYWc&t=8s(2024)

  32. [32]

    Hedrick, How I Brew Tasty Coffee: My Espresso Ob- session Observed,https://youtu.be/eUc3lvF8mHM?si= cPBzma9w1T2XEnJi(2024)

    L. Hedrick, How I Brew Tasty Coffee: My Espresso Ob- session Observed,https://youtu.be/eUc3lvF8mHM?si= cPBzma9w1T2XEnJi(2024)

  33. [33]

    Gagné,The physics of espresso(Scott Rao, 2025)

    J. Gagné,The physics of espresso(Scott Rao, 2025)

  34. [34]

    A. J. Ladd and P. Szymczak, Reactive flows in porous media: Challenges in theoretical and numerical methods, Annual review of chemical and biomolecular engineering 12, 543 (2021)

  35. [35]

    Illy,Espresso coffee: the science of quality(Elsevier, 2004)

    A. Illy,Espresso coffee: the science of quality(Elsevier, 2004)

  36. [36]

    D. R. Hewitt, J. S. Nijjer, M. G. Worster, and J. A. Neufeld, Flow-induced compaction of a deformable porous medium, Phys. Rev. E93, 023116 (2016)

  37. [37]

    Guerriero and S

    V. Guerriero and S. Mazzoli, Theory of effective stress in soil and rock and implications for fracturing processes: a review, Geosciences11, 119 (2021)

  38. [38]

    L. S. Bennethum, M. Murad, and J. Cushman, Modified darcy’s law, terzaghi’s effective stress principle and fick’s law for swelling clay soils, Computers and Geotechnics 20, 245 (1997)

  39. [39]

    P. C. Carman, Fluid flow through granular beds, Trans. Inst. Chem. Eng. London15, 150 (1937)

  40. [40]

    W. L. McCabe, J. C. Smith, and P. Harriott,Unit oper- ations of chemical engineering, Vol. 5 (McGraw-hill New York, 1967)

  41. [41]

    C. Mo, L. Navarini, F. S. Liverani, and M. Ellero, Modeling swelling effects during coffee extraction with smoothed particle hydrodynamics, Physics of Fluids34, 10.1063/5.0086897 (2022)

  42. [42]

    Tien and B

    C. Tien and B. V. Ramarao, Can filter cake porosity be estimated based on the kozeny–carman equation?, Pow- der Technology237, 233–240 (2013)

  43. [43]

    Z. R. Lindsey, J. R. Williams, J. S. Burgess, N. T. Moore, and P. M. Splichal, Caffeine content in filter coffee brews as a function of degree of roast and extraction yield, Sci- entific Reports14, 10.1038/s41598-024-80385-3 (2024)

  44. [44]

    H. G. Severin and B. Lindemann, Elasticity of coffee beans: A novel approach to understanding the roasting process, J. Food Eng.383, 112212 (2024)

  45. [45]

    Zanoni, E.Pagliarini, andC

    B. Zanoni, E.Pagliarini, andC. Peri, Modelling the aque- ous extraction of soluble substances from ground roasted coffee, Journal of the Science of Food and Agriculture58, 275–279 (1992)

  46. [46]

    Spiro, Modelling the aqueous extraction of soluble substances from ground roast coffee, Journal of the Sci- ence of Food and Agriculture61, 371–372 (1993)

    M. Spiro, Modelling the aqueous extraction of soluble substances from ground roast coffee, Journal of the Sci- ence of Food and Agriculture61, 371–372 (1993)

  47. [47]

    Waszkiewicz, F

    R. Waszkiewicz, F. Myck, Ł. Białas, M. Puciata- Mroczynska, M. Dzikowski, P. Szymczak, and M. Lisicki, Under pressure: Poroelastic regulation of flow in espresso brewing, GitHub repository (2025)

  48. [48]

    Waszkiewicz, F

    R. Waszkiewicz, F. Myck, Ł. Białas, M. Puciata- Mroczynska, M. Dzikowski, P. Szymczak, and M. Lisicki, Under pressure: Poroelastic regulation of flow in espresso brewing, Zenodo (2025)