Studying Acoustic Behavior of BFRP Laminated Composite in Dual-Chamber Muffler Application Using Deep Learning Algorithm

Over the last two decades, several experimental and numerical studies have been performed in order to investigate the acoustic behavior of different muffler materials. However, there is a problem in which it is necessary to perform large, important, time-consuming calculations particularly if the muffler was made from advanced materials such as composite materials. Therefore, this work focused on developing the concept of the indirect dual-chamber muffler made from a basalt fiber reinforced polymer (BFRP) laminated composite, which is a monitoring system that uses a deep learning algorithm to predict the acoustic behavior of the muffler material in order to save effort and time on muffler design optimization. Two types of deep neural networks (DNNs) architectures are developed in Python. The first DNN is called a recurrent neural network with long short-term memory blocks (RNN-LSTM), where the other is called a convolutional neural network (CNN). First, a dual-chamber laminated composite muffler (DCLCM) model is developed in MATLAB to provide the acoustic behavior datasets of mufflers such as acoustic transmission loss (TL) and the power transmission coefficient (PTC). The model training parameters are optimized by using Bayesian genetic algorithms (BGA) optimization. The acoustic results from the proposed method are compared with available experimental results in literature, thus validating the accuracy and reliability of the proposed technique. The results indicate that the present approach is efficient and significantly reduced the time and effort to select the muffler material and optimal design, where both models CNN and RNN-LSTM achieved accuracy above 90% on the test and validation dataset. This work will reinforce the mufflers’ industrials, and its design may one day be equipped with deep learning based algorithms.


Introduction
A muffler is considered as a common device for passive noise elimination. This device is widely used in various equipment and machines which are connected with ducts or exhaust systems, such as diesel engines, HVAC system, compressors, blowers and fans for ventilation [1]. The laminated composite muffler is also used in military equipment and machineries such as tanks, heavy, and light cannons, due to the unique thermo-acoustic and mechanical characteristics of composite materials (see Figure 1). However, mufflers are often limited by space availability. Consequently, there has been rising interest in designing mufflers to optimize acoustic characteristics. There are several optimized TL methods for the optimization of muffler shape under space constraints that have been developed by several researchers over the last two decades [2][3][4]. Numerous studies have also been and machineries such as tanks, heavy, and light cannons, due to the unique thermo-acoustic and mechanical characteristics of composite materials (see Figure 1). However, mufflers are often limited by space availability. Consequently, there has been rising interest in designing mufflers to optimize acoustic characteristics. There are several optimized TL methods for the optimization of muffler shape under space constraints that have been developed by several researchers over the last two decades [2][3][4]. Numerous studies have also been carried out by many researchers on mufflers to estimate the noise attenuation performance. These studies have suggested various evaluation methods for the actual noise reduction in mufflers when an optimally designed muffler is mounted on a duct [5][6][7].

Related Work
Since the focus of this research is on the application of laminated composite mufflers, the literature on acoustic behavior of composite materials had to also be investigated. Presently, research on the acoustic behavior of composite materials is becoming increasingly mature. This is largely due to good dynamic and mechanical properties, especially the high strength/lightweight properties, of these types of materials. Numerous studies have been carried out on composite materials to estimate the acoustic behavior using both experimental [8,9] and numerical prediction methods [10,11], while the conclusions have shown that the acoustic behavior of composite laminated mufflers need to be further investigated in order to develop a better understanding of the acoustic performance of these types of mufflers.

Acoustic Transmission Loss (TL)
Transmission loss (TL) generally describes the cumulative degradations in the acoustic power transmission of a waveform as a wave propagates out a source, or as it propagates through a given area or through a particular type of structure such as mufflers. It is a frequently used terminology in acoustics and is measured in decibels. TL measurements are very important in the industry of acoustic devices such as mufflers and sonars.
To improve the TL accuracy of mufflers, several simulation studies have been reported on how to obtain an optimized shape for mufflers.
Numerical methods such as the finite element method (FEM) [12] and the boundary element method (BEM) [13,14] are suitable for predicting the TL performance of reactive

Related Work
Since the focus of this research is on the application of laminated composite mufflers, the literature on acoustic behavior of composite materials had to also be investigated. Presently, research on the acoustic behavior of composite materials is becoming increasingly mature. This is largely due to good dynamic and mechanical properties, especially the high strength/lightweight properties, of these types of materials. Numerous studies have been carried out on composite materials to estimate the acoustic behavior using both experimental [8,9] and numerical prediction methods [10,11], while the conclusions have shown that the acoustic behavior of composite laminated mufflers need to be further investigated in order to develop a better understanding of the acoustic performance of these types of mufflers.

Acoustic Transmission Loss (TL)
Transmission loss (TL) generally describes the cumulative degradations in the acoustic power transmission of a waveform as a wave propagates out a source, or as it propagates through a given area or through a particular type of structure such as mufflers. It is a frequently used terminology in acoustics and is measured in decibels. TL measurements are very important in the industry of acoustic devices such as mufflers and sonars.
To improve the TL accuracy of mufflers, several simulation studies have been reported on how to obtain an optimized shape for mufflers.
Numerical methods such as the finite element method (FEM) [12] and the boundary element method (BEM) [13,14] are suitable for predicting the TL performance of reactive and dissipative mufflers and are not limited to the geometry method of mufflers. The radiation impedance form of an unflanged circular pipe was developed by Levine and Schwinger [15], and another form of radiation impedance procedure was suggested for an infinitely large flanged circular pipe by Norris and Sheng [16]. For a radiation impedance formula, Silva et al. [17] presented approximate mathematical expressions for a complicated end shape. The radiation impedance and the reflection coefficient relationship for a circular pipe was presented by Polack et al. [18]. The radiation impedance measurement methods were summarized and calculated for several boundary conditions by Dalmont et al. [19,20].

Mufflers Design Optimization
In most muffler design problems for sound reduction, the muffler has been designed in two separate parts, the main pipe or duct and the muffler part. The TL in the muffler has been widely used for evaluating the muffler performance during the optimization process of the muffler size [21][22][23]. However, the optimal muffler performance of the sound attenuation may be degraded during the installation process due to the long distance between the main pipe and the radiation impedance. The length of the main pipe mounted on the muffler is usually longer than the length of the outlet of a muffler. These problems motivated the authors to find an effective method to optimize the muffler size in order to improve the evaluation of the sound attenuation performance in a muffler. Lee et al. [24] studied developing the multi-chamber muffler with dissipative elements to reduce TL of noise. They used FEM and experimental methods to design a muffler that consists of dual expansion connected in a series of flow systems that effectively optimize sound in the speech interference range [25]. Lee et al. [26] proposed a multi-objective topology optimization problem to maximize the TL at the target frequency while minimizing the voltage drop. Shi et al. [27] studied experimentally and theoretically the propagation of sound wave in a periodic array pipe mufflers with micro-holes. They found that the periodic structure affected the performance of the muffler with micro-holes. Du et al. [28] studied experimentally the acoustic performance of water muffler based on a Kevlar-reinforced rubber tube for reducing the noise by optimizing vibration damping and the hydrodynamic noise reduction method.

Deep Neural Networks (DNNs)
The artificial neuronal network (ANN) is a mathematical model of informationprocessing established by imitating the structure and function of the brain neuronal network. This method implies repeated learning and training of known information, and includes a scheme to gradually adjust the learning weights of neurons to imitate the relationship between the input and output data. Compared to the traditional method, the ANN technique has obvious advantages in the processing of blurred data, random data and nonlinear data, and is especially suitable for complex large-scale structures and unclear information systems.
DNN is an ANN architecture with multiple hidden layers. Depth generally refers to the amount of hidden layers in the structure of the DNN. Learning is a cognitive process of the unknown to the known, corresponding to the input and output process in the structure of the DNN. Deep learning strengthens the learning ability of feature and makes the prediction more accurately through the hierarchical expression of the input information and the transformation of present layer feature by a previous layer.

Contribution
Generally, mufflers reduce the sound pressure generated by sound-generating sources such as motors, fans, etc. existing in systems of a vehicle exhaust or the home ventilation by generating dissipative sound waves caused by geometric discontinuity by means of acoustic impedance difference. Therefore, internal partitions should be optimally placed to improve the noise attenuation performance in the main noise frequency range when the outer size of the muffler is limited. Muffler material is one of the important parameters in sound insulation and absorption. The most methods used for evaluating the acoustical performance of a muffler (TL) are the FEM, traditional laboratory methods, four-pole transfer matrix, and the three-point methods. The presented work utilizes two types of ANN architectures which are called RNN-LSTM and CNN. These two types of ANNs are developed in Python in order to predict the acoustic behavior of a DCLCM manufactured from a BFRP laminated composite. Figure 2 shows a flowchart of the proposed method. The PTC and the TL of an acoustic muffler will be calculated from the exact solution of the governing acoustic equations of the muffler model. After the essential characteristics of the muffler have been extracted, the muffler characteristics are fed to an ANN as labeled data. The inputs and outputs weights of the network are calculated. Then, we compute the derivative of the error for weights by the backpropagation (BP) algorithm. If the error target is not acceptable "case no", the muffler characteristics data are fed again to an ANN to train and calculate new weights to compute the new error and so on until the error target is acceptable "case yes". The training and model parameters were optimized by using BGA Optimization. The proposed work will be discussed more in depth in the following sections.  Figure 3 shows the geometry of a circular DCLCM considered in this study. An acoustic DCLCM was constructed from a BFRP laminated composite, the staking distribution to three symmetrically plies is [0/90°/0]s, and the thickness of each ply is 5 mm. The BFRP Physical and mechanical properties are shown in Table 1. The geometric information used for the calculations in this work is summarized in Table 2. The values of the muffler thickness, B, to the total chamber length, L, are listed in Table 2. The range of frequencies used for PTCs calculation of muffler is [0-3 KHz]. The muffler length, B, is optimized in order to maximize the acoustic TL. At the right end of the muffler, an anechoic termination is assumed.  In this work, to develop a monitoring strategy for the acoustic behavior performance of DCLCM, a novel deep learning algorithm is utilized. The muffler acoustic performance is calculated for a wide range of geometric data of the DCLCM using commercial software. An RNN-LSTM and a CNN are trained, evaluated, and developed via Python. BGA optimization is used to optimize the training and the parameter selection for the model. To the best of the authors' knowledge, the methodology presented in this paper and the results they obtained is an original contribution to scientific research for monitoring an acoustic behavior in composite muffler materials.

Case Study
To achieve the goals presented above, this work is arranged as follows: Section 2 presents system modeling and discusses a mathematical model including case studies for composite muffler systems. Results and discussion are presented in Section 3, including an algorithm established and proposed utilizing a deep neural network. Conclusions are presented in Section 4, and a brief summary as well as a plan of future work are reported in Section 5. Figure 3 shows the geometry of a circular DCLCM considered in this study. An acoustic DCLCM was constructed from a BFRP laminated composite, the staking distribution to three symmetrically plies is [0/90 • /0] s , and the thickness of each ply is 5 mm. The BFRP Physical and mechanical properties are shown in Table 1. The geometric information used for the calculations in this work is summarized in Table 2. The values of the muffler thickness, B, to the total chamber length, L, are listed in Table 2. The range of frequencies used for PTCs calculation of muffler is [0-3 KHz]. The muffler length, B, is optimized in order to maximize the acoustic TL. At the right end of the muffler, an anechoic termination is assumed.    Figure 4 presents the mathematical model of DCLCM which consists of pipes and two expansion chamber pipes being identified. As shown in th = is the area of the straight pipe, and = is the area of the muffler    Table 1 shows the physical and mechanical properties of the BFRP, where ρ is the material density, E 11 , E 22 are elastic modulus in the '1'and '2' directions, respectively, G 12 , G 21 are the shear modulus, υ 21 , υ 12 are the Poisson's ratio of the transverse strain in the directions '1' and '2' caused by the normal stress in the directions '2' and '1', respectively.

Case Study
2.1.1. Basic Acoustic Equations of the Dual-Chamber Muffler Figure 4 presents the mathematical model of DCLCM which consists of three straight pipes and two expansion chamber pipes being identified. As shown in the figure, S 1 = S 3 = S 5 is the area of the straight pipe, and S 2 = S 4 is the area of the muffler. Eight points were chosen to represent the flow condition inside the muffler (pt 1~p t 8 ). At each two consecutive points, the standing pressure P s and reflected P r pressure are similar e.g., P s 1 + P r 1 = P s 2 + P r 2 ; in addition, the standing U s volume velocity and reflected U r volume velocity are similar e.g., U s 1 + U r 1 = U s 2 + U r 2 . The continuity algorithm of pressure and volume velocity on continuity junctions numbers (1 − 1), (1 − 2), (2 − 1), (2 − 2) are applied to compute junctions' transmission coefficients A 11 , B 11 , A 12 , B 12 , A 21 , B 21 , A 22 , B 22 , by solving the Helmholtz Equation (1) [30], and applying the junction boundary condition (B.C) at each junction: where is the junction number, = is wavenumber ratio, = ⁄ is the wavenumber, is the angular frequency, = 330 [ ] ⁄ is the sound speed in air, and is the bending wavenumber of the laminated composite pipe.

Acoustic Properties of Composite Laminated Muffler
The laminated muffler studied in this paper is made of BFRP composites in general. Due to the low strength and stiffness characteristics in the transverse direction of laminate, it does not consist only of the unidirectional lamina. Thus, some laminas in most laminates are placed at an angle. Therefore, the stress-strain relationship for an angle lamina must be developed. The new axis is called local axes in the 1-2 coordinate system, where direction 1 is parallel to the fibers, and direction 2 is perpendicular to the fibers. The angle between the local axes in the 1-2 coordinate system and the global axes in the x-y coordinate system is the fiber angle . The plane stress transformed reduced stiffness coefficients For the equal size chamber systems shown in Figure 4, the three-point method was used for describing acoustical properties via absorbing materials assuming notation of the plane wave propagation since the transverse system dimensions are smaller than one wavelength over the frequency range of interest. A standing pressure propagating wave amplitude A j and a reflected pressure propagating wave amplitude B j are assumed. The standing P s and reflected P r pressure wave can be expressed in Equations (2) and (3) respectively as: where j is the junction number, is wavenumber ratio, k = ω/c is the wavenumber, ω is the angular frequency, c = 330[m/s] is the sound speed in air, and k b is the bending wavenumber of the laminated composite pipe.

Acoustic Properties of Composite Laminated Muffler
The laminated muffler studied in this paper is made of BFRP composites in general. Due to the low strength and stiffness characteristics in the transverse direction of laminate, it does not consist only of the unidirectional lamina. Thus, some laminas in most laminates are placed at an angle. Therefore, the stress-strain relationship for an angle lamina must be developed. The new axis is called local axes in the 1-2 coordinate system, where direction 1 is parallel to the fibers, and direction 2 is perpendicular to the fibers. The angle between the local axes in the 1-2 coordinate system and the global axes in the x-y coordinate system is the fiber angle θ. The plane stress transformed reduced stiffness coefficients Q ij of the lamina can be expressed in terms of the engineering notations (see Table 1) as: Q 11 = Q 11 cos 4 θ + 2(Q 12 + 2Q 66 ) sin 2 θ cos 2 θ + Q 22 sin 4 θ (5) The bending stiffness D ij can be calculated from: Using the above analysis of composite material, it can be found that the muffler performance can be greatly affected by the laminating plies staking and fiber angle θ. When the muffler is excited by sound waves, there is a relationship between the sound pressure inside such a DCLCM and the normal vibration velocity [31].

Acoustic Transmission Loss (TL)
In mufflers, the theoretical definition of TL is the logarithmic ratio of incident to power transmission for the case of reflection-free terminations. This can be expressed in terms of sound pressure by solving Equation (10): where k b is the bending wavenumber of the laminated composite pipe: where m = S 2 S 1 is the area ratio. Figure 5 presents the algorithm flow chart of the Equation (10) solution technique to find the acoustic pressure P j wave in each junction (j) of the DCLCM model (see Figure 4). By applying algorithms of continuity of pressure (P) in Equation (12) and volume velocity (U) in Equation (13) and substituting with boundary conditions (B.C) at each junction (j) in the DCLCM model based on the flow condition inside the muffler (pt 1~p t 8 ) and section area S j (see Figure 4): where j is the junction number, and U denotes the axial acoustic velocity which can be obtained by the momentum equation: where ρ is the fluid density, and i = √ −1 is the imaginary unit. Then, the eigenfunctions' orthogonal properties are used [32], where L 1 = L 2 = L and S 1 = S 3 = S 5 , S 2 = S 4 , and the junction boundary condition (B.C) is applied. Consequently, the resultant sound TL of DCLCM model resulted in being:

Artificial Neural Networks
ANN is the most used algorithm between different artificial intelligence (AI) algorithms for the advanced nonlinear problems solution [33,34]. Each individual network of ANN consists of the number of computational nodes, and each node is used for processing the inputs and transferring the input calculation results to output connections. Each node output may be an input to another node or more. Weights and biases are used to scale and bias the outputs, respectively, e.g., in the function y = mx + b, y and x refer to output and new output, respectively, m is the weight, and b is the bias. Some networks are activated to determine the output with the type of function being linear or non-linear by adding activation functions. The rectified linear unit (ReLU) is the activation function most used in deep learning [35].
Before adjusting the weights and biases in ANN, the model does not work well, i.e., the ANN model is not trained. The learning of NN can be done automatically from raw data as a hierarchical feature representation [36,37] or can be trained via case study examples. In our models designed in this work, we used supervised learning to train, by comparing the training data and model's predictions to actual data. The training accuracy can be improved via updating the trainable parameters to optimize the error between prediction and actual data.
ANN can be determined depending on the type of input data, such as the case study in this work; the ANN models are used to analyze the time-series data.
In this work, the two major ANN model types for analyzing the time series data that were discussed are the RNN-LSTM and CNN.

Artificial Neural Networks
ANN is the most used algorithm between different artificial intelligence (AI) algo rithms for the advanced nonlinear problems solution [33,34]. Each individual network o ANN consists of the number of computational nodes, and each node is used for processin the inputs and transferring the input calculation results to output connections. Each nod output may be an input to another node or more. Weights and biases are used to scale an bias the outputs, respectively, e.g., in the function y = mx + b, y and x refer to output an new output, respectively, m is the weight, and b is the bias. Some networks are activate to determine the output with the type of function being linear or non-linear by addin activation functions. The rectified linear unit (ReLU) is the activation function most use in deep learning [35].
Before adjusting the weights and biases in ANN, the model does not work well, i.e the ANN model is not trained. The learning of NN can be done automatically from raw data as a hierarchical feature representation [36,37] or can be trained via case study exam ples. In our models designed in this work, we used supervised learning to train, by com paring the training data and model's predictions to actual data. The training accuracy ca be improved via updating the trainable parameters to optimize the error between predic tion and actual data.
ANN can be determined depending on the type of input data, such as the case stud in this work; the ANN models are used to analyze the time-series data.
In this work, the two major ANN model types for analyzing the time series data tha were discussed are the RNN-LSTM and CNN.
CNNs are used to analyze the groups of data such as the time series, images, sen tences, sound recordings, etc. Weight matrices in CNNs are applied as kernels or filters t

Convolutional Neural Network (CNN)
CNNs are used to analyze the groups of data such as the time series, images, sentences, sound recordings, etc. Weight matrices in CNNs are applied as kernels or filters to extract the features [38,39].
As shown in Figure 6, a proposed CNN model typically consists of feature extraction through a stack of layers on the input layer such as convolution, activation and poling, and classification through fully connected layers for outputting the scores for each class. Each layer is responsible for different functions and uses the result from the previous layer as the input.
Equation (16) describes the operation of proposed CNN, and this process is defined as: where x l j is the ith output map in layer l; x l−1 i is the ith output map in layer l−1; w l ij is the weight; b l j is the bias; f (·) is a nonlinear function that is applied component-wise.
As shown in Figure 6, a proposed CNN model typically consists of feature extraction through a stack of layers on the input layer such as convolution, activation and poling, and classification through fully connected layers for outputting the scores for each class. Each layer is responsible for different functions and uses the result from the previous layer as the input. Equation (16) describes the operation of proposed CNN, and this process is defined as: where is the ith output map in layer ; −1 is the ith output map in layer −1; is the weight; is the bias; (•) is a nonlinear function that is applied component-wise.

Recurrent Neural Network with Long Short-Term Memory Blocks (RNN-LSTM)
RNNs deeply analyze the time series data via applying the feedback loops to original ANN [40]. The biggest disadvantage in RNNs is known as the vanishing gradient problem, where, during the backpropagation process, the error signal used to train the network exponentially decreases the further you travel backwards in RNN, thus sometimes using computational nodes known as LSTM to relieve this problem, as shown in Figure 7. The data feature extraction is done from the first layers of ANN. These layers are responsible for extracting significant information from the input data [41].
LSTM is a special type of RNN with gating mechanism and memory cells, which greatly improves the performance of RNNs. There are three types of gates within each LSTM cell: input gate, forget gate, and output gate, and these gates define the state of each memory cell by using sigmoid as the activation function to cause information to be transmitted selectively. The memory cell that retains the long-term state is the key architecture of each LSTM cell. The internal architecture of a single LSTM cell is shown in Figure  8.

Recurrent Neural Network with Long Short-Term Memory Blocks (RNN-LSTM)
RNNs deeply analyze the time series data via applying the feedback loops to original ANN [40]. The biggest disadvantage in RNNs is known as the vanishing gradient problem, where, during the backpropagation process, the error signal used to train the network exponentially decreases the further you travel backwards in RNN, thus sometimes using computational nodes known as LSTM to relieve this problem, as shown in Figure 7. The data feature extraction is done from the first layers of ANN. These layers are responsible for extracting significant information from the input data [41].  LSTM is a special type of RNN with gating mechanism and memory cells, which greatly improves the performance of RNNs. There are three types of gates within each LSTM cell: input gate, forget gate, and output gate, and these gates define the state of each memory cell by using sigmoid as the activation function to cause information to be transmitted selectively. The memory cell that retains the long-term state c t is the key architecture of each LSTM cell. The internal architecture of a single LSTM cell is shown in Figure 8.  We can describe the operation of the three gates presented in Figure 8 in Equations (17)- (19). Equations (20) where , , , and represent the weight matrices of LSTM; , , , and denote the bias vector of LSTM; , , and are forget gate, input gate, and output gate vectors at time t; −1 and ́ mean, respectively, the previous cell state and a new candidate value. ( ) and ℎ( ) are utilized as the activation functions, as shown below: We can describe the operation of the three gates presented in Figure 8 in Equations (17)- (19). Equations (20)- (22) suggest the cell states c t and the hidden state h t of each LSTM unit at time t: where W f , W i , W c , and W o represent the weight matrices of LSTM; b f , b i , b c , and b o denote the bias vector of LSTM; f t , i t , and o t are forget gate, input gate, and output gate vectors at time t; c t−1 and c t mean, respectively, the previous cell state and a new candidate value. σ(z) and tanh(z) are utilized as the activation functions, as shown below:

Bayesian Genetic Algorithms (BGA) Optimization
Generally, there are several strategies for BGA optimization modeling of objective functions (f ), such as Gaussian processes [42][43][44], random forests [41], and tree-structured Parzen estimators [45,46]. Figure 9 presents the comparison between the Grid search and BGA Optimization tuning method for tuning the model's hyperparameters. As shown in the figure, the yellow dot refers to the assessment of the model in each method; note that the grid search method may be searched in rough space, but the BGA optimization methods can test any possible combination in space and intelligently suggest combinations to obtain optimal solutions with fewer evaluations. In this research, we applied the BGA techniques for quick, intelligent tuning. functions (f), such as Gaussian processes [42][43][44], random forests [41], and tree-structured Parzen estimators [45,46]. Figure 9 presents the comparison between the Grid search and BGA Optimization tuning method for tuning the model's hyperparameters. As shown in the figure, the yellow dot refers to the assessment of the model in each method; note that the grid search method may be searched in rough space, but the BGA optimization methods can test any possible combination in space and intelligently suggest combinations to obtain optimal solutions with fewer evaluations. In this research, we applied the BGA techniques for quick, intelligent tuning.

The Work Description
With the desire to explore alternative indirect Monitoring Acoustic Behavior of DCLCM frameworks and inspired by deep learning, this work discusses the framework of deep learning that analyzes TL data of DCLCM to optimize the acoustic TL of the DCLCM. To evaluate the idea of the presented work, the work was divided into the following sections: (1) Data Collecting: The ability and accuracy of the designed ANN are based on the volume of feeding data to the ANN; this accuracy will be improved when there are more data for training. In the present work, the data were computed via MATLAB from a DCLCM model. This dataset was used for training, validation, and test for the ANN. (2) The Established Algorithm: from the data collected in the previous step, via Python, an RNN-LSTM and CNN are developed and then, using BGA, we are training parameters and tuning models. Figure 10 presents the flowchart of framework for the proposed ANNs used in this work. As shown in Figure 10, the weight coefficients between the processing neurons are adjusted to reach some desired goal during the training process. The neurons are interconnected via feed-forward considered links for the basic neuron computation. The multiplication results of the neuron input and the connection weight between input layer and hidden layer are obtained. Then, a bias value is added to them using Equation (16). The results of calculation processes are subjected to an activation function. Activation function generates neuron output results. Similar operations can be applied to the output layer. The activation function determines the input/output behaviors of the network. The error,

The Work Description
With the desire to explore alternative indirect Monitoring Acoustic Behavior of DCLCM frameworks and inspired by deep learning, this work discusses the framework of deep learning that analyzes TL data of DCLCM to optimize the acoustic TL of the DCLCM. To evaluate the idea of the presented work, the work was divided into the following sections: (1) Data Collecting: The ability and accuracy of the designed ANN are based on the volume of feeding data to the ANN; this accuracy will be improved when there are more data for training. In the present work, the data were computed via MATLAB from a DCLCM model. This dataset was used for training, validation, and test for the ANN. (2) The Established Algorithm: from the data collected in the previous step, via Python, an RNN-LSTM and CNN are developed and then, using BGA, we are training parameters and tuning models. Figure 10 presents the flowchart of framework for the proposed ANNs used in this work. As shown in Figure 10, the weight coefficients between the processing neurons are adjusted to reach some desired goal during the training process. The neurons are interconnected via feed-forward considered links for the basic neuron computation. The multiplication results of the neuron input and the connection weight between input layer and hidden layer are obtained. Then, a bias value is added to them using Equation (16). The results of calculation processes are subjected to an activation function. Activation function generates neuron output results. Similar operations can be applied to the output layer. The activation function determines the input/output behaviors of the network. The error, known as the difference between the ANN output and desired output, is calculated. Finally, the new weight values are obtained. This method must be repeated until the error is acceptable. This training method is continued for all the data in the training dataset. known as the difference between the ANN output and desired output, is nally, the new weight values are obtained. This method must be repeated is acceptable. This training method is continued for all the data in the train

Validation of the Proposed Method
In this subsection, a convergence investigation is carried out for method. The acoustic TL at all octave Band Center Frequencies (OBCFs) and compared with available experimental results in literature. Table 3 pres gence and comparison study for a muffler structure optimization design to noise in the pumping system by Liu et al. [47]. The muffler is silenced at th pump pipeline system, the muffler material is commercial steel and the str fler is shown in Figure 11.

Validation of the Proposed Method
In this subsection, a convergence investigation is carried out for the proposed method. The acoustic TL at all octave Band Center Frequencies (OBCFs) are calculated and compared with available experimental results in literature. Table 3 presents a convergence and comparison study for a muffler structure optimization design to eliminate the noise in the pumping system by Liu et al. [47]. The muffler is silenced at the outlet of the pump pipeline system, the muffler material is commercial steel and the structure of muffler is shown in Figure 11.

Validation of the Proposed Method
In this subsection, a convergence investigation is carried out method. The acoustic TL at all octave Band Center Frequencies (OBC and compared with available experimental results in literature. Table 3 gence and comparison study for a muffler structure optimization desi noise in the pumping system by Liu et al. [47]. The muffler is silenced pump pipeline system, the muffler material is commercial steel and th fler is shown in Figure 11.  The experimental test device used by Liu consists of a booster pump, motor, valve, pressure transmitter, flowmeter, muffler, data acquisition system, etc. They studied the noise reduction performance by using three values for the expansion angle of the flow channel (θ) which are 120 • , 145 • , and 160 • . They found that, when the extension angle was 145 • , the muffler had the best sound attenuation effect. In this validation, the measured data of the experiment for muffler with extension angle 145 • were compared with the data of the proposed method simulation data. As shown in Table 3 and Figure 12, we can see the acoustic TL at all OBCFs that were obtained by test of the muffler had the same trend as that of the proposed method calculation.

Data Collecting
In this work, the dimensions of the DCLCM model and all S/C ratios (B/L) that were used in the most calculations in the analytical section are included in Table 2

Data Collecting
In this work, the dimensions of the DCLCM model and all S/C ratios (B L ⁄ ) that were used in the most calculations in the analytical section are included in Table 2. A MATLAB model for the DCLCM representation as shown in Figure 3 was run at various S/C ratios to generate simulated examples of DCLCM system response. The simulation solves the system acoustic equations to find acoustic pressure in each junction of the DCLCM for every pressure and volume velocity at each junction in the DCLCM system. The simulation inputs are extracted from below Figures 13 and 14 and their accompanying derivations are presented in Equations (9)- (30).
The acoustic TL properties for the proposed muffler are illustrated in Figures 12 and  13 (for B/L = 0, 1.0, 0.2, 0.4, 0.6, 0.8, 1.0), where high attenuation is evident over a wide frequency range. Figure 13 presents the acoustic TL distribution over all OBCFs.
The DCLCM behavior is presented in Figure 13, where the acoustic TL distribution changes over the frequency range with the change in muffler dimensions. Therefore, we find that the expansion chamber without a muffler (S/C ratio equal zero) has equal domes, and the higher the S/C ratio, the more it changes the domes to an unequal shape, and the the first dome has the smaller amplitude and frequency band than another one at each muffler geometrical configuration in S/C ratios presented in Table 2. The rate of rise and fall of domes increases the acoustic TL with the increase in frequency so that it reaches a peak at a certain frequency that then decreases and so on. In addition, the S/C ratio ( / ) in the muffler geometry has an effective effect in sound attenuation of TL value, and we find that the increases in the S/C ratio ( / ) have an effect on widening of the second dome and tend to cover two domes or more of the acoustic TL for the lower S/C ratios. The cut-off frequency on which the behavior of acoustic TL operates is 3000 Hz. At all S/C ratios ( / ) presented in Table 2, the muffler effect on the TL value has almost sequentially vanished at frequencies 130, 615, 825, 1440, 1780, 2270, and 2600 Hz, respectively. Figure 14 shows the acoustic TL distribution with respect to S/C ratio over one completed period of frequency from 0 to 825 Hz. As shown in Figure 14, at stationary frequency, the value of acoustic TL gain changes with the increases of S/C ratio ( / ), which indicates that the acoustic TL has a high sensitivity to S/C ratio ( / ) changes. Consequently, we can consider that the S/C ratio is the main parameter of frequency shifting, and that it is the powerful key in the acoustic TL gain.

The Established Algorithm
The establishing, training, and evaluating of a CNN and or RNN were done in Python software, and each neural network was developed via the same steps as listed in Algorithm 1 below based on steps discussed in Figure 10.

ANN Models Development
In this work, as shown in Figure 15, the deep learning neural network used mainly has an input layer of the muffler characteristics data, and three 1D convolutions' full connection (FC) layers for training each layer are 56 × 128, 28 × 256, 14 × 512, respectively, for two sub-sampling LSTM layers, each layer is 14 × 512, 7 × 512, respectively,  Figure 13 presents the acoustic TL distribution over all OBCFs.
The DCLCM behavior is presented in Figure 13, where the acoustic TL distribution changes over the frequency range with the change in muffler dimensions. Therefore, we find that the expansion chamber without a muffler (S/C ratio equal zero) has equal domes, and the higher the S/C ratio, the more it changes the domes to an unequal shape, and the the first dome has the smaller amplitude and frequency band than another one at each muffler geometrical configuration in S/C ratios presented in Table 2. The rate of rise and fall of domes increases the acoustic TL with the increase in frequency so that it reaches a peak at a certain frequency that then decreases and so on. In addition, the S/C ratio (B/L) in the muffler geometry has an effective effect in sound attenuation of TL value, and we find that the increases in the S/C ratio (B/L) have an effect on widening of the second dome and tend to cover two domes or more of the acoustic TL for the lower S/C ratios. The cut-off frequency on which the behavior of acoustic TL operates is 3000 Hz. At all S/C ratios (B/L) presented in Table 2, the muffler effect on the TL value has almost sequentially vanished at frequencies 130, 615, 825, 1440, 1780, 2270, and 2600 Hz, respectively. Figure 14 shows the acoustic TL distribution with respect to S/C ratio over one completed period of frequency from 0 to 825 Hz. As shown in Figure 14, at stationary frequency, the value of acoustic TL gain changes with the increases of S/C ratio (B/L), which indicates that the acoustic TL has a high sensitivity to S/C ratio (B/L) changes. Consequently, we can consider that the S/C ratio is the main parameter of frequency shifting, and that it is the powerful key in the acoustic TL gain.

The Established Algorithm
The establishing, training, and evaluating of a CNN and or RNN were done in Python software, and each neural network was developed via the same steps as listed in Algorithm 1 below based on steps discussed in Figure 10. output: score of CNN or RNN trained model on test dataset to predict TL for various S/C ratio (B/L) 4: let f be the feature set 3d matrix 5: for i in dataset do 6: let f i be the feature set matrix of sample I 7: for j in i do 8: v j ← vectorize (j, w) 9: append v j to f i 10: append f i to f 11: f train , f test , l train , l test ← split feature set and prediction into train subset and test subset 12: M ← CNN (f train , l train ) or RNN (f train , l train ) 13: score ← evaluate (I, l test , M) 14: return score 15: end for 16: end for

ANN Models Development
In this work, as shown in Figure 15, the deep learning neural network used mainly has an input layer of the muffler characteristics data, and three 1D convolutions' full connection (FC) layers for training each layer are 56 × 128, 28 × 256, 14 × 512, respectively, for two sub-sampling LSTM layers, each layer is 14 × 512, 7 × 512, respectively, and then all the resultant 2D arrays from pooled feature maps are converted into a single long continuous linear vector in a flattening layer that has 25, 088 elements in one linear vector, and a softmax layer as the activation function in the output layer that predicts a multinomial probability distribution of muffler characteristics such as TL datasets. In this work, both models of RNN and CNN have similar configurations, except the region close to the input layer for the feature extraction maps. The cost and accuracy of both models are evaluated, and the test set accuracy is also evaluated before and after training. The BGA technique was used to optimize the training speed and accuracy, respectively. Figure 16 shows the final performances for both models. As shown in Figure 16 presenting accuracy curves during the training of training and validation datasets, both models CNN and RNN achieve accuracy above 90% on test and validation datasets, saving sufficient time. Furthermore, we also calculate the training and testing time of RNN-LSTM and CNN as shown in Table 4. vector, and a softmax layer as the activation function in the output layer that predicts a multinomial probability distribution of muffler characteristics such as TL datasets. In thi work, both models of RNN and CNN have similar configurations, except the region close to the input layer for the feature extraction maps. The cost and accuracy of both model are evaluated, and the test set accuracy is also evaluated before and after training. Th BGA technique was used to optimize the training speed and accuracy, respectively.  the CNNs' weaknesses, but, in this work, this problem was not observed. As a result, for further understanding the weaknesses of the CNN, more work must be done using a more complex dataset. Simply, if directly compared against fundamental equations of the convolution layer, we can see the contrast in complexity. If the previous layer inputs have a large number of filters, the CNN model still easier and faster to train and more simple than the RNN-LSTM model .   In addition, the performance of the CNN model is better than the RNN-LSTM model in classification ability. The additional computations are what distinguish the performance of CNN over RNN-LSTM. These computations are for feeding the hidden layer from the previous step in order to provide long-range contextual information into the next step. This means that the RNN-LSTM is incorrectly based on the built-up memory. On the other hand, the CNN can train and classify input data quickly and accurately because the CNN extracts the features within windows of time including time-series data.

Performance Evaluation of ANN Models
The drawbacks of the CNN in the application of Monitoring Acoustic Behavior of DCLCM are the same as any applications of time-based series. The CNN performance is based on the dataset size and quality, but, in the presented work, the dataset is small and clean from noisy. In addition, in general, the poor random overfitting problem is one of the CNNs' weaknesses, but, in this work, this problem was not observed. As a result, for further understanding the weaknesses of the CNN, more work must be done using a more complex dataset. Simply, if directly compared against fundamental equations of the convolution layer, we can see the contrast in complexity. If the previous layer inputs have a large number of filters, the CNN model still easier and faster to train and more simple than the RNN-LSTM model.

Acoustic DCLCM Geometry Design Optimization
For selecting the optimum DCLCM Geometry Design, a genetic optimization method is used. A genetic algorithm (GA) is a stochastic global search and optimization method that mimics the metaphor of natural biological evolution. This process leads to the evolution of populations of individuals that are better suited to their environment than the individuals from which they were created, just as in natural adaptation. Figure 17 illustrates the flow chart of a genetic optimization method applied in this paper. In this case, the DCLCM junction locations as a string of integer numbers are encoded analogous to the genetic code on a DNA string; accordingly, the DCLCM acoustic outputs of TL are propagated over different OBCFs, and it is analogously identical to the behavior of the breeding population for a number of individuals, each characterized by its DNA, and each individual is determined according to some fitness function. Here, it is preferable to raise individuals with high fitness (high OBCFs) in the breeding process, so that useful genes tend to propagate through the generations (High OBCFs) and detrimental genes disappear (low OBCFs), just like with the range of acoustic frequencies generated by DCLCM. Usually, a GA code was written using MATLAB.
behavior of the breeding population for a number of individuals, each characterized by its DNA, and each individual is determined according to some fitness function. Here, it is preferable to raise individuals with high fitness (high OBCFs) in the breeding process, so that useful genes tend to propagate through the generations (High OBCFs) and detrimental genes disappear (low OBCFs), just like with the range of acoustic frequencies generated by DCLCM. Usually, a GA code was written using MATLAB.
As the muffler development depends on optimizing its geometry design for high performance, and because TL is an essential characteristic of the muffler, as well as the S/C ratio ( ⁄ ) having a high sensitivity to acoustic TL, the ( ⁄ ) is therefore a powerful parameter in the gain of acoustic TL. Actually, the ( ⁄ ) is the key parameter for muffler development. In this section, we will use the generated simulated acoustic output of deep learning analysis of TL for acoustic TL maximization. The optimal value of the dimension (B) is planned and carried out. From the final performances for both models of ANNs used in this work, the performance of the CNN model is better than the RNN-LSTM model; therefore, we will use the acoustic TL output of the CNN model in the optimization processes. It should be noted that the derivation processes of the maximum value of the acoustic TL were achieved using MATLAB software. Figures 18 and 19 show the results of these derivations in order to pick out the suitable dimension (B) which maximizes the value of TL. Figure 19 shows the relationship between the values of TL and frequency range of DCLCM with optimal value of the S/C ratio (B/L), compared with the system without a As the muffler development depends on optimizing its geometry design for high performance, and because TL is an essential characteristic of the muffler, as well as the S/C ratio (B/L) having a high sensitivity to acoustic TL, the (B/L) is therefore a powerful parameter in the gain of acoustic TL. Actually, the (B/L) is the key parameter for muffler development.
In this section, we will use the generated simulated acoustic output of deep learning analysis of TL for acoustic TL maximization. The optimal value of the dimension (B) is planned and carried out. From the final performances for both models of ANNs used in this work, the performance of the CNN model is better than the RNN-LSTM model; therefore, we will use the acoustic TL output of the CNN model in the optimization processes. It should be noted that the derivation processes of the maximum value of the acoustic TL were achieved using MATLAB software. Figures 18 and 19 show the results of these derivations in order to pick out the suitable dimension (B) which maximizes the value of TL. muffler. As shown in Figure 19, the maximum value of TL can be obtained at the DCLCM resonance frequency (612 Hz), which is estimated analytically. Therefore, the maximum attenuation value will be 60 dB. It is possible to conclude that the optimization method is the better way to obtain the maximum value of TL corresponding to the defined dimensions range of DLCS. As a result of the DLCS dimensions optimization, the optimal value of S/C ratio (B/L) is 0.106. This value is corresponding to the higher-level sound attenuation of TL value.

Conclusions
This work studies the solutions of drawbacks of experimental and numerical works to study an acoustic behavior of different mufflers' materials to improve its performance, where it is necessary to perform large, important, time-consuming calculations, particularly if the muffler is made from advanced materials such as composite materials. All of these problems were resolved with modern methods in this study by using a deep learning algorithm to predict the acoustic behavior of the muffler material in order to save effort and time on muffler design optimization. The acoustic behavior of DCLCM made from a BFRP laminated composite was predicted in Python through two types of DNN muffler. As shown in Figure 19, the maximum value of TL can be obtained at the DCLCM resonance frequency (612 Hz), which is estimated analytically. Therefore, the maximum attenuation value will be 60 dB. It is possible to conclude that the optimization method is the better way to obtain the maximum value of TL corresponding to the defined dimensions range of DLCS. As a result of the DLCS dimensions optimization, the optimal value of S/C ratio (B/L) is 0.106. This value is corresponding to the higher-level sound attenuation of TL value.

Conclusions
This work studies the solutions of drawbacks of experimental and numerical works to study an acoustic behavior of different mufflers' materials to improve its performance, where it is necessary to perform large, important, time-consuming calculations, particularly if the muffler is made from advanced materials such as composite materials. All of these problems were resolved with modern methods in this study by using a deep learning algorithm to predict the acoustic behavior of the muffler material in order to save effort and time on muffler design optimization. The acoustic behavior of DCLCM made from a BFRP laminated composite was predicted in Python through two types of DNN Figure 19. TL of DCLCM with optimal value of S/C ratio (B/L), compared with the system without muffler B/L = 0. Figure 19 shows the relationship between the values of TL and frequency range of DCLCM with optimal value of the S/C ratio (B/L), compared with the system without a muffler. As shown in Figure 19, the maximum value of TL can be obtained at the DCLCM resonance frequency (612 Hz), which is estimated analytically. Therefore, the maximum attenuation value will be 60 dB. It is possible to conclude that the optimization method is the better way to obtain the maximum value of TL corresponding to the defined dimensions range of DLCS. As a result of the DLCS dimensions optimization, the optimal value of S/C ratio (B/L) is 0.106. This value is corresponding to the higher-level sound attenuation of TL value.

Conclusions
This work studies the solutions of drawbacks of experimental and numerical works to study an acoustic behavior of different mufflers' materials to improve its performance, where it is necessary to perform large, important, time-consuming calculations, particularly if the muffler is made from advanced materials such as composite materials. All of these problems were resolved with modern methods in this study by using a deep learning algorithm to predict the acoustic behavior of the muffler material in order to save effort and time on muffler design optimization. The acoustic behavior of DCLCM made from a BFRP laminated composite was predicted in Python through two types of DNN architec-tures which are RNN-LSTM and CNN. A comparison between these two architectures is presented in terms of the speed of training and accuracy of predictive data. We found that both the CNN and RNN-LSTM models achieved accuracy above 90% on test and validation datasets. In addition, the performance of the CNN model is better than the RNN-LSTM model in acoustic monitoring ability. The acoustic parameters of a transmission loss (TL) and the power transmission coefficient (PTC) were computed by solving the exact solution of the governing acoustic equations of the muffler model in MATLAB. The model training parameters are optimized by using Bayesian genetic algorithms (BGA) optimization. The acoustic TL output of the CNN model defined that the optimum muffler length is 0.0212 m with a fixed pipe radius and encountered the maximum acoustic TL of 60 dB in the resonance frequency of DCLCM (612 Hz). This work will reinforce the mufflers' industrials in the future, and its design may one day be equipped with deep learning based algorithms.

Future Work
It will be interesting in the future to apply the proposed methodology in this work with some special boundary conditions and change the composite laminate stacking and filer angle or to add some special structures such as honeycomb materials between the consecutive stacking. Furthermore, a realistic experiment in the standing wave tube may make the results more convincing and useful. I hope that the study can not only give a basic reference to engineering applications but also inspire the researchers who are interested in composite laminates' acoustic properties.