Learning

Yacht hydrodynamics modeling using machine learning

Estimate yacht residuary resistance before committing to a hull design

This hydrodynamic surrogate predicts residuary resistance from five hull-form coefficients and the Froude number. It provides naval architects with a fast screening model for comparing candidate geometries and speeds during early-stage design.

22tested hull forms
308experimental cases
0.947recalculated testing R²
3.114recalculated testing RMSE

Residuary resistance influences the power required to reach a target speed and is therefore an important quantity during preliminary hull design. Physical towing-tank campaigns and high-fidelity simulations are valuable but costly to repeat for every candidate. A trained surrogate can provide rapid estimates inside an experimentally represented design space and help engineers decide which alternatives deserve detailed analysis.

1. Industrial challenge

This is an approximation problem: the model maps hull geometry and speed to the continuous response resistance, defined in the dataset as residuary resistance per unit weight of displacement.

Screen hull alternativesCompare candidate forms before allocating towing-tank or CFD resources.
Explore speed sensitivityQuantify how resistance changes with Froude number for a defined hull geometry.
Support early powering studiesUse rapid resistance estimates as one input to preliminary performance and propulsion decisions.

Potential users include naval architects, hydrodynamicists, yacht designers, shipyards, marine engineering consultancies and research teams.

Naval architectureHydrodynamicsYacht designShipyardsMarine R&D
Model role. This is a data-driven hydrodynamic surrogate for early design-space exploration. It complements, rather than replaces, towing-tank measurements, CFD, powering analysis and class or safety assessment.

2. Data set

The Delft Yacht Hydrodynamics dataset contains 308 experiments performed at the Delft Ship Hydromechanics Laboratory. It covers 22 hull forms derived from a parent form related to the Standfast 43, with 14 Froude numbers per hull form.

VariableEngineering meaningRoleRange
center_of_buoyancyLongitudinal centre of buoyancyInput-5.00 to 0.00
prismatic_coefficientPrismatic coefficientInput0.53 to 0.60
length_displacementLength–displacement ratioInput4.34 to 5.14
beam_draught_ratioBeam–draught ratioInput2.81 to 5.35
length_beam_ratioLength–beam ratioInput2.73 to 3.64
froude_numberDimensionless speed parameterInput0.125 to 0.450
resistanceResiduary resistance per unit weight of displacementTarget0.01 to 62.42

There are no missing values. The configured split uses 186 samples for training, 61 for selection and 61 for testing.

Distribution of yacht residuary resistance
Target distribution. Most experiments lie in the lower resistance range, with a long upper tail.
Correlations between yacht design inputs and residuary resistance
Input–target relationships. Froude number is expected to be the dominant driver, while hull coefficients modify the response.
Validation note. The example uses a random row split. Because every hull form is evaluated at several speeds, rows from the same hull may appear in both training and testing. For a stronger estimate of generalization to new designs, reserve complete hull forms for testing.

3. Model

The exported neural network receives six dimensionless inputs and returns one continuous resistance estimate. Mean-and-standard-deviation scaling feeds three tanh hidden neurons; a linear output is then unscaled to the resistance units used by the dataset. The resulting 6–3–1 model contains 25 trainable parameters and applies no output bounding.

Exported yacht resistance neural network with six inputs, three hidden neurons and one output

4. Training strategy

The model minimizes normalized squared error with L2 regularization (weight 0.01) and the quasi-Newton method. Training stopped after 62 epochs because the loss improvement fell below the configured threshold. The final recorded training and selection errors are 0.017 NSE and 0.007 NSE.

Updated quasi-Newton training and selection error history for the yacht resistance model

5. Model selection

The current exported deployment model uses three hidden neurons. This is the same architecture described above, evaluated in the testing section and implemented in the browser calculator and downloadable Python model.

For production use, architecture selection should be repeated with complete hull forms held out. A random row split can reward models that interpolate the speed curve of a hull geometry already represented during training.

6. Testing analysis

The exported Python model was independently recalculated on the 61 observations marked as testing in the supplied Neural Designer project. R² measures explained variation, while MAE and RMSE express error in the resistance units used by the dataset.

Testing samplesMAERMSEMaximum absolute error
610.94741.8813.11416.849

The model captures the dominant speed–resistance trend, but the largest errors occur at the high-resistance end. This matters operationally because those cases can drive powering and feasibility decisions.

Updated predicted versus measured yacht residuary resistance chart
Validation boundary. These metrics describe a random row holdout from one related hull family. They are not evidence of equal accuracy for an unseen hull geometry.

7. Model deployment

The surrogate can be integrated into an early-stage design workflow to screen geometries, generate speed–resistance curves and identify candidates for higher-fidelity analysis.

Hull coefficients and target speed
Range and feasibility checks
Hydrodynamic surrogate
Resistance estimate and design decision

Try the hydrodynamic surrogate model

Enter a hull geometry and Froude number within the experimental ranges. The browser evaluates the compact analytical expression published with this example.

-5.00 to 0.00
0.53 to 0.60
4.34 to 5.14
2.81 to 5.35
2.73 to 3.64
0.125 to 0.450
Predicted residuary resistance

Demonstration surrogate — not a certified hull-design or powering tool.

Response optimization

The updated run minimizes predicted resistance while constraining the longitudinal centre of buoyancy between -3 and -2. It returns the following feasible point inside the declared bounds:

QuantityOptimized value
Centre of buoyancy-2.341
Prismatic coefficient0.564
Length–displacement ratio4.943
Beam–draught ratio4.902
Length–beam ratio2.898
Froude number0.152
Predicted resistance0.008
Engineering interpretation. Because speed and nearly every hull coefficient are free, minimizing resistance naturally selects a low Froude number. A decision-grade optimization should hold the required mission speed, displacement, stability, capacity and geometric feasibility as constraints.

Directional response

The updated directional output varies Froude number around a reference hull defined by centre of buoyancy -2.3, prismatic coefficient 0.53, length–displacement ratio 4.76, beam–draught ratio 3.68 and length–beam ratio 3.16.

Updated predicted yacht resistance as a function of Froude number

Download and reproduce

The deployment package contains the current Python and browser models. The project package contains the Neural Designer project and its CSV dataset.

8. Scope and limitations

  • The dataset represents 22 related hull forms derived from one parent family, not the full yacht design space.
  • The model is valid only inside the coefficient and Froude-number ranges shown above.
  • The response covers residuary resistance per unit weight of displacement; total powering also requires other resistance components and propulsion-system assumptions.
  • The random row split does not prove generalization to unseen hull forms.
  • Feasible design constraints, stability, seakeeping, structural requirements and regulatory criteria are outside this surrogate.
  • Final decisions should be confirmed with naval-architecture analysis, CFD, model testing or full-scale evidence as appropriate.

References

  • Gerritsma, J., Onnink, R., & Versluis, A. (1981). Yacht Hydrodynamics dataset, UCI Machine Learning Repository.
  • Ortigosa, I., López, R., & García, J. (2007). A neural networks approach to residuary resistance of sailing yachts prediction. Proceedings of MARINE 2007.