_{1}

^{*}

The purpose of this letter is to present the corrected expressions for the four-pole transmission matrix for a duct with a linear temperature gradient and an exponential temperature profile, described in Sujith [1]. The corrected equations are used in the analyses of a duct that is driven by a piston at one end and a rigid termination at the other end and the gas has a linear and exponential temperature gradients. The acoustic pressure and particle velocity along the duct are calculated and the theoretical results are compared with predictions using finite element analysis.

The four-pole, or transmission line method is a useful theoretical tool for the acoustic analysis of duct systems with plane waves [2,3].

Sujith [

Consider a duct filled with a gas that has linear temperature profile given by

(1)

where is the gradient of the temperature profile, and is the temperature at. The speed of sound and density of the gas vary with temperature as [

(2)

(3)

The pressure and acoustic particle velocities at the ends of the duct are related by the four-pole transmission matrix as

Where p_{i}_{ }and u_{i} are the acoustic pressure and acoustic particle velocity at the ends of the duct, respectively. The equations presented by Sujith [

The elements of the four-pole transmission matrix are

where the specific gas constant is and the constant is

Equations (7) and (8) presented here are the corrected versions of Equations (14) and (15) in Sujith [

Note that if one were to define a constant temperature profile in the duct, such that and, the terms and in Equations (7) and (8) equate to, and Equation (10) equates to zero which causes numerical difficulties. Instead, an approximation to a constant temperature profile can be made by specifying a small temperature difference between the ends of the duct, say 0.1˚C. It can be shown numerically that this will approximate the transmission matrix for a duct with a constant temperature profile and the same cross-sectional areas at the inlet and outlet of the duct given by [

Sujith [

The equation for the acoustic pressure along the duct is given by

The equation for the acoustic particle velocity is given by

Equation (15) presented here is the corrected version of Equation (20) in Sujith [

where

Equations (17) to (20) presented here are the corrected versions of Equations (21) to (24) in Sujith [

A finite element model of a duct with a piston excitation at one end and rigid termination at the other end was created using the finite element analysis software Ansys, release 14.5. A capability introduced in release 14.5 is the ability to define acoustic elements that have temperatures defined at nodes. In previous releases of the software, the speed of sound and density of the gas had to be defined using a set of material properties. A model of duct with temperature variations could only be created using small duct segments with constant material definitions in each segment. Hence, an impedance discontinuity would have been created at the interface between two duct segments with dissimilar material properties, and would have caused acoustic reflections due to the impedance mismatch at the interface. The new capability enables one to define a temperature gradient across an element so that there is no impedance discontinuity.

The process for conducting an acoustic finite element analysis where there are variations in the temperature of the gas involves several steps as follows:

1) A solid model is created that defines the geometry of the system.

2) The solid model is meshed with thermal elements (SOLID70).

3) The temperature boundary conditions are applied for each region.

4) A static thermal analysis is conducted to calculate the temperature distribution throughout the duct network.

5) The temperatures at each node are stored in an array.

6) The thermal elements are replaced with acoustic elements (FLUID30).

7) The values of temperatures stored in the array are used to define the temperature at each node.

8) Anechoic boundary conditions are set at the duct inlet and outlet (using the MAPDL command SF, IMPD, INF).

9) The acoustic velocity at the duct inlet is defined (using the MAPDL command SF, SHLD, velocity).

10) A harmonic analysis is conducted over the analysis frequency range.

11) The sound pressure levels at the duct inlet, outlet, and the entrance and closed end of the QWT are calculated.