Abstract
Duct noise control is an important practical problem. This paper explores the vibro-acoustic mechanism of duct noise through the membrane in tandem. Validity and feasibility of the proposed analytical model is demonstrated by comparing existing simplified models and full direct aeroacoustic simulation solved with the CE/SE model. It is shown that the coupling effect between two membranes in tandem is not negligible to predict system response. Moreover, introduction of multimembranes is important or even the only efficient way to apply this passive control method in practice.
1. Introduction
How to obtain better duct noise control method has long been a long-term studying problem in industrial field [1–3]. It is known that two traditional passive control methods (i.e., expansion-chamber design and duct lining) have their own disadvantage. Specifically, there exist periodic pass bands for the expansion-chamber case, while the duct lining method is not efficient in the low frequency range. In 2001, Huang [4] creatively proposed a new passive control method theoretically based on wave reflection and structural damping of a flexible membrane inspired by the investigation on flexible tube in human respiratory systems. Furthermore, a series of systematic research studies have done by Huang et al., such as investigation on detailed physical mechanism of flexible membrane [5], drum-like design to avoid break-out noise [6], and practical experimental test work. Besides, other researchers also conduct many valuable works recently. Du et al. [8] have discussed the influence of different boundary constraints on sound transmission loss by a modified improved Fourier series method. Yu et al. [9] replaced ordinary membrane with an electro-active membrane (i.e., dielectric elastomer) and extended this passive control method into adaptive-passive control strategy. Recently, Fan et al. [10–12] conducted systematic numerical analysis on Huang’s flexible membrane model from the viewpoint of aeroacoustic-structural interaction based on the conservation element and solution element (CE/SE) method; they found that the presence of flow with the incident wave in the flow direction mainly weakens the reflection.
For the practical application of flexible membrane, the request of high tension force may be one of major barriers since larger tension force can produce better sound insulation effect [4]. To solve this problem, our idea is to introduce another membrane in tandem. Both an analytical model and CE/SE method are built to explore the role of additional membrane in consideration of complete coupling effect between two membranes. In detail, the frequency domain model and time domain model are described in Section 2 and Section 3, respectively. Validation of the proposed model and corresponding discussion are given in Section 4. Conclusions are summarized in Section 5.
2. Frequency Domain Model
As shown in Figure 1, the first membrane is assumed to be hinged at , where is the dimensionless membrane length. An incident plane sound wave is assumed coming from left to right with unit amplitude:
Governing dynamic equations of the membranes 1 and 2 are written aswhere , Mi, Di, and Ti (i = 1, 2) are, respectively, displacement, mass ratio, dimensionless damping coefficient, and tension force. Z11 and Z22 are local radiation impedance structure itself. Z21 and Z12 are mutual radiation impedance which reflect the coupling effect in this kind of the structural acoustic system. The standard Galerkin method is employed to solve equations above. If velocity is defined, then Vi can be expanded in terms of sine transform using two local dimensionless variable and ; then, both and span from 0 to 1 for membrane 1 and membrane 2:
After sine transform, equations (2) can be reduced to a set of linear equations:where , , and are wave speed of these two membranes, respectively, and is the sine transform of incident wave, i.e., and . Detailed expressions of local radiation impedance can be found in Huang’s paper and not shown here for brevity. And, mutual modal radiation impedance will be given in the following.
2.1. Mutual Modal Radiation Impedance
The starting point is the calculation of radiated sound pressure by vibrating membrane. Formula [4] is used as follows:where H is the Heaviside function, and are coordinates of exciting point and response point, is the channel modal function, is the Kronecker delta, and is the complex modal wave speed. Then, mutual modal radiation impedance Z21 and Z12 can be expressed aswhere channel modal function has degenerated to since exciting and response point located at y = 0 and . It is noted that and always exist once the present coordinate system is built; then, the Heaviside function can be easy to calculate. The rest of the derivation process is routine and not be repeated here.
2.2. Sound Transmission Loss
Once the linear equation (4) are solved, according to equation (5), the expanded velocity terms can be used to predict radiated sound pressure , which consists of two parts from membrane 1 and membrane 2; then, the transmission loss is defined as
3. Time Domain Model
An aeroacoustic model based on direct aeroacoustic simulation (DAS) approach is also employed in the present study. The aeroacoustic problem is governed by the two-dimensional compressible Navier–Stokes equations together with ideal gas law for calorically perfect gas. The normalized Navier–Stokes equations without the source can be written in the strong conservation form as follows:where , and
is the density of fluid, and are the velocities in x and y direction, respectively, is the time, stress , and , total energy , pressure , heat flux , , the specific heat ratio , Mach number , where is the duct mean flow velocity, and , specific gas constant for . Choosing a reference length , a reference density , a reference acoustic speed , Reynolds number , and Prandtl number , the Navier–Stokes equations are solved by the CE/SE-based DAS solver. Detailed implementation progress can be found in Fan and Leung’s another work [10].
4. Numerical Results and Discussion
4.1. Validation of Model
Firstly, Huang’s model [4] is employed to validate present vibro-acoustic model preliminary. In the following discussions, the same parameters are chosen as the default. The height of duct , mass per unit length of membrane , tension force , the length of membrane 1 and 2, and space is . Since Huang’s analytical model refers to one membrane setup, the second membrane is neglected in real calculation to approximate Huang’s result. As shown in Figure 2, good agreement between the present degenerative model and Huang model is achieved, which validate present frequency domain model primarily.
The significant STL reduction at frequency f = 0.36 is attributed to the critical frequency phenomenon where wavenumber of the flexural wave is equal to acoustic wavenumber of surrounding fluid. If frequency is higher than critical frequency, supersonic radiation sound wave cannot provide effective wave reflection; then, most of sound wave energy will propagate to downstream except consumed part because of structural damping.
4.2. Frequency Domain Model vs. Time Domain Model
Furthermore, the present frequency domain model is compared to the time domain method to further illustrate the validity of the model. Here, two acoustic response point at and are chosen to reflect the sound insulation effect of these two membranes, as shown in Figures 3 and 4, respectively. Generally, it is shown that excellent agreement between present two models can be achieved including oscillations’ peaks and troughs on the curves, which validate the present analytical model and the numerical model effectively again. The quantitative difference between them is mainly attributed to the structural damping model of the membrane. Specifically, the modal damping model in the frequency domain model cannot be realized completely in the time domain model. More detailed explanation can be found in [10–12].
4.3. Frequency Domain Model vs. Transfer Matrix Method
For the multimembrane cases, another straightforward way is to use the transfer matrix method, which neglects the coupling effect between different membranes, as shown in Choy and Huang’s work [13]. Therefore, it is significant to discuss the coupling effect of membranes in tandem through the comparison between the transfer matrix method and the present coupled model. Intuitively, the coupling effect varies with the distance of two membranes. As shown in Figure 5, the discrepancy between two models decrease with increasing spacing when incident dimensionless frequency f = 0.2, and the maximum gap (about 4 dB) happens when the gap is zero, i.e., the membranes are arranged continuously. By the way, the coupling effect for the whole system is embodied in the mutual modal radiation impedance terms Z12 and Z21 from the mathematical point of view.
Furthermore, the variation of STL with incident sound frequency is shown in Figure 6 for gap Lm = 0. It is found that there is no difference between these two models in low frequency range, which means coupling effect can be neglected in this situation. In the middle frequency range smaller than critical frequency, the discrepancy remains the same and 5 dB approximately. Notably, the quantitative deviation varies with specific structural and fluid parameters (i.e., heavy water case), and present numerical discussion highlights the coupling radiation effect owing to two membranes adequately.
4.4. Membrane in Tandem vs. Single Membrane
In contrast to the traditional one-membrane design (L = 20), the predictable advantage of multimembranes (L1 = Lm = L2 = 10) is larger transmission loss and lower requirement of tension force to maintain certain sound insulation effect if total length of membranes is fixed. As shown in Figure 7, STL of the membranes in tandem is higher than the single membrane case in the low frequency range (i.e., f < 0.24). However, in the frequency range between f = 0.24 and critical frequency region, the single membrane setup exhibits better sound insulation property than the membranes in tandem factually, which shows that it may be invalid to divide one complete membrane into two pieces for better sound insulation.
5. Conclusions
A completed coupling frequency domain model is developed to investigate the sound transmission loss property of the membrane in tandem. The validity and feasibility of the model is verified by comparing Huang’s classical single membrane model and the authors’ time domain model. The present model is capable of capturing the coupling effect between double membranes and helpful for further parameters’ discussion and practical design.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (11502110), Open Fund of the State Key Laboratory for Strength and Vibration of Mechanical Structures, Xi’an Jiaotong University (nos. SV2018-KF-01 and SV2018-KF-22), and Open Project of Space Structure and Mechanism Technology Laboratory of China Aerospace Science and Technology Group Co., Ltd. The third and fourth authors are grateful to the support from the Research Grants Council of the Government of Hong Kong Special Administrative Region, under Grant no. A-PolyU 503/15.