Abstract

Two WJ-3 anemometers placed at the same height on the top of an architectural engineering building in Wenzhou University are used to determine the wind speed of Typhoon Morakot during its landing in real time. This study aims to explore Typhoon Morakot’s wind field characteristics, including mean wind speed, probability density distribution of fluctuating wind speed, power spectral density, correlation analysis, and coherence, on the basis of data measured by the two anemometers. Results show that the probability density distribution of the fluctuating wind speed of the typhoon follows the Gaussian distribution, and the measured cross-power spectrum of fluctuating wind speed is in good agreement with the modified Karman spectrum. The autocorrelation decreases with the increase in time interval (τ). The longitudinal autocorrelation coefficient decays rapidly with the increase in τ, and the lateral autocorrelation coefficient decays at an unchanged rate. The exponential attenuation coefficients of the longitudinal and transverse fluctuating wind speeds increase with the increase in the mean wind speed, and their mean values are 10.86 and 15.33, respectively. The change trends of the coherence coefficients of the two wind speed components with the mean wind speed are the same. The measured coherence coefficients of the two wind speed components are in good agreement with the exponential function.

1. Introduction

Statistics show that the loss caused by wind disasters ranks first amongst losses arising from all kinds of natural disasters [1]. The collapse of houses under strong wind conditions is the primary cause of property damage and casualties [2]. Studying wind disaster mechanisms and wind field characteristics is important for improving the antiwind performance of structures and reducing property damage and casualties [35]. Scholars at home and abroad have conducted extensive research on near-ground wind to master the mechanism of wind disasters. So far, several methods could be used to measure wind characteristics, which are mostly field test [68], numerical simulation [9, 10], and wind tunnel test [1113]. Nowadays, these methods are applied in many studies.

Since bridges become more and more complicated with the increasing long spans, the operability and controllability of numerical simulation are more and more prominent. Computational fluid dynamics (CFD) is an effective method used to analyse wind characteristics [14]; the simulation method is rapidly developing and widely used in wind characteristics [15, 16]. CFD simulation was carried out by Baghaei Daemei et al. to reduce drag coefficient performance of aerodynamic modification approaches. The results showed that modification of rounded corners was able to cause a reduction in the drag coefficient [17]. The accuracy of numerical simulation results depended on whether the numerical inflow boundary conditions match the mean and turbulent inflow conditions in the CFD method [18]. Li et al. carried out a numerical simulation for the wind characteristics of bridge site with complex terrain and provided a basis for the determination of design wind speed of the bridge [19]. Wind tunnel test and CFD simulation were conducted to study the wind characteristics at a bridge site in mountainous terrain, and the results of these two methods were compared [20]. Zhang et al. studied the influence of thermal effects on the wind characteristics at a bridge in a deep-cutting gorge with high-altitude and temperature difference in FLUENT. Results showed that the surface temperature distribution of the bridge was uneven under different temperatures [21]. Taking the surface conditions of a hill as an object, Tamura et al. used the large-eddy simulation (LES) method to study the turbulent boundary-layer flows over a hill and studied the influence of vegetation on turbulence statistics [14]. In addition, many scholars investigated wind characteristics and wind pressure via wind tunnel tests. A mountain terrain model with a scale of 1 : 2000 was built, and wind characteristics were carried out using wind tunnel tests [22]. By a combination of wind tunnel tests and numerical simulation, the profiles of wind speed and turbulence intensity varied with terrains were proved [23]. Li et al. built a terrain model with a scale of 1 : 1000 and studied wind characteristics with the method of wind tunnel tests, and the results showed that the incoming wind direction affected the wind characteristics [24]. Besides, to measure the unsteady aerodynamic force on a galloping prism and galloping response, HAPB (hybrid aeroelastic-pressure balance) wind tunnel tests were performed by Chen et al. [25]. The wind tunnel test of a tapered model showed that wind pressure distributions and unsteady force coefficients have different characteristics at low and high wind speeds and found a novel and unique phenomenon-termed partial reattachment [26, 27]. In these previously mentioned studies, most of these produced good results. However, as the most direct method, field measurements are essential for understanding the wind characteristics of a bridge site and are adopted by many scholars.

Comparing the numerical simulation methods, field measurements have been regarded as a reliable method to provide accurate evaluations of wind effects [28]. Huang et al. [29] conducted a comparative analysis of the near-ground fluctuating power spectrum of Pudong, Shanghai, under typhoon and monsoon action. The analysis results showed that the change trends of the attenuation coefficients of the coherence index of the fluctuating wind component under typhoon and monsoon action are consistent with the change in the mean wind speed. Cao et al. [30] analysed data synchronously measured by several anemometers at the same position. Results showed that when the wind speed is low, turbulence intensity decreases with the increase in mean wind speed. Turbulence intensity ceases to change with further increments in wind speed, and the measured power spectrum of fluctuating wind speed is in good agreement with the Karman spectrum only at a low frequency. Miyata et al. [31], using the measured data of Nos. 9807 and 9918 typhoons, analysed the power spectral density of longitudinal fluctuating wind speed and spatial coherence. Li et al. [32] explored the wind characteristic parameters at the site of Sutong Bridge by using measured data. They discovered that the ratio of longitudinal-to-transverse turbulence intensity is larger than the empirical value. Pang et al. [33] studied near-ground wind characteristics in Pudong, Shanghai, on the basis of a 3D ultrasonic wind speed sample measured at a height of 20 m and discovered that turbulence intensity and gust factor are high near the ground. Li et al. [34] explored the wind load characteristics of low buildings and the wind field characteristics of a strong typhoon before and after landing. The results showed that turbulence intensity and gust factor decrease significantly after typhoon landing, and the power spectral density of the fluctuating wind component is consistent with the Von Karman empirical spectrum. However, existing research results are against the wind field characteristics of low buildings. Most high-rise buildings in the coastal cities of China frequently suffer from typhoon disasters. Research on the field characteristics of high-rise buildings is rare, which is unfavourable for the wind resistance design of coastal structures.

In this study, the wind speed of Typhoon Morakot was obtained in real time for a duration of 13 h by using the topmost part of a construction engineering building in Wenzhou University as the observation site. The typhoon’s wind field characteristics, including mean wind speed, turbulence intensity, gust factor, peak factor, fluctuating wind speed power spectral density, and coherence, were analysed.

2. Typhon Morakot and Measurement Procedure

Typhoon No. 8 named “Morakot” was generated on August 4, 2009, in the Western Pacific Ocean and intensified into a strong typhoon the next day. The typhoon landed in Hualian, Taiwan, on August 7 with a maximum wind of Level 13 (40 m/s). On August 9 at 16 : 00, the typhoon landed in Fujian, Mainland China, with the maximum wind in the centre being Level 12 (33 m/s). Then, the typhoon was reduced into a tropical storm on the night of August 9. This typhoon was characterised by high intensity, long influence time, wide range, and low movement speed. The typhoon’s path is shown in Figure 1.

The measured site was at the top of a construction engineering building in Wenzhou University. The site was surrounded by hills in the east and south and open plains in the west and north. The building was about 8 m high and was in the suburb of the city. The roughness was classified as Class III in accordance with China’s regulations [35]. Two WJ-3 anemometers were set on the top of the construction engineering building at the same height. The anemometers were installed on 9 m high straight poles to avoid disturbance in the wind field. The transverse distance between the two straight poles was 17 m, and the distance from the anemometer to the ground was 30 m. The two anemometers collected wind speed synchronously at a sampling frequency of 20 Hz. The arrangement of the anemometers is shown in Figure 2 [36]. The specific parameters of the anemometers are shown in Table 1.

3. Data Processing Method

The real-time sampling period was from 22 : 00 on August 9, 2009, to 12 : 00 on August 11, 2009, and involved a total of 42 hours. The sample was divided by using 10 min as the basic time interval in accordance with regulations [35].

3.1. Mean Wind Speed and Direction

The measured sample includes two time series, namely, wind speed and wind direction. Wind speed has two components in the horizontal direction, namely, and . The expression is as follows:where u(t) and refer to wind speed series and wind direction series, respectively.

With 10 min as the basic time interval, mean wind speed U and horizontal wind angle ϕ were calculated via vector decomposition. The expression is as follows:where and refer to the mean wind speeds of ux(t) and uy(t), respectively, at 10 min.

3.2. Power Spectral Density of Fluctuating Wind Speed
3.2.1. Auto-Power Spectrum

The turbulence power spectrum represents the distribution of turbulence energy in the frequency domain and can describe the characteristics of fluctuating wind accurately. Its general expression can be written as follows [37]:where refers to Monin coordinate and A, B, , and refer to four parameters.

Previous field measurements and wind tunnel test results have indicated that the Von Karman spectrum can accurately reflect the characteristics of fluctuating wind. The functional expression of the power spectral density of longitudinal fluctuating wind is as follows [38]:where refers to the mean wind speed; refers to the integral scale of longitudinal turbulence; refers to the variance of longitudinal fluctuating wind speed; n is the frequency of longitudinal fluctuating wind; and c is 4.2065.

The functional expression of the power spectral density of horizontal and vertical wind can be obtained based on the hypothesis of isotropic turbulence as follows:

For isotropic turbulence, . The final expression of the Karman spectrum is obtained as follows:where is the variance of fluctuating wind speed component along i direction.

3.2.2. Cross-Power Spectrum

Research on the cross-power spectrum of the fluctuating wind speed of near-ground wind is rare. Kaimal et al. [39] derived the expression of the cross-power spectrum of longitudinal and vertical fluctuating wind speed by using well-conditioned experimental wind data:where refers to the cross-power spectrum, refers to friction speed, and f refers to Moin coordinates. Equation (7) is modified to obtain the expression of the cross-power spectrum of longitudinal and transverse fluctuating wind speed.where is equal to . Parameters a, b, and c are obtained by fitting and processing the sample at a 10 min time interval based on equation (8).

3.3. Autocorrelation Analysis

The dependence between the instantaneous value of a signal at one time and the instantaneous value at another time is called autocorrelation; it reflects the correlation between the values of the same sequence at different times. If X(t) refers to a time series, then the expression of the correlation function is as follows:

If X(t) refers to a stationary stochastic process, thenwhere refers to the autocorrelation function and refers to the time difference.

3.4. Coherence Analysis

The correlation between signals in the frequency domain is expressed by coherence. Davenport [40] proposed an empirical expression of the coherence coefficient of the fluctuating wind speed component in the form of an exponential function as follows:where f refers to the frequency; U refers to the mean wind speed; Cy and Cz refer to the index attenuation coefficients in horizontal and vertical directions, respectively; Δy and Δz refer to the horizontal and vertical distances between signal measuring points, respectively; and Coh refers to the coherence coefficient.

In this study, if Δy = 0, then equation (11) can be simplified as follows:

4. Results and Discussion

4.1. Mean Wind Speed and Angle of Wind Direction

In this study, the due north is 0° of the angle of wind direction, and the clockwise direction is positive. Figures 3(a) and 3(b) show the changing curves of the mean wind speed and angle of wind direction of two measuring points on the east and west, respectively, over 10 min. The mean wind speed and angle of wind direction of the two measuring points varying with time were consistent. From 22 : 00 on August 9, 2009, the mean wind speed increased with time and peaked at 14 : 00 on August 10. The maximum wind speeds at the eastern and western measuring points were 11.75 and 10.99 m/s, respectively. The mean wind speed decreased with time. The angle of wind direction varied greatly at 14 : 00. Figure 3(a) shows that the wind speed at that moment varied significantly, indicating that the typhoon gradually approached the measuring point, and the angle of wind direction gradually increased. Comparison of the wind speed and wind direction at both measuring points revealed that the mean wind speed of the eastern measuring point was slightly larger than that of the western measuring point. The difference at the peak value was 0.76 m/s. The angle of wind direction in the eastern measuring point was larger than that in the western measuring point, and the maximum difference was 22.87°.

4.2. Probability Density Distribution

The probability density function of fluctuating wind speed usually follows Gaussian distribution. Figure 4 shows the probability density distribution diagrams of longitudinal and transverse fluctuating wind speeds at measuring points 1 and 2 to determine whether the fluctuating wind speed follows Gaussian distribution under typhoon action. The Gaussian function distribution curve was obtained via moment estimation. The results showed that the probability density distribution of fluctuating wind speed in longitudinal and transverse directions approximately followed the Gaussian distribution, indicating that the probability density distribution of the fluctuating wind speed of this typhoon followed the Gaussian distribution.

4.3. Auto-Power Spectrum

Figure 5 shows power spectrum curves of fluctuating wind speed at measuring points 1 and 2. The Von Karman empirical spectrum is also provided for comparison. The results show that the power spectrum curve of longitudinal fluctuating wind speed coincided with the Von Karman empirical spectrum better than the power spectrum curve of transverse fluctuating wind speed did. The power spectrum curves of longitudinal and transverse fluctuating wind speeds coincided with the Von Karman empirical spectrum at a low frequency. However, at a high frequency, the measured power spectrum was smaller than the Von Karman empirical spectrum. Overall, the fluctuating wind power spectrum of this typhoon was in good agreement with the Von Karman empirical spectrum, which shows that the Von Karman empirical spectrum can be used to describe the fluctuating wind power spectrum change rule of this typhoon.

4.4. Cross-Power Spectrum

Figure 6 shows the cross-power spectrum of fluctuating wind speed in longitudinal and transverse directions at measuring points 1 and 2. A fitting analysis was conducted on the measured data on the basis of equation (7). The fitting parameters (a, b, and c) of the Karman correction spectrum at measuring point 1 were 6.36, 34.75, and 2.21. The fitting parameters (a, b, and c) of the Karman correction spectrum at measuring point 2 were 4.47, 16.59, and 2.48. The fitting parameters of the Karman correction spectrum were unaffected by the position of the measuring point. The corresponding Karman correction spectrum is provided in this paper for comparison. The results show that the measured cross-power spectrum of fluctuating wind speed was in good agreement with the Karman correction spectrum; thus, the Karman correction spectrum is also in good agreement with the variation law of the cross-power spectrum of fluctuating wind speed of this typhoon.

4.5. Correlation Analysis

Figure 7 shows the autocorrelation curves of longitudinal and transverse fluctuating wind speeds at measuring points 1 and 2. The variation trends of the autocorrelation coefficients of longitudinal and transverse fluctuating wind speeds were the same. The autocorrelation coefficients of longitudinal fluctuating wind speed at the two measuring points differed obviously. The longitudinal autocorrelation coefficient at measuring point 2 was larger than that in the transverse direction. The autocorrelation coefficient decreased with the increase in τ. Notably, the longitudinal autocorrelation coefficient decayed rapidly with the increase in τ, whereas the transverse autocorrelation coefficient decayed at a constant rate.

4.6. Coherence

The coherence of fluctuating wind speed components was also analysed. Figure 8 shows the change in the exponential attenuation coefficient of longitudinal and transverse fluctuating wind speed components with the average wind speed. The exponential attenuation coefficients of longitudinal and transverse fluctuating wind speed components increased with the increase in the mean wind speed, and the distribution became increasingly discrete. The mean values of the exponential attenuation coefficients of longitudinal and transverse fluctuating wind speed components were 10.86 and 15.33, respectively, which are larger than the values in the research of Cao et al. [30] and Hui et al. [41]. Figure 9 shows the variations of the coherence coefficients of the longitudinal and transverse fluctuating wind speed components with the mean wind speed. The variations of the coherence coefficients of the two components with the mean wind speed were consistent. With the increase in the mean wind speed, the coherence coefficient decreased gradually and eventually stabilised at 0.1.

The coherence coefficients of the longitudinal and transverse fluctuating wind speed components were compared with the fitting curves of exponential functions. Figures 10 and 11 show the coherence coefficient distribution of the longitudinal and transverse fluctuating wind speed components in different wind speed periods and the corresponding exponential function fitting curves. The measured results of the coherence coefficients of the longitudinal and transverse fluctuating wind speed components were in good agreement with the exponential function. This finding indicates that equation (12) is suitable for describing the coherence coefficients of the fluctuating wind speed components of this typhoon at eastern and western measuring points.

5. Conclusions

The characteristics of the wind field of Typhoon Morakot, including mean wind speed, turbulence intensity, gust factor, power spectrum density, and coherence, were analysed by using data measured by two anemometers placed on the top of a construction engineering building in Wenzhou University (30 cm high from the ground). The following conclusions were obtained:(1)The probability density of fluctuating wind speed in longitudinal and transverse directions approximately followed the Gaussian distribution, indicating that the probability density distribution of the fluctuating wind speed of this typhoon exhibited a Gaussian distribution.(2)The power spectrum curve of longitudinal fluctuating wind speed coincided with the Von Karman empirical spectrum better than that of transverse fluctuating wind speed did. The power spectrum curves of longitudinal and transverse fluctuating wind speeds coincided with the Von Karman empirical spectrum well at a low frequency. However, at a high frequency, the measured power spectrum was smaller than the Von Karman empirical spectrum. In general, the fluctuating wind power spectrum of this typhoon was consistent with the Von Karman empirical spectrum. Thus, the Von Karman empirical spectrum can be used to describe the fluctuating wind power spectrum change rule of this typhoon.(3)A fitting analysis of the measured data was conducted based on equation (10). a, b, and c fitting parameters of the Karman correction spectrum in the east measuring point were 6.36, 34.75, and 2.21, respectively; a, b, and c fitting parameters of the Karman correction spectrum in the east measuring point were 6.21, 32.57, and 2.42, respectively. The fitting parameters of the Karman correction spectrum were unaffected by the position of the measuring point. The corresponding Karman correction spectrum was provided in this paper for comparison. The measured cross-power spectrum of fluctuating wind speed was in good agreement with the Von Karman correction spectrum. This result indicates that the Von Karman correction spectrum is in good agreement with the variation law of the cross-power spectrum of fluctuating wind speed of this typhoon.(4)The variation trends of the autocorrelation coefficients of longitudinal and transverse fluctuating wind speeds were the same. The autocorrelation coefficients of longitudinal fluctuating wind speed at the two measuring points differed obviously. The longitudinal autocorrelation coefficient at measuring point 2 was larger than that in the transverse direction. The autocorrelation coefficient decreased with the increase in τ. Notably, the longitudinal autocorrelation coefficient decayed rapidly with the increase in τ, whereas the transverse autocorrelation coefficient decayed at a constant rate.(5)The exponential attenuation coefficients of longitudinal and transverse fluctuating wind speed components increased with the increase in the mean wind speed, and the distribution became increasingly discrete. The mean values of the exponential attenuation coefficients of longitudinal and transverse fluctuating wind speed components were 10.86 and 15.33, respectively, which are larger than the findings of Cao and Hui. The variations of the coherence coefficients of the longitudinal and transverse fluctuating wind speed components with the mean wind speed were consistent. With the increase in the mean wind speed, the coherence coefficient decreased gradually and eventually stabilised at 0.1. The measured results of the coherence coefficients of the longitudinal and transverse fluctuating wind components were consistent with the exponential function.(6)The fluctuating wind field based on the stationary theory can be simulated by the fluctuating wind spectrum model, but the simulation of the nonstationary wind field based on the nonstationary evolution power spectrum is still difficult and needs further analyses and research studies.

Data Availability

Some or all data, models, or code generated or used during the study are available from the corresponding author by request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The paper was supported by Zhejiang Province Natural Science Foundation Project (LY19E080022), Natural Science Foundation of China (51508419 and 51678455), and Zhejiang Provincial Department of Education Project (The title of the Project: Wind Field Characteristic Monitoring and Wind Load Monitoring Research on Coastal Low-rise Buildings, and the number: Y202147409).