Theory
The run stores the scattering matrix of the muffler and of its reference pipe: for every frequency, the waves leaving every port for a unit wave arriving at every port, for all modes. A source is described by the wave it sends, , and by the reflection coefficient of its plane mode for the wave coming back. The wave arriving at the inlet is then , where is the wave leaving the inlet, and the closed loop of source, model and outlet termination is solved exactly from the stored . The muffler and the reference pipe get the same source, so their radiated powers form the insertion loss .
An anechoic source absorbs everything that comes back. A constant-velocity source (a rigid piston) reflects it fully in phase, a constant-pressure source in opposite phase. A table gives the measured normalised source impedance , interpolated linearly between its rows and held at its end values outside them. With a reflecting source the reference pipe's power depends on its length (Reference pipe length).
Equations
| where | meaning | unit |
|---|---|---|
| Plane-mode reflection coefficient of the source | - | |
| Source impedance normalised by of the inlet gas | - |
| where | meaning | unit |
|---|---|---|
| Waves leaving the ports | Pa | |
| Stored scattering matrix (all ports, all modes) | - | |
| Reflection matrices of the source (inlets) and of the outlet terminations | - | |
| Waves sent by the source(s) | Pa |
Default and rule
- Default
- Anechoic (matched)
- Rule
Anechoic.
- Source
The algebra of the program on its finite-element scattering matrices. It was checked against the plane-wave four-pole insertion loss of Munjal (1987) (section 2.10.1, eqs. 2.58-2.65, pp. 55-56, with the uniform-tube matrix eq. 2.134, p. 77) on an expansion chamber for every source and radiation choice; that definition refers to the source connected straight to the radiation impedance, so the comparison subtracts the same expression for the reference pipe.
Valid range anechoic | constant velocity | constant pressure | table of z_s(f)
No solve
Changing the source re-evaluates the stored matrices in well under a second.
Table format
Three numbers per line: frequency in Hz, real and imaginary part of the normalised source impedance.