Important Principles in Wind Simulation

Technical Article on the Topic Structural Analysis Using Dlubal Software

  • Knowledge Base

Technical Article

RWIND 2 is a program for generating wind loads based on CFD (Computational Fluid Dynamics). The wind flow numerical simulation is generated around any building, including irregular or unique geometry types, to determine the wind loads on surfaces and members. RWIND 2 can be integrated with RFEM/RSTAB for the structural analysis and design or as a stand-alone application.

RFEM and RWIND are used to generate a tensile membrane structure model so that wind simulation can be started, along with the implementation of important criteria. RWIND is a powerful tool for creating wind loads on general structures and complicated forms. CFD solver is an OpenFOAM® software package (version 17.10), which gives very good results and is a widely used tool for CFD simulations. The numerical solver is steady-state for incompressible, turbulent flow, using the SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm.

The wind loads are regulated by specific standards, such as EN 1991-1-4, ASCE/SEI 7-16, or NBC 2015. RFEM is well-equipped technical software for creating tensile membrane structures; it considers nonlinear form-finding analysis for prestressed double-curvature surfaces. Image 02 presents a numerical airflow flowchart and FEM modeling for performing a verification example regarding tensile membrane structures.

In order to determine partial differential equations numerically, all differential expressions (space and time derivatives) need to be discretized. There is a wide range of discretization methods with different approaches in terms of accuracy, stability, and convergence. Generally, the order of the discretization illustrates how accurate the numerical simulation is when compared to the solutions of the original non-discretized equations. The first-order numerical discretization basically generates better convergence than the second-order scheme. In the current study, second-order discretization is used. Also, when using the second-order numerical scheme, we recommend increasing the maximum number of iterations to achieve better convergence (Image 03).

Verification example:

To verify the process for wind simulation, a double curvature model is developed as shown in [1] [2], and the results are investigated. The scale of the slight curve model is selected as 1/25, which is the same as the experimental model in reference [1], illustrating a hypar roof 10 m by 10 m by 1.25 m. The slight curve is considered for verifying an example with an angle θ=45o. The pretension force for the real scale was applied 2.5 kN/m on the surface, and mechanical properties such as Young's modulus and Poisson's ratio are defined as Ex=1000 kN/m. Ey=800 kN/m. Gxy=100 kN/m, vxy=0.20. Image 04 illustrates the geometry of the double curvature model. The input information and wind velocity input for CFD simulation are shown in Image 05.

Wind Tunnel Dimension:

It is important to note that the wind tunnel dimension can produce errors if the tunnel size is smaller than the standard type. The following picture shows a standard dimension of a wind tunnel [3]. Also, the results are sensitive to mesh sizes, so the calculation should perform for at least three different mesh numbers, and when the results are close enough to the previous stage, grid independence is achieved (Image 06).

A side view of the mesh generation of the model is presented in Image 07; as can be seen, the algorithm of mesh refinement is employed at a close distance to the model surface.

Computational grid study:

The results in CFD simulation are sensitive to mesh size, so grid independence should be performed for at least three different numbers of mesh elements. Here are the results of the Cp-value on the roof center line; as can be seen, there are very slight differences that show the results of wind simulation become independent of the third grid (Image 08).

Enhanced Wall Function:

RWIND uses the Blended Wall Function (BWF), also known as Enhanced Wall Function (EWF), which shows much better performance than the Standard Wall Function (SWF). Thus, generally, you can be sure that you will receive accurate results for the wide range of y+ numbers. The BWF is not a symmetric function. We have a solid black line (as shown in image 09), which is the Direct Numerical Simulation (DNS) that we are trying to reproduce. What you can see straight away is that the EWF is definitely much closer to the DNS data than the SWF. Thus, in the buffer region between a y+ of 5 and 30, the EWF is definitely much more accurate compared to the SWF. This is why EWF is often recommended in CFD simulations [4].

The mean pressure diagram and contour for the two turbulence models are shown in Images 10 and 11. Based on the center line of the hypar roof, the diagram of Cpe distribution is plotted and compared with the experimental wind tunnel test. The Cp value can be driven by the following equation, where P is the wind pressure at the measured point, Pref the reference pressure (atmosphere pressure), ρ is the air density, and Uref the reference velocity equal to 15.3 m/s.

Cp value

Cp=P-Pref12ρUref2

As shown here, the K-omega turbulence model shows better performance in predicting the wind pressure coefficient; in the current study, the K-omega turbulence model captures the effects of vortex shedding in high gradient negative wind pressure better than the K-epsilon models. We recommend using this turbulence model as a more accurate option in wind-structure interactions.

Author

Mahyar Kazemian

Mahyar Kazemian

Product Engineering

Mr. Kazemian is responsible for providing technical articles, product engineering and marketing for Dlubal Software, specially RWIND 2.

Keywords

Wind simulation Wind load Computational fluid dynamics Tensile membrane structure Mahyar Kazemian

Reference

[1]   Colliers, J., et al., Prototyping of thin shell wind tunnel models to facilitate experimental wind load analysis on curved canopy structures. Journal of Wind Engineering and Industrial Aerodynamics, 2019. 188: p. 308-322.
[2]   Rizzo, F., et al., Wind action evaluation on tension roofs of hyperbolic paraboloid shape. Engineering structures, 2011. 33(2): p. 445-461.
[3]   Zhang, C., et al., Wind pressure coefficients for buildings with air curtains. Journal of Wind Engineering and Industrial Aerodynamics, 2020. 205: p. 104265.
[4]   What are Wall Functions and How do they work? https://www.youtube.com/watch?v=h5OiFpu0L4M

Downloads

Write Comment...

Write Comment...

  • Views 591x
  • Updated 18 January 2023

Contact Us

Contact Dlubal

Do you have any further questions or need advice? Contact us via phone, email, chat or forum or search the FAQ page, available 24/7.

+49 (0) 9673 9203-0

[email protected]

Event Invitation

International Mass Timber Conference

Conference 27 March 2023 - 29 March 2023

Form-Finding and Calculation of Membrane Structures in RFEM 6

Form-Finding and Calculation of Membrane Structures in RFEM 6

Webinar 25 August 2022 2:00 PM - 3:00 PM CEST

CFD Wind Simulation with RWIND 2

CFD Wind Simulation with RWIND 2 (USA)

Webinar 29 June 2022 2:00 PM - 3:00 PM EDT

Tensile Membrane Structure Design in RFEM 6

Tensile Membrane Structure Design in RFEM 6 (USA)

Webinar 17 March 2022 2:00 PM - 3:00 PM EDT

Wind Load Generation Based on CFD in RWIND 2

Wind Load Generation Based on CFD in RWIND 2

Webinar 4 March 2022 10:00 AM - 11:00 AM CET

Form-Finding and Calculation of Membrane Structures in RFEM 6

Form-Finding and Calculation of Membrane Structures in RFEM 6

Webinar 3 February 2022 2:00 PM - 3:00 PM CET

Online Training | English

Eurocode 5 | Timber Structures According to DIN EN 1995-1-1

Online Training 23 September 2021 8:30 AM - 12:30 PM CEST

Online Training | English

Eurocode 3 | Steel Structures According to DIN EN 1993-1-1

Online Training 25 August 2021 8:30 AM - 12:30 PM CEST

Online Training | English

RFEM for Students | USA

Online Training 11 August 2021 1:00 PM - 4:00 PM EDT

Online Training | English

RFEM | Structural Dynamics and Seismic Design According to EC 8

Online Training 11 August 2021 8:30 AM - 12:30 PM CEST

Online Training | English

Eurocode 2 | Concrete Structures According to DIN EN 1992-1-1

Online Training 29 July 2021 8:30 AM - 12:30 PM CEST

Online Training | English

RFEM | Basics

Online Training 13 July 2021 9:00 AM - 1:00 PM CEST

Online Training | English

RFEM | Basics | USA

Online Training 17 June 2021 9:00 AM - 1:00 PM EDT

Online Training | English

RFEM for Students | Part 3

Online Training 15 June 2021 2:00 PM - 4:30 PM CEST

Glass Design with Dlubal Software

Glass Design with Dlubal Software

Webinar 8 June 2021 2:00 PM - 2:45 PM CEST

Online Training | English

RFEM | Structural Dynamics and Seismic Design According to EC 8

Online Training 2 June 2021 8:30 AM - 12:30 PM CEST

Online Training | English

Eurocode 5 | Timber Structures According to DIN EN 1995-1-1

Online Training 20 May 2021 8:30 AM - 12:30 PM CEST

Online Training | English

RFEM for Students | Part 2

Online Training 17 May 2021 2:00 PM - 4:30 PM CEST

Blast Time History Analysis in RFEM

Blast Time History Analysis in RFEM

Webinar 13 May 2021 2:00 PM - 3:00 PM EDT

RWIND

Stand-Alone Program

RWIND is a program (digital wind tunnel) for the numerical simulation of wind flows around any building geometries with determination of the wind loads on their surfaces. It can be used as a stand-alone application or together with RFEM and RSTAB for complete structural analysis and design.

Price of First License
3,450.00 EUR
RFEM 5
RFEM

Main Program

Structural engineering software for finite element analysis (FEA) of planar and spatial structural systems consisting of plates, walls, shells, members (beams), solids, and contact elements

Price of First License
4,350.00 EUR