WeirNet: A Large-Scale 3D CFD Benchmark for Geometric Surrogate Modeling of Piano Key Weirs
Authors: Lisa Lüddecke, Michael Hohmann, Sebastian Eilermann, Jan Tillmann-Mumm, Pezhman Pourabdollah, Mario Oertel, Oliver Niggemann
Organizations: Hydraulic Engineering Section, Civil Engineering, Helmut-Schmidt-University, Hamburg, Germany · Institute for Artificial Intelligence, Helmut-Schmidt-University, Hamburg, Germany · Professorship of Computer Science in Mechanical Engineering, Helmut-Schmidt-University, Hamburg, Germany · Federal Waterways Engineering and Research Institute, Karlsruhe, Germany · Chair for Automation Technology, Helmut-Schmidt-University, Hamburg, Germany
Reliable prediction of hydraulic performance is challenging for Piano Key Weir (PKW) design because discharge capacity depends on three-dimensional geometry and operating conditions. Surrogate models can accelerate hydraulic-structure design, but progress is limited by scarce large, well-documented datasets that jointly capture geometric variation, operating conditions, and functional performance. This study presents WeirNet, a large 3D CFD benchmark dataset for geometric surrogate modeling of PKWs. WeirNet contains 3,794 parametric, feasibility-constrained rectangular and trapezoidal PKW geometries, each scheduled at 19 discharge conditions using a consistent free-surface OpenFOAM workflow, resulting in 71,387 completed simulations that form the benchmark and with complete discharge coefficient labels. The dataset is released as multiple modalities compact parametric descriptors, watertight surface meshes and high-resolution point clouds together with standardized tasks and in-distribution and out-of-distribution splits. Representative surrogate families are benchmarked for discharge coefficient prediction. Tree-based regressors on parametric descriptors achieve the best overall accuracy, while point- and mesh-based models remain competitive and offer parameterization-agnostic inference. All surrogates evaluate in milliseconds per sample, providing orders-of-magnitude speedups over CFD runtimes. Out-of-distribution results identify geometry shift as the dominant failure mode compared to unseen discharge values, and data-efficiency experiments show diminishing returns beyond roughly 60% of the training data. By publicly releasing the dataset together with simulation setups and evaluation pipelines, WeirNet establishes a reproducible framework for data-driven hydraulic modeling and enables faster exploration of PKW designs during the early stages of hydraulic planning.
Figures & tables
Figure 1: Main geometric parameters of a rectangular PKW in plan view (left) and side view (right) defined by Pralong et al. [ 56 ] with Q indicating the direction of flow.
Parameters
Definition
B
Upstream-downstream length of the PKW (total weir length)
Bo
Upstream (outlet key) overhang crest length
Bi
Downstream (inlet key) overhang crest length
Bb
Base length
Pi
Height of the inlet entrance measured from the PKW crest (including possible parapet walls)
Po
Height of the outlet entrance measured from the PKW crest (including possible parapet walls)
Table 1: Nomenclature of the parameters of the PKW geometry by Pralong et al. [ 56 ] .
Figure 2: Exemplary results of a CFD simulation of water over a PKW geometry.
Study
Size
Ngeom
SpD
Dim.
Struct. type
Rep.
Outputs
Open
[sims/cases]
[–]
[–]
(P/M/PC)
Studies on PKW discharge coefficients
Bansal et al. [ 4 ]
162
9
18
3D
PKW key width ratio
✓/ – / –
cD , Energy dissipation (exp.)
–
Iqbal & Ghani. [ 30 ]
88
11
8
3D
PKW key width ratio & slope ratio
✓/ – / –
cD (exp.)
–
Studies on Weir geometry datasets
Eilermann et al. [ 14 ]
15
15
1
2D
Weir geometries
✓/ – / –
h(x) , free surface, cD
✓
Table 2: Comparison of existing studies and datasets on hydraulic and multiphase flow processes with geometric variation. Size counts CFD simulations or experimental cases, Ngeom denotes the number of distinct geometries where reported and SpD is the approximate number of simulated operating points (e.g., discharges or boundary conditions) per geometry. P/M/PC indicate availability of parametric descriptors (P), surface or volume meshes (M), and point clouds (PC). The Outputs column summarizes the primary simulated or measured quantities used for surrogate modeling.
Figure 3: Exemplary geometries from the dataset showing variations in key widths, overhang lengths and wall thickness.
Parameters
Definition
Wi,u
Inlet key width (from sidewall to sidewall) on the upstream side *1
Wi,d
Inlet key width (from sidewall to sidewall) on the downstream side *2
Wo,u
Outlet key width (from sidewall to sidewall) on the upstream side *2
Wo,d
Outlet key width (from sidewall to sidewall) on the downstream side *1
RB,i
Ratio of the downstream (inlet key) overhang crest length to base length
RB,o
Ratio of the upstream (outlet key) overhang crest length to base length
Table 3: Differentiation of the key widths for trapezoidal PKWs and extended nomenclature of the sidewall thickness for trapezoidal PKWs.
Figure 4: Detailed presentation of geometric relationships between inlet and outlet key widths for TPKWs with Q indicating the direction of flow.
Figure 5: Additional defined parameters at the PKW sidewall. Single parameters are shown on the left, while the corresponding reference locations are indicated on the right.
Parameters
Min
Max
Constant
W
1000 mm
P
330 mm
Bb
0.33⋅P
1.67⋅P
RB,i
0.25
1
Nu
3
Ts
0.015⋅P
0.18⋅P
Table 4: Limitations for the generation of PKW geometries via parametric design.
Parameter
Computational Relationship
Wu
W/Nu
Bi
RB,i⋅Bb
Bo
RB,o⋅Bb
B
Bb+Bi+Bo
α
tan−1(2⋅(B−Ts)Wi,u−Wi,d)
Ts,2
Ts/cos(α)
Table 5: Derived parameters in the parametric model
Q [l/s]
50
80
120
160
200
55
90
130
170
225
60
100
140
180
250
70
110
150
190
Table 6: Simulated discharges Q for all geometries.
Figure 6: Example of flow behavior over an exemplary PKW geometry in a scaled physical model test with flow direction from left to right.
Figure 7: Calculated discharge coefficients from CFD simulations of all investigated geometries for a representative exemplary selection of discharges simulated for the dataset.
Figure 8: Relative deviation between numerical and experimental results of upstream flow depths used for model validation.
Figure 9: Pearson correlation of selected weir parameters.
Figure 10: Pairwise correlation between the features and discharge coefficient.
Model
MSE (↓)
R2 (↑)
MAE (↓)
Max AE (↓)
Time/epoch (↓)
Inference Time (↓)
#Params
PointNet [ 57 ]
15.30
98.96
9.06
1.67
10 s
0.49 ms
5,435,009
Mesh-GCN [ 36 ]
241.60
64.33
36.88
2.99
24 s
0.90 ms
100,481
RegDGCNN [ 16 ]
97.40
83.11
22.01
1.90
1054 s
53.15 ms
4,482,657
Table 7: In-distribution surrogate modeling results for predicting cD at given Q using geometric representation data (Task 1). In addition, the Inference Time for batch size 1 as well as the total number of parameters (#Params) per model are reported. Reported errors are scaled: MSE×10−5 , R2×10−2 , MAE×10−3 , and MaxAE×10−1 .
Model
MSE (↓)
R2 (↑)
MAE (↓)
Max AE (↓)
RandomForest [ 9 ]
5.20
99.67
3.94
2.90
XGBoost [ 10 ]
6.80
99.58
5.53
1.44
LightGBM [ 34 ]
7.50
99.54
6.00
1.05
GradientBoosting [ 21 ]
119.50
92.63
19.61
3.54
Table 8: In-distribution surrogate modeling results with regression models for predicting cD at given Q using parametric data (Task 1). Reported errors are scaled: MSE ×10−5 , R2×10−2 , MAE ×10−3 , and Max MAE ×10−1 .
Figure 11: Parity plot of predicted versus CFD discharge coefficients cD on the in-distribution test set (Task 1).
Figure 12: SHAP summary plot for predictions based on all parametric models.
Figure 13: Distribution of the discharge coefficiency cD (left), discharge Q (middle), and inclination angle α (right).
Model
Split
MSE (↓)
R2 (↑)
MAE (↓)
Max MAE (↓)
RandomForest [ 9 ]
α≤2∘
1628.40
-19.43
94.34
5.27
XGBoost [ 10 ]
α≤2∘
1350.50
0.95
86.13
4.78
LightGBM [ 34 ]
α≤2∘
1242.40
8.88
83.35
4.80
GradientBoosting [ 21 ]
α≤2∘
1582.60
-16.07
94.28
5.19
RandomForest [ 9 ]
α3∘–5∘
87.20
87.59
19.11
4.93
XGBoost [ 10 ]
α3∘–5∘
37.10
94.72
14.53
1.57
Table 9: Out-of-distribution performance ( MSE , R2 , MAE , and MaxAE ) for all surrogate models on in-distribution (ID), OOD-Geom, and OOD-Head splits. Reported errors are scaled: MSE ×10−5 , R2×10−2 , MAE ×10−3 , and Max MAE ×10−1 .
Figure 14: Data-efficiency of surrogate models on WeirNet.
ETSIAE-UPM-School of Aeronautics, Universidad Politécnica de Madrid, Plaza Cardenal Cisneros 3, E-28040 Madrid, Spain · Institute for Cross-Disciplinary Physics and Complex Systems (IFISC, CSIC-UIB), 07122 Palma de Mallorca (Spain) · Center for Computational Simulation, Universidad Politécnica de Madrid, Campus de Montegancedo, Boadilla del Monte, 28660 Madrid, Spain