## 2 Materials and Methods  
Two junctions, shown schematically in Figure 1, were manufactured. A T-junction and a Y-junction were considered, in order to allow the analysis of the effect of the angle of the side branch. An internal diameter of 51 mm was used in all the branches of the junctions.  
While several formalisms may be used for the representation of the transient response of a system, the most intuitive one for the present case of a junction is that based on the consideration of wave components, so that the junction is actually regarded as a multi-port. In this framework, for a junction such as that represented in Figure 2, one has three excitations and three responses, and writing the relations between them directly in matrix form in the frequency domain, one has:  
\[\left[\begin{array}{c}B_{1}\\ B_{2}\\ B_{3}\end{array}\right]=\left[\begin{array}{ccc}R_{1}&T_{21}&T_{31}\\ T_{12}&R_{2}&T_{32}\\ T_{13}&T_{23}&R_{3}\end{array}\right]\left[\begin{array}{c}A_{1}\\ A_{2}\\ A_{3}\end{array}\right]. \tag{1}\]  
where, as indicated in Figure 2, \(A_{i}\) represents the wave component moving towards the junction in port \(i\) and \(B_{j}\) represents the wave component moving away from the junction in port \(j\). Regarding  
Figure 1: Junctions considered in the study.  
the matrix elements, \(R_{i}\) denotes the reflection coefficient as seen from port \(i\) whereas \(T_{ij}\) denotes the transmission coefficient between ports \(i\) and \(j\). All these magnitudes are functions of the frequency \(f\).  
In this way, one has a reflection coefficient for each of the pipes arriving at the junction, and transmission coefficients for all the possible transmission paths, indicated by the corresponding subscripts. The experimental setup and the corresponding measurement procedure are described in detail in [32; 33], and here only a brief overview is given in Appendix A.  
Two different modeling approaches were evaluated: a staggered mesh finite volume method and, as a reference, a more conventional pressure loss-based model. The staggered mesh finite volume method used is described in detail in [34], where a flux-corrected transport (FCT) technique was used in order to suppress the spurious oscillations that these numerical methods exhibit in the vicinity of discontinuities in the flow variables. It was found that satisfactory results were obtained through the combination of dissipation via damping together with the phoenical form of the anti-diffusion term. As an alternative, the momentum diffusion term described in [35] was also used. A brief summary is given in Appendix B.  
The pressure loss-based model uses a conventional one-dimensional finite volume model with a collocated mesh, derived from the code available in [36]. The junction is modeled as a small volume with three connections to which different pressure loss coefficients are assigned, and in which the mass and energy conservation equations are solved. For the connection to the duct where the incident pressure pulse propagates it has been assumed that the total pressure loss is zero, whereas for the other two connections their corresponding pressure loss coefficients are computed following the expressions given in [15; 16]. The procedure is summarized in Appendix C.