Abstract
Boundary element solutions of the Helmholtz integral equation for the sound attenuation of acoustically lined curved pipes are presented. The curved pipes are in the form of cylindrical and toroidal coordinates. The main objective of the study is a comparative analysis of sound attenuation characteristics of lined toroidal pipes. Although toroidal pipes with circular cross-section are used extensively in industrial and agricultural applications, they have been taken little attention in literature. One reason may be the impossibility of theoretical solution since wave equation is not separable in toroidal coordinates. A numerical method should be utilized in the solution; in fact, not many researchers are familiar with toroidal coordinates on the contrary of other coordinate systems, such as cylindrical coordinates. The use of toroidal coordinates and generation of a pipe in these coordinates are illustrated herein. Sound waves in the pipes are produced by the vibrating inlet surface and resulting sound pressure distribution on the interior surfaces of side walls and anechoic exit are computed. Sound transmission loss of the pipes is represented by the ratio of inlet to exit pressures, as sound attenuation spectra. The results of cylindrical and toroidal pipes with the same cross-sectional area are compared. The sound attenuation performance of lined toroidal pipes is compared in terms of bend sharpness as a measure of curvature. The comparison is extended to cover the hard-walled and lined toroidal pipes to conclude their common and/or different acoustic behaviours.
Keywords
1. Introduction
The sound propagation in pipes and ducts is a prominent subject in both theoretical and applied acoustics. Straight pipe problems have become a classical topic; however, curved pipes have been relatively less endeavoured due to their complex geometries (Cummings, 1974; Rostafinski, 1974; Cabelli, 1980; Keefe and Benade, 1983; Firth and Fahy, 1984; Furnell and Bies, 1989; Campos and Serrão, 2003; Dequand et al., 2003; Kim et al., 2003; Bakkali et al., 2014; Beauvais et al., 2021). A basic problem in applied acoustics is to attenuate the sound output of pipes. Although a curved shape causes attenuation at certain level (Sarıgül, 1999), the most efficient application incontestably is to line the pipe with an acoustic material. The combination of acoustic lining with curvature magnifies the sound attenuation and complexity of the problem. In literature, there are a few distinguished studies on this field (Grigor’yan, 1970; Myers and Mungur, 1976; Ko and Ho, 1977; Ko, 1979; Rostafinski, 1982; Felix and Pagneux, 2004; Yang, 2019).
The wave equation is separable only in very few coordinates, and toroidal coordinates are not one of them. Since the theoretical solution for toroidal pipes has not been accomplished, a reliable numerical technique is required. It has been shown that the boundary element method (BEM) that bases on the interior surface Helmholtz integral equation is a perfect numerical tool to perform the solution in these coordinates (Sarıgül, 1999). The solution for hard-walled pipes is relatively straightforward. However, the modelling of lined pipes includes a lot of questions and steps because they are explicit to many alternatives. Therefore, a validation study is needed to determine the most suitable BEM modelling that exposes closest results to the analytical solutions. This pre-study has been carried out for the 45°-cylindrical pipe with rectangular cross-section and the best BEM configurations have been determined (Sarıgül, 2019).
This study is mainly devoted to lined toroidal pipes. Cylindrical pipe is solved only for the comparison of its acoustic behaviour with that of the toroidal pipe. Lined toroidal pipes are examined extensively to put forward their curvature and sound attenuation relation. The paper finalizes with a comparative analysis between sound attenuation characteristics of lined and hard-walled toroidal pipes.
2. Formulation of the problem
If the domain of interest is the interior acoustic field of a body with a boundary
Before the numerical solution of the Helmholtz integral in equation (1), application of the boundary conditions is necessary. In general, if the normal velocity amplitude of the boundary is denoted by
Helmholtz integral equation (1) relates the sound pressure amplitude 1. The surface integral of the general body in equation (1) is represented by the numerical summation of the surface integrals of the discretized boundary elements. Boundary elements are quadrilateral and quadratic and have eight nodes. 2. Coordinates, pressures and normal velocities of all points on the elements are represented in terms of nodal values by using the same shape functions in accordance with the isoparametric element formulation. 3. Surface integrals of the boundary elements including singularity are handled by using the Gaussian quadrature technique. Sixteen collocation points are utilized for each element. 4. At the end of this procedure, a set of complex algebraic equations is obtained
Here,
The numerical solution is performed by an ‘in-house’ computer code which has the capability of generating equation (5) for arbitrarily shaped three-dimensional bodies possessing any type of boundary conditions. This code has been used to solve different interior and exterior problems (Dokumacı and Sarıgül, 1995; Sarıgül, 1999, 2019; Sarıgül and Kıral, 1999; Sarıgül and Seçgin, 2004; Seçgin and Sarıgül, 2010; Sarıgül et al., 2014).
3. Toroidal coordinates
Toroidal coordinates are three-dimensional orthogonal coordinate system that results from rotating two-dimensional bipolar coordinate system about the vertical axis that separates its two foci (Spiegel, 1968; Arfken, 1970). Figure 1 shows half of a torus with foci F1 and F2 located at distances in the Representation of toroidal coordinates on a half torus. Distances ρ1 and ρ2 from a point on the torus to two foci.


The transformation from toroidal to Cartesian coordinates is given as (Spiegel, 1968; Arfken, 1970)
Some notes on toroidal coordinates may be gathered as follows:
1. Coordinate v is specific to a given torus and the same for all points on this torus.
2. Coordinate
3. Different boundary points on any circular section have different u coordinates.
4. For any torus, focal distance e and v coordinate may be calculated by equating the transformations from polar and toroidal coordinates to Cartesian coordinates.
-Writing x coordinate of point A on the right section given in Figure 3 ( Right circular section of the half torus.

-Writing
Solving equations (9) and (10) simultaneously, focal distance
Taking bend sharpness of the torus
5.
For
6. Toroidal coordinates of other boundary points located on the subsequent sections are given as follows:
-
-
-
4. Sound attenuation in pipes
No fluid flow is present in the pipes that have an entrance section vibrating with constant normal velocity amplitude like a piston to produce planar sound waves. The exit section behaves as if an anechoic surface representing no reflection condition. This section possesses an acoustic admittance of Geometry of 45°-cylindrical pipe with rectangular cross-section.

Since the cross-sectional area of the pipes is constant, sound reduction characteristics are represented by Sound Attenuation Level, that is measured by the ratio of the inlet and exit pressures
Sound pressure is constant throughout the inlet surface; however, nodal pressures of the exit have quite different values. The average value is used for the exit pressure as performed by Sarıgül (1999, 2019).
In the present study, pipes are handled in three headings: (1) 90°-cylindrical pipe with rectangular cross-section. (2) 90°-toroidal pipe with circular cross-section with a definite curvature. (3) 90°-toroidal pipes with circular cross-section and varying curvatures.
For all pipes, the boundary element modelling features proposed for lined cylindrical pipes by Sarıgül (2019) are used. Basically, this paper focuses on the third heading. The numerical solution of toroidal pipes is essential since they have no analytical solution on the contrary of cylindrical coordinates. In the study, 90° pipes with a shape of quarter torus are examined due to their common use in industry as right-angle bends. The paper basically aims to obtain the relation between the bend sharpness and the sound attenuation characteristics of lined toroidal pipes. This relation has been put forward for the same hard-walled pipes by Sarıgül (1999). The question to be answered here is that ‘the trend is similar or some new characteristics are gained with lining?’
4.1. 90°-cylindrical pipe with rectangular cross-section
A 90°-cylindrical pipe with rectangular cross-section is shown in Figure 5. The discretization is the same as the model with 384 elements and 1186 nodes applied for the 45°-pipe by Sarıgül (2019). In this model, side walls have 256 elements and 800 nodes; each of the inlet and outlet surfaces has 64 elements and 193 nodes. 90°-pipe model is constructed by extending the side elements of the 45° pipe. Also, in the analyses, the same numerical parameters are used ( 90°-cylindrical pipe with rectangular cross-section.
Using the same configuration of the selected model by Sarıgül (2019) (A3), the sound attenuation spectrum of this pipe is obtained and presented in Figure 6 with the spectra of the 45° pipe. Spectrum of the 90° pipe continues smoothly and effectively up to almost Sound attenuation spectra of cylindrical pipes with rectangular cross-section (- 45° analytical, - - 45° BEM and 
4.2. 90°-toroidal pipe with circular cross-section with a definite curvature
A 90°-toroidal pipe with circular cross-section is seen in Figure 7 and its geometry is given in Figure 8. Here, 90°-toroidal pipe with circular cross-section. Geometry of 90°-toroidal pipe with circular cross-section.

The pipe is modelled by using 256 elements and 800 nodes on the side walls again and 60 elements and 205 nodes on each of the inlet and outlet surfaces. Therefore, 376 elements and 1210 nodes are used in the numerical implementation. The discretization of inlet–outlet surfaces is shown in Figure 9. The attenuation spectrum of this pipe is presented in Figure 10 comparatively with that of the 90°-cylindrical pipe with rectangular cross-section. No apparent difference is observed between the two spectra in trend, in the position and level of the maximum attenuation. However, the cylindrical pipe exposes smoother variation throughout the entire spectrum. Although the discretization sensitivity of both types is quite similar, the spectrum of the toroidal pipe has some imperfections. These fluctuations may be considered as small for lower frequencies, whereas large for higher frequencies. Therefore, this spectrum is extended up to Discretization of the inlet–outlet surfaces of 90°-toroidal pipe with circular cross-section (60 elements and Sound attenuation spectra of 90°-pipes (

4.3. 90°-toroidal pipes with circular cross-section and different curvatures
These pipes have the same general shape as that of the second heading. In this section, the curvature Final discretization of the inlet–outlet surfaces of 90°-toroidal pipe with circular cross-section (128 elements and 
In the generation of these geometries, two parameters are fixed and equal to those of the 90°-cylindrical pipe and toroidal pipe with 1. Median arc radius 2. Relative depth of the lining (lining core depth/pipe radius)
The variation of pipe parameters with bend sharpness.
The sound attenuation spectra with respect to non-dimensional frequency Comparative sound attenuation spectra of toroidal pipes with different bend sharpness The variation of sound attenuation with bend sharpness for lined toroidal pipes.
The comparison among the spectra in Figure 12 puts forward the following features of lined toroidal pipes: 1. As the bend sharpness increases, the magnitude of the peak attenuation decreases. 2. This decrease is rapid for pipes with smaller bend sharpness and slow for larger ones. 3. The maximum attenuation occurs at higher non-dimensional frequencies for pipes with smaller bend sharpness. 4. As the bend sharpness increases, these frequencies shift to smaller values, rapidly at first, slowly afterwards. 5. Sound attenuation spectra of pipes with small bend sharpness have a smooth distribution. The attenuation capacity of these pipes also has a large frequency range from zero to ten. 6. As the bend sharpness increases, the upper limit of attenuation range decreases towards smaller frequencies.
Although Table 2 shows the aforementioned trends apparently, the peak attenuation values computed for pipes with high bend sharpness somewhat violate the decrease trend. This is due to the computational inaccuracy regarding to the increase in boundary element size.
Specifically, since in the present study, increase in bend sharpness
Table 3 includes the numerical values of the sound attenuation spectra for the hard-walled pipes presented by Sarıgül (1999). In that study, the variation in bend sharpness has been provided again by changing 1. Addition of lining greatly increases the attenuation capability of pipes both in level and range. 2. In the hard-walled pipes, as the bend sharpness increases, the peak attenuation level increases on the contrary of lined pipes. This means pipes with wider radii have better sound attenuation. This is a conflicting characteristic between lined and hard-walled pipes. 3. In the hard-walled pipes, the position of peak attenuation shifts to higher kr values, as the bend sharpness increases. For the lined pipes, there isn’t a uniform trend related to kr; however, there is an obvious opposite trend in terms of 4. For the hard-walled pipes, attenuation range increases and shifts towards the higher kr values, as the bend sharpness increases. There is no specific trend for lined pipes in terms of kr; however, an opposite trend is present with regard to 5. It is seen that the relation between pipe dimension and frequency is more sensitively represented by the non-dimensional frequency The variation of sound attenuation with bend sharpness for hard-walled toroidal pipes in Sarıgül (1999).
5. Conclusion
The present study had the objective of examining sound attenuation characteristics of lined curved pipes, specifically 90°-toroidal pipes. The study put forward basic acoustical characteristics of lined toroidal pipes, the effects of dimension on sound attenuation and the difference of lined and hard-walled pipes in terms of sound attenuation behaviour. A brief comparative analysis on the effect of cross-sectional shape was also comprised. In order to expose these characteristics comparatively, solutions of a 90°-cylindrical pipe and many 90°-toroidal pipes with different dimensions leading to following inferences were accomplished. 1. The sound attenuation performance of toroidal pipes is almost the same as that of the cylindrical pipes with the same cross-sectional area. That is, there is almost no effect of the cross-sectional shape on the sound attenuation. 2. As the pipes get narrower, the sound attenuation performance of lined toroidal pipes increases both in level and range. Lining of a narrow pipe with an acoustical material would yield a better sound attenuation performance. 3. In fact, the attenuation in widest pipes lined with fibrous sheet mounted on a locally reacting core with impervious backing is almost 30 dB and may not be considered as low. However, the frequency range of this performance is quite limited and attention should be paid in practical applications. 4. Different lining materials may be applied to examine the effect of the impedance model. This application may alter the level of attenuation; however, the general behaviour of toroidal pipes would not be much affected.
