# Automotive Muffler Simulation

## Introduction

As a car engine runs, a significant amount of noise is generated due to the combustion. An automotive muffler is an essential device that attenuates the noise exhaustion from the engine.

This study is to simulate the sound propagation within three different kinds of mufflers, and compare their performance varying from the frequency to . The mufflers analyzed here are a classical muffler, a perforated muffler and a padded muffler. The sound pressure distribution and the transmission loss spectrum will be calculated to show the damping behavior of these types of mufflers.

The symbols and corresponding units used throughout this tutorial are summarized in the Nomenclature section.

Please refer to the information provided in "Acoustics in the Frequency Domain" for more general theoretical background for acoustics.

Load the finite element package and set the \$HistoryLength to 0.

## Frequency Domain Acoustic Model

To study the performance of a muffler over a range of frequencies, an acoustic model is built in the frequency domain based on the source-free Helmholtz partial differential equation (PDE): Define the variables and select the parameters for a frequency domain acoustics model.
Air is used as the sound medium.
Define material parameters for the acoustics PDE.

## Domain

The classical acoustic muffler consists of three parts: an inlet pipe, an outlet pipe and expansion chambers. Note that the inlet and outlet pipes are assumed to have infinite extension, and are modeled with a radiation boundary condition and an absorbing boundary condition, respectively.

There are two common ways to further improve the muffler performance. The first method is to use perforated pipes within the muffler, and the second method is to cover the walls of the muffler with sound absorptive pads. Each technique has its advantages and disadvantages, and will be explained in the later section.

Specify the parameters of the geometry and the location of inlet/outlet.

## Boundary Conditions

There are three types of the boundary conditions involved in this example. To model the incoming sound wave at the inlet a radiation boundary condition is used.

This radiation boundary condition is very general. However, in the 3D case it may use up more memory that is available. To reduce the memory requirement, we observe that the inlet is at . Therefor the boundary normal is in the direction . This boundary unit normal can be specified as a parameter and will simplify the boundary condition.

Specify the incoming sound amplitude , and set a radiation boundary condition at the inlet.

To model the outgoing wave at the outlet an absorbing boundary condition is specified.

Set an absorbing boundary condition at the outlet.

On the remaining walls of the muffler the default sound hard boundary conditions are used. Since the sound hard boundary condition is the default boundary condition nothing needs to be further specified.

## Model 1: Classical Muffler

As a benchmark we first consider a classical muffler without any perforation or sound absorbing pads.

Import an externally generated boundary mesh of the acoustic muffler.

In acoustics simulations the wavelength of a sound wave needs to be resolved by a sufficiently fine mesh in order to get an accurate numerical solution. Here we set the max edge length to 12 nodes per , which means that there will be at least twelve elements per wavelength in each direction of the wave propagation.

Set the mesh element size according to the maximum frequency of interest.
Generate the full mesh of the domain.
Setup the PDE.

To study the behavior of the muffler over the frequency range , the PDE model is solved repetitively with ParametricNDSolveValue.

Solve the PDE varying from the frequency to with an increment of .

## Model 2: Muffler with Perforation

Next, we consider a muffler with perforation. The perforated pipes allow the sound pressure wave to scatter out in various directions within the expansion chamber.

Import a mesh of the perforated acoustic muffler.
Generate the full mesh of the domain.

In this scenario the PDE is the same as in the previous case.

Setup the PDE.

Note that ruining the following simulation will require hardware that has more than 8GB memory.

Check the available amount of memory.
Solve the PDE varying from the frequency to with an increment of .

## Model 3: Muffler with Sound absorptive Pads

Last, we consider a muffler with sound absorption materials padded on the inner walls. Common sound absorptive materials are made of fiberglass or steel wool. Depending on the density and the mean radius of the material , the flow resistivity values typically lie in range of to .

In the following model a fiberglass pad with the flow resistivity and the thickness is used.

Define the flow resistivity and the thickness of the sound absorptive pad.

A technique to model the sound propagation within absorptive materials is called "Equivalent Fluid Method" . The idea is to characterize the sound attenuation by inserting the effective density and the effective speed of sound in the original Helmholtz equation (3). and are complex valued and can be calculated with an empirical formula derived by Delany and Bazley . For the low-frequency range where the non-dimensional parameter , Kirby and Cummings  presented the following more accurate approximations:   Here and are the complex wave number and complex impedance, respectively.

Calculate the parameter based on the frequency range of interest: .

Since the parameter over the entire spectrum, the equation (6) should be used to calculate the complex parameters and .

Calculate the complex valued speed of sound and the complex valued density .
Define material parameter rules for the padded muffler.
Insert the material parameters into the model.
Solve the PDE varying from the frequency to with an increment of .

## Post-Processing and Visualization

### Sound Pressure Distribution

To see the general damping behavior of acoustic waves, the sound propagation is visualized within the classical muffler. The harmonic wave relation (7) is applied to transform the solution into the time domain: Find the maximum value of the sound pressure at , and set up a legend bar and ContourPlot options.
Visualize the sound propagation within the classical muffler.

See this note about improving the visual quality of the animation.

Within the inlet pipe and the first expansion chamber, the sound wave appears to be fixed in space but it simply oscillates in time. This type of wave is known as the standing wave, and often forms several local extrema of the sound pressure.

This excess sound pressure is called "Back-pressure" , which will restrict the gas flow within the muffler and reduce the power output of the car engine. To see how the perforated and padded mufflers improve this issue, it is better to visualize the amplitude distribution of the sound pressure waves.

Visualize the sound distribution at frequency .

Both the perforated and padded mufflers improve the effect of back-pressure that is seen in the classical muffler. Because the sound pressure wave is allowed to dissipate within the muffler pipes, the perforated muffler created the lowest back-pressure among the three cases.

To further investigate this results at different frequencies, the index of the solution data pfunTable can be altered.

### Transmission Loss

Another important quantity to measure the performance of an acoustic muffler is the transmission loss (TL), which is defined as the ratio between the incident and the transmitted sound power. The formulas for the transmission loss, sound power and sound intensity are given by: Calculate the sound intensity and the sound power at the inlet pipe.
Calculate the sound intensity and the sound power at the outlet pipe for three types of the muffler.
Calculate the and visualize transmission loss spectrum.

For the classical muffler, the transmission loss curve appears in a wavelike pattern with few dips. These dips correspond to the resonance frequencies of the muffler, where the sound transmission is enhanced and the muffler becomes less effective. The padded muffler, however, not only smooths out these dips but also increases the overall transmission loss values, especially for the higher frequencies, making it a more effective muffler in that range. Note that although the perforated muffler has the best effect on the back-pressure issue, it has generally smaller transmission loss values over the spectrum except for the frequency around .

More advanced mufflers can be constructed by combining the features of absorptive pads and perforation . If there is a targeting frequency range where the noise is to be controlled, the dimension of each expansion chamber can also be customized to meet the functional requirements.

### Unmuffled versus Muffled Engine Noise

To hear the damping effect of the acoustic muffler, an audio file of unmuffled engine noise is used as the sound input.

Extract the time-domain sound pressure data from the audio file.
Get the data length and the sampling rate of the sound input.

In order to process the sound signal, a Fourier transform is used to convert the input noise from the time to the frequency domain.

Convert the noise signal by the Fourier transform.

As explained in a separate tutorial, the discrete Fourier transform generates Fourier coefficients with the symmetric property that . For each pair of the Fourier coefficients and , the corresponding frequency component is .

Set the table of frequency components.

Each frequency component is then processed by the frequency response function, which is defined as the ratio between the input and the output sound amplitude: .

Get the frequency response function at the muffler outlet.
Define a function to process each frequency component and its symmetric component .
Convert the processed signal back into the time domain, and generate the audio of the muffled noise.
Compare the amplitude spectrum of the unmuffled and muffled engine noise.

Based on the output audio and the amplitude spectrum, the engine noise has been significantly damped by all of the mufflers. For the low frequency range which the noise is mostly composed of, the padded muffler shows the best sound attenuation.

## Nomenclature

 Symbol Description Unit ρ density [kg/m3] c speed of sound [m/s] p sound pressure [Pa] p sound amplitude [Pa] ω sound wave angular frequency [rad/s] f sound wave frequency [Hz] X position vector [m] R radius of the expansion chamber [m] L length of the expansion chamber [m] r radius of the inlet/outlet pipe [m] l length of the inlet/outlet pipe [m] Rf flow resistivity [N·s/m3] λ eigenvalue [rad2·m/kg] TL transmission loss N/A W sound power [W] I sound intensity [W/m2]

## References

1.  Potente, D., General Design Principles for an Automotive Muffler. Proceedings of ACOUSTICS, Busselton, Western Australia: 2005, pp. 153-158.

2.  Kirby, R. and Cummings, A., Bulk acoustic properties of rigid fibrous absortives extended to low frequencies. Proceedings of Euro-noise 95, Lyon, France: 1995, pp. 835-840.

3.  Delany, E. and Bazley, N., Acoustical properties of fibrous materials. Applied Acoustics 3: 1970, pp. 105-116.

4.  Fahy, F., Foundations of Engineering Acoustics. Academic Press: 2001.

5.  Cox, J. and D'Antonio, P., Acoustic absorbers and diffusers: Theory, design, and application. London: Spon Press: 2004.