Abstract

Nonaxisymmetric wave propagation in an inviscid fluid with a pipeline shear flow is investigated. Mathematical equation is deduced from the conservations of mass and momentum, leading to a second-order differential equation in terms of the acoustic pressure. Meanwhile a general boundary condition is formulated to cover different types of wall configurations. A semianalytical method based on the Fourier-Bessel theory is provided to transform the differential equation to algebraic equations. Numerical analysis of phase velocity and wave attenuation in water is addressed in the laminar and turbulent flow. Meanwhile comparison among different kinds of boundary condition is given. In the end, the measurement performance of an ultrasonic flow meter is demonstrated.

1. Introduction

Wave propagation in an inviscid fluid with a shear mean flow has been studied mostly by numerical and semianalytical methods [1, 2] which is of considerable interest in industrial applications such as noise damping [3, 4] and ultrasonic flow meter [5, 6], to mention a few. On the neglect of the fluid viscosity and thermal conductivity, Pridmore-Brown [7] firstly established the governing equation of the transverse modes in a shear mean flow confined by a rectangular duct. Mungur and Plumblee [8] derived the corresponding convected wave equation in the case of a circular pipeline.

Following the contribution of Mungur and Plumblee [8], various aspects of wave propagation in the circular pipeline containing the shear flow were studied by many researchers, such as Gogate and Munjal [9], Ko [10], Eversman [11], Agarwal and Bull [12], and Pagneux and Froelich [13]. Boucheron et al. [2] made a comprehensive review of axisymmetric acoustic wave propagating in the rigid-walled pipeline. In the case of a lined-walled pipeline where the effect of acoustic impedance of the wall was considered, theoretical and numerical contributions were made by Rienstra and Darau [3], Rienstra [14], Smith et al. [15], Brambley et al. [16], and Buske et al. [17]. Some experiment data can be found in the work of Marx et al. [18].

In the application of an ultrasonic flow measurement, particular considerations are given to improve the measurement performance. Due to the complexity of wave propagation in the pipeline flow, various compensation methods in the engineering aspect have been provided, which was comprehensively summarized by Chen et al. [5]. On the other hand, theoretical analysis based on the wave equation was conducted by Lechner [19], Willatzen [2024], and Chen et al. [6]. Based on the axisymmetric acoustic wave, the performance for the low Mach flow measurement was analyzed in the abovementioned paper. Although the effect of the acoustic impedance at the wall on the measurement performance was analyzed by Willatzen [23], the features of wave absorption at the wall were not addressed.

The present paper concentrates on the nonaxisymmetric acoustic wave propagating in the inviscid fluid with the presence of a shear mean flow confined by a circular pipeline. Particular attention is given to the influence of the acoustic impedance on the wave propagation and measurement performance. Mathematical formulation starts from the conservations of mass and momentum which leads to a second-order differential equation in terms of the acoustic pressure. The boundary condition is then discussed under the assumptions of rigid wall, free-boundary wall (pressure-release wall) [25, 26], and lined wall. In the previous research, the authors introduced a general solution based on the Fourier-Bessel theory [27, 28] to study the axisymmetric acoustic wave propagating along the rigid-walled pipeline [29] and heating pipeline [6, 30]. Such method was extended to handle the nonaxisymmetric wave propagation. Numerical analysis of phase velocity and wave attenuation is addressed in the laminar and turbulent flow. Furthermore, the measurement performance of an ultrasonic flow meter based on the wave propagation is demonstrated.

2. Mathematical Formulation

In this section, the convected equation of the nonaxisymmetric acoustic wave propagating in the inviscid fluid flowing through a circular pipeline is mathematically deduced as shown in Figure 1. The flow profile is assumed to be shear and denoted by . The radial, circumferential, and axial coordinates are expressed by , , and , respectively. The fluid density, velocity, and pressure are denoted by , , and , respectively. represents the pipeline inner radius. In what follows, the thermoviscous effect is not taken into account, and then the entropy perturbation during wave propagation can be taken to be zero.

2.1. Governing Equation

In the presence of wave propagation, acoustic disturbances impose small fluctuating terms on the fluid’s ambient variables which can be expressed by , , and . If these fluctuations are assumed to be linear, according to the work of Mungur and Plumblee [8] and Chen et al. [31], the governing equation of the acoustic wave can be expressed by where and represent the adiabatic sound speed and the local Mach number, respectively. is the normalized radial coordinate through .

The fluctuating quantities can be taken in the form of , where is the angular frequency, the dimensionless axial wavenumber, the free-space wavenumber, and is the circumferential mode which is constrained to be integral due to the periodic wave propagation. The radial component of the acoustic velocity may be expressed by Through expressing the acoustic pressure as , where is the radial component of the acoustic pressure, (1) can be simplified to while the radial acoustic velocity is

Obviously, the nonaxisymmetric acoustic wave propagating in the inviscid fluid with the existence of the pipeline shear flow can be governed by a second-order differential equation with respect to the acoustic pressure. To solve the deduced equation, the boundary condition should be discussed.

2.2. Boundary Condition

On the ignorance of the fluid viscosity and thermal conductivity, the boundary condition can be expressed by the standard Ingard-Myers condition [3, 16]: where denotes the acoustic impedance at the wall and is assumed constant in this paper. Due to the vanishing of the shear flow velocity at the wall, one obtains . Then substituting (4) into the boundary condition yields

According to the work of Rienstra and Darau [3], Brambley et al. [16], and Marx et al. [18], the constant acoustic impedance can be modeled as an enhanced Helmholtz resonator where is the lining depth, is the damping resistance, and is the inertance. In the case of a rigid wall ( ), (6) is simplified to In the situation of a free-boundary wall ( ), (6) is transformed to

As a result, the boundary condition can be generally expressed by where is independent of the radial coordinate. Specifically , , and represent the rigid wall, free-boundary wall, and lined wall, respectively.

3. Solution Based on Fourier-Bessel Theory

In this section, a solution based on the Fourier-Bessel theory, which is mathematically complete and orthogonal in the Lebesgue space [27, 28], is proposed to solve (3) with the boundary condition (10).

Mathematically, the bounded continuous function in the domain may be expressed in terms of a Fourier-Bessel sequence where denotes the th-order Bessel function of the first kind and corresponds to the circumferential mode. Substitution into (10) leads to which can determine the coefficients . Meanwhile, the orthogonality [27, 28] of the Fourier-Bessel sequence may be represented by where is the Kronecker delta function. If the acoustic pressure is given, the corresponding coefficients can be obtained by which shows that is independent of the radial coordinate.

Obviously, the Bessel functions are solutions of cylindrical waves without flow, which are not the eigenfunctions with shear flow. While the eigenfunctions without flow ( ) are orthogonal and complete, the eigenfunctions with shear flow may be not. As a result, (11) is a valid orthogonal representation but is not an eigenfunction decomposition. It specifies pseudo-modes that are the approximate representation of the acoustic waves with shear in terms of acoustic waves without shear. The calculation of the coefficients and pseudo-eigenvalues minimizes the mean square error of the representation. As a result, substituting (11) into (3) yields In the above equation, the properties of Bessel function are used [27, 28]. Similar deductions can be found by the authors’ previous papers [6, 29, 30]. Multiplying and then integrating in the domain result in where

Obviously, given specific local Mach profile , acoustic angular frequency , pipeline radius , fluid density , and the orthogonal Bessel sequence with the corresponding eigenvalues , being the roots of (12), expanding (16) gives a set of homogeneous algebraic equations with the unknown parameters and the dimensionless axial wavenumber . It is obvious that if the number of equations expanded from (16) is , the number of coefficients is . Clearly, the problem is homogeneous and the wave amplitude can be multiplied by an arbitrary constant.

While the specific coefficients cannot be solved from (16), the dimensionless axial wavenumber can be obtained. Furthermore, the ratio of all to one of them such as (if ) can be calculated. As we are interested in the wavenumber , by taking the coefficients as unknown variables with (the subscript “T” means the operation of matrix transpose), (16) can be written as where is a matrix whose element is a function of the dimensionless axial wavenumber. From (14), it can be learned that the coefficients do not vanish simultaneously, leading to the condition of . The existence condition of a nontrivial solution to (16) is then the vanishment of the corresponding determinant which can be solved by various numerical methods. A Newton-related method can be found in the authors’ papers [6, 29, 30].

4. Numerical Analysis of Wave Propagation

In this section, wave propagation in water is numerically studied. The corresponding constant parameters are , , , and . Following the work of Rienstra and Darau [3], the parameters of acoustic impedance are assumed to be , , and . The laminar flow profile can be described by a Poisson’s equation while the turbulent flow profile can be represented as the seventh-root law where is the cross-section average Mach number. Due to the acoustic impedance of the lined wall, the acoustic energy is absorbed at the wall. The following analysis focuses on the phase velocity ( ) and attenuation coefficient ( ) where the subscripts “ ” and “ ” denote the real and imaginary components of the axial wavenumber, respectively. To get a dimensionless description of the phase velocity, the relative phase velocity is used in the following discussion. Furthermore, the number of Bessel functions is set to and a convincible calculation convergence of the axial wavenumber can be reached [6, 29, 30]. In what follows, the mean Mach number is assumed to be in the interval . Readers may calculate the axial wavenumber for any value of the mean Mach number using presented method.

4.1. Phase Velocity
4.1.1. Comparison among Different Kinds of Boundary Condition

Figure 2 represents the relative phase velocity under the three different wall configurations (explained in Section 2.2) in the laminar flow, while Figure 3 shows the corresponding relative phase velocity in the turbulent flow. Numerical comparison is imposed on the relative phase velocity of the first four axisymmetric ( ) radial pseudo-modes as a function of the average Mach number. In the legends, ‘‘down/up’’ means wave propagation in the downstream/upstream direction while ‘‘free/gen/rigid’’ represents the free-boundary/lined/rigid wall configuration. The first four axisymmetric radial pseudo-modes are denoted by “1,” “2,” “3,” and “4,” respectively.

When the fluid is at rest, the relative phase velocity of each axisymmetric acoustic mode in the rigid wall turns out to be the largest among the three types of the boundary condition. The relative phase velocity in the lined-walled pipeline shows the smallest. In the case of the high modes, the difference of the relative phase velocity among the three boundary types becomes more distinct. Physically speaking, the complicated refraction and reflection at the wall lead to the existence of the high modes (“2,” “3,” and “4”). With a higher mode, the corresponding wave confronts more complex reflections and refractions. The wall has influence on the acoustic wave and then becomes more obvious, leading to a slower propagation speed in the lined wall. However, the wall influence on the fundamental mode (“1”) is slight.

In the case of a moving fluid, the relative phase velocity increases with the average Mach number in the downstream propagation but decreases in the upstream propagation. In the turbulent flow, the relative phase velocity in the downstream propagation is semisymmetric to that in the upstream propagation with respect to the case of the static fluid. However, the symmetric relationship is altered in the laminar flow due to the shear effect. Such phenomenon is more obvious in the fundamental mode wave than the other higher acoustic modes.

Physically speaking, with the increase of the average flow Mach number, the convective flow (uniform flow) accelerates the downstream propagation as the propagation direction is along the flow direction but decelerates the upstream propagation. The consequence is that the relative phase velocity in the downstream and upstream propagation due to the flow convection shows to be symmetric with respect to the static fluid case. Such symmetric relationship can be found in Figure 3 as the turbulent flow profile is similar to the convective flow profile.

In the laminar flow, the phenomenon becomes more complicated as the shear effect alters the abovementioned symmetric relationship. From Figure 2(a), it can be inferred that the shear effect channels the fundamental mode (“1”) into a narrow layer near the wall in the downstream propagation. In the upstream propagation, the acoustic wave is channeled into a narrow layer near the pipeline center [7]. As the local Mach number near the wall is smaller than that near the pipeline center, the contribution of the flow profile on the propagation speed is slight in the downstream propagation. In the upstream propagation, the increased local Mach number near the center heavily decelerates the propagation speed. For the high acoustic modes (“2,” “3,” and “4”), the shear effect channels the acoustic wave into a narrow layer next to the centerline. As a result, the local Mach number near the pipeline center accelerates the downstream propagation but decelerates the upstream propagation. Such phenomenon can be seen by comparing the relative phase velocity between Figures 2 and 3. It should be noticed that the shear effect exists in both the laminar and turbulent flow. However, such influence is more significant in the laminar flow as the fluid vorticity is larger for the parabolic velocity profile (20) than for the logarithmic (21) flow profile.

It should be noticed that the shear effect on wave propagation confined by the free-boundary wall is less obvious. On the other hand, the existence of acoustic impedance at the wall prefers to regulate the wave propagation to a narrow layer near the pipeline wall. As a result, in the downstream propagation, the local Mach number near the pipeline wall slightly accelerates the propagation speed compared with that near the pipeline center. In the upstream propagation, the local flow profile near the wall weakens the decreasing rate of the relative phase velocity due to the flow convection. Such phenomenon can be found in both the laminar and turbulent flow as shown in Figures 2 and 3.

4.1.2. Comparison among Different Circumferential Modes

Figure 4 compares the relative phase velocity of the first four axial pseudo-modes (“1,” “2,” “3,” and “4”) among the three different circumferential modes ( , , and ) in the laminar flow confined by the lined-walled pipeline, while Figure 5 shows the corresponding comparisons in the turbulent flow.

In the static fluid ( ), the relative phase velocity of a higher circumferential mode is larger. In the case of a higher axial mode, the difference of the relative phase velocity among the three circumferential modes becomes more visible. In the presence of a flow, the flow convection speeds up the phase velocity of the acoustic wave propagating in the downstream direction but slows down the relative phase velocity in the upstream propagation. The influence of flow convection on each acoustic mode can be explicitly found in the case of turbulent flow in Figure 5 as the core profile of the turbulent flow is nearly the same with the uniform flow profile.

In the laminar flow (Figure 4), the tendency becomes more complicated compared with that in the turbulent flow as shown in Figure 5. For the first axial pseudo-mode (Figure 4(a)), the shear effect regulates each circumferential mode into a narrow layer next to the wall in the downstream propagation; meanwhile, it channels each circumferential mode into a narrow layer next to the pipeline center in the upstream propagation. By comparing the tendencies of the three circumferential modes as functions of the average Mach number, it can be learned that the shear effect shows the most significant influence on the axisymmetric mode ( ). In the case of higher axial modes (Figures 4, 5(b), 5(c), and 5(d)), the influence of shear flow behaves similarly as shown in Figure 2. Comparison among the three higher axial modes reveals that the influence of wall impedance tends to channel the wave mode propagate near the wall. Such tendency is more obvious for a higher axial mode wave as shown in Figures 4(d) and 5(d).

4.2. Acoustic Attenuation
4.2.1. Comparison between Laminar and Turbulent Flow

In this subsection, the attenuation coefficient due to the acoustic impedance at the wall is investigated. Specifically, Figure 6 compares the attenuation coefficient of the first four axial modes between the laminar (“lam”) and turbulent (“tur”) flow profiles in both the downstream (Figure 6(a)) and upstream (Figure 6(b)) propagations. The circumferential wave mode is assumed to be ; then the feature of the axisymmetric acoustic wave is addressed.

In the absence of a moving fluid ( in Figures 6(a) and 6(b)), a higher axial mode leads to a larger attenuation coefficient. As shown in Section 4.1.1, the relative phase velocity of a higher axial mode is larger which corresponds to shorter wavelength, steeper waveform, and larger attenuation. In the presence of a flow, the behavior of attenuation coefficient shows significant differences between the first axial mode and the other three axial modes.

For the first axial mode propagating in the downstream direction (Figure 6(a)), the attenuation coefficient becomes large in both the laminar and turbulent flows with the ascent of the average Mach number. Meanwhile, the ascending rate of attenuation coefficient in the laminar flow is quicker than that in the turbulent flow. As described in Figure 3, the shear effect regulates the first axial mode to a narrow layer near the wall in the downstream propagation. Then, the wave damping due to the acoustic impedance becomes larger in the presence of a higher Mach number. As the shear effect in the laminar flow is more noticeable, the wave attenuation in the laminar flow behaves severer. In the upstream propagation (Figure 6(b)), the attenuation coefficient slows down against the average Mach number in the laminar and turbulent flow. The attenuation descent is quicker in the laminar flow than that in turbulent flow. As the shear effect channels the acoustic wave in the upstream propagation into a narrow layer next to the pipeline center, the wave damping at the wall becomes slight. As the shear effect in the laminar flow is more significant than that in the turbulent flow, the decrease ratio of wave attenuation becomes quicker in the laminar flow.

Comparison of attenuation coefficients among the three high modes propagating in the downstream direction (Figure 6(a)) shows that the flow convection eases the wave attenuation, but the shear flow complicates the tendency of attenuation coefficients. In the upstream propagation, the flow convection slows down the propagation speed; then the occurrence probability of the wave damping at the wall increases, leading to an ascending of the attenuation coefficient. On the other hand, the shear effect regulates the high acoustic wave in the upstream propagation to a narrow layer near the pipeline center, leading to a decreasing of the wave damping at the wall. As the shear effect in the laminar flow is stronger than that in the turbulent flow, the average Mach number corresponding to the decrease of wave attenuation is smaller in the laminar flow. It should be noticed that the shear influence on a higher acoustic mode is less significant as the higher acoustic mode leads to a more complex propagation path.

4.2.2. Comparison among Different Circumferential Modes

Figure 7 compares the attenuation coefficients of the first four axial wave modes among the first three circumferential modes ( , , and ) in the laminar flow while Figure 8 displays the corresponding comparison in the turbulent flow.

In the absence of a shear flow, the attenuation coefficient of a higher circumferential mode is larger for each axial mode. The difference of the attenuation coefficients among the three circumferential modes becomes more visible in the case of a higher axial mode. However, such tendency is not suitable for the first axial mode.

In the presence of a shear flow, the attenuation coefficient changes due to the existence of the shear and convective effects. Specifically, the influences of shear and convective actions on the other two circumferential modes ( and ) are similar to the case of the axisymmetric mode ( ). It should be noticed that the shear effect on the attenuation coefficient of the axisymmetric acoustic mode shows more sensitivity compared with the case of the other two circumferential modes.

5. Performance of Ultrasonic Flow Measurement

In this section, the measurement performance of an ultrasonic flow meter based on the wave propagation is analyzed theoretically. To describe the measurement performance explicitly, the measurement error is firstly defined [5, 19, 22]. Specifically, when a plane wave propagates in an inviscid fluid with a uniform flow ( ), travelling a distance in the downstream direction requires while the consumed time in the upstream direction is Subtracting (23) from (24) gives

Under the assumption of , one obtains and the volume flow ( ) then can be obtained by

For a specific wave mode ( ), from (4), traveling a distance requires Then we have and the corresponding measured volume flow is Thus, the measurement error of a specific wave mode can be expressed by

It can be inferred that the measurement error represents the symmetric characteristics of the real component of the axial wavenumber in the downstream and upstream propagations with respect to the case of the static fluid.

5.1. Comparison among Different Kinds of Boundary Condition

Figure 9 presents the difference of the measurement errors of the first four axial pseudo-modes among the three different boundary configurations (free-boundary/lined/rigid wall) in the laminar flow, while Figure 10 details the corresponding comparison in the turbulent flow. The circumferential mode is assumed to be ; thus the feature of the axisymmetric acoustic wave is analyzed. Obviously, the measurement error in the lined-walled pipeline is the smallest in the laminar and turbulent flow. Furthermore, the measurement error of each acoustic mode is higher in the laminar flow than that in the turbulent flow.

In the turbulent flow (Figure 10), the measurement error of each axial mode grows up semilinearly with the mean Mach number whatever the boundary condition is. As the influence of flow convection on the relative phase velocity in the downstream and upstream wave propagations is semisymmetric with respect to the case of the static fluid (Figure 3), the linearity of the measurement error in terms of the average Mach number is maintained. For a higher axial mode, the variation range of the measurement error with respect to the average Mach number becomes shorter. A possible interpretation may be that the complex reflection and refraction at the fluid-wall interface for the high acoustic mode result in a complicated propagation path in the pipeline. Then, the inconsistent effect of flow profile on the acoustic wave propagating in downstream and upstream directions is compensated. Such tendency becomes more obvious in the laminar flow as shown in Figure 9.

Figure 9 reveals that a higher axial mode has a slower variation of the measurement error. Comparing with the measurement error in the turbulent flow shows that the shear effect sharply enlarges the measurement error. The complex propagation path of a higher axial mode due to the reflection and refraction compensates the shear effect on the measurement error. The mechanism of the acoustic impedance decreases the measurement error compared with the cases of the rigid and free-boundary configurations.

With the increase of the mean Mach number, the shear effect on the measurement error begins to play the dominant role which worsens the measurement performance. As a result, the measurement error in a lined-walled pipeline firstly decreases against the average Mach number but finally increases. The complex reflection and refraction ease the shear effect and then postpone the increase ratio of the measurement error. In the literature, Lechner [19] and Willatzen [2024] did similar work in the case of very small Mach number. Their work reveals that the measurement error turns to be insensitive to the Mach number. From Figure 9, it can be learned that their contributions are consistent with the present study.

5.2. Comparison among Different Circumferential Modes

Figure 11 compares the measurement errors of the first four axial modes among the three circumferential modes in the laminar flow confined by the lined-walled pipeline and Figure 12 reveals the corresponding comparison in the turbulent flow.

Clearly, the shear and convective effects on the measurement performance of the two nonaxisymmetric modes ( and ) are similar to the case of the axisymmetric mode ( ). Comparison among the three circumferential modes reveals that the measurement error of the axisymmetric mode is the most sensitive. For each axial mode, the measurement error of a lower circumferential mode is larger. However, such tendency is not suitable when the Mach number is large.

In the end, it can be learned that the effect of flow convection on the phase velocity in the downstream and upstream propagations keeps a symmetric way which can be inferred from Figures 3 and 5. In the presence of a shear flow, the inconsistency of phase velocity between the downstream and upstream propagations is obvious as shown by Figures 2 and 4. The shear effect severely alters the influence of flow convection. It should be pointed out that in many industrial applications, the actual flow Mach number may be small. Then the measurement error of each mode in the laminar and turbulent flow changes semilinearly with respect to the average Mach number. Furthermore, with the increase of the flow Mach number, the turbulent flow profile seems to be a better candidate to describe the actual condition than the laminar flow.

6. Conclusion

In this paper, the nonaxisymmetric wave propagation in the inviscid fluid with a shear mean flow confined by a circular pipeline is studied while the acoustic wave is assumed to be isentropic and linear. The mathematical formulation is deduced into the form of a second-order differential equation with respect to the acoustic pressure. A general boundary condition is presented to handle three different configurations with free-boundary, rigid, and lined walls respectively. Using the Fourier-Bessel theory, a method is provided to numerically calculate the axial wavenumber, which describes the phase velocity and wave attenuation. Numerical analysis reveals the following.(1)The phase velocity of wave propagation in the rigid-walled pipeline turns out to be the largest among the three types of boundary condition while the phase velocity in the lined-walled pipeline is the smallest.(2)Wave propagation with a higher circumferential mode has a larger phase velocity. In the case of a higher axial mode, the distinction of the phase velocity among the circumferential modes is more explicit.(3)The effect of flow convection accelerates the wave propagation in the downstream direction and decelerates the wave propagation in the upstream direction. Meanwhile, the corresponding attenuation coefficient increases in the upstream propagation and decreases in the downstream propagation. However, the shear effect greatly alters the influence of flow convection on both phase velocity and attenuation coefficient.(4)Acoustic wave with a higher circumferential mode results in severer wave attenuation. Furthermore, a high circumferential mode postpones the shear effect on the wave attenuation.

In the end, the measurement performance of an ultrasonic flow meter is numerically investigated. The effects of circumferential mode, boundary condition, and shear profile are parametrically demonstrated.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work described in this paper is funded by the National Natural Science Foundation of China (nos. 91216201 and 51205403). The authors gratefully acknowledge the funding.