REVIEW 3 major objections 4 minor 16 references
Data-driven pressure field prediction for ships in regular sea states
T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read A convolutional-autoencoder/regression surrogate trained on time-averaged CFD pressure fields predicts the full wave-pressure field around a free-floating container ship from three sea-state parameters, yielding a held-out mean absolute…
desk verdict Solid engineering study of a CAE/NN surrogate for ship wave pressures, but the routing claim is unsupported by the ship-fixed-only validation and the time-averaging penalty in inertial-frame drag. 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 mechanism is the two-block CAE/NN surrogate. The convolutional autoencoder performs nonlinear order reduction: it squeezes a 128x64x96 pressure field down to ten latent numbers and reconstructs the field from them, capturing spatial dependencies with relatively few parameters. The fully connected regression network (200-200-10) maps the three sea-state control parameters into that same latent space, closing the loop from weather data to pressure field. Training uses a simultaneous loss that combines reconstruction error and prediction error on the full field, with the latent-matching term set to zero; interior ship cells are filled by iterated boundary means (Int.-Mean) so the convolutional filters near the hull are not corrupted by zeros.
What would settle it
Run the surrogate in the inertial frame by combining its predicted ship-fixed pressure field with the CFD pitch-angle history (or a separately predicted pitch angle), and compare the global-coordinate mean absolute drag error against the 0.94e-5 benchmark for instantaneous integration; if the error stays at the time-averaged level of 3.42e-5 rather than dropping, the surrogate does not deliver routing-grade resistance in earth-fixed coordinates.
Extended reading notes
Core claim
The central claim is that time-averaged, ship-fixed pressure deviations from calm water lie on a low-dimensional manifold parameterized by wave height, wave length and encounter angle, and that a CAE/NN pair can learn that manifold from 285 CFD-derived fields (150 simulated, 135 mirrored). The encoder reduces each field to ten latent variables, the regression network predicts the latent code from the three sea-state inputs, and the decoder outputs the full pressure field, which is then projected onto the hull and integrated to give added drag. The paper reports a held-out mean absolute error of 1.77e-5 for the ship-fixed drag coefficient with the finer meta-grid and iterative interior averaging, with head-sea cases reproduced well and the beam-sea short-wave case e1-3m as a notable outlier attributed to sparse data coverage. The authors evaluate only in the ship-fixed coordinate system, because for the inertial frame their own transfer analysis shows time-averaging introduces a larger error (3.42e-5) than integrating instantaneous fields (0.94e-5).
Load-bearing premise
The load-bearing premise is that a time-averaged pressure field written in the ship's own moving frame contains enough information to predict the added resistance a routing system would act on; the paper only tests that premise in the ship-fixed frame, and its own data-transfer check shows that averaging raises the global-coordinate drag error from 0.94e-5 to 3.42e-5.
Editorial extensions
If this is right
- A routing system could query the surrogate for many wave scenarios in the time it takes one CFD run to advance a few time steps, making added-resistance-aware route optimization feasible in practice.
- Because the machine-learning stage adds only about twice the error of the pure interpolation/integration pipeline, the accuracy ceiling is currently set by data transfer to the structured meta-grid rather than by the autoencoder itself.
- Head-sea conditions (encounter angles near 180 degrees) are predicted with the lowest errors, so initial deployment is most credible for head-sea routing.
- The time-averaged surrogate is validated only for ship-fixed drag; a global-coordinate routing target would need a pitch-angle estimate or an instantaneous-field variant.
- A latent dimension of ten already sits at the accuracy/cost sweet spot compared with dimensions five and twenty, so the approach does not need a large latent space to capture the wave-response manifold.
Reading between the lines
- The strongest unstated consequence is negative: if routing requires resistance in the earth-fixed frame, the surrogate in its current form is incomplete, and the missing piece is quantified by the paper's own data (3.42e-5 for time-averaged versus 0.94e-5 for instantaneous integration).
- A natural test of the method's ceiling is whether adding training cases near beam seas or explicitly encoding the mirror symmetry between beta and 360 degrees minus beta shrinks the outlier error of case e1-3m; the paper's data-coverage argument predicts it would.
- The same CAE/NN pipeline should transfer to other hulls and to other field quantities such as air concentration or local wave elevation, provided a similarly smooth latent manifold exists; the paper's broader claim about complex multiphase flow fields is extrapolation rather than demonstration.
- For an operational routing system, coupling this pressure-field surrogate with a cheap pitch-angle regression could recover most of the global-frame accuracy gap, since pitch enters the inertial-coordinate force integration in a nearly linear way.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes an offline/online surrogate-modeling pipeline for predicting wave-induced added resistance of a free-floating KRISO container ship. A convolutional autoencoder reduces time-averaged, ship-fixed CFD pressure fields to a low-dimensional latent space, and a fully connected neural network maps the three sea-state parameters (wave height, encounter angle, wave length) to that latent space; the decoded pressure fields are then projected onto the hull and integrated to obtain added-resistance coefficients. The manuscript evaluates the pipeline for 150 CFD-simulated regular sea states (augmented to 285 cases by mirroring) on two structured meta-grids, reports data-transfer errors separately from machine-learning errors, and presents per-case drag and pressure-field errors on 18 held-out test cases. The stated headline result is a mean ship-fixed added-drag error of 1.77e-5 for the best configuration on meta-grid B with interior cells filled by iterative averaging (Section 5.2, Table 5).
Significance. If the claimed accuracy carried over to the routing-relevant quantity, the work would be a useful engineering contribution: it demonstrates a complete chain from CFD data through projection, dimension reduction, regression, and force integration, with per-case error tables, a careful separation of data-transfer versus machine-learning errors, openly available training data, and honest reporting of outlier cases such as the beam-sea case e1-3m. The data-transfer verification in Section 5.1 is a particular strength because it isolates the interpolation error before any ML is introduced. However, the central routing claim in the abstract and introduction requires added resistance in the inertial/global frame, and the manuscript deliberately evaluates the surrogate only in the ship-fixed frame; Section 5.1 shows that time-averaging degrades the global-frame drag error substantially. That gap, together with unresolved issues around test-set selection and mirrored-case leakage, means the paper currently supports a narrower claim than the one it advertises.
major comments (3)
- [§5.2, Table 5 vs. §5.1] The headline accuracy of E'_D = 1.77e-5 (Table 5) is computed in the ship-fixed coordinate system, while the routing motivation requires the inertial-frame added resistance, i.e., the x1-component in Eq. (1). Section 5.1 reports that transferring time-averaged pressure fields to meta-grid B and integrating with an averaged pitch angle yields a mean global-frame drag error of 3.42e-5, compared with 0.94e-5 when instantaneous fields are integrated; the paper explicitly states that only the ship-fixed system will be used to judge the surrogate. Because the ML surrogate is trained on time-averaged fields and no model of the instantaneous pitch angle is described, the route-relevant global-frame error is at least the time-averaging error of 3.42e-5 unless the ML contribution compensates, for which no evidence is provided. The authors should either evaluate the full surrogate in the inertial frame or clearly restrict the claims to ship-fixed added resistance.
- [§4.1, Dataset split and mirrored test cases] The test set in Table 5 contains nine mirrored cases marked 'm' (e.g., b2-3m, d1-2m, e1-3m), which are exact reflections of simulated cases at the complementary encounter angle. The paper does not state whether the non-mirrored counterpart of each mirrored test case was excluded from the training and validation sets. If any such counterpart appears in training, the corresponding test errors are not independent because the mirrored pressure field is a deterministic transformation of the original. The authors should describe the random split explicitly and either exclude counterpart cases or report the test errors separately for mirrored and non-mirrored cases.
- [§5.2, Table 5, configuration selection] The reported mean error of 1.77e-5 corresponds to the configuration (meta-grid B, Int.-Mean, simultaneous training) that performs best on the same 18-case test table used for evaluation. Selecting the best of four configurations on the test set introduces selection bias, so the reported accuracy is optimistic relative to what would be obtained on a truly unseen configuration choice. The authors should either use a separate hold-out set for configuration selection or present the selection rule and its uncertainty transparently.
minor comments (4)
- [§1, Introduction] The sentence 'Thanks to the the now widely available hardware and software' contains a duplicated article; it should read 'Thanks to the now widely available hardware and software'.
- [§5.2, Table 5 caption and text] The symbol E'_D is introduced in the Table 5 caption but the text in Section 5.2 defines E_D and E_p; the notation should be unified so the ship-fixed drag error has a single name.
- [Appendix A, Table 6] Some numeric entries contain stray characters, e.g., '7 .765' in the c1-2 row, '4 .160' in the f3-2 row, and '2 .467' in the j4-1 row; these should be corrected.
- [References] The citation in Section 4.3 contains a stray question mark ('(?Wu et al. [2021], Schwarz et al. [2025b])'), and the reference to Muzaferija et al. contains a typo ('Syposium' instead of 'Symposium').
Circularity Check
No significant circularity: the reported surrogate accuracy is evaluated on held-out sea states, not on fitted values or equations that reduce to their own inputs.
full rationale
The derivation chain is self-contained. The CAE/NN surrogate is trained on time-averaged pressure fields from 249 CFD cases and evaluated on 18 held-out parameter combinations (Section 4.1, Table 5); the reported mean drag error of 1.77e-5 is an out-of-sample error rather than a training residual or a fitted quantity renamed as a prediction. The pressure-to-force integration path is verified separately against CFD in Section 5.1, and the paper explicitly quantifies the global-frame penalty of time averaging (3.42e-5 versus 0.94e-5 relative to instantaneous integration). The later decision to evaluate the full surrogate only in ship-fixed coordinates is an acknowledged limitation of the claimed routing applicability, not a circular reduction. Self-citations to the authors' own CFD solver, motion model, Hilbert-window selection, and CAE/NN architecture choices provide methods and prior design context, but none of these citations defines the measured test metric or logically forces the reported outcome. No equation is defined in terms of its own target, and no fitted parameter is relabeled as a prediction. The malformed citation marker '(?Wu et al. [2021], Schwarz et al. [2025b])' in Section 4.3 is noted; it is a formatting/reference issue and does not carry any circularity load.
Assumptions & free parameters
free parameters (9)
- Latent space dimension =
10
- Loss weights nu_reconst, nu_predict, nu_latent =
1, 1, 0
- Mini-batch size =
2 for CAE and simultaneous training, 128 for regression
- Learning rate =
1e-3 (Adam)
- Meta-grid B resolution =
128x64x96 cells
- Inverse-distance interpolation neighbor count =
5
- Early stopping patience =
50 epochs
- LeakyReLU slope =
0.01
- Virtual spring stiffness per simulation
assumptions (8)
- domain assumption RANS k-omega MSST and VoF accurately represent wave-induced pressure fields
- domain assumption Added resistance is pressure-dominated and friction changes are negligible
- ad hoc to paper 4DoF motion with roll and yaw suppressed, plus a virtual spring, gives representative vessel response
- domain assumption Time-averaged pressure fields in ship-fixed coordinates are sufficient for surrogate modeling and routing
- ad hoc to paper Mirroring pressure fields across the x'2 axis exactly represents beta+180 scenarios
- domain assumption Linear wave theory generates the incident waves
- domain assumption Interpolation to the structured meta-grid preserves the pressure information needed for drag integration
- ad hoc to paper The random train/validation/test split does not place mirrored twins in different splits in a way that leaks information
Cite this review
Pith. "Pith review of Data-driven pressure field prediction for ships in regular sea states." pith.science (2026). https://pith.science/paper/FRYAGJLR
@misc{pith2026250506014,
author = {Pith},
title = {Pith review of: Data-driven pressure field prediction for ships in regular sea states},
year = {2026},
howpublished = {\url{https://pith.science/paper/FRYAGJLR}},
note = {Machine review of arXiv:2505.06014}
}
read the original abstract
Merchant shipping is responsible for more than 90% of the global trade and has a significant environmental impact, accounting for over 2% of global greenhouse gas emissions. Therefore, fuel-saving measures are becoming increasingly important in reducing the ecological footprint and increasing the fuel efficiency of maritime transport. Routing optimization systems, which require a rapid prediction of ambient-dependent fuel consumption, represent an essential pillar here, e.g. to reduce added resistances due to seaways and/or wind. The paper aims to predict the added resistance due to seaways. In contrast to conventional methods the goal is achieved by surrogate modeling of the entire pressure fields. To this end, an online/offline-procedure is applied to an exemplarily free-floating container vessel. The online approach to be trained consists of two building blocks, namely a convolutional autoencoder (CAE)-based order reduction step and a neural network-based (NN) regression step that links the reduced space of the autoencoder with three control parameters that describe the sea state (wave height/steepness, encounter angle and wave length). Training data is obtained from time-averaged values for simulating instantaneous ship motion and pressure fields. During the offline phase, the combined CAE/NN is trained to capture the time-averaged pressure fields for a variety of sea-state conditions. During ship operation (online phase), the surrogate model predicts the three-dimensional pressure fields in response to sea state conditions, projects the pressure fields onto the ship hull, and integrates the corresponding resistances to guide the route. The evaluation of the method shows promising results for the different building blocks and the concept could therefore represent an attractive approach for cost-effective surrogate modeling of complex multiphase flow fields.
Figures
Figures from the paper (9 more)
Reference graph
Works this paper leans on
-
[2]
doi:10.1063/5.0020526. Y. Hou, H. Li, H. Chen, W. Wei, J. Wang, and Y. Huang. A novel deep U-Net-LSTM framework for time- sequenced hydrodynamics prediction of the SUBOFF AFF-8. Engineering Applications of Computational Fluid Mechanics, 16(1):630–645,
-
[4]
doi:10.1007/s00158-022-03338-2. 14 Table 6: Overview of all local drag coefficients based on CFD simulations (∆ C′ D) as well as reconstructed values from time averaged (∆C′ D,a) and instantaneous flow data (∆C′ D,i) for two meta-grids A and B. The mean drag error E′ D :=||C′ D−C′ D,a,i||1 is given at the end of the table. Run ∆C′ D/10−5 ∆C′ D,a/10−5 ∆C′ ...
-
[5]
doi:https://doi.org/10.1016/j.ast.2022.107629
ISSN 1270-9638. doi:https://doi.org/10.1016/j.ast.2022.107629. URL https://www.sciencedirect.com/science/article/ 16 pii/S1270963822003030. S. Liu, A. Papanikolaou, and G. Zaraphonitis. Prediction of added resistance of ships in waves.Ocean Engineering, 38(4):641–650,
-
[14]
doi:https://doi.org/10.1016/j.compfluid.2018.07.021
ISSN 00457930. doi:https://doi.org/10.1016/j.compfluid.2018.07.021. S. V¨ olkner, J. Brunswig, and T. Rung. Analysis of Non-Conservative Interpolation Techniques in Overset Grid Finite-Volume Methods. Computers & Fluids , 148:39–55,
-
[15]
doi:10.1016/j.compfluid.2017.02.010. P. Wu, S. Gong, K. Pan, F. Qiu, W. Feng, and C. Pain. Reduced order model using convolutional auto-encoder with self-attention. Physics of Fluids , 33(7):077107,
-
[1994]
doi:10.2514/3.12149. B. Mi and W. Cheng. Intelligent aerodynamic modelling method for steady/unsteady flow fields of airfoils driven by flow field images based on modified u-net neural network. Engineering Applications of Computational Fluid Mechanics, 19(1):2440075,
-
[2008]
Accessed: 2023-09-14. A. Solera-Rico, C. Sanmiguel Vila, M. G´ omez-L´ opez, Y. Wang, A. Almashjary, S. Dawson, , and R. Vinuesa.β- variational autoencoders and transformers for reduced-order modelling of fluid flows. Nature Communications, 15, 02
work page 2023
-
[2009]
H. Schwarz, P. P. Lin, J.-P. M. Zemke, and T. Rung. Disentangled latent spaces for reduced order models using deterministic autoencoders. arXiv preprint arXiv:2502.14679 , 2025a. URL https://arxiv.org/abs/2502. 14679. H. Schwarz, M. ¨Uberr¨ uck, J.-P. M. Zemke, and T. Rung. Machine learning based prediction of ditching loads. AIAA Journal, 63(5):1835–1854...
Show all 16 references
-
[2017]
doi:10.1080/09377255.2017.1348654. X. Luo-Theilen and T. Rung. Numerical Analysis of the Installation Procedures of Offshore Structures. Ocean Engineering, 179:116–127,
2017
-
[2019]
doi:10.1016/j.oceaneng.2019.03.004. A. L. Maas, A. Y. Hannun, and A. Y. Ng. Rectifier nonlinearities improve neural network acoustic models. In Proceedings of the 30th International Conference on Machine Learning , volume 28,
2019 doi
-
[2020]
doi:10.1063/5.0012906
ISSN 1070-6631. doi:10.1063/5.0012906. F. Chollet et al. Keras. https://keras.io,
-
[2021]
doi:10.1063/5.0051155. 17
-
[2022]
URL https://doi.org/10.1080/ 19942060.2022.2030802
doi:10.1080/19942060.2022.2030802. URL https://doi.org/10.1080/ 19942060.2022.2030802. D. P. Kingma and J. Ba. Adam: A method for stochastic optimization. In Proceedings of the 3rd International Conference on Learning Representations (ICLR),
2022
-
[2023]
doi:10.1063/5.0156963. J. L. Lumley. The structure of inhomogeneous turbulent flows. Atmospheric turbulence and radio wave propagation, pages 166–178,
-
[2024]
doi:10.1038/s41467-024-45578-4. R. Swischuk, L. Mainini, B. Peherstorfer, and K. Willcox. Projection-based model reduction: Formu- lations for physics-based machine learning. Computers & Fluids , 179:704–717,
-
[2025]
URL https://doi.org/10.1080/ 19942060.2024.2440075
doi:10.1080/19942060.2024.2440075. URL https://doi.org/10.1080/ 19942060.2024.2440075. M. Milano and P. Koumoutsakos. Neural network modeling for near wall turbulent flow.Journal of Computational Physics, 182(1):1–26,
2024
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.