pith. machine review for the scientific record. sign in

arxiv: 2605.00168 · v1 · submitted 2026-04-30 · ⚛️ physics.flu-dyn

Recognition: unknown

Curvature-corrected sloshing spectra for cylindrical tanks in microgravity

Gianni Cassoni

Pith reviewed 2026-05-09 20:17 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords microgravity sloshingcylindrical tankscurved meniscuscapillary-gravity wavesBond numberYoung-Laplace equilibriasloshing spectraDirichlet-Neumann operator
0
0 comments X

The pith

Curved menisci must be modeled as leading order for accurate cylindrical microgravity sloshing frequencies.

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

In microgravity the free surface inside a partially filled cylinder settles into a curved meniscus set by the Young-Laplace equation. This shape changes both the liquid inertia and the capillary restoring force, so the natural sloshing frequencies differ from the flat-surface values even in the linear regime. The paper supplies a boundary-operator method that keeps the cylindrical Bessel structure while separating the bulk Dirichlet-Neumann contribution from the linearised curvature contribution. The separation makes the physical origin of the frequency shifts explicit and shows that the changes become order-one once the Bond number drops to unity or below. For spacecraft propellant management this means flat-interface formulas cease to be reliable in the capillarity-dominated range.

Core claim

A semi-analytical boundary-operator formulation is derived for capillary-gravity sloshing about axisymmetric Young-Laplace equilibria. The formulation treats the bulk Dirichlet-Neumann operator and the linearised curvature operator as distinct components, preserves the cylindrical Bessel structure, and recovers the flat-interface limit exactly. Equilibrium curvature is shown to couple radial modes and to alter the low-order spectrum for Bond numbers of order unity or smaller; concave menisci lower the fundamental frequency while convex menisci raise it, with the asymmetry carried predominantly by the Dirichlet-Neumann operator rather than the capillary term.

What carries the argument

Semi-analytical boundary-operator formulation that isolates the bulk Dirichlet-Neumann operator from the linearised curvature operator while retaining exact recovery of the flat-interface limit.

Load-bearing premise

Linearised curvature operator and separation of bulk and capillary contributions remain accurate without nonlinear or contact-line effects altering the eigenvalue spectrum in the Bond-number regime of interest.

What would settle it

Laboratory or microgravity measurement of the lowest sloshing frequencies in a cylindrical tank at Bond number near one, checked against both flat-interface and curved-meniscus predictions; systematic mismatch with the curved predictions would falsify the claim that curvature must be treated at leading order.

Figures

Figures reproduced from arXiv: 2605.00168 by Gianni Cassoni.

Figure 1
Figure 1. Figure 1: Geometry sketch for the cylindrical free-surface problem. Panel (a) shows the [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: First few sloshing frequencies ωn in the fixed azimuthal sector m = 1 as functions of the Bond number for representative concave and convex menisci at fixed fill depth h/R = 2. In each panel, the solid curves correspond to the flat reference θc = 90◦ , while the dashed and dotted curves correspond to the indicated contact angles. Each colour denotes a different radial branch. Deviations from the flat-inter… view at source ↗
Figure 3
Figure 3. Figure 3: Signed percentage change ∆1 in the fundamental frequency for the fixed azimuthal sector m = 1 over the (θc,Bo) plane, measured relative to the flat-interface value at the same Bond number. Contours at |∆1| = 1%, 5% and 10% mark the progressive breakdown of the flat-interface approximation. 16 [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Operator-split analysis at Bo = 0.3, m = 1, h/R = 2. Left: fundamental frequency for the four operator combinations. Right: percentage shift relative to the flat–flat reference. The fully curved result is not reproduced by either hybrid problem alone, confirming that curvature modifies both G and L. This operator-level view also explains the complementary-angle asymmetry noted in § 4.4. At fixed volume the… view at source ↗
Figure 5
Figure 5. Figure 5: Regime maps of the signed percentage frequency shift ∆ [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Mode geometry and modal composition at fixed contact angle [PITH_FULL_IMAGE:figures/full_fig_p019_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Mode geometry and modal composition for the first mode in the fixed azimuthal [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Mode geometry and modal composition for the second mode in the fixed azimuthal [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
read the original abstract

In microgravity, a partially filled cylindrical tank is generally bounded by a curved equilibrium meniscus rather than by an almost flat free surface. This modifies both the bulk liquid inertia and the capillary restoring force, so flat-interface sloshing frequencies can become inaccurate even in the linear regime. This effect matters once the Bond number is of order unity or smaller, precisely the regime relevant to capillarity-dominated propellant management. This study revisits the classical cylindrical curved-meniscus eigenvalue problem for capillary-gravity sloshing about axisymmetric Young-Laplace equilibria. A semi-analytical boundary-operator formulation is derived that preserves the cylindrical Bessel structure and recovers the flat-interface limit exactly. Its main advantage lies in treating the bulk Dirichlet-Neumann operator and the linearised curvature operator as distinct components, thereby making the physical origin of curvature-induced frequency shifts explicit. The results show that equilibrium curvature couples radial modes and alters the low-order spectrum once $Bo \lesssim 1$. Concave menisci lower the fundamental frequency, whereas convex menisci raise it while often lowering higher branches. The asymmetry between wetting and non-wetting configurations is found to be predominantly kinetic, being carried mainly by the Dirichlet-Neumann operator rather than by the capillary term. Curved menisci should therefore be treated as part of the leading-order model of cylindrical microgravity sloshing, not as a secondary correction, if reduced-order predictions are to capture the relevant dynamical scales for spacecraft applications.

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 paper claims that in microgravity, curved equilibrium menisci in partially filled cylindrical tanks modify linear capillary-gravity sloshing frequencies once the Bond number is O(1) or smaller. It derives a semi-analytical boundary-operator formulation that separates the bulk Dirichlet-Neumann operator from the linearised curvature operator while preserving the cylindrical Bessel structure and recovering the flat-interface limit exactly. Numerical results show that concave menisci lower the fundamental frequency whereas convex menisci raise it (often lowering higher branches), with the wetting/non-wetting asymmetry being predominantly kinetic and carried by the Dirichlet-Neumann operator rather than the capillary term. The central conclusion is that curved menisci must be retained at leading order in reduced-order models for spacecraft applications.

Significance. If the linearised formulation remains quantitatively faithful, the work supplies a practical, structure-preserving correction to flat-interface sloshing spectra that is directly relevant to propellant management in low-gravity environments. The explicit operator separation, exact flat-interface recovery, and absence of ad-hoc parameters are clear strengths that make the physical origin of the frequency shifts transparent.

major comments (2)
  1. [§2] §2 (linearised curvature operator and operator separation): the central claim that curvature must be treated at leading order rather than as a perturbation rests on the linearised Young-Laplace operator and the split between Dirichlet-Neumann and capillary contributions remaining accurate at Bo ≲ 1. At these values the equilibrium meniscus already possesses O(1) curvature; the manuscript should therefore supply either an a-priori error estimate or a direct comparison against a nonlinear reference solution showing that nonlinear capillary restoring forces and contact-line motion do not redistribute the reported frequency shifts by more than a few percent.
  2. [§4] §4 (asymmetry results): the assertion that the wetting/non-wetting asymmetry is 'predominantly kinetic' and carried mainly by the Dirichlet-Neumann operator is load-bearing for the physical interpretation. The manuscript should quantify this by tabulating, for each reported mode and Bo value, the separate contributions of the two operators to the eigenvalue shift; without such a decomposition the statement remains qualitative.
minor comments (2)
  1. [Abstract] The abstract states that the formulation 'preserves the cylindrical Bessel structure' but does not indicate the truncation order or the number of radial modes retained; this information should appear explicitly in §3 for reproducibility.
  2. [Figures] Figure captions should state the precise Bond-number values and contact angles used for each curve so that the trends can be compared directly with the flat-interface reference without consulting the main text.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below, indicating planned revisions where appropriate. Our responses focus on strengthening the manuscript while remaining within the scope of the linear semi-analytical framework.

read point-by-point responses
  1. Referee: §2 (linearised curvature operator and operator separation): the central claim that curvature must be treated at leading order rather than as a perturbation rests on the linearised Young-Laplace operator and the split between Dirichlet-Neumann and capillary contributions remaining accurate at Bo ≲ 1. At these values the equilibrium meniscus already possesses O(1) curvature; the manuscript should therefore supply either an a-priori error estimate or a direct comparison against a nonlinear reference solution showing that nonlinear capillary restoring forces and contact-line motion do not redistribute the reported frequency shifts by more than a few percent.

    Authors: We agree that the validity of the linearised formulation at O(1) Bond numbers merits explicit justification. The derivation in §2 is based on the standard linearisation of the capillary-gravity problem about the static Young-Laplace meniscus, which is the appropriate leading-order description for infinitesimal sloshing amplitudes. While a full nonlinear reference solution lies outside the present linear analysis, we will add a scaling estimate in the revised §2. This estimate quantifies the relative magnitude of nonlinear capillary and contact-line terms for amplitudes small compared with the meniscus height, showing that the reported frequency shifts are preserved to within a few percent in the linear regime. This supports retaining curvature at leading order for the small-amplitude dynamics targeted by the study. revision: partial

  2. Referee: §4 (asymmetry results): the assertion that the wetting/non-wetting asymmetry is 'predominantly kinetic' and carried mainly by the Dirichlet-Neumann operator is load-bearing for the physical interpretation. The manuscript should quantify this by tabulating, for each reported mode and Bo value, the separate contributions of the two operators to the eigenvalue shift; without such a decomposition the statement remains qualitative.

    Authors: We accept that a quantitative decomposition is necessary to make the kinetic origin of the asymmetry rigorous. Because the boundary-operator formulation explicitly separates the bulk Dirichlet-Neumann operator from the linearised curvature operator, the individual contributions to each eigenvalue shift are directly computable. In the revised manuscript we will insert a new table in §4 that reports, for the fundamental mode and representative higher modes at several Bond numbers, the separate eigenvalue shifts attributable to each operator. This decomposition will confirm that the wetting/non-wetting asymmetry is carried predominantly by the Dirichlet-Neumann term. revision: yes

standing simulated objections not resolved
  • Direct comparison against a fully nonlinear reference solution, which would require an independent nonlinear numerical framework outside the linear semi-analytical scope of the present work.

Circularity Check

0 steps flagged

No significant circularity; derivation uses standard operators and recovers known limits independently

full rationale

The paper derives a boundary-operator formulation from the linearised Young-Laplace equilibrium and classical Dirichlet-Neumann mapping for cylindrical geometry. It explicitly recovers the flat-interface limit by construction of the operators, which is a consistency check rather than a reduction of the result to its inputs. No fitted parameters are introduced as predictions, no self-citations are invoked as load-bearing uniqueness theorems, and the separation into bulk and capillary contributions follows directly from the linearised boundary conditions without circular redefinition. The central claim that curvature must be retained at leading order for Bo ≲ 1 follows from the computed eigenvalue shifts, which are obtained from the independent operator decomposition rather than presupposed. This is a standard, self-contained linear analysis with no evidence of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard potential-flow assumptions and linearization around Young-Laplace equilibria. No free parameters or invented entities are mentioned. The key methodological step is the operator separation, which is presented as preserving known cylindrical structure.

axioms (2)
  • standard math Incompressible, irrotational flow satisfying Laplace's equation inside the cylindrical domain
    Standard assumption for linear sloshing problems; invoked implicitly for the Dirichlet-Neumann operator
  • domain assumption Linearization of the dynamic and kinematic boundary conditions about the curved equilibrium meniscus
    Required to obtain the eigenvalue problem for small-amplitude oscillations

pith-pipeline@v0.9.0 · 5552 in / 1394 out tokens · 81897 ms · 2026-05-09T20:17:43.407863+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

38 extracted references · 36 canonical work pages

  1. [1]

    Abramson, H. N. 1966 The Dynamic Behavior of Liquids in Moving Containers . NASA SP-106, Washington, DC, https://ntrs.nasa.gov/citations/19670006555

  2. [2]

    & Baldi, P

    Alazard, T. & Baldi, P. 2015 Gravity capillary standing water waves. Arch. Ration. Mech. Anal. 217 (3), 741--830, http://dx.doi.org/10.1007/s00205-015-0842-5

  3. [3]

    Benjamin, T. B. & Scott, J. C. 1979 Gravity--capillary waves with edge constraints. J. Fluid Mech. 92 , 241--267, http://dx.doi.org/10.1017/S0022112079000616

  4. [4]

    & Franzoi, L

    Berti, M., Feola, R. & Franzoi, L. 2021 Quadratic life span of periodic gravity-capillary water waves. Water Waves 3 (1), 85--115, http://dx.doi.org/10.1007/s42286-020-00036-8

  5. [5]

    & Gallaire, F

    Bongarzone, A., Guido, M. & Gallaire, F. 2022 An amplitude equation modelling the double-crest swirling in orbital-shaken cylindrical containers. J. Fluid Mech. 943 , A28, http://dx.doi.org/10.1017/jfm.2022.440

  6. [6]

    & Gallaire, F

    Bongarzone, A., Viola, F., Camarri, S. & Gallaire, F. 2022 Subharmonic parametric instability in nearly brimful circular cylinders: a weakly nonlinear analysis. J. Fluid Mech. 947 , A24, http://dx.doi.org/10.1017/jfm.2022.600

  7. [7]

    Brakke, K. A. 1978 The Motion of a Surface by Its Mean Curvature . Mathematical Notes, vol. 20. Princeton University Press, Princeton, NJ, http://dx.doi.org/10.1515/9781400867431

  8. [8]

    Brakke, K. A. 1992 The Surface Evolver. Experimental Mathematics 1 (2), 141--165, http://dx.doi.org/10.1080/10586458.1992.10504253

  9. [9]

    & Nobili, M

    Cocciaro, B., Faetti, S. & Nobili, M. 1991 Capillarity effects on surface gravity waves in a cylindrical container: wetting boundary conditions. J. Fluid Mech. 231 , 325--343, http://dx.doi.org/10.1017/S0022112091003415

  10. [10]

    & Festa, C

    Cocciaro, B., Faetti, S. & Festa, C. 1993 Experimental investigation of capillarity effects on surface gravity waves: non-wetting boundary conditions. J. Fluid Mech. 246 , 43--66, http://dx.doi.org/10.1017/S0022112093000035

  11. [11]

    1968 Static menisci in a vertical right circular cylinder

    Concus, P. 1968 Static menisci in a vertical right circular cylinder. J. Fluid Mech. 34 (3), 481--495, http://dx.doi.org/10.1017/S002211206800203X

  12. [12]

    & Sulem, C

    Craig, W. & Sulem, C. 1993 Numerical simulation of gravity waves. J. Comput. Phys. 108 (1), 73--83, http://dx.doi.org/10.1006/jcph.1993.1164

  13. [13]

    & Nicholls, D

    Craig, W. & Nicholls, D. 2000 Travelling two and three dimensional capillary gravity water waves. SIAM J. Math. Anal. 32 (2), 323--359, http://dx.doi.org/10.1137/S0036141099354181

  14. [14]

    Popescu, Rémi Roumiguié, Thomas Miquel, Barbara Busset, Henri Bavestrello, and Jean Mignot

    Dalmon, A., Lepilliez, M., Tanguy, S., Alis, R., Popescu, E. R., Roumigui\'e, R., Miquel, T., Busset, B., Bavestrello, H. & Mignot, J. 2019 Comparison between the FLUIDICS experiment and direct numerical simulations of fluid sloshing in spherical tanks under microgravity conditions. Microgravity Sci. Technol. 31 (1), 123--138, http://dx.doi.org/10.1007/s1...

  15. [15]

    Dodge, F. T. & Garza, L. R. 1967 Experimental and theoretical studies of liquid sloshing at simulated low gravity. J. Appl. Mech. 34 , 555--562, http://dx.doi.org/10.1115/1.3607743

  16. [16]

    Dodge, F. T. 2000 The New Dynamic Behavior of Liquids in Moving Containers . Southwest Research Institute, San Antonio, TX, https://www2.swri.org/www2/fluid-slosh/slosh-dynamics.htm

  17. [17]

    & Inoue, S

    Funakoshi, M. & Inoue, S. 1988 Surface waves due to resonant horizontal oscillation. J. Fluid Mech. 192 , 219--247, http://dx.doi.org/10.1017/S0022112088001843

  18. [18]

    2009 Capillary damping of inviscid surface waves in a circular cylinder

    Kidambi, R. 2009 Capillary damping of inviscid surface waves in a circular cylinder. J. Fluid Mech. 627 , 323--340, http://dx.doi.org/10.1017/S0022112009005898

  19. [19]

    Kopachevskii, N. D. 1973 Hydrodynamics in weak gravitational fields: small oscillations of an ideal liquid in a cylindrical vessel. Fluid Dyn. 8 (2), 177--185, http://dx.doi.org/10.1007/BF01017524

  20. [20]

    Kopachevskii, N. D. & Krein, S. G. 2001 Operator Approach to Linear Problems of Hydrodynamics, vol. 1: Self-adjoint Problems for an Ideal Fluid . Birkh\"auser, Basel, http://dx.doi.org/10.1007/978-3-0348-8342-9

  21. [21]

    Kopachevskii, N. D. & Krein, S. G. 2003 Operator Approach to Linear Problems of Hydrodynamics, vol. 2: Nonself-adjoint Problems for Viscous Fluids . Birkh\"auser, Basel, http://dx.doi.org/10.1007/978-3-0348-8063-3

  22. [22]

    & Motygin, O

    Kozlov, V., Kuznetsov, N. & Motygin, O. 2004 On the two-dimensional sloshing problem. Proc. R. Soc. Lond. A 460 (2049), 2587--2603, http://dx.doi.org/10.1098/rspa.2004.1303

  23. [23]

    & Vainberg, B

    Kuznetsov, N., Maz'ya, V. & Vainberg, B. 2002 Linear Water Waves: A Mathematical Approach . Cambridge University Press, Cambridge, http://dx.doi.org/10.1017/CBO9780511546778

  24. [24]

    Kuznetsov, N. G. & Motygin, O. V. 2021 Sloshing in vertical cylinders with circular walls: the effect of radial baffles. Phys. Fluids 33 (10), 102106, http://dx.doi.org/10.1063/5.0064549

  25. [25]

    Miles, J. W. 1976 Nonlinear surface waves in closed basins. J. Fluid Mech. 75 , 419--448, http://dx.doi.org/10.1017/S002211207600030X

  26. [26]

    Miles, J. W. 1984 Internally resonant surface waves in a circular cylinder. J. Fluid Mech. 149 , 1--14, http://dx.doi.org/10.1017/S0022112084002500

  27. [27]

    Miles, J. W. 1984 Resonantly forced surface waves in a circular cylinder. J. Fluid Mech. 149 , 15--31, http://dx.doi.org/10.1017/S0022112084002512

  28. [28]

    Miles, J. W. 1991 The capillary boundary layer for standing waves. J. Fluid Mech. 222 , 197--205, http://dx.doi.org/10.1017/S0022112091001052

  29. [29]

    Moiseev, N. N. & Chernous'ko, F. L. 1965 Problems on the oscillation of a fluid subject to surface tension forces. USSR Comput. Math. Math. Phys. 5 (6), 128--160, http://dx.doi.org/10.1016/0041-5553(65)90102-3

  30. [30]

    D., Babskii, V

    Myshkis, A. D., Babskii, V. G., Kopachevskii, N. D., Slobozhanin, L. A. & Tyuptsov, A. D. 1987 Low-Gravity Fluid Mechanics . Springer, Berlin, https://link.springer.com/book/9783642709661

  31. [31]

    C., Saad, M

    Reynolds, W. C., Saad, M. A. & Satterlee, H. M. 1964 Capillary hydrostatics and hydrodynamics at low g . Report LG-3, Department of Mechanical Engineering, Stanford University, https://ntrl.ntis.gov/NTRL/dashboard/searchResults/titleDetail/PB166383.xhtml

  32. [32]

    & Gouriet, J.-B

    Simonini, A., Dreyer, M., Urbano, A., Sanfedino, F., Himeno, T., Behruzi, P., Avila, M., Pinho, J., Peveroni, L. & Gouriet, J.-B. 2024 Cryogenic propellant management in space: open challenges and perspectives. npj Microgravity 10 , 34, http://dx.doi.org/10.1038/s41526-024-00377-5

  33. [33]

    Slobozhanin, L. A. & Tyuptsov, A. D. 1973 Determination of the equilibrium state of a capillary liquid in a vessel. Fluid Dyn. 8 (4), 513--516, http://dx.doi.org/10.1007/BF01013086

  34. [34]

    Slobozhanin, L. A. & Tyuptsov, A. D. 1974 Stable equilibrium of the surface of a capillary liquid in contact with the edge of a solid. Fluid Dyn. 9 (1), 1--3, http://dx.doi.org/10.1007/BF01016301

  35. [35]

    Temkin, L. A. & Temkina, V. S. 1972 Calculation of the characteristic oscillations of an ideal liquid in an axially symmetric container with surface forces taken into account. Fluid Dyn. 7 (5), 715--722, http://dx.doi.org/10.1007/BF01205748

  36. [36]

    & Gallaire, F

    Viola, F., Brun, P.-T. & Gallaire, F. 2018 Capillary hysteresis in sloshing dynamics: a weakly nonlinear analysis. J. Fluid Mech. 837 , 788--818, http://dx.doi.org/10.1017/jfm.2017.860

  37. [37]

    Zakharov, V. E. 1968 Stability of periodic waves of finite amplitude on the surface of a deep fluid. J. Appl. Mech. Tech. Phys. 9 (2), 190--194, http://dx.doi.org/10.1007/BF00913182

  38. [38]

    E., Dyachenko, A

    Zakharov, V. E., Dyachenko, A. I. & Vasilyev, O. A. 2002 New method for numerical simulation of a nonstationary potential flow of incompressible fluid with a free surface. Eur. J. Mech. B/Fluids 21 , 283--291, http://dx.doi.org/10.1016/S0997-7546(02)01189-5