60x
006036
2025-07-22

H.1.5.1. Group 1: Force and Moment Coefficient

User Story

In this example, we are going to calculate the average force and moment coefficient for Tokyo Polytechnic University (TPU) as an experimental example (M0/S0), such as those applicable to the structural design process based on WTG-Merkblatt-M3.

This part of the results belongs to Group 1 according to Figure 2.2 in WTG-Merkblatt-M3:

  • G2: Absolute values with medium accuracy requirements. The area of application can include parameters or preliminary studies when later investigations with higher accuracy are planned (e.g., wind tunnel examination of class G3).
  • R2: Solitary, all relevant wind directions with sufficiently fine directional resolution.
  • Z1: Statistical mean values, provided they concern stationary flow processes where fluctuations (e.g., due to approaching flow turbulence) can be captured sufficiently by other measures.
  • S1: Static effects. It is sufficient to represent the structural model with the necessary mechanical detail, but without mass and damping properties.

Description

By investigating the aerodynamic behavior across different planes, this comparison aims to provide a more comprehensive understanding of the force distribution. Such insights are critical for evaluating the model's overall stability, aerodynamic performance, and structural response under flow conditions. The average velocity and turbulence intensity inflow profiles reproduced in the wind tunnel correspond to those of Terrain Category IV according to the standard of the Architectural Institute of Japan (AIJ). The exponent of 0.25 cited in some figures was classified as a typographical error upon inquiry. The assumed model is shown in Image 1:

The assumptions underlying the analysis and simulations are summarized in Table 1, providing a clear overview of the parameters and conditions considered throughout the study.:

Table 1: Input data of the 3D rectangular model

Parameter Symbol Value Unit
Reference Wind Velocity UH 11 m/s
Roof Height Href 0.4 m
Profile Exponent α 0.25 -
Terrain Category - IV -
Air Density – RWIND ρ 1.25 kg/m³
Turbulence Model – RWIND RANS K-Omega - -
Kinematic Viscosity – RWIND ν 1.5×10⁻⁵ m²/s
Scheme Order – RWIND Second - -
Residual Target Value – RWIND 10⁻⁴ - -
Residual Type – RWIND Pressure - -
Minimum Number of Iterations – RWIND 800 - -
Boundary Layer – RWIND NL 10 -
Type of Wall Function – RWIND Enhanced / Blended - -

Computational Mesh Study

Image 2 summarizes a mesh sensitivity study, showing that as the mesh density increases from 20% to 55%, the force coefficient (Cf) increases from 0.95 to 1.05 and then stabilizes. This indicates that the simulation reaches mesh independence at 55%, meaning that further refinement no longer changes the results significantly. The image highlights the importance of mesh quality in ensuring reliable wind simulation outcomes.

Also, the computational mesh study needs to be performed according to the following link:

Image 3 presents a diagram comparing RWIND results to experimental data, concerning the plotting of normalized wind velocity profiles according to height. The results show very good agreement.

WTG-Merkblatt Accuracy Requirement

The WTG-Merkblatt M3 provides two key methods for validating simulation results. The Hit Rate Method evaluates how many of the simulated values Pi correctly match the reference values Oi within a defined tolerance, using a binary classification approach (hit or miss). This approach assesses the reliability of the simulation by calculating a hit rate q, similar to confidence functions used in reliability theory. In contrast, the Normalized Mean Squared Error (e2) method offers a more detailed accuracy assessment by quantifying the average squared deviation between simulated and reference values, normalized to account for scale differences. Together, these methods provide both qualitative and quantitative measures for simulation validation.

Results and Discussion

Table 3 shows very good agreement for the validation metric between the RWIND results and the WTG reference data for normalized velocity values. All the deviations are in an acceptable range (below 10%), so the hit rate is obtained q = 100%, and the normalized mean square error e2 = 0.00001 is extremely low. These results confirm that RWIND accurately reproduces the reference wind profile and meets strict validation criteria.

Table 3: Validation Metric for Velocity Profile

WTG - u / uref RWIND - u / uref Deviation (%) n
0.394 0.390 1.045 1.00
0.478 0.473 1.093 1.00
0.566 0.564 0.291 1.00
0.629 0.629 0.081 1.00
0.675 0.674 0.183 1.00
0.713 0.712 0.072 1.00
0.750 0.753 0.422 1.00
0.784 0.783 0.098 1.00
0.822 0.819 0.362 1.00
0.869 0.865 0.463 1.00
0.897 0.901 0.496 1.00
0.939 0.940 0.113 1.00
1.010 1.005 0.493 1.00
1.065 1.070 0.540 1.00
1.123 1.124 0.044 1.00
1.161 1.162 0.121 1.00
1.195 1.195 0.051 1.00
1.237 1.240 0.263 1.00
1.266 1.270 0.291 1.00
1.299 1.302 0.223 1.00

Image 4 illustrates how surface zones are defined in RWIND to accurately calculate the force coefficients on each face of a building model exposed to wind. The 3D model, with a simple rectangular shape, has each face front, back, left, right, and top assigned a distinct color-coded zone. These zones enable RWIND to separately analyze and calculate the aerodynamic forces and pressure distribution independently for each surface. Instead of placing numerous point probes, we chose to define surface zones, as RWIND can directly determine the force coefficient for each defined zone.

In the center of Image 4, the "Edit Zone" window is open for Zone No. 2 – Front as an example. This detailed output includes key aerodynamic parameters. The projected area of this zone in the wind direction is 0.04 m². The calculated drag force on this surface is 1.7 N, and the corresponding force coefficient (Cx) is 0.56. In the following formula, the wind force coefficient (aerodynamic force coefficient) for a specific pressure zone is calculated as an example. This value can then be compared to the value displayed under the Info tab in the Edit Model Data dialog.

Image 5 shows the post-processing results of a wind simulation conducted for a high-rise building using the k-omega turbulence model. On the right side of the image, a comparison table displays the force coefficients obtained from the CFD simulation alongside benchmark values from the WTG guideline. The front surface has a force coefficient of 0.562, compared to the WTG reference of 0.540, which corresponds to a deviation of 4.07%. The rear surface, which is subjected to negative pressure due to wake effects, shows a slightly larger deviation of 6.09%. The side surfaces also exhibit minor differences, with deviations ranging between 3.68% and 5.38%. Notably, the total force coefficient in the wind direction representing the global drag effect shows only a 1.95% deviation, highlighting the accuracy and reliability of the CFD simulation in capturing the overall aerodynamic behavior of the structure.

The actual force values in newtons for each building face are displayed below this table. The front surface is subjected to a positive drag force of 1.7 N, while the rear, right, and left surfaces experience negative forces due to suction and side drag. The total net force in the wind direction amounts to 3.2 N, matching the global force output stated in the mesh summary. The simulation used a reference wind speed of 11 m/s and a face area of 0.04 m².

Table 4 shows that RWIND accurately replicates the reference WTG mean force coefficients on all surfaces, with deviations ranging from 1.95% to 6.09% and q=100% as hit rate. The low normalized mean square error e2=0.0015 confirms strong agreement between simulation and measurements, meeting the validation standards effectively.

Table 4: Validation Metric for Mean Force Coefficients Between WTG and RWIND

Force Coefficient Cf,mean – WTG Cf – RWIND – k-omega Deviation (%) n
Front surface 0.540 0.562 4.07 1.00
Rear surface –0.528 –0.496 6.09 1.00
Right side surface –0.829 –0.860 3.68 1.00
Left side surface –0.847 –0.893 5.38 1.00
Total force in wind direction 1.070 1.058 1.95 1.00
Total force in crosswind direction 0.018 0.020 10.00 1.00
Moment in wind direction 0.582 0.520 10.65 1.00
Moment in crosswind direction 0.010 0.0099 1.00 1.00