Issue 
Acta Acust.
Volume 8, 2024



Article Number  7  
Number of page(s)  11  
Section  Aeroacoustics  
DOI  https://doi.org/10.1051/aacus/2023062  
Published online  01 February 2024 
Scientific Article
Experimental and numerical observation of flowacoustics feedback phenomena due to two diaphragms in tandem inserted in a rectangular duct
^{1}
Université de Technologie de Compiègne, Laboratoire Roberval, BP 20529, 60205 Compiègne Cedex, France
^{2}
ESIGroup, 39 Chemin Des Ramassiers, ZAC Des Ramassiers Les Ailes de l’Europe, 31170 Colomiers CEDEX, France
^{3}
Naval Group, Technopôle de la mer, 199 Avenue PierreGilles de Gennes, 83190, Ollioules, France
^{4}
Université de Sherbrooke, 2500 boulevard de l'Université, Sherbrooke, Canada
^{*} Corresponding author: s.de.reboul@gmail.com
Received:
26
April
2023
Accepted:
6
November
2023
The proximity of obstacles in a ventilation network can result in a significant increase of the noise production. In this paper, a tandem diaphragm inserted in a rigid duct with rectangular cross section is investigated and compared to a single diaphragm obstruction. The disturbed flow which reaches the second obstacle causes an amplification of the radiated noise which goes beyond a simple doubling of the power. Experimental results, carried with various spacing between the two diaphragms and bulk velocity, allows identifying other aeroacoustic source mechanisms which does not appear with a single diaphragm. Depending on the configuration, the broadband level can increase by more than 10 dB and feedback phenomena can appear at low frequency. Explanations and illustrations of those phenomena are given thanks to compressiblefluid simulations.
Key words: Duct acoustics / Flowinduced noise / Direct noise computation / Feedback loop
© The Author(s), Published by EDP Sciences, 2024
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The Heating Ventilation and Air Conditioning (HVAC) system of a car has to provide airflow to ensure the passengers comfort regarding the temperature inside the vehicle cabin without damaging the acoustic environment. To obtain a satisfactory indoor environment, the acoustical design of a ventilation system is as important as its thermal design and becomes even more significant since the development of the electrical and hybrid electric vehicles. In a HVAC, the noise is caused by the complicated acoustic and turbulent flow interactions with multiple discontinuities located in a compact housing duct at a low mean velocity interacting mainly with a moving element (the blower) but also with fixed elements such as the heat exchanger, the mixing and distributing flaps, the filter, bends, constrictions and finally with the outlet grille before radiating into the passenger compartment. Because of the compactness of the HVAC housing required by the car manufacturers, high amplification and whistling can sometimes be produced [1, 2] and there is a need for an appropriate description of the physical phenomenon which causes these amplifications. The prediction of the aeroacoustic sound power spectra (SWL) of single elements in ducts [3, 4] such as spoilers, bends, orifices is usually based on the NelsonMorfey’s theory [5]. This formulation assumes a dipole source distribution resulting of the drag fluctuating forces arising from the turbulent flow near the element and that the ratio between these forces and the static drag force depends on the Strouhal number. This approach is very convenient from an industrial point of view because the SWL can be easily predicted for a wide range of configurations thanks to similarity considerations [6]. The prediction of noise due to the presence of multiple elements in air ducts have also been the subject of research in the past. The theory, which is also based on Nelson and Morfey’s formulation, relies on the crossspectra between the fluctuating forces between each element [7, 8]. These predictive models are however based on the fluctuations of the aerodynamic forces on the obstacles and the scattering of acoustic waves as well as possible resonance effects are usually ignored in the analysis. Depending on the situation, these resonances can be purely acoustical [9] or can result from a more complex flowacoustic phenomena whereby a vortex shedding mechanism and the resonant acoustic frequencies are locked together [10]. Although this is not the place for a complete survey on this difficult matter, one can cite the work of Ziada and Lafon [11] in the field of industrial aeroacoustics. The authors identified three different mechanisms which generate upstream feedback of disturbance: (1) structural vibrations, (2) fluidresonant mechanism triggered when the frequency of vortex shedding becomes close to that of an acoustic resonator, and (3) the fluiddynamic mechanism caused by flow impingement on a downstream object.
The aim of this paper is twofold: to show that mechanism (3) can be responsible for the production of high level sound power around specific frequencies because of the closeness of the HVAC elements and also to give all necessary details allowing to capture these fluidacoustic interaction correctly via the use of a commercial Computational Fluid Dynamics (CFD) software. The academic study consists of two identical diaphragms with rectangular openings inserted in a duct of rectangular crosssection and separated by three different distances which are comparable with the transverse dimensions of the duct. It is interesting to note that similar configurations, mainly conducted with circular geometries both for the duct and the orifice, have been the subject of many papers in the recent years and we can cite the work of Sovardi et al. [12], Karban and Schram [13, 14], Sack et al. [15, 16], and also [17]. In a somewhat more industrial context, other research work investigated the ability of different CFD software to predict broadband noise and the presence of tonal noise, see [18–20] for instance. Both experimental and numerical aspects are presented here in a rectangular duct. The article is organized as follows: in Section 2, experimental results are shown in terms of power spectra for the sound power radiated downstream by the obstacle which consists either of a single diaphragm or two diaphragms in tandem. It is shown that the latter configuration can result in a significant increase of the noise production with the presence of peaks which are typical of flowacoustics resonance or acoustic feedbackloop phenomena. This is evidenced by matching the frequencies of the peaks with a classical feedbackloop equation as found for cavity or airfoil tonal noise. In Section 3, numerical calculations using LES simulations are performed, analyzed and compared with experimental results. The modal pattern observed at the resonant peaks, signature of large coherent structures formed downstream of the diaphragms, is a further evidence of the flowacoustics resonance. The paper ends with a conclusion and some discussion.
2 Experimental results
2.1 Experimental setup and measurement procedure
Although the detailed description of the test set up can be found in [21], it is worth recalling its main characteristics. The flow generation part is composed of a fan, mufflers, a plenum chamber and a hot wire air flow meter. Two mufflers connected with an elbow located upstream of the blower are designed to reduce fan and flow noise. The air passes through a flow tranquilizing box which minimizes the turbulence and a flow nozzle equipped with a differential pressure gauge to measure the volumetric flow rate and sensors have been used to measure the humidity as well as the temperature of the environment. The measurement part (visible in Fig. 1) includes a total of 96 flush mounted microphones located upstream and downstream the test section. Figure 1 also shows two anechoic terminations placed on both sides of the duct.
Figure 1 Sketch of the AeroAcoustics test rig. 
Acoustic measurements are based on the 2N port method [22]. The formulation takes the form of a matrix system (here ports 1/2 correspond to the left/right part of the test section)
where is the scattering matrix and the source vectors and contain the amplitudes of the duct acoustic modes propagating downstream (+) and upstream (−).
Pressure signals are recorded by a 96 channels LMS Scadas acquisition with a bandwidth of 6400 Hz and 4096 spectral lines. The measurement procedure is performed using an averaging process over 600 realizations to achieve statistical convergence (the acquisition time is about 2 min long for each configuration). The acoustic scattering matrix of the obstructed duct and the impedance of the surrounding environment are first identified by measuring the acoustic response to given external excitations. The test rig was designed to perform modal decomposition up to the 8th propagating mode (N = 8). Therefore, the computation of the scattering matrix must require at least 2N configurations which must be linearly independent. This is achieved by using 10 loudspeaker positions both upstream and downstream the test section (see Fig. 1). In the case of a single diaphragm inserted in the duct, measured scattering matrix coefficients have been validated by comparison with numerical simulations [21].
In a second step, the aeroacoustic noise produced by the flowobstacle interaction can be measured. In order to obtain a formulation for both random and periodic types of signals, the sound power from values of the source vectors is obtained via the modal crossspectrum matrix [15, 22]. The formulation extracts the source strength using crossspectrum densities between the four different sections S_{1}, S_{4}, S_{7}, S_{8} according to the procedure first described in [21] and also applied in [23] (see Appendix B in [23] for all necessary details). The procedure allows us to suppress the hydrodynamic fluctuations, particularly on the downstream side where the flow is expected to be turbulent.
The convection of the propagating waves due to the flow can be neglected here and the sound power radiated downstream by the obstacle is estimated via
where k_{mn} are the longitudinal wavenumbers associated with mode (m, n), i.e.
where k = ω/c_{0} (note the summation is restricted to propagating modes only, i.e. k_{mn} must be a real quantity). Here, are the amplitudes of the duct acoustic modes from the source vector. The dimensions of the crosssection of the duct are a = 20 cm and b = 10 cm. The geometry of the diaphragm inserted in the test section consists of a constriction of rectangular shape with dimension 10 cm × 5 cm × 0.8 cm centered in the rectangular duct (Fig. 2). In the case of two identical diaphragms in tandem, three separation distances with L = 6 cm, 13 cm, and 21 cm are considered.
Figure 2 Two diaphragms in tandem and dimensions of a single diaphragm. 
2.2 Case of a single diaphragm and scaling law for the fundamental mode
The hypothesis used in Nelson and Morphey paper is that the root mean square fluctuating drag force on the obstacle is proportional to the steady drag force at a certain frequency band and the expression for the sound power carried by the fundamental mode travelling in one direction is given by
In the case of a diaphragm, we follow the recent development of Kårekull et al. [24] who revisited the NelsonMorfey semiempirical scaling law by suggesting that the dynamic force is assumed to scale with the momentum flux: where U_{d} is a characteristic velocity taken as the velocity in the vena contracta evaluated from the formula presented by Durrieu et al. [25] and given by U_{d} = U/σ and
The key result in the analysis is to view the source strength spectrum K^{2}(St) as function of the Strouhal number
Here, the characteristic length is chosen to be which corresponds to the hydraulic diameter of the rectangular duct. Experimental results are scaled following the previous model and are given Figure 3 for a diaphragm whose area orifice is A_{o} = 0.005 m^{2} and a vena contracta ratio σ = 0.16 (evaluated using Eq. (5)). A good collapse is found for all velocities over the majority of the frequency range (up to around 3000 Hz). Note that these results are restricted to the contribution of the plane wave mode only and the change in inclination which appears when St ≈ 2.2 correspond to the cutoff frequency of the first transverse mode (around 860 Hz) of the duct. Generalization to the multimodal context could be made following the original work [5] but this is not investigated here. It should be noted that if results are somewhat consistent with experimental investigations reviewed in [4], they are in excellent agreement with those reported in Figure 14 in [5]. In particular below cutoff where it is shown that the curve of 20log(K(St)) exhibits a decay of 20 dB per decade with increasing Strouhal number.
Figure 3 Dimensionless source strength spectra for a single diaphragm (only the contribution of the plane wav mode is considered). 
2.3 Case of two diaphragms in tandem
The prediction of noise due to the presence of multiple elements in air ducts have been the subject of research in the past and we can cite [7, 8]. The theory relies on the interaction factor which is a function of frequency, distance between the various duct elements and geometries. Depending on the crossspectra between the fluctuating forces on both diaphragms, it is anticipated that the radiated acoustic power should increase by 3 dB up to a maximum of 6 dB (assuming that forces on both diaphragms are identical). A series of measurements has been conducted for three separation distances. Results are shown in Figure 4 where the sound power levels are compared with the case of a single diaphragm for a mean flow velocity U = 7 m/s.
Figure 4 Total sound power measured downstream for two diaphragms in tandem with U = 7 m/s. 
Two main observations can be made: first, the broadband noise is increased by 10–20 dB due to impinging turbulence at the downstream orifice. The increase in power decreases with frequency and also when the distance separation is larger. These results are in line with those published in [16] in the case of two diaphragms inserted in a circular duct (see Fig. 14 in [16]) where it is shown that the increase in power becomes marginal when the distance exceeds at least 10 times the hydraulic diameter. Second, a series of consecutive and clearly distinguishable peaks at low frequencies below the first cutoff frequency. Two additional peaks can be identified, the first one appears at the first cutoff frequency of the duct and the second one, around 1500 Hz, corresponds to an acoustic resonance due to the cavity made by the two diaphragms. The modulation of the sound power spectra due to resonance between the orifice plates is also observed and discussed in [16]. The peaks at low frequencies are due to high amplitude acoustic waves and this suggests the presence of flowacoustic feedback mechanisms. This is further investigated in Figure 5 where the effect of the flow velocity is now shown. Again, a series of low frequency peaks can be identified, except when the flow velocity is too low. At this stage, it is instructive to report these resonant frequencies with respect to the mode number n corresponding to the nth identified peak in Figures 4 and 5. This is illustrated in Figure 6. For each configuration, it is observed that peak frequencies behaves like an affine function with respect to mode number. This is typical of the Rossiter’s formula [26, 27] which originates from the study of feedback mechanisms of flows past open cavities. The frequential peaks depends upon the flow velocity and the separation distance as
Figure 5 Total sound power measured downstream for two diaphragms in tandem with L = 13 cm. 
Figure 6 Peak frequencies for all configurations. 
The physical interpretation for this formula is that these tones are the result of a feedback effect produced by the vortices impingement on the second diaphragm causing acoustic waves propagating upstream and triggering a new set of vortices. Here n is the number of vortices located between both diaphragms and γ is an empirical constant. In the original paper [26], the convection velocity of the vortices is usually estimated as a fraction of the external upstream flow velocity U_{∞} and U_{c} = κU_{∞} and in our situation, it is more relevant to consider the flow velocity in the vena contracta so we put U_{∞} = U_{d}. In Table 1 are reported the values of the empirical constants γ and κ which yield a fairly good collapse with the experimental data. Despite the fact that the present configuration notably differs from the more classical cavity flow, the values are somewhat in line with those reported in the scientific literature and formula (7) is representative of the flow mechanisms. We may notice that the empirical constant γ is found to be positive here (as in the experimental data in [28]) though most studies report a negative value [27].
Empirical constants of Rossiter’s model.
3 Numerical results
The aim of this section is to offer a better understanding of the flow oscillations and the radiated acoustic power by means of numerical simulations. It also shows some useful guidelines which permits to capture the physical mechanisms correctly taking place in the duct.
3.1 Numerical setup and parameters
The numerical results are obtained for the separation distance L = 13 cm with an inlet velocity of U = 7 m/s. To prevent acoustic waves to bounce off the boundaries, two Large Eddy Simulations have been tested: the first one, called LES1, utilizes acoustic damping zones at both ends. These zones must stretch over a sufficient distance in order to cover a full wavelength. The lowest frequency of interest is 200 Hz which means that the damping zone is 1.7 m long. In a second setup, named LES2, the acoustic damping is used only in the upstream direction and the nonreflective boundary condition, via a waveTransmissive pressure boundary condition, is prescribed downstream. In both configurations, the same computational domain has been used and the duct length on both sides of the two diaphragms is 3 m long. All surfaces, which includes the duct walls and the obstacles are considered perfectly rigid and a no slip wall condition is applied. The boundary conditions at both end of the duct are: a constant velocity is imposed at the inlet and a constant pressure is imposed at the outlet. Note that planes of symmetry have not been exploited here and the full domain has been simulated.
The computation has been performed with the compressible solver rhoPimpleFoam from OpenFOAM CFD software. In order to assess the convergence, two different meshes are tested. The first mesh G1 consists of a structured hexahedral grid containing 14 million cells. The size of the elements is chosen to ensure that the propagation of acoustic waves is well simulated over a sufficient large distance which must include the position of the virtual microphones, placed on the duct wall as in the experimental setup. A simple rule of thumb shows there are at least 30 cells per wavelength at 1 kHz. Cells are much finer in the vicinity of the shear layer, and are stretched away from the obstacles and in the damping zones. The size of the smallest cells around the edges of the diaphragm where vortices are created is 5 × 10^{−4} m. A finer mesh, G2, is built with the same refinement strategy and contains 33 millions hexahedral cells, the smallest cell size around the obstacle is 3.25 × 10^{−4} m. Figure 7 shows a crosssectional closeup view of both meshes in the plane x_{3} = 0 (the axes are defined in Fig. 8).
Figure 7 Crosssectional view (x_{3} = 0 plane) of the mesh G1 (top) and G2 (bottom), centered around the edges of the two diaphragms. 
Figure 8 Instantaneous velocity and pressure fields in the planes x_{3} = 0 (top) and x_{2} = 0 (bottom). 
LES computations are performed with a given initial state obtained with a RANS k–ω SST (Shear Stress transport) computation. The WALE subgridscale model is used, and second order schemes both in space and time are employed, using respectively linear and backwards schemes. Since the time scheme is implicit, the maximum Courant number is around 5 for both meshes and this corresponds to a time step of 3 × 10^{−5} s and 2.3 × 10^{−5} s respectively. In all results shown here, the simulation time is around 1 s. As a comparison, i.e. the flowthrough time, the amount of time required for the fluid to cross the gap between the obstacles is around 0.05 s.
3.2 Convergence check
Results are postprocessed after a transient period of 0.54 s. This value, which corresponds to 11 flowthrough times, has been determined using the procedure proposed by Mockett et al. [29] to detect the presence of an initial transient content. This procedure is based on the hypothesis that the initial transient of a signal may be reflected through a distortion of its mean or standard deviation. In Figure 9 (top) is shown the time history of the drag force exerted on the second obstacle F_{2}(t). To illustrate the presence of this initial transient, Figure 9 (bottom) shows the evolution of the product of a statistical random error on the mean μ and standard deviation σ of F_{2}, as the start time t_{0} is shifted. Then, the detection of the minimum of this product indicates the end of the initial transient t_{t}.
Figure 9 Time history of the force exerted on the second obstacle (top), and variation of the product with t_{0}, illustrating the detection of t_{t} (bottom). 
3.3 Aerodynamic predictions
Figure 8 shows the instantaneous velocity (left) and pressure (right) fields in the planes x_{3} = 0 and x_{2} = 0. While the geometry and boundary conditions possess symmetry properties, the computed flow is not, as shown in the planar streamlines representation in Figure 10 (left). These streamlines have been obtained from a velocity field averaged on a period of Δt ≈ 0.3 s (approximately 6 flowthrough times). Furthermore, these planar representations do not give insight on the normal velocity component: the 3D streamlines (Fig. 10 (right)) show a strong swirl movement inside the cavity. Isosurfaces of the Qcriterion (Q = 2 × 10^{7}) of the instantaneous flow colored by velocity magnitude is shown in Figure 11. In the side view (Fig. 11 (left)), ring vortices (flagged with black arrows) are clearly identifiable after the first obstacle. Another type of large coherent structure is more visible when the Qcriterion is computed from the mean velocity, see Figure 11 (right) with Q = 2 × 10^{6}. These streamwise coherent structures, flagged with red arrows and clearly seen in Figure 11 (right), seem to break apart the ring vortices as they travel past the first obstacle.
Figure 10 Planar (left) and three dimensional (right) streamlines. 
Figure 11 Qcriterion, top view (left) and axial view (right). 
In order to distinguish the flow features at different frequencies, Figure 12 shows the instantaneous dilatation field, defined as ρ_{0}∂_{t}ρ, filtered at different frequencies. The top picture shows the unfiltered field, and the pictures below show the same field filtered at the peak frequencies identified in Figure 13. The filtered fields exhibit a noticeable pattern as a succession of white and black spots (modal structure), which are present both in the cavity and downstream of the second diaphragm. It needs to be pointed out that (i) the same filtering procedure at another arbitrary frequency would lead to a uniform grey zone in the cavity with no particular pattern and (ii) the fields shown here correspond to LES2 results and using LES1 results would lead to the same observations. From these observations, it emerges that each frequency f_{n} is associated with a number n of vortices in the cavity, and
where λ_{a} = U_{c}/f_{n} is the aerodynamic wavelength, that is the spacing between two vortices, and this is in line with the empirical equation (7).
Figure 12 Dilatation field at an arbitrary time step: unfiltered field (top) and filtered at the mode frequency f_{n}. 
Figure 13 Acoustic power radiated downstream (computed with Eq. (2) where modal amplitudes are extracted using the 2N port formulation): influence of the mesh (left) and influence of the boundary condition (right). 
3.4 Acoustic predictions
In the compressible simulation, virtual microphones have been placed in the computational domain at the same location, i.e. at sections S_{5}, S_{6}, S_{7}, S_{8} as in the experimental setup (see Fig. 1). Amplitudes of the duct acoustic modes are extracted using the 2N port method as described earlier. Contrary to the experimental procedure, radiation conditions are assumed to be exact at both ends of the duct and no acoustic wave travel toward the test section (i.e. ). Thus the acoustic power radiated downstream is straightforward without the need for computing the scattering matrix. Note that this greatly simplifies the analysis and we can refer to a complete computational procedure developed in [15] which includes both the aeroacoustic source strength and the scattering matrix of the induct components.
The sound power spectra obtained with LES1 for both meshes are compared with the experimental one in Figure 13 (left). The grey curves correspond to the numerical results: the solid one has been obtained with the coarsest mesh G1 and the dotted one with the finest mesh G2. Both numerical results are seen to agree, and this implies that the mesh convergence has been reached. Furthermore, the numerical results show similar trend with the experimental ones, which suggests that the feedback phenomenon is indeed captured by the simulation. However, two main differences with the experimental results may be observed: (i) the numerical results exhibit an overall underestimation of the acoustic level and (ii) peak frequencies do not coincide perfectly, especially the second and third peak.
The influence of the downstream boundary condition is shown in Figure 13 (right). While both LES computed results show the same trend, the acoustic power level computed with LES2 is higher by a few dB up to 600 Hz, and peaks seem slightly sharper with a better agreement with experimental data. Even though the exit section is located approximately 3 m downstream, the numerical treatment of the radiation condition can have a discernible effect on the radiated sound power.
In order to highlight the importance of the role played by the acoustic waves on the feedbackloop mechanism, an incompressible computation (LES3) is carried out with a mesh similar to G1 and by using the same schemes and subgridscale model. Here the length of the simulated domain is reduced to remove the damping zones and the incompressible solver pimpleFoam is used. Since acoustic waves can not be captured by the numerical simulation, the approach utilized earlier can not be followed. Instead, the formulation of Nelson and Morfey [5], which allows one to compute the acoustic power from the sole knowledge of the pressure fluctuations on the obstacles can be used. Since all peaks appear below the first cutoff frequency, the expression for the radiated acoustic power can be simplified with the plane wave hypothesis, as in [30]. As two obstacles, assumed acoustically compact, are present here, three different terms appear in the final expression: the first two terms relate to the drag fluctuations on both obstacles and the last one takes into account the crossspectra between the fluctuating forces (see a derivation of the formula in Appendix)
In the present case, most of the noise is generated by the vortices impinging on the second obstacle and the corresponding term is predominant over the whole bandwidth. This is illustrated in Figure 14 (left) where the three components of (9) are plotted separately. The contribution of the first obstacle increases with frequency but remains below the one of the second obstacle. The correlation term is particularly high at the peak frequencies and rapidly decreases with frequency. It is worth observing that the total radiated power, computed with LES2 compressible simulation, is in good agreement with results of Figure 13 (recall that Fig. 14 is obtained with equation (9) while equation (2) is used in Fig. 13)). This means that formula (9) is very reliable, at least for the plane wave mode, and the source term which normally appears in the Lighthill’s acoustic analogy can be discarded in the integral formulation (A2). The reasons for this is that solid walls, which includes the duct and the obstacles, act as perfect sound reflectors and are therefore responsible for most of the radiated power [31].
Figure 14 SWL for the plane wave mode only computed via (9) from the fluctuating forces exerted on the two diaphragms. Contribution of the different terms of (9) computed with LES2 (left). Comparison between LES2 (compressible) and LES3 (incompressible) simulations (right). 
Results obtained with and without compressibility effects (LES2 and LES3 respectively) are now shown in Figure 14 (right). It can be seen that, although broadband noise spectral density are in good agreement, no peaks are identifiable if the flow is assumed incompressible. This point is not as trivial as it may appears since, depending on the configuration, incompressible flow solvers are sometimes sufficient to capture flow oscillation effects, see for instance [32]. Finally, it is worth noting that other numerical setups have also been tried and they all failed to capture the acoustic feedback loop. In particular when symmetries of the geometry are used and only one fourth of the domain is simulated (the fact that symmetry of the flow is broken is confirmed by the simulation in Fig. 10 for instance) and also when limited schemes such as the vanLeer and linearUpwind are used.
4 Conclusions
In this paper, the fluidacoustic interactions caused by flow impingement on two identical diaphragms with rectangular openings inserted in a duct of rectangular crosssection are investigated both experimentally and numerically. The former relies on the the 2N port method, while the latter compares results for the first time from both compressible and incompressible large eddy simulations with the opensource software OpenFOAM. Good overall agreement is achieved between all approaches, providing new insight in the noise mechanisms.
It is first shown that the presence of the second diaphragm can increase the broadband noise by 10–20 dB with the presence of additional peaks at low frequencies. In fact, most of the noise is generated by the vortices impinging on the second obstacle and the corresponding term in the equation of the acoustic power solely obtained from the pressure fluctuations on the diaphragms is predominant over the whole bandwidth.
The location of the peak frequencies are then observed to fit well with Rossiter’s formula which originates from the study of acoustic feedbackloop mechanisms of flows past open cavities. This phenomenon is highlighted using Direct Noise Computations. In particular, the dilatation field filtered at the resonant frequencies shows a distinct pattern corresponding to the Rossiter modes. It is also shown that, although broadband noise spectral density are in good agreement, no resonance peaks are identifiable if the flow is assumed incompressible. Although the Rossiter model has been shown to be relevant in order to give some useful interpretation of the feedback mechanisms, the authors of the present paper are aware that the studied configuration departs noticeably from the more academic case of a flow past a rectangular cavity and there is room for further studies. In particular a more thorough analysis could be conducted and this could be done in the spirit of [33] for instance.
Acknowledgments
The computations were made on the supercomputers Cedar from Simon Fraser University and Niagara from University of Toronto, both parts of Compute Canada’s national platform of Advanced Research Computing resources.
Conflict of interest
The authors declare that they have no conflicts of interest in relation to this article.
Data availability statement
Data are available on request from the authors.
Appendix
Lighthill’s equation for the pressure fluctuations in a turbulent flow takes the form of the inhomogeneous Helmholtz equation
Under the hypothesis of an isentropic lowMach number flow at high Reynolds number, sound sources due to density, viscous stress and entropy changes are neglected and Lighthill’s stress tensor becomes T_{ij} ≈ ρ_{0}u_{i}u_{j}, where scalar quantities u_{i} can be provided either by compressible or incompressible flow simulations. Note that (A1) is written in the frequency domain with the convention e^{−iωt}. Now, using the Taylored Green’s function G_{d} for the rectangular duct [30], the acoustic pressure radiated away has the explicit closed form expression:
where p(y, ω) stands for the pressure fluctuations on the rigid obstacle(s) Γ and n_{y} is the normal unit vector to the surface. Because we are interested in the low frequency propagation, i.e. below the first cutoff duct acoustic mode, the first term of the modal series is predominant [30] and
The theory of Nelson and Morfey [5] relies on the assumption that the volume integral in (A2) is negligible compared to the surface integral. This yields (we use the same notation as in (4)):
where
is the cross power spectral density for the pressure fluctuations. The formula can be simplified in the presence of two acoustically compact obstacles Γ_{j}(j = 1,2) separated by distance L and by calling
the drag force on each obstacle, we end up with formula (9) or (4) in the presence of a single obstacle.
References
 S. Guérin, E. Thomy, M.C.M. Wright: Aeroacoustics of automotive vents. Journal of Sound and Vibration 285, 4 (2005) 859–875. [CrossRef] [Google Scholar]
 D. Kurniawan, E. Rogers: Investigation of airflow induced whistle noise by HVAC control doors utilizing a “vshape” rubber seal. SAE Technical Paper 2011011615, 2011. [Google Scholar]
 D. Oldham, A. Ukpoho: A pressurebased technique for predicting regenerated noise levels in ventilation systems. Journal of Sound and Vibration 140, 2 (1990) 259–272. [CrossRef] [Google Scholar]
 O. Kårekull, G. Efraimsson, M. Åbom: Prediction model of flow duct constriction noise. Applied Acoustics 82 (2014) 45–52. [CrossRef] [Google Scholar]
 P. Nelson, C. Morfey: Aerodynamic sound production in low speed flow ducts. Journal of Sound and Vibration 79, 2 (1981) 263–289. [CrossRef] [Google Scholar]
 D.J. Oldham, D.C. Waddington: The prediction of airflowgenerated noise in ducts from considerations of similarity, Journal of Sound and Vibration 248 (2001) 780–787. [CrossRef] [Google Scholar]
 N. Han, C.M. Mak, Prediction of flowgenerated noise produced by acoustic and aerodynamic interactions of multiple induct elements. Applied Acoustics 69, 6 (2008) 566–573. [CrossRef] [Google Scholar]
 C. Cai, C.M. Mak: Generalized flowgenerated noise prediction method for multiple elements in air ducts. Applied Acoustics 135 (2018) 136–141. [CrossRef] [Google Scholar]
 J. Kreuzinger, F. Schwertfirm, M. Hartmann, N. Peller: Analysis of resonance phenomena caused by obstacles in HVAC exhaust nozzles using a combined CFDCAA approach, in: 19th AIAA/CEAS Aeroacoustics Conference, Berlin, Germany, May 27–29, 2013. [Google Scholar]
 R. Parker, S.A.T. Stoneman: The excitation and consequences of acoustic resonances in enclosed fluid flow around solid bodies. Proceedings of the Institution of Mechanical Engineers, Part C: Mechanical Engineering Science 203, 1 (1989) 9–19. [CrossRef] [Google Scholar]
 S. Ziada, P. Lafon: Flowexcited acoustic resonance excitation mechanism, design guidelines, and counter measures. Applied Mechanics Reviews 66, 1 (2013) 010802 [Google Scholar]
 C. Sovardi, Y. Aurégan, W. Polifke: Parametric LES/SI based aeroacoustic characterization of tandem orifices in low Mach number flows, Acta Acustica united with Acustica 102, 5 (2016) 793–803. [CrossRef] [Google Scholar]
 U. Karban, C. Schram: Modal identification of aeroacoustic systems using passive and active approaches. Journal of the Acoustical Society of America 142, 6 (2017) 3804–3812. [CrossRef] [PubMed] [Google Scholar]
 U. Karban, G. Ogus, K. Kucukcoskun, C.F. Schram, C. Sovardi, W. Polifke: Noise produced by a tandem diaphragm: Experimental and numerical investigations. in: 20th AIAA/CEAS Aeroacoustics Conference, Atlanta, GA, 16–20 June, 2014. [Google Scholar]
 S. Sack, M. Shur, M. Åbom, M. Strelets, A. Travin: Numerical education of active multiport data for induct obstructions. Journal of Sound and Vibration 411 (2017) 328–345. [CrossRef] [Google Scholar]
 S. Sack, M. Åbom: Investigation of orifice aeroacoustics by means of multiport methods. Journal of Sound and Vibration 407 (2017) 32–45. [CrossRef] [Google Scholar]
 F. Li, P. Wang, Y. Liu: Unsteady flow behaviors and noise generation mechanisms of tandem orifices in a circular duct. Physics of Fluids 35, 1 (2023) 015142. [CrossRef] [Google Scholar]
 P. MartinezLera, K. Kucukcoskun, M. Shur, A. Travin, M. Tournour: Hybrid aeroacoustic computations for flows in ducts with single and tandem diaphragms, in: 22nd AIAA/CEAS Aeroacoustics Conference Lyon, France, 30 May–1 June, 2016. [Google Scholar]
 A. Sengissen, B. Caruelle, P. Souchotte, E. Jondeau, T. Poinsot: LES of noise induced by flow through a double diaphragm system, in: 15th AIAA/CEAS Aeroacoustics Conference (30th AIAA Aeroacoustics Conference), Miami, Florida, 11 May 2009–13 May 2009. [Google Scholar]
 F. Mendonca, A. Read, S. Caro, K. Debatin, B. Caruelle: Aeroacoustic simulation of double diaphragm orifices in an aircraft climate control system, in: 11th AIAA/CEAS Aeroacoustics Conference, Monterey, California, 23 May 2005 – 25 May 2005. [Google Scholar]
 H. Trabelsi, N. Zerbib, J.M. Ville, F. Foucart: Passive and active acoustic properties of a diaphragm at low Mach number: Experimental procedure and numerical simulation. European Journal of Computational Mechanics 20 (2011) 49–71. [CrossRef] [Google Scholar]
 J. Lavrentjev, M. Åbom: Characterization of fluid machines as acoustic multiport sources. Journal of Sound and Vibration 197, 1 (1996) 1–16. [CrossRef] [Google Scholar]
 F. Hugues: Modelling the vibrations generated by turbulent flows in ducts. PhD Thesis, Université de Technologie de Compiène, 2018. Available at https://theses.hal.science/tel02094391/file/These_Florian_Hugues.pdf. [Google Scholar]
 O. Kårekull, G. Efraimsson, M. Åbom: Revisiting the Nelson–Morfey scaling law for flow noise from duct constrictions. Journal of Sound and Vibration 357 (2015) 233–244. [CrossRef] [Google Scholar]
 P. Durrieu, G. Hofmans, G. Ajello, R. Boot, Y. Aurégan, A. Hirschberg, M.C.A.M. Peters: Quasisteady aeroacoustic response of orifices. The Journal of the Acoustical Society of America 110 (2001) 1859–1872. [CrossRef] [PubMed] [Google Scholar]
 J. Rossiter: Windtunnel experiments on the flow over rectangular cavities at subsonic and transonic speeds, volume 3438 of Aeronautical research council report and memorandum, Ministry of Aviation, United Kingdom, 1964. [Google Scholar]
 X. Gloerfelt, C. Bogey, C. Bailly: Cavity noise [lecture notes], Arts et Métiers ParisTech, 2007, pp. 262–272. [Google Scholar]
 V. Sarohia: Experimental investigation of oscillations in flows over shallow cavities, AIAA Journal 15, 7 (1977) 984–991. [CrossRef] [Google Scholar]
 C. Mockett, T. Knacke, F. Thiele: Detection of initial transient and estimation of statistical error in timeresolved turbulent flow data, in: 8th International Symposium on Engineering Turbulence Modelling and Measurements (ETMM8), Marseille, France, 9–11 June 2010. [Google Scholar]
 N. Papaxanthos, E. PerreyDebain, S. Bennouna, B. Ouedraogo, S. Moreau, J. Ville: Pressurebased integral formulations of Lighthill–Curle’s analogy for internal aeroacoustics at low mach numbers. Journal of Sound and Vibration 393 (2017) 176–186. [CrossRef] [Google Scholar]
 N. Papaxanthos, E. PerreyDebain: Integral formulations for the prediction of low Mach number flow noise with noncompact solid surfaces, in: 22nd AIAA/CEAS Aeroacoustics Conference, Lyon, France, 30 May–1 June 2016. [Google Scholar]
 V. Suponitsky, E. Avital, M. Gaster: On threedimensionality and control of incompressible cavity flow, Physics of Fluids 17, 10 (2005) 104103. [CrossRef] [Google Scholar]
 C.W. Rowley, T. Colonius, A.J. Basu: On selfsustained oscillations in twodimensional compressible flow over rectangular cavities. Journal of Fluid Mechanics 455 (2002) 315–346. [CrossRef] [Google Scholar]
Cite this article as: de Reboul S. PerreyDebain E. Ville JM. Zerbib N. Hugues F. & Moreau S. 2024. Experimental and numerical observation of flowacoustics feedback phenomena due to two diaphragms in tandem inserted in a rectangular duct. Acta Acustica, 8, 7.
All Tables
All Figures
Figure 1 Sketch of the AeroAcoustics test rig. 

In the text 
Figure 2 Two diaphragms in tandem and dimensions of a single diaphragm. 

In the text 
Figure 3 Dimensionless source strength spectra for a single diaphragm (only the contribution of the plane wav mode is considered). 

In the text 
Figure 4 Total sound power measured downstream for two diaphragms in tandem with U = 7 m/s. 

In the text 
Figure 5 Total sound power measured downstream for two diaphragms in tandem with L = 13 cm. 

In the text 
Figure 6 Peak frequencies for all configurations. 

In the text 
Figure 7 Crosssectional view (x_{3} = 0 plane) of the mesh G1 (top) and G2 (bottom), centered around the edges of the two diaphragms. 

In the text 
Figure 8 Instantaneous velocity and pressure fields in the planes x_{3} = 0 (top) and x_{2} = 0 (bottom). 

In the text 
Figure 9 Time history of the force exerted on the second obstacle (top), and variation of the product with t_{0}, illustrating the detection of t_{t} (bottom). 

In the text 
Figure 10 Planar (left) and three dimensional (right) streamlines. 

In the text 
Figure 11 Qcriterion, top view (left) and axial view (right). 

In the text 
Figure 12 Dilatation field at an arbitrary time step: unfiltered field (top) and filtered at the mode frequency f_{n}. 

In the text 
Figure 13 Acoustic power radiated downstream (computed with Eq. (2) where modal amplitudes are extracted using the 2N port formulation): influence of the mesh (left) and influence of the boundary condition (right). 

In the text 
Figure 14 SWL for the plane wave mode only computed via (9) from the fluctuating forces exerted on the two diaphragms. Contribution of the different terms of (9) computed with LES2 (left). Comparison between LES2 (compressible) and LES3 (incompressible) simulations (right). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.