Extreme-Point Symmetric Mode Decomposition to Define the Turbulence Characteristics of a Flume Flow

Shi, Zhi-Peng; He, Ting-Ting; Zhang, Gen-Guang

doi:https://doi.org/10.1155/2021/5595918

Mathematical Problems in Engineering

On this page

Abstract Introduction Results and Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 5595918 | https://doi.org/10.1155/2021/5595918

Extreme-Point Symmetric Mode Decomposition to Define the Turbulence Characteristics of a Flume Flow

Zhi-Peng Shi,¹Ting-Ting He,¹and Gen-Guang Zhang²

Academic Editor: Adnan Maqsood

Received01 Feb 2021

Accepted29 Jun 2021

Published09 Jul 2021

Abstract

Turbulence is a key feature of solid-liquid two-phase flows, and the pulsating velocity is the basis for calculating turbulence characteristics. In general, the method of mathematical expectation is used to calculate pulsating velocity. However, this method does not reflect the fluctuating state of the instantaneous velocity. Therefore, the method of extreme-point symmetric mode decomposition (ESMD) is adopted to calculate pulsating velocity and turbulence characteristics. The ESMD involves two stages, namely, modal decomposition and time-frequency analysis. The optimal adaptive global mean (AGM), which is the result of modal decomposition, can accurately reflect the fluctuation state of the instantaneous velocity, and the theory of the pulsating velocity defined on this basis is reasonable. Moreover, the flow pattern and turbulence behaviour of a two-phase flow can be predicted using the calculated turbulence characteristics. The method is used to analyse the pulsating velocity of the flume, and its rationality in theoretically predicting the turbulence behaviour of flume flows is demonstrated.

1. Introduction

The suspension of solid particles in turbulent flows can be widely observed in nature, such as in the case of fine sediments in fluvial and estuarine environments [1–3]. Owing to the frequent occurrence of such flows, the complex phenomena pertaining to the interaction of solid particles and turbulence must be examined. Many experimental and computational studies related to solid particles in open-channel turbulent flows have been reported, for example, studies pertaining to sediment transportation and the measurement of the pulsating velocity of liquid-solid two-phase flows [4–7].

At present, the use of dynamic response signals to analyse turbulence is a key research topic pertaining to two-phase flows. In this context, the turbulence characteristics must be tracked based on the acquired nonstationary velocity response signals. In particular, due to the time-varying process of the flume flow, the properties for the turbulence characteristics must be tracked based on nonstationary response signals acquired through acoustic Doppler velocimetry (ADV). Pulsating velocity is a key to identify turbulence characteristics; however, data measured using ADV pertain to the instantaneous velocity. In general, the mathematical mean is usually applied to determine the pulsating velocity [8]. This process of data analysis is known as time-frequency analysis (TFA), which can be conducted using fast Fourier transform (FFT) or wavelet transform (WT). However, FFT is applicable only for linear stationary signals, and prespecified basis functions are required to apply WT. Consequently, FFT and WT cannot be applied to effectively process nonstationary and nonlinear signals [9]. Therefore, an extreme-point symmetric mode decomposition (ESMD) method was proposed to solve the above problems through an optimal global mean curve in an adaptive way and determine the instantaneous frequency and amplitude in a direct way [10]. The ESMD was mostly used in climate research, such as the spatiotemporal variations of nonlinear trends of precipitation [11] or the temporal and spatial patterns of drought in the Yellow River [12]. Ren et al. [13] studied the atmospheric turbulence structure and intermittency during heavy haze pollution in the Beijing area with the method of Hilbert–Huang transform. The ESMD is an improvement of Hilbert–Huang transform. This study was aimed at using the ESMD to extract the properties of nonstationary and transient signals. Subsequently, the turbulence characteristics of a flume flow were analysed based on ESMD.

2. Model Description and Mathematical Method

2.1. Experimental Model

A flume made of plexiglass with a length, width, and height of 232 cm, 16 cm, and 15 cm, respectively, was used. To ensure that the flume flow exhibited a steady state, the ratio of the flume width to its height was approximately 1, and the height of the sections between sections 2 and 12 was 26.5 cm (Figure 1) [14]. The flume was divided into 13 sections. The region between sections 2 and 12 was the key research area, and the distance between these two sections was 15 cm. Sand ( = 39.63 μm) with a concentration of 1 kg/m³ was applied. The grid and gate valve switch were used to realize a fully developed turbulent flow and stable liquid-solid two-phase flow. Acoustic Doppler velocimetry (Vector, Nortek China Ltd., China) was performed at a sample frequency of 25 Hz to test the flow velocity.

2.2. Extreme-Point Symmetric Mode Decomposition (ESMD)

ESMD is an innovative algorithm based on the Hilbert–Huang transform (HHT). ESMD involves two parts: modal decomposition, in which several intrinsic model functions (IMFs) and an optimal adaptive global mean (AGM) can be generated, and time-frequency analysis pertains to the total energy conversion [9, 10, 12].

2.2.1. Algorithm of Modal Decomposition

The purpose of modal decomposition is to yield a series of intrinsic model functions (IMFs) and an optimal adaptive mean curve. The process flow is as follows: Step 1: find all the extreme points of the input data and record the points as sequentially Step 2: connect the adjacent points with lines and record the midpoints of the lines as Step 3: add the left and right boundary midpoints and , respectively. Step 4: construct interpolation lines using the midpoints obtained previously and calculate the mean curves as Step 5: if or the number of sifting times reaches the maximum value , the first empirical mode is decomposed. If not, the abovementioned steps are repeated. Step 6: repeat the abovementioned steps on the residual and obtain the other modes until the final value with only a certain number of extreme points is obtained Step 7: ensure that is within the integer interval and that a series of decomposition modes are obtained through the previous steps. Subsequently, calculate the variance and plot its curve with , where and denote the relative standard deviation of and the standard deviation of the original data, respectively. Step 8: find the maximum number of sifting times corresponding to the minimum variance ratio in the interval . Subsequently, repeat the abovementioned steps to output the decomposition results. The last residual is the optimal adaptive mean curve.

2.2.2. Algorithm of Instantaneous Frequency

When applying FFT and WT, integral operations must be performed on the data, and the original discrete signals must be converted into analytical functions for processing; however, this processing is limited in terms of mathematical theory. The direct interpolation algorithm uses the average frequency over the local period as the interpolation point to generate a smooth curve. It is assumed that the adjustment of the frequency between the adjacent local periods is gradual. This assumption is consistent with the basic methods of data analysis. The direct interpolation algorithm is used to perform time-frequency analysis. The process flow is as follows.

If the discrete form of the intrinsic model functions (IMFs) is , the algorithm to determine the instantaneous frequency involves the following steps: Step 1: determine all of the quasiextreme points of each intrinsic model function that satisfies the following conditions (1). Subsequently, enumerate all the points as in set Step 2: define the interpolation frequency points, considering Step 3: add the boundary points by using linear interpolation. For the left boundary, if , . Otherwise, use linear interpolation to define the boundary. Identify the right boundary in the same manner as the left boundary. Step 4: perform cubic spline interpolation to obtain the interpolation curve by using the time-frequency points. The instantaneous frequency curve can be defined as

3. Results and Discussion

3.1. Signal Acquisition and Analysis

In the experiments, 13 sections were defined in the flume model, as shown in Figure 1. A total of 13 points corresponding to all the sections (sections 2–12) from the water surface to the bed were identified, and the instantaneous velocity on the 143 points was examined through the acoustic Doppler velocimetry technique. In addition, an independent coordinate system was established, whose X, Y, and Z directions were parallel to the flow direction, parallel to the water depth direction, and perpendicular to the flume width, respectively. The X, Y, and Z directions corresponded to the u, , and velocity components, respectively. The velocities at a distance of 3 cm from the riverbed at sections 3, 7, and 11 were considered as the data samples for the analysis. Figure 2 shows the instantaneous velocity sample for the considered location.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 2 shows that the instantaneous velocities fluctuate significantly although they exhibit certain regularity in the flow and water depth directions, as shown in Figures 2(d), 2(e), 2(g), and 2(h). The regularity of the instantaneous velocity at section 7 was the highest; the corresponding regularity at section 11 was low and that at section 3 was the lowest. These phenomena can be attributed to the flow pattern of the section. At section 3, similar to jet flows, a strong turbulence occurs, and the flow has a large amount of kinetic energy; therefore, the instantaneous velocity exhibits a higher randomness (Figures 2(a)–2(c)). Section 7 lies in the middle of the flume, and its velocity is the most representative of the main flume flow characteristics. Although section 11 lies at the end of the main flume, the instantaneous velocity of the section still exhibits certain regularity. This regularity, similar to the Sin function, indicates that the liquid-solid two-phase flow exhibits a high fusion and uniformity. In addition, along the flow direction, the value of the instantaneous velocity gradually decreases (Figures 2(a), 2(d), and 2(g)), which is consistent with the flow characteristics of an actual open-channel flow.

3.2. Optimal Adaptive Global Mean

As described in Section 3.l, the flume exhibits typical pulsating characteristics. An actual open channel is considerably influenced by the pulsating characteristics, such as the vibration and ability to carry sediments. The instantaneous velocity consists of two parts: the mean velocity and pulsating velocity. The pulsating velocity is a key-pulsating characteristic and reflects the difference between the instantaneous and mean velocities. The conventional approach to calculate the mean velocity corresponds to the mathematical expectation, which is defined as method 1 herein [8, 15, 16]. If the instantaneous velocity is , the mean velocity can be obtained using equation (3) (method 1).where is the mean velocity, is the instantaneous velocity, and is the number of points for which the instantaneous velocity is considered.

Although method 1 represents the simplest method to calculate the mean velocity, it is not the optimal method. Specifically, method 1 cannot reflect the variation trend of the instantaneous velocity over time. The method involving the use of the optimal adaptive global mean is defined as method 2. Method 2 is a part of the ESMD technique and corresponds to the result of modal decomposition. The specific algorithm is described in Section 2.2. Figure 3 shows the mean velocities calculated using methods 1 and 2.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

As shown in Figure 3, the mean velocities calculated using methods 1 and 2 are considerably different. The mean velocity calculated using method 1 is a fixed value and does not reflect the fluctuations in the instantaneous velocity. However, the mean velocity calculated using method 2 can accurately reflect the fluctuation in the instantaneous velocity. Moreover, the fluctuation law of the mean velocity obtained using method 2 is similar to the law of the Sin function curve, as shown in Figures 3(a), 3(b), 3(d), and 3(e). This law is consistent with that of the instantaneous velocity described in Section 3.1. Therefore, the optimal adaptive global mean velocity (method 2) can more accurately reflect the instantaneous velocity, regardless of the direction and value. The pulsating velocity can be defined in a more precise manner under a more accurate mean velocity.

3.3. Pulsating Velocity

Figure 4 compares the pulsating velocities calculated using methods 1 and 2. The pulsating velocity calculated using method 1 is equal to the difference between the instantaneous velocity and the mathematical expectation of the velocity. In contrast, the pulsating velocity calculated using method 2 is equal to the difference between the instantaneous velocity and the optimal adaptive global mean velocity. Figures 4(a)–4(y) show the comparison of the pulsating velocities obtained using the two different methods over the entire time period. Figures 4(b), 4(c), 4(e), 4(f), 4(h), 4(i), 4(k), 4(l), 4(n), 4(o), 4(q), 4(r), 4(t), 4(u), 4(w), 4(x), 4(z), and 4(a1) show the segments of the maximum difference in the pulsating velocity corresponding to Figures 4(a)–4(y).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

(p)

(q)

(r)

(s)

(t)

(u)

(v)

(w)

(x)

(y)

(z)

(a1)

As shown in Figure 4, the pulsating velocities obtained using the two methods are consistent in terms of the velocity magnitude, fluctuation range, and fluctuation trend (Figures 4(a), 4(d), 4(g), 4(j), 4(m), 4(p), 4(s), 4(v), and 4(y)). This finding indicates that method 2 is an effective approach to calculate the pulsating velocity. Nevertheless, the pulsating velocities obtained using the two methods involve certain differences, which are more notable in the middle and later segments of the curves, especially in the later segments (Figures 4(c), 4(i), 4(k), 4(l), 4(o), 4(r), 4(u), and 4(x)). To further examine the differences in the data obtained using the two methods, the SPSS software was used to conduct a paired sample test (Table 1). As indicated in Table 1, except for pair 4, the value for all the pairs was less than 0.05, indicating that each pair of data involved significant differences.

3.4. Turbulent Characteristics

Turbulent characteristics are key features of flume turbulence. The relative turbulence intensity (TI), Reynolds shear stress (), and turbulent kinetic energy (TKE) are the most commonly used indexes in practice [15, 17]. The parameters can be calculated as follows:where , , and denote the relative turbulence intensities in directions, respectively; , , and denote the different components of the pulsating velocity; is the average velocity at a given section; and is the flow density. Figure 5 compares the turbulent characteristics determined using methods 1 and 2.

(a)

(b)

(c)

As shown in Figure 5, the three turbulence characteristics calculated using the two methods are highly consistent in terms of magnitude and range. This finding validates that the ESMD method (method 2) is accurate and can be used to analyse the turbulence characteristics. Furthermore, the decomposed characteristics shown in the figure clearly highlight that (1) along the direction of the flow (u3, u7, and u11), the turbulence intensity does not change considerably; the Reynolds shear stress and the turbulent kinetic energy first increase and later decrease, with a considerable change in the turbulent kinetic energy. (2) Along the direction of the water depth (, and ), due to the flume shape, a drop flow occurs (section 3), and the water surface tends to rise after section 11. Subsequently, the turbulence intensity and Reynolds shear stress exhibit negative values. This finding indicates that a strong turbulence occurs at section 7, in which all the three turbulence characteristics exhibit a large magnitude. Vowinckel et al. [18] also observed the same trend and indicated that in the case, the particles altered the flow field, thereby substantially enhancing the turbulence intensity and modifying the size of the coherent particle structure. (3) Along the direction of the flume width (, and ), the turbulence intensity and Reynolds shear stress did not change considerably, thereby indicating that the flume width did not considerably influence the turbulence.

4. Conclusion

In this study, to analyse the feasibility of the ESMD method in determining the turbulence characteristics in flume flow, the traditional mathematical expectation method and ESMD method were used to analyse the turbulence characteristics of a solid-liquid flow in a flume. The results demonstrated that the method can be efficiently used to perform the turbulence analysis of two-phase flows in a flume. The following conclusions were derived from the results:(1)In the context of the traditional mathematical expectation approach to calculate the pulsating velocity, the ESMD method is a new and powerful alternative method to perform the turbulent analysis of solid-liquid flows. The theoretical basis of the ESMD method is correct, and the obtained results are more consistent with the turbulence behaviour in the flume.(2)The optimal adaptive global mean (AGM), which is the result of modal decomposition, can suitably reflect the overall tendency of the instantaneous velocity, along with the turbulence response at different times(3)The turbulence characteristics calculated based on the ESMD method can effectively reflect the turbulence behaviour in the flume which is important to further understand and predict the flow pattern, the recirculation zone or the scour area, and so on.

Nomenclature

:	Median diameter of a solid particle
:	Instantaneous frequency
:	Mean velocity
:	Instantaneous velocity
:	Number of points for which the instantaneous velocity is considered
:	Relative turbulence intensity in the direction
:	Relative turbulence intensity in the direction
:	Relative turbulence intensity in the direction
, , :	Different components of the pulsating velocity
:	Average velocity for a given section
:	Flow density
:	Reynolds shear stress
TKE:	Turbulent kinetic energy.

Data Availability

The data of Figures 2–5 used to support the findings of this study are available from the Jiangsu Vocational Institute of Architectural Technology (szp821002@126.com).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study benefited from the valuable and critical suggestions provided by the referees and the editor. This work was supported by the National Natural Science Foundation of China (51879227), Science and Technology Project of Department of Housing and Urban-Rural Development of Jiangsu Province (2019ZD082), Qinglan Project in the University of Jiangsu Province (2021), and Natural Science Project of Jiangsu Vocational Institute of Architectural Technology (JYJBZX19-05).

References

M. Fettweis, M. Baeye, F. Francken et al., “Monitoring the effects of disposal of fine sediments from maintenance dredging on suspended particulate matter concentration in the Belgian nearshore area (Southern North Sea),” Marine Pollution Bulletin, vol. 62, no. 2, pp. 258–269, 2011.
View at: Publisher Site | Google Scholar
G. Wang, S. Zhou, J. B. Joshi, G. J. Jameson, and G. M. Evans, “An energy model on particle detachment in the turbulent field,” Minerals Engineering, vol. 69, pp. 165–169, 2014.
View at: Publisher Site | Google Scholar
J. Chen, H. Li, X. Lv, and Q. Zhu, “A structure-based drag model for the simulation of Geldart A and B particles in turbulent fluidized beds,” Powder Technology, vol. 274, pp. 112–122, 2015.
View at: Publisher Site | Google Scholar
K. Sun, L. Lu, H. Jiang, and H. H. Jin, “Experimental study of solid particle deposition in 90° ventilated bends of rectangular cross section with turbulent flow,” Aerosol Science and Technology, vol. 47, no. 2, pp. 115–124, 2012.
View at: Google Scholar
X. Liu and M. ASCE, “New near-wall treatment for suspended sediment transport simulations with high Reynolds number turbulence models,” Journal of Hydraulic Engineering, vol. 140, no. 3, pp. 333–339, 2014.
View at: Publisher Site | Google Scholar
N. S. Chen, “Vegetation effects on open channel flows and sediment transport,” Journal of Hydroelectric Engineering, vol. 37, no. 9, pp. 1–7, 2018, in Chinese.
View at: Google Scholar
L. Zhong, T. Jiang, J. M. Zhang, and J. J. Liu, “PIV experimental study on velocity fluctuations of turbulent flow in open channel,” Advanced Engineering Sciences, vol. 51, no. 4, pp. 84–93, 2019, in Chinese.
View at: Google Scholar
C. Venier, J. Figueiredo da Silva, S. J. Mclelland, R. W. Duck, and S. Lanzoni, “Experimental investigation of the impact of macroalgal mats on flow dynamics and sediment stability in shallow tidal areas,” Estuarine, Coastal and Shelf Science, vol. 112, pp. 52–60, 2012.
View at: Publisher Site | Google Scholar
X. Liu, Y. Tang, Z. Lu, H. Huang, X. Tong, and J. Ma, “ESMD-based stability analysis in the progressive collapse of a building model: a case study of a reinforced concrete frame-shear wall model,” Measurement, vol. 120, pp. 34–42, 2018.
View at: Publisher Site | Google Scholar
J. L. Wang and Z. J. Li, “Extreme-point symmetric mode decomposition method for data analysis,” Advance in Adaptive Data Analysis, vol. 5, no. 3, Article ID 1350015, 2013.
View at: Publisher Site | Google Scholar
Y. Qin, B. Li, Z. Chen, Y. Chen, and L. Lian, “Spatio-temporal variations of nonlinear trends of precipitation over an arid region of northwest China according to the extreme-point symmetric mode decomposition method,” International Journal of Climatology, vol. 38, no. 5, pp. 2239–2249, 2017.
View at: Publisher Site | Google Scholar
F. Wang, Z. Wang, H. Yang, and Y. Zhao, “Study of the temporal and spatial patterns of drought in the Yellow River basin based on SPEI,” Science China Earth Sciences, vol. 61, no. 8, pp. 1098–1111, 2018.
View at: Publisher Site | Google Scholar
Y. Ren, H. Zhang, W. Wei, X. Cai, Y. Song, and L. Kang, “A study on atmospheric turbulence structure and intermittency during heavy haze pollution in the Beijing area,” Science China Earth Sciences, vol. 62, no. 12, pp. 2058–2068, 2019, in Chinese.
View at: Publisher Site | Google Scholar
K. M. O. Håkansson, M. Kvick, F. Lundell, L. Prahl Wittberg, and L. D. Söderberg, “Measurement of width and intensity of particle streaks in turbulent flows,” Experiments in Fluids, vol. 54, no. 6, p. 1555, 2013.
View at: Publisher Site | Google Scholar
E. Izadinia, M. Heidarpour, and A. J. Schleiss, “Investigation of turbulence flow and sediment entrainment around a bridge pier,” Stochastic Environmental Research and Risk Assessment, vol. 27, no. 6, pp. 1303–1314, 2013.
View at: Publisher Site | Google Scholar
D. Cava, A. Donateo, and D. Contini, “Combined stationarity index for the estimation of turbulent fluxes of scalars and particles in the atmospheric surface layer,” Agricultural and Forest Meteorology, vol. 194, pp. 88–103, 2014.
View at: Publisher Site | Google Scholar
M. Bizhani, F. E. R. Corredor, and E. Kuru, “An experimental study of turbulent non-Newtonian Fluid flow in concentric Annuli using particle image velocimetry technique,” Flow, Turbulence and Combustion, vol. 94, no. 3, pp. 527–554, 2015.
View at: Publisher Site | Google Scholar
B. Vowinckel, T. Kempe, and J. Fröhlich, “Fluid-particle interaction in turbulent open channel flow with fully-resolved mobile beds,” Advances in Water Resources, vol. 72, pp. 32–44, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Zhi-Peng Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies