Connected-tank study · damping surface · conduit sizing

Damping & Sizing

A connected-tank roll device is decided by two numbers: where its exchange mode sits, and how heavily that mode is damped. Both are now measured across a grid rather than inferred from a single case, and both feed a closed-form model that predicts the conduit area which maximises roll reduction. Geometry is the source-neutral dm1528 reference; no client data appears here.

DispositionDesign-grade for conduit sizing on this geometry family. The inertia law predicts the exchange period to ±3 % across a four-fold conduit-area range and 25–70 % fill; the loss law reproduces 35 measured cases to 8.6 % rms; the tank moment is closed form to 4.3 % rms for periods of 20 s and above.

The exchange period is 14.31 s, not the response peak

Earlier work on this capability read the exchange period off the amplitude peak of a forced-roll sweep, at about 23 s. That is not the natural period. Phase settles it without ambiguity: the lag through the exchange mode crosses 90° exactly at resonance whatever the damping, and it does so at 14.31 s. The amplitude peak is a different quantity — it is loss-controlled, and it moves.

The distinction is not academic. An effective conduit length calibrated to the response peak absorbs a damping-induced shift into an inertia parameter, and will not transfer to another conduit area or fill.

14.31 sexchange period, from phase
n = 2.14conduit-area exponent of the loss
8.6%loss law, rms over 35 cases
4–6 m²conduit area at maximum roll reduction

Equivalent damping over period and roll amplitude

Twenty-one forced-roll CFD cases, one per cell, identical geometry, mesh and solver settings throughout. Damping is not a property of the tank alone: it rises steeply with roll amplitude at every period, which is what a quadratic conduit loss does and a linear-viscous one cannot.

roll 10 deg, period 10 s: zeta_eq = 0.711roll 10 deg, period 13 s: zeta_eq = 0.682roll 10 deg, period 16 s: zeta_eq = 0.705roll 10 deg, period 20 s: zeta_eq = 0.669roll 10 deg, period 23 s: zeta_eq = 0.625roll 10 deg, period 26 s: zeta_eq = 0.568roll 10 deg, period 40 s: zeta_eq = 0.37410°roll 5 deg, period 10 s: zeta_eq = 0.446roll 5 deg, period 13 s: zeta_eq = 0.587roll 5 deg, period 16 s: zeta_eq = 0.527roll 5 deg, period 20 s: zeta_eq = 0.441roll 5 deg, period 23 s: zeta_eq = 0.389roll 5 deg, period 26 s: zeta_eq = 0.326roll 5 deg, period 40 s: zeta_eq = 0.188roll 2.5 deg, period 10 s: zeta_eq = 0.258roll 2.5 deg, period 13 s: zeta_eq = 0.437roll 2.5 deg, period 16 s: zeta_eq = 0.306roll 2.5 deg, period 20 s: zeta_eq = 0.269roll 2.5 deg, period 23 s: zeta_eq = 0.216roll 2.5 deg, period 26 s: zeta_eq = 0.176roll 2.5 deg, period 40 s: zeta_eq = 0.1002.5°10131620232640forcing period (s)0.20.40.6equivalent damping ζroll amplitude →
  • roll 2.5°
  • roll 5°
  • roll 10°
Equivalent damping ratio ζ of the exchange mode. Each ridge is one roll amplitude; each marker is one CFD case. Hover a marker for its value.
Equivalent damping ζ over the measured grid — one CFD case per cell, 21 cells
roll10 s13 s16 s20 s23 s26 s40 s
2.5°0.2580.4370.3060.2690.2160.1760.100
0.4460.5870.5270.4410.3890.3260.188
10°0.7110.6820.7050.6690.6250.5680.374

ζ is a linearisation of a quadratic loss and is valid at the amplitude it was measured at — it is not a material constant. The 13 s column is the least reliable: it sits closest to resonance, where the extraction is numerically ill-conditioned.

The natural period does not move; the response peak does

A single-amplitude sweep cannot separate these two. The grid can, and the separation is decisive: the exchange period is an inertial property and holds to under 2 % across a four-fold change in roll amplitude, while the amplitude peak doubles over the same range.

Natural period against amplitude peak, by roll amplitude
roll amplitudenatural period (90° phase)amplitude peakpeak amplification
2.5°14.15 s20 s×1.539
14.31 s23 s×1.269
10°14.41 s40 s×1.099

The 10° peak sits at the edge of the tested period range, so its amplification is a lower bound.

Predictions made before the runs, not fitted after

Two numbers are fitted, both at one conduit area and one fill: an effective conduit length, and a loss coefficient. Everything else is geometry. The cases below were then placed where the resulting predictions are falsifiable — at the predicted resonance of conduit areas and fills the model had never seen.

Exchange period: predicted before the run, measured after
casepredictedmeasurederrorphase lag at the predicted period
conduit area 3.4 m²19.73 s19.27 s+2.4 %88.3°
conduit area 13.5 m²10.64 s10.97 s-3.0 %95.0°
fill 25 %13.95 s13.45 s+3.7 %85.4°
fill 70 %14.59 s14.38 s+1.4 %88.3°

The loss law needed correcting by this test, and the correction is instructive. Derived as a form drag it scales as 1/Ac, which reproduced the calibration area but failed in opposite directions either side of it — under by 55 % at the smaller conduit, over by 108 % at the larger. Fitting the exponent rather than assuming it gives n = 2.14, and rms error over 35 cases at three conduit areas falls from 36.5 % to 8.6 %.

Tank roll moment, closed form against CFD (5° roll, validated band)
periodmeasuredclosed formerror
20 s9.103 MN·m8.570 MN·m-5.9 %
22 s9.498 MN·m9.040 MN·m-4.8 %
23 s9.606 MN·m9.181 MN·m-4.4 %
24 s9.673 MN·m9.279 MN·m-4.1 %
26 s9.683 MN·m9.342 MN·m-3.5 %
40 s9.206 MN·m9.075 MN·m-1.4 %
60 s8.807 MN·m8.763 MN·m-0.5 %

The moment adds two static terms the redistribution estimate omits — the weight moment of the whole fluid mass about the roll axis, and the within-leg free surface. Below about 20 s the closed form degrades, because the legs’ own sloshing modes begin to participate.

Roll reduction against conduit area

With every link validated, the chain runs from geometry to roll response, and the design question becomes answerable: for a given hull, which conduit area returns the most roll reduction?

superseded sizing 13–22 m²0%25%50%75%optimum: A_c 5.0 m2, 50.8 % reduction (roll 18 s, hull zeta 0.02)optimum: A_c 6.0 m2, 25.2 % reduction (roll 18 s, hull zeta 0.05)optimum: A_c 4.5 m2, 46.1 % reduction (roll 21 s, hull zeta 0.02)optimum: A_c 5.0 m2, 22.6 % reduction (roll 21 s, hull zeta 0.05)251015202530conduit area Ac (m²)peak roll reduction
  • roll period 18 s
  • roll period 21 s
  • solid = hull ζ 0.02 · dashed = hull ζ 0.05
Peak roll reduction against conduit area, tank authority α = 0.05. Filled markers are the optimum for each hull. Hover any optimum for its values.
Optimum conduit area by hull
roll periodhull ζoptimum Actank periodtuning ratiopeak reduction
18 s0.025.0 m²16.47 s0.91550.8 % extrapolated
18 s0.056.0 m²15.15 s0.84125.2 % extrapolated
21 s0.024.5 m²17.30 s0.82446.1 %
21 s0.055.0 m²16.47 s0.78422.6 %

Two results matter more than the peak value. The optimum is deliberately detuned — the tank period sits below the roll period, which is correct for a heavily damped absorber. And the optimum is broad: 5 and 6 m² are within one percent of each other, so the as-built routing uncertainty that worries a tuned device is survivable here.

Sensitivity to conduit area (roll 18 s, hull ζ = 0.02)
conduit areatank periodpeak reduction
3.0 m²20.94 s33.3 %
4.0 m²18.28 s43.1 %
5.0 m²16.47 s50.8 %
6.0 m²15.15 s50.0 %
8.0 m²13.31 s43.1 %
10.0 m²12.07 s33.9 %
13.5 m²10.64 s20.2 %
18.0 m²9.48 s9.8 %
22.0 m²8.79 s5.1 %
25.0 m²8.39 s3.0 %

Hull damping roughly halves the benefit at every conduit area. It is a property of the vessel, not of the tank, and it moves the answer by more than any sizing choice does.

What this does not establish

Every case and derived value on this page is published in the immutable release — digest 1034bb4efc5d9b39…, pinned dataset revision 51ba5ddbca8d…. That release grew from 24 to 57 cases and 133 to 441 derived metrics to carry this work, and all 47 declared source files resolve on the pinned revision.

Three of the 33 new cases are flagged, not clean. fill-f25-t12-a5, fill-f25-t13p95-a5 and grid-t10-a10 exceed the declared 2 % exchange cycle-change limit (2.25 %, 2.30 % and 3.24 %). They are published with status accepted_with_exception and counted in the release’s dispositions table rather than dropped, because they are real runs and because where they fall is informative: the two lowest-fill cases and the largest-amplitude shortest-period cell are the most nonlinear conditions in the set, and the slowest to settle. Statistics that assume every cell is equally converged should exclude them.

No time series were published. The samples table stands at 9,933 rows against a declared limit of 10,000, so adding series would breach the release’s own limit; the raw histories stay in the pinned private source.

Generated by scripts/generate_sloshing_damping_page.py, which regenerates this page from the reviewed analysis files — no value here is hand-transcribed, and --check fails the build if the page drifts from its sources. Sources: