Next Article in Journal
Efficiency-Enhancing Horizontal Mergers in Spatial Competition with Network Externalities
Previous Article in Journal
Feed-Forward Neural Networks Training with Hybrid Taguchi Vortex Search Algorithm for Transmission Line Fault Classification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heuristic Approaches Based on Modified Three-Parameter Model for Inverse Acoustic Characterisation of Sintered Metal Fibre Materials

1
School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China
2
Unmanned System Research Institute, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2022, 10(18), 3264; https://doi.org/10.3390/math10183264
Submission received: 18 August 2022 / Revised: 2 September 2022 / Accepted: 5 September 2022 / Published: 8 September 2022
(This article belongs to the Section Engineering Mathematics)

Abstract

:
Modelling of sound propagation in porous media generally requires the knowledge of several transport properties of the materials. In this study, a three-parameter analytical model that links microstructure properties of sintered metal fibre materials and non-acoustical parameters of the JCAL model is used and modified, and two heuristic approaches based on the established model for inverse acoustic characterisation of fibrous metal felts are developed. The geometric microstructure of sintered fibrous metals is simplified to derive the relationship between pores and fibre diameters. The new set of transport parameters in the modified three-parameter model can cover two controllable parameters during the fabrication process of fibrous metals. With two known transport parameters, six sintered specimens are characterised using a deterministic algorithm, and a satisfactory result is achieved in fitting the normalised surface impedance measured by an acoustic measurement system. Moreover, the forward evaluation shows that our modified three-parameter theoretical model is capable of yielding accurate results for the sintered metal fibre materials. A numerical investigation of the complete inverse acoustic characterisation of fibrous metals by a global non-deterministic algorithm indicates that inversion from two porous material properties is preferable to the normalised surface impedance.

1. Introduction

Sintered metal fibre felts, typically fabricated using micron-sized FeCrAl fibres, are a new engineering material. Compared with conventional fibrous materials, e.g., glass fibre and polyester, sintered metal fibre materials possess not only excellent sound absorption ability, but also exhibit many functional properties: good mechanical properties, high thermal stability, easy processability, and long-life, etc. [1,2,3]. Therefore, sintered metal fibre materials are particularly suitable for applying noise reduction and vibration control under extreme conditions, such as high pressure and temperature [4].
In order to give a physical basis to the description of the macroscopic acoustic behaviour of porous media, a variety of semi-phenomenological models involving several physical parameters have been presented, such as the Delany–Bazley–Miki model (one parameter) [5,6], the Johnson–Champoux–Allard (JCA) model (five parameters) [7,8,9], the Johnson–Champoux–Allard–Lafarge (JCAL) model (six parameters) [7,8,10], the Johnson–Champoux–Allard–Pride–Lafarge (JCAPL) model (eight parameters) [7,8,10,11] and the Biot’s theory [12,13]. Apart from the direct use of these models to predict sound propagation in porous media [14,15,16,17,18,19,20], the inverse acoustic characterisation from the acoustical properties of materials is also of great interest for its practical values in engineering application [21], e.g., inverse design of fibrous metals for achieving high absorption at desired frequencies. Hence, the modelling of sound propagation in sintered metal fibre materials requires the known set of several transport properties, such as porosity, permeability, and characteristic length [22]. To estimate these properties, many methods for characterisation of porous media have been developed, three main of which include: direct, indirect, and inverse methods. Direct methods are based on various experiments, which can directly obtain a set of properties, e.g., porosity, static air flow resistivity, tortuosity, two characteristic lengths, and two permeabilities [23,24,25,26]. Such approaches provide the most accurate and reliable measured parameters among the three methods while suffering the disadvantage of the huge cost and the complex experiment procedures. For indirect methods, based on equivalent fluid models assume, the equivalent density and equivalent bulk modulus of porous media obtained from several different methods of impedance tube measurement are used to analyse the transport properties of porous materials [27,28,29,30].
As an alternative solution to another two methods for obtaining the physical parameters from materials, inverse methods transform the physical problem into a mathematical optimisation problem, which directly estimates transport parameters from acoustical properties measured by impedance tube. Previous existing studies on the inverse characterisation of porous media can be divided into two categories: deterministic or non-deterministic parametric inversion methods and statistical inversion methods. The first relies on local or global optimisation algorithms in fitting the curves of acoustical properties from impedance tube or ultrasonic measurements Atalla and Panneton [31] proposed an inverse characterisation method based on normalised surface impedance and yielded three of five JCA transport parameters by using a differential evolution algorithm Dauchez and Yvars [32] employed the constraint satisfaction problem approach to an inverse procedure for recovering all five parameters of the JCA model from the density and bulk modulus of the equivalent fluid. Furthermore, Dossi et al. [33] proposed an inverse calculation of the elastic parameters of polyurethane foams by fitting measured specific acoustic impedance from samples with different thicknesses. The statistical inversion methods, capable of taking into account knowledge of uncertainty of measurement process and prior information, have been developed to estimate uncertainty for the pore parameter values Chazot et al. [34] introduced a Bayesian approach for evaluating eight properties of Biot’s model with a measurement in a standing wave tube and returning probabilistic data such as the confidence interval of all parameters Niskanen et al. [35] and Cuenca et al. [36] proposed a deterministic and statistical framework for inverse characterisation, in which the minimum obtained by the least-squares method was validated using the information given by the statistical inversion and uncertainties in the estimations could be quantified.
Most existing works on the inverse calculation of transport parameters employed the JCA/JCAL or Biot-JCA model [12,13]. The first two are based on equivalent fluid models assume, while the latter one considers the solid phase of porous media as an elastic material. For sintered metal fibre felts, the frame of this material can be assumed motionless in all frequency ranges due to its quite low decoupling frequency and barely visible resonance behaviour. Thus, the equivalent fluid description for fibrous metals is employed in this study. Some researchers suggested that the classical parametrisation of the acoustical behaviour of certain types of porous media is somewhat redundant [36,37]. Additionally, all the above-mentioned methods can suffer from the accuracy of the experiment and the ill condition of the inverse problem, respectively.
For that purpose, in this study, we develop two heuristic approaches using a modified three-parameter analytical model to investigate the inverse acoustic characterisation of sintered metal fibre materials. First, the three-parameter analytical model for the non-acoustical properties based on the JCAL model is provided and modified based on the average pore size approximation according to the microstructure characteristics of the fibrous metals. Then, the normalised surface impedance of six samples is measured by a classic two-microphone impedance tube measurement. The proposed modified model is validated through the curve fitting with a single parameter estimation of measured data using a local optimisation algorithm. The theoretical model established is employed to examine the difference between the values recovered from surface impedance and those from two porous material properties with a global optimisation algorithm in the numerical investigation.
This paper is organised as follows. First, Section 2 presents the three-parameter analytical model and derives the geometric relationship between pores and fibre diameters according to the simplified microstructure of sintered fibrous metals. The measurement set-up based on a two-microphone configuration is outlined in Section 3. Section 4 describes the observation model of measurements along with the two representative heuristic approaches. Next, Section 5 presents the verification of the theoretical model and numerical investigation with two heuristic methods. Section 6 concludes this study and discusses related future work aspects.

2. Acoustic Modelling of the Sintered Fibrous Metals

2.1. Three-Parameter Analytical Model

The well-established Johnson–Champoux–Allard–Lafarge (JCAL) model for sound propagation in porous media proposed by Champoux, Allard, and Lafarge relies on six transport parameters: porosity ϕ , flow resistivity σ , tortuosity α , viscous characteristic length Λ , thermal characteristic length Λ and static thermal permeability k 0 . It has been shown that the popular six-parameter model can be simplified, and the use of three parameters in the JCAL model is enough for some types of porous media. In our study, a three-parameter analytical model based on the JCAL model proposed by Horoshenkov et al. [37] is adopted. The choice of using this model instead of others is motivated by the following three reasons:
1.
This theoretical model links three pore parameters of porous media to six non-acoustical parameters in the JCAL model, which lay the groundwork for the study of using customisable morphological parameters of real fibrous metals replacement to the transport parameters that are difficult to measure;
2.
The JCAL model, which introduces corrections to the bulk modulus thermal behaviour at low frequencies that is not captured by the JCA model, its robustness has been validated by many studies.
3.
Halving the number of variables can significantly reduce the complexity of the inverse characterisation problem, thus reducing uncertainty and consumption of time in computation during the inversion procedure.
As previously discussed, with the skeleton of sintered metal fibre materials being far heavier than air, the frame can be almost motionless for large ranges of acoustical frequencies. In other words, fibrous metals can be assumed as an equivalent fluid with an equivalent effective density ρ e q ( ω ) and equivalent bulk modulus K e q ( ω ) used to describe the macroscopic acoustical behaviour of the material. Based on the six transport parameters presented previously, the expressions of ρ e q ( ω ) and K e q used for the JCAL model are given by [7,8,10]:
ρ e q ( ω ) = α ρ 0 ϕ 1 + σ ϕ j ω ρ 0 α 1 + j 4 α 2 η ρ 0 ω σ 2 Λ 2 ϕ 2
K e q ( ω ) = γ P 0 / ϕ γ ( γ 1 ) 1 j ϕ κ k 0 C p ρ 0 ω 1 + j 4 k 0 2 C p ρ 0 ω κ Λ 2 ϕ 2 1
where j = 1 , ρ 0 is the density of air, α is the tortuosity, ϕ is the porosity, σ is the flow resistivity, η is the dynamic viscosity of the fluid, γ is the specific heat ratio, k 0 is the static thermal permeability, Λ is the viscous characteristic length, Λ is the thermal characteristic length, P 0 is the quiescent pressure, κ is the thermal conductivity, C p is the specific heat.
The shape of pores in the majority of realistic porous materials is variable, and the size normally obeys a particular statistical distribution. For sintered metal fibre material, based on the prior assumptions of log-normal distribution of pore size, the five of six transport parameters in Equations (1) and (2) can be expressed via three pore parameters using the three-parameter model: the porosity ϕ , average pore size s ¯ and standard deviation in the pore size σ s , which are given, respectively, as follows [37]:
Λ = s ¯ e 5 / 2 σ s log 2 2 ,
Λ = s ¯ e 3 / 2 σ s log 2 2 ,
k 0 = s ¯ 2 ϕ 8 α e 6 σ s log 2 2 ,
α = e 4 σ s log 2 2 ,
σ = 8 η α s ¯ 2 ϕ e 6 σ s log 2 2 .
Under the assumption of an equivalent fluid model, the specific characteristic impedance Z c and the complex wave number k c of fibrous metal materials can be calculated as [22]:
Z c = ρ e q ( w ) K e q ( w ) ρ 0 c 0 ,
k c = w ρ e q ( w ) K e q ( w ) .
where c 0 is the sound velocity in the air.

2.2. Average Pore Size Approximation

As a kind of fibrous material, the commonly used non-acoustical parameters to characterise the sintered metal fibre felts are porosity ϕ and fibre diameter d. These two parameters can be precisely controlled during the fabrication process of fibrous metals. This section aims to link porosity ϕ and fibre diameter d to the average pore size s ¯ that is not a controllable parameter during the manufacturing process. The purpose is to obtain the quantitative relation among three geometrical parameters.
The 3D digital model of a fibrous metal felt reconstructed by Micro-computed tomography (micro-CT) [38] is displayed in Figure 1. As shown in Figure 1a, the composition of the real fibrous metals is twisted metal fibres which are connected by the sintered points after the sintering process. The front view of the slice image in the red square plane is illustrated in Figure 1a is provided in Figure 2a. A thin slice thickness along the x axis is made, thus, the fibre orientation can be observed intuitively in this 3D slice image. Fibres are randomly positioned and horizontal obliquely to the image plane. Combined with Figure 1b,c and Figure 2a, it can be seen that the intersecting fibres lie in the x-y plane parallel to the surface of fibrous metal felts. In order to build relations between geometrical parameters, a cell method proposed by Tarnow [39], Tarnow [40] is used. A Voronoi-based segmentation is performed on the 2D slice image after the binary process, as shown in Figure 2b. Each oval fibre cross-section is surrounded by a so-called Voronoi polygon. It is constructed by drawing perpendicular bisectors passing through the midpoints between each centre and the nearest neighbour centres. These lines intersect and form the Voronoi polygons. The distance between one centre and the nearest neighbour centre can be regarded as the pore size. It has been demonstrated that any Voronoi polygon is necessarily convex. Since sound normally incidents on the surface of felts, we consider a simplified model that the fibres are all parallel and regularly placed and in the same diameter, as shown schematically in Figure 2c, the sound wave incident perpendicularly passes all cylindrical fibres, i.e., parallel to the z axis. In this case, the irregular Voronoi polygons displayed in Figure 2b can be approximated to regular hexagons. The distance between the centre of one cylindrical fibre and any nearest neighbour fibre is equal to the average pore size s ¯ . For simplicity, each ideal hexagonal Voronoi region is approximated by an inscribed circle having the approximate area. The porosity of one unit cell in the parallel cylindrical fibres array is identical to that of the macroscopical homogeneous sintered metal fibrous material. Then the equivalent porosity of the randomly distributed fibres can be given by
ϕ = π ( s ¯ / 2 ) 2 π ( d / 2 ) 2 π ( s ¯ / 2 ) 2 ,
where the porosity ϕ and the fibre diameter d are assumed known, thus the average pore size s ¯ can be eventually represented as
s ¯ = d 1 ϕ .
By substituting Equation (11) into Equations (3)–(7), the transport parameters of the JCAL model can be calculated by a new set of parameters in modified three-parameter model, i.e., the porosity ϕ , fibre diameter d and standard deviation in the pore size σ s , as follows:
Λ = d e 5 / 2 σ s log 2 2 1 ϕ ,
Λ = d e 3 / 2 σ s log 2 2 1 ϕ ,
k 0 = d 2 ϕ 8 α ( 1 ϕ ) e 6 σ s log 2 2 ,
α = e 4 σ s log 2 2 ,
σ = 8 η α ( 1 ϕ ) d 2 ϕ e 6 σ s log 2 2 .

3. Measurement Configuration

The impedance tube is an acoustic measurement system for determining the sound absorption coefficient, the sound reflection coefficient, and the surface impedance of materials in the laboratory on the basis of the transfer function method according to ISO 10534-2 Standard [41].
For testing the surface impedance of our prepared fibrous metal felt samples, a custom-built impedance tube with a 50 m m × 50 m m square cross-section is designed. The data-acquisition module (m+p VibPilot) system includes two microphones (G.R.A.S. 46BD) inserted into a hole in the tube. As shown in Figure 3a, s is the distance between the microphones, and l is the distance between the surface of the felt being tested and the nearest microphone, which are equal to 50 m m . The measurements are made in the frequency range of 200 Hz to 3150 Hz .
Overall, the sound pressure in the tube produced by the loudspeaker is measured by the microphones, and the sound reflection coefficient R is calculated via the transfer function method. Based on the obtained reflection coefficient R, the sound absorption coefficient α and surface impedance Z s can be calculated as [22]:
α = 1 | R | 2 ,
Z s = 1 + R 1 R Z 0 .
where Z 0 = ρ 0 c 0 is the impedance of air.
With the rigid frame assumption, the characteristic impedance Z c and the complex wave number k c are the fundamental acoustic properties of the rigid backing fibrous materials that can completely characterise their acoustical behaviour. In this case, the measured surface impedance Z s can also be expressed by Z c and k c as [22]:
Z s = j Z c cot ( k c t ) ,
where t is the test sample thickness corresponding to the measured surface impedance Z s .

4. Inverse Characterisation

4.1. Observation Model

The observational model is an error model representing the additive random noise and the uncertainty of the measurements using the impedance tube. In general, the experimental observational model can be described as:
y = h ( θ , ξ ) + ϵ ,
where y = [ Z s m e a s ] or [ Z c m e a s , k c m e a s ] corresponding to the measurement procedures, h ( θ , ξ ) is the modified three-parameter model based on the JCAL model derived in Section 2, θ is a set of unknown transport parameters [ ϕ , d , σ s ] , ξ consists of the known environmental conditions (e.g., specific heat ratio γ , quiescent pressure P 0 ) and the measured thickness of felt samples t, and ϵ is the unknown factor may result from media inhomogeneity, measurement error and so on. Here, we consider the white Gaussian noise from the electric signal as the only source of the unknown factor ϵ .

4.2. Heuristic Approaches

For the purpose of estimating θ from a given y, two representative heuristic algorithms are performed. Based on the number of unknown transport parameters, inverse problems can be grouped into two categories: single-parameter and multi-parameter inversion.
In general, the porosity ϕ and fibre diameter d are the given parameters for one sintered metal fibre felt. For our modified three-parameter model, this means that only one parameter is unknown, which is the standard deviation in the pore size σ s . Therefore, the minimisation problem reduces to a single parameter estimation, which can be solved by a local optimisation algorithm.
As a non-gradient method, the Nelder–-Mead Method has been successfully implemented in some similar nonlinear optimisation schemes [42]. Here, the standard deviation σ s is recovered from the surface impedance Z s . The objective function is formed as:
min θ F ( ω ; θ ) = min θ 1 N f f = f l f u Z s m e a s ( ω ) Z s m o d e l ( ω ; σ s )
where ω is the angular frequency, N f is the number of frequency points, f l is the lower frequency, and f u is the upper frequency.
During the inversion, the frequency range f [ 500 , 2500 ] at intervals of 50 Hz , the standard deviation σ s [ 0 , 1 ] , the starting point is set to 0.5, and the maximum number of iterations is set to 500.
In the case of totally unknown parameters of one fibrous metal sample, a global optimisation algorithm for inversion is proposed. The genetic algorithm (GA) is a non-gradient method inspired by the process of natural selection, which is the class of evolutionary algorithms [43]. Unlike the Nelder–-Mead simplex algorithm, which can only find a local optimum close to the starting point, however, GA is capable of finding generally good global solutions, and no starting point is required, which makes GA easier to implement.
In this case, to study the effect of the choice of initial acoustical properties for inversion, [ ϕ , d , σ s ] are recovered from the surface impedance [ Z s m e a s ] and the fundamental acoustic properties [ Z c m e a s , k c m e a s ] , respectively.
The flow of GA is illustrated in Figure 4. Three transport parameters [ ϕ , d , σ s ] (summarised in Table 1) are encoded as chromosomes and optimised by GA. The GA optimisation process starts from a population of individuals randomly generated in constraints, which is called initialisation. The next step is to generate another two generation populations of solutions through crossover and mutation, respectively. The crossover probability is basically set to 0.95, and the mutation probability is 0.1 for two inversion problems. To avoid producing a solution that is not physically reasonable, constraints are also employed in the process of crossover and mutation. A portion of the existing population is selected by rating the fitness of each solution to breed the next generation. Elitist selection is introduced here to guarantee that the solution quality will not decrease from one generation to the next [44]. The fitness function for recovering [ ϕ , d , σ s ] from [ Z s m e a s ] is analogous to Equation (21), which can be given by
min θ F ( ω ; θ ) = min θ 1 N f f = f l f u Z s m e a s ( ω ) Z s m o d e l ( ω ; σ s ) Z s m e a s ( ω )
for inversion from [ Z c m e a s , k c m e a s ] , the fitness function can be written as:
min θ F ( ω ; θ ) = min θ 1 N f f = f l f u Z c m e a s ( ω ) Z c m o d e l ( ω ; ϕ , d , σ s ) Z c m e a s ( ω ) + k c m e a s ( ω ) k c m o d e l ( ω ; ϕ , d , σ s ) k c m e a s ( ω )
Considering the characteristic impedance and the complex wave number may have different contributions to the fitness function, the fitness function is formed as a sum of the error percentage of Z c and k c to make sure two complex properties can be minimised simultaneously and equally.
To increase the chance of searching the global minimum, the initial population size and the maximum number of iterations are both set to 200, and each inversion procedure will run several times due to the initial population of GA being randomly generated.

5. Results and Discussion

5.1. Theoretical Model Verification

The proposed theoretical model for sintered metal fibre materials presented in Section 2 has been validated with a set of independent experiment data. Six fibrous metal samples with 50 m m × 50 m m square cross-section are used to test the sound absorption property using the acoustic measurement system performed in Section 3. The porosity ϕ and fibre diameter d of each sintered specimen are listed in Table 2. The thickness is averaged from five sample points located around and in the centre of the specimen by using a spiral micrometer. These six sintered specimens, generally, can be grouped into two categories: samples 1–3 with gradient porosity and sample 4–6 with gradient thickness. Furthermore, the only variable between samples 2 and 5 is the fibre diameter if the slight difference in thickness is neglected. Note that samples 4–6 are cut from one complete fibrous metal material.
Based on the prior assumption of the log-normal distribution of pore size, the standard deviation in the pore size σ s is the only unknown parameter that is difficult to measure directly. For this sake, the Nelder–-Mead simplex algorithm stated in Section 4.2 is performed to fit the measured surface impedance Z s curves of six felt samples. The fitting of experimental measurements by minimising the objective function (Equation (21)) is illustrated in Figure 5. Following the inverse characterisation results of six felt samples, it can be seen that the real and imaginary parts of the measured Z s are generally fitted well by the algorithm using the theoretical model. The optimal parameter σ s recovered from the inversion process for each sample is summarised in Table 3.
In order to avoid the possibility of mistaking a local minimum for a global minimum by the algorithm, the evolution of the objective function minima in the function of the standard deviation σ s related to all samples is displayed in Figure 6. Visually, only one minimum exists in the objective function for σ s values ranging from 0 to 1. Therefore, the transport parameter σ s recovered from experimental Z s by the Nelder–-Mead simplex algorithm can be confirmed as a global minimum.
As previously derived in the theoretical model, the five of six transport parameters in the JCAL model can be expressed via three pore parameters, using Equations (12)–(16). The complete set of non-acoustical parameters in the JCAL model for these samples, calculated by two known physical parameters [ ϕ , d ] and the identified parameter σ s based on the modified three-parameter model, are summarised in Table 4. The estimated values of five parameters seem sensible in their ranges, respectively, compared with those of fibrous felts reported in other literature. By comparing the transport parameters of each sintered felt sample in Table 3 and Table 4, the standard deviation σ s of samples 1–3 increases gradually with the decrease in the porosity. For the comparison between samples 2 and 5, fibre diameter d is inversely proportional to the standard deviation σ s . As for samples 4–6, a good coincidence between samples 4 and 5 is observed, while sample 6 has apparent discrepancies with non-acoustical parameters in the other two samples. This may attribute to inhomogeneity across the complete sintered metal fibre materials incurred by manufacturing defects. In summary, the value of standard deviation σ s range from 0.4 to 0.5 and be insensitive to the changes of the other two parameters. It is noted that the tortuosity of a fibrous metal (defined in Equation (7)) with the log-normal distributed pore size depends only on the standard deviation σ s .
Moreover, a forward evaluation for the sound absorption coefficient using the JCAL model, which is calculated by six non-acoustical parameter [ ϕ , σ , α , Λ , Λ , k 0 ] listed in Table 4, are compared with the measured data, as shown in Figure 7. The predicted sound absorption coefficient matches well with the measured curve for each felt sample.
Since the fit and predicted curves are both in good agreement with the measured normalised surface impedance and sound absorption coefficient, respectively, the prior assumption on the log-normal distribution of pore size in these fibrous metal felts is verified. Overall, results have shown that our modified three-parameter theoretical model is capable of yielding accurate results for the sintered metal fibre materials.

5.2. Numerical Investigation

In a further analysis, the theoretical acoustical data generated from the established modified model is used to study the difference between two inverse problems: recovering the transport parameters from surface impedance Z s and two porous material properties [ Z c , k c ] , respectively. According to the Equations (8), (9) and (19), surface impedance Z s can be expressed by two porous material properties [ Z c , k c ] with the known material thickness. For the purpose of analysing the complexity of two inverse problems, numerical parameters are adopted. In order to take both generality and particularity into consideration, three instead of one set of randomly generated transport parameters (listed in Table 5) and the corresponding acoustical data are selected. The GA presented in Section 4.2 is applied for characterising these three generated samples owing to its power for multi-parameter identification.
Due to the nature of the non-deterministic algorithm, the GA is continuously run three times with the initial population that is randomly generated within the range values in Table 1. The optimal parameters obtained in each run for two inverse problems are summarised in Table 6, and the relative errors between the theoretical and optimal parameters are given together for intuitive understanding. The best fitness values of each run are also provided as a reference. The iterative process of GA takes approximately 2.8   s on a Intel ® Core TM i7-9700 processor for each run.
As shown in Table 6, the same transport parameter obtained in any two trials for the same input parameters differs from each other. Meanwhile, the variation trend of the error relative to the theoretical value is also irregular. Actually, as a typical population-based global search meta-heuristic, the GA is capable of identifying near-optimal solutions among the local optimums but cannot guarantee that a global optimal solution can be found. In comparison, local search-based strategies such as the Newton–-Raphson Method and Nelder–-Mead simplex algorithm (presented in Section 4.2) are dependent on initial guesses and easily plunge into the local optimum. This is an important factor in the choice of GA to compare the complexity of two inverse problems.
For ease of comparison, the mean relative errors of transport parameters identified from two input parameters are illustrated in Figure 8 for three generated samples, respectively. In terms of comparing the effect of input parameters upon the inversion of transport parameters in one generated sample, inversion from [ Z c , k c ] shows lower errors than the Z s . For these two inverse problems, variations in the input parameters lead to differences in the fitness function. The best fitness value in each trial is very close to the optimal solution ( F ( ω ; θ ) = 0 ) while always being unequal at our limited attempts, which means the fitness functions may have multiple near-optimal solutions that are sufficiently close to the optimum. The optimisation algorithm seems to be easier to get stuck at a locally optimal value in inverse characterisation from Z s than from [ Z c , k c ] . The nearest optimal solution searched by the GA appeared at trial #1 in inversion from [ Z c , k c ] of sample A, whose porosity with 0.001% error can be regarded as an exact value. Meanwhile, it is noted that the relative error of porosity is found to be generally lower than the other two parameters in comparing each transport parameter. Additionally, the identified results recovered from two input parameters of sample C are visibly inferior to the other two samples, especially for the standard deviation σ s .
For further investigating the causes of error in searching the global minima, sample C is selected as a typical one to examine the difference in characteristics of fitness function between using Z s and [ Z c , k c ] as input parameters to determine three transport parameters [ ϕ , d , σ s ] . As illustrated in Figure 9, Figure 10 and Figure 11, the three-dimensional (3D) profiles of the fitness function are performed by fixing one of three parameters and varying the other two. Due to the fact that the objective and near-optimal solutions of functions distribute in a narrow region close to 0, which could hardly be identified on a 3D surface, the logarithmic transformation is performed for two fitness functions to amplify the differences in number and characteristics of local minimums. Moreover, a two-dimensional (2D) cross section is plotted at L o g ( F i t n e s s ) = 2 , i.e., fitness value equals to 0.01. The section derived from this plane intersecting the curved surface is the so-called contour, which is plotted on the 2D projection plane of function. According to the Equations (22) and (23), fitness functions of GA are measured by the normalised mean error (NME). Thus, the 3D surface under the drawn 2D cross-section is composed of the solutions with the NME of ± 1 % . Correspondingly, the local and global minimums are located in the area composed of axes and contour lines on the 2D projection plane.
Considering the mean relative error of identified porosity in sample C is significantly lower than the other two parameters, a detailed analysis is performed with a theoretical value and two fixed porosity values selected from ( 0.70 , 0.88 ) obtained in six trials, as shown in Figure 9. The range values of fibre diameter are also narrowed to ( 5 , 25 ) based on the results of six trials. It can be clearly seen that every sampling surface has at least one valley, which means there are an amount of equal and unequal local optimums for both two fitness functions. The depth of a valley measures the degree of closeness between the local and global minimum. Note that the global minima are both F ( ω ; θ ) = 0 for two inversion problems, which means the logarithm of global minima actually goes to negative infinity as sampling space increases. A combination of prominence and the number of valleys directly leads to the difficulty in achieving global optimisation. The lots of comb-like valleys in Figure 9a–c and e would lead the search algorithm to converge to a wrong optimal point and hard to jump from the local optimal position. Conversely, global optimisation techniques could generally avoid the smooth valleys displayed in Figure 9d,f.
Visually, fitness function using Z s as input parameters shows a complex search space which has numerous minimum “traps” in Figure 9c,e, while function with [ Z c , k c ] shows a relatively smooth search space with fewer minimum “traps” in Figure 9d,f. This explains the reason why the relative errors of the identified porosity recovered from Z s are apparently higher than those from [ Z c , k c ] in general, even though part of the identified parameters inversion from Z s is a little better than those from [ Z c , k c ] .
The surfaces presented in Figure 9a,b show the common feature that both have lots of comb-like valleys. In other words, there are multiple sets of near-optimal solutions [ d , σ s ] around the theoretical value of porosity, which interprets the cause of significant errors in the fibre diameter d and standard deviation σ s recovered from both two inversion problem.
Further, another two 3D profiles of the fitness function are also provided by letting d and σ s to be constant, respectively, as shown in Figure 10 and Figure 11. In comparison of (a) and (b) in these two figures, overall, the valleys in the fitness function using [ Z c , k c ] as input parameters have fewer numbers and more concentrated distribution than those in the function with Z s . Especially, only one minimum, i.e., optimal solution, can be found in Figure 11b. This in turn explains the intrinsic cause of the larger error in searching the standard deviation. Overall, one can conclude from the comparison of Figure 9, Figure 10 and Figure 11 horizontally that the fixing of one value in [ d , σ s ] can largely improve the optimisation quality and efficiency, or the difficulty in determining fibre diameter and standard deviation from Z s or [ Z c , k c ] is obviously higher than in determining porosity. Actually, global optimisation of the fitness function is processed on all three parameters simultaneously, which means the search space is far more complicated than the 3D surface shown in Figure 9, Figure 10 and Figure 11. The combination of these factors results in the complexity of inverse characterisation from acoustical data and deviations of objective variables to a different degree. In general, the multiple execution results of GA and the detailed investigation of sample C for two inverse problems both indicate that inversion from two porous material properties [ Z c , k c ] is preferable to surface impedance Z s .

6. Conclusions

In this study, two heuristic approaches are, respectively, performed to validate the proposed modified three-parameter model by impedance tube measurement and investigate the effects of two input acoustical parameters on the inverse characterisation of sintered metal fibre materials. A new set of transport parameters in the modified three-parameter model based on the JCAL model for the acoustical properties of porous media is provided. A satisfactory result is achieved in fitting the normalised surface impedance measured by an acoustic measurement system for six sintered specimens with varying porosity, fibre diameter, and thickness. With the established modified model, a numerical investigation is performed to study the difference between recovering all three transport parameters from surface impedance and from two porous material properties. The several execution results of GA and further analysis of two fitness functions indicate that the search space in inversion from surface impedance is more complex than that from the characteristic impedance and the complex wave number.
Differing from the previous works on direct, indirect, or inverse methods based on the five or six-parameter acoustical model, our modified three-parameter model and inverse characterisation method realise effective and accurate inverse acoustic characterisation of sintered metal fibre materials. The accuracy of our modified three-parameter theoretical model is validated by experiments. Based on the proposed modified three-parameter model, the sintered metal fibre materials with two known customised microstructure parameters can be accurately characterised using a local optimisation algorithm. Complete inverse acoustic characterisation of fibrous metals from the characteristic impedance and the complex wave number is proved to be more accurate than from the normalised surface impedance by a global optimisation algorithm. Our future work will concentrate on developing new methods to improve accuracy for inverse characterising various porous materials based on this modified three-parameter model.

Author Contributions

Conceptualization, B.P. and J.Z.; methodology, B.P.; software, B.P.; validation, B.P., X.S. and T.Z.; formal analysis, B.P.; investigation, B.P.; resources, T.Z. and J.Z.; data curation, B.P.; writing—original draft preparation, B.P.; writing—review and editing, J.Z. and T.Z.; visualization, B.P.; supervision, J.Z., H.X., D.S. and T.Z.; project administration, J.Z., H.X., D.S. and T.Z.; funding acquisition, J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China [grant number 12072277]; and the Natural Science Foundation of China [No. 12072271]; and the Fundamental Research Funds for the Central Universities of China [grant numbers G2022KY0608].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank Ke Li and Tiqiao Xiao and Yanling Xue, of Shanghai Synchrotron Radiation Facility (SSRF) in Shanghai, for helpful discussions on X-ray tomography.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, T.; He, D.; Chen, C.; Zhao, C.; Fang, D.; Wang, X. The multi-functionality of ultra-light porous metals and their applications. Adv. Mech. 2006, 36, 517–535. [Google Scholar]
  2. Huiping, T.; Jianzhong, W.; Qingbo, A.; Hao, Z. Effect of pore structure on performance of porous metal fiber materials. Rare Met. Mater. Eng. 2015, 44, 1821–1826. [Google Scholar] [CrossRef]
  3. Sun, W.; Pan, B.; Song, X.; Xiao, H.; Zhou, J.; Sui, D. A novel sound absorber design of nanofibrous composite porous material. Mater. Des. 2022, 214, 110418. [Google Scholar] [CrossRef]
  4. Sun, F.; Chen, H.; Wu, J.; Feng, K. Sound absorbing characteristics of fibrous metal materials at high temperatures. Appl. Acoust. 2010, 71, 221–235. [Google Scholar] [CrossRef]
  5. Delany, M.; Bazley, E. Acoustical properties of fibrous absorbent materials. Appl. Acoust. 1970, 3, 105–116. [Google Scholar] [CrossRef]
  6. Miki, Y. Acoustical properties of porous materials-Modifications of Delany-Bazley models. J. Acoust. Soc. Jpn. (E) 1990, 11, 19–24. [Google Scholar] [CrossRef]
  7. Johnson, D.L.; Koplik, J.; Dashen, R. Theory of dynamic permeability and tortuosity in fluid-saturated porous media. J. Fluid Mech. 1987, 176, 379–402. [Google Scholar] [CrossRef]
  8. Champoux, Y.; Allard, J.F. Dynamic tortuosity and bulk modulus in air-saturated porous media. J. Appl. Phys. 1991, 70, 1975–1979. [Google Scholar] [CrossRef]
  9. Allard, J.F.; Champoux, Y. New empirical equations for sound propagation in rigid frame fibrous materials. J. Acoust. Soc. Am. 1992, 91, 3346–3353. [Google Scholar] [CrossRef]
  10. Lafarge, D.; Lemarinier, P.; Allard, J.F.; Tarnow, V. Dynamic compressibility of air in porous structures at audible frequencies. J. Acoust. Soc. Am. 1997, 102, 1995–2006. [Google Scholar] [CrossRef]
  11. Pride, S.R.; Morgan, F.D.; Gangi, A.F. Drag forces of porous-medium acoustics. Phys. Rev. B 1993, 47, 4964. [Google Scholar] [CrossRef]
  12. Biot, M.A. Theory of elastic waves in a fluid-saturated porous solid. 1. Low frequency range. J. Acoust. Soc. Am. 1956, 28, 168–178. [Google Scholar] [CrossRef]
  13. Biot, M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range. J. Acoust. Soc. Am. 1956, 28, 179–191. [Google Scholar] [CrossRef]
  14. Yuan, T.; Song, X.; Xu, J.; Pan, B.; Sui, D.; Xiao, H.; Zhou, J. Tunable acoustic composite metasurface based porous material for broadband sound absorption. Compos. Struct. 2022, 116014. [Google Scholar] [CrossRef]
  15. Zarastvand, M.; Asadijafari, M.; Talebitooti, R. Improvement of the low-frequency sound insulation of the poroelastic aerospace constructions considering Pasternak elastic foundation. Aerosp. Sci. Technol. 2021, 112, 106620. [Google Scholar] [CrossRef]
  16. Rahmatnezhad, K.; Zarastvand, M.; Talebitooti, R. Mechanism study and power transmission feature of acoustically stimulated and thermally loaded composite shell structures with double curvature. Compos. Struct. 2021, 276, 114557. [Google Scholar] [CrossRef]
  17. Zarastvand, M.; Ghassabi, M.; Talebitooti, R. Prediction of acoustic wave transmission features of the multilayered plate constructions: A review. J. Sandw. Struct. Mater. 2022, 24, 218–293. [Google Scholar] [CrossRef]
  18. Khayat, M.; Baghlani, A.; Najafgholipour, M. The propagation of uncertainty in the geometrically nonlinear responses of smart sandwich porous cylindrical shells reinforced with graphene platelets. Compos. Struct. 2021, 258, 113209. [Google Scholar] [CrossRef]
  19. Khayat, M.; Baghlani, A.; Najafgholipour, M.A. The effect of uncertainty sources on the dynamic instability of CNT-reinforced porous cylindrical shells integrated with piezoelectric layers under electro-mechanical loadings. Compos. Struct. 2021, 273, 114336. [Google Scholar] [CrossRef]
  20. Khayat, M.; Baghlani, A.; Dehghan, S.M.; Najafgholipour, M.A. Geometrically nonlinear dynamic analysis of functionally graded porous partially fluid-filled cylindrical shells subjected to exponential loads. J. Vib. Control. 2022, 28, 758–772. [Google Scholar] [CrossRef]
  21. Horoshenkov, K.V. A review of acoustical methods for porous material characterisation. Int. J. Acoust. Vib 2017, 22, 92–103. [Google Scholar] [CrossRef]
  22. Allard, J.; Atalla, N. Propagation of Sound in Porous Media: Modelling Sound Absorbing Materials 2e; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  23. Leclaire, P.; Kelders, L.; Lauriks, W.; Melon, M.; Brown, N.; Castagnede, B. Determination of the viscous and thermal characteristic lengths of plastic foams by ultrasonic measurements in helium and air. J. Appl. Phys. 1996, 80, 2009–2012. [Google Scholar] [CrossRef]
  24. Leclaire, P.; Kelders, L.; Lauriks, W.; Glorieux, C.; Thoen, J. Determination of the viscous characteristic length in air-filled porous materials by ultrasonic attenuation measurements. J. Acoust. Soc. Am. 1996, 99, 1944–1948. [Google Scholar] [CrossRef]
  25. Leclaire, P.; Umnova, O.; Horoshenkov, K.; Maillet, L. Porosity measurement by comparison of air volumes. Rev. Sci. Instruments 2003, 74, 1366–1370. [Google Scholar] [CrossRef]
  26. Moussatov, A.; Ayrault, C.; Castagnède, B. Porous material characterization–ultrasonic method for estimation of tortuosity and characteristic length using a barometric chamber. Ultrasonics 2001, 39, 195–202. [Google Scholar] [CrossRef]
  27. Panneton, R.; Olny, X. Acoustical determination of the parameters governing viscous dissipation in porous media. J. Acoust. Soc. Am. 2006, 119, 2027–2040. [Google Scholar] [CrossRef] [PubMed]
  28. Olny, X.; Panneton, R. Acoustical determination of the parameters governing thermal dissipation in porous media. J. Acoust. Soc. Am. 2008, 123, 814–824. [Google Scholar] [CrossRef] [PubMed]
  29. Bonfiglio, P.; Pompoli, F. Inversion problems for determining physical parameters of porous materials: Overview and comparison between different methods. Acta Acust. United Acust. 2013, 99, 341–351. [Google Scholar] [CrossRef]
  30. Groby, J.P.; Ogam, E.; De Ryck, L.; Sebaa, N.; Lauriks, W. Analytical method for the ultrasonic characterization of homogeneous rigid porous materials from transmitted and reflected coefficients. J. Acoust. Soc. Am. 2010, 127, 764–772. [Google Scholar] [CrossRef]
  31. Atalla, Y.; Panneton, R. Inverse Acoustical Characterization of Open Cell Porous Media Using Impedance Tube Measurements. Can. Acoust. Acoust. Can. 2005, 33. [Google Scholar]
  32. Dauchez, N.; Yvars, P.A. Inverse Method for Porous Material Characterization Using the Constraint Satisfaction Problem Approach; Acoustics: Bengaluru, India, 2012. [Google Scholar]
  33. Dossi, M.; Brennan, M.; Moesen, M.; Vandenbroeck, J.; Huo, L. An Inverse Method to Determine Acoustic Parameters of Polyurethane Foams. Inter-Noise -Noise-Con Congr. Conf. Proc. 2019, 259, 1170–1181. [Google Scholar]
  34. Chazot, J.D.; Zhang, E.; Antoni, J. Acoustical and mechanical characterization of poroelastic materials using a Bayesian approach. J. Acoust. Soc. Am. 2012, 131, 4584–4595. [Google Scholar] [CrossRef]
  35. Niskanen, M.; Groby, J.P.; Duclos, A.; Dazel, O.; Le Roux, J.C.; Poulain, N.; Huttunen, T.; Lähivaara, T. Deterministic and Statistical Characterization of Rigid Frame Porous Materials from Impedance Tube Measurements. J. Acoust. Soc. Am. 2017, 142, 2407–2418. [Google Scholar] [CrossRef]
  36. Cuenca, J.; Göransson, P.; De Ryck, L.; Lähivaara, T. Deterministic and Statistical Methods for the Characterisation of Poroelastic Media from Multi-Observation Sound Absorption Measurements. Mech. Syst. Signal Process. 2022, 163, 108186. [Google Scholar] [CrossRef]
  37. Horoshenkov, K.V.; Hurrell, A.; Groby, J.P. A Three-Parameter Analytical Model for the Acoustical Properties of Porous Media. J. Acoust. Soc. Am. 2019, 145, 2512–2517. [Google Scholar] [CrossRef]
  38. Xie, H.-L.; Deng, B.; Du, G.-H.; Fu, Y.-N. Latest advances of X-ray imaging and biomedical applications beamline at SSRF. Nucl. Sci. Tech. 2015, 26, 20102-020102. [Google Scholar]
  39. Tarnow, V. Airflow resistivity of models of fibrous acoustic materials. J. Acoust. Soc. Am. 1996, 100, 3706–3713. [Google Scholar] [CrossRef]
  40. Tarnow, V. Calculation of the dynamic air flow resistivity of fiber materials. J. Acoust. Soc. Am. 1997, 102, 1680–1688. [Google Scholar] [CrossRef]
  41. ISO 10534-2:1998; Acoustics–Determination of Sound Absorption Coefficient And Impedance in Impedance Tubes—Part 2: Transfer-Function Method. International Organization for Standardization: Geneva, Switzerland, 1998.
  42. Nelder, J.A.; Mead, R. A simplex method for function minimization. Comput. J. 1965, 7, 308–313. [Google Scholar] [CrossRef]
  43. Holland, J.H. Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  44. Baluja, S.; Caruana, R. Removing the genetics from the standard genetic algorithm. In Machine Learning Proceedings 1995; Elsevier: Amsterdam, The Netherlands, 1995; pp. 38–46. [Google Scholar]
Figure 1. 3D geometrical model of sintered metal fibre materials based on micro-CT image reconstruction. The red square plane in (a) denotes a slice of the micro CT scan. (b) Top view of the model in x-y plane. (c) Right view of the model in x-z plane.
Figure 1. 3D geometrical model of sintered metal fibre materials based on micro-CT image reconstruction. The red square plane in (a) denotes a slice of the micro CT scan. (b) Top view of the model in x-y plane. (c) Right view of the model in x-z plane.
Mathematics 10 03264 g001
Figure 2. (a) Front view of the 3D micro CT slice image at the position shown in the isometric view of the 3D model. (b) Voronoi diagram of the 2D slice image. (c) Array of parallel cylindrical fibres with ideal hexagonal Voronoi regions.
Figure 2. (a) Front view of the 3D micro CT slice image at the position shown in the isometric view of the 3D model. (b) Voronoi diagram of the 2D slice image. (c) Array of parallel cylindrical fibres with ideal hexagonal Voronoi regions.
Mathematics 10 03264 g002
Figure 3. Acoustic measurement system used for measuring the normalised surface impedance of the prepared sintered specimens. (a) Schematic of the classic two-microphone impedance tube. (b) Realistic facility used in our test.
Figure 3. Acoustic measurement system used for measuring the normalised surface impedance of the prepared sintered specimens. (a) Schematic of the classic two-microphone impedance tube. (b) Realistic facility used in our test.
Mathematics 10 03264 g003
Figure 4. Outline of the genetic algorithm.
Figure 4. Outline of the genetic algorithm.
Mathematics 10 03264 g004
Figure 5. Fitting for the measured normalised surface impedance of six felt samples, (a) Sample 1, (b) Sample 2, (c) Sample 3, (d) Sample 4, (e) Sample 5, (f) Sample 6.
Figure 5. Fitting for the measured normalised surface impedance of six felt samples, (a) Sample 1, (b) Sample 2, (c) Sample 3, (d) Sample 4, (e) Sample 5, (f) Sample 6.
Mathematics 10 03264 g005
Figure 6. Minimum of the objective function for six inversion processes, (a) Sample 1, (b) Sample 2, (c) Sample 3, (d) Sample 4, (e) Sample 5, (f) Sample 6.
Figure 6. Minimum of the objective function for six inversion processes, (a) Sample 1, (b) Sample 2, (c) Sample 3, (d) Sample 4, (e) Sample 5, (f) Sample 6.
Mathematics 10 03264 g006
Figure 7. Comparison of the measured sound absorption coefficient and theoretical one calculated by six non-acoustical parameters based on the JCAL model for six sintered felt samples, (a) Sample 1, (b) Sample 2, (c) Sample 3, (d) Sample 4, (e) Sample 5, (f) Sample 6.
Figure 7. Comparison of the measured sound absorption coefficient and theoretical one calculated by six non-acoustical parameters based on the JCAL model for six sintered felt samples, (a) Sample 1, (b) Sample 2, (c) Sample 3, (d) Sample 4, (e) Sample 5, (f) Sample 6.
Mathematics 10 03264 g007
Figure 8. Comparison of the mean relative error of transport parameters identified from the surface impedance Z s and two acoustic properties [ Z c , k c ] for three generated samples, (a) Sample A, (b) Sample B, (c) Sample A.
Figure 8. Comparison of the mean relative error of transport parameters identified from the surface impedance Z s and two acoustic properties [ Z c , k c ] for three generated samples, (a) Sample A, (b) Sample B, (c) Sample A.
Mathematics 10 03264 g008
Figure 9. 3D profile of the logarithm of fitness function versus the fibre diameter d and standard deviation σ s with the fixed porosity ϕ , and 2D projection along d σ s plane with contour plot at L o g ( F i t n e s s ) = 2 . (a,c,e) Varying different porosity to investigate inverse characterisation from the surface impedance Z s . (b,d,f) Varying different porosity to investigate inverse characterisation from two porous material properties [ Z c , k c ] .
Figure 9. 3D profile of the logarithm of fitness function versus the fibre diameter d and standard deviation σ s with the fixed porosity ϕ , and 2D projection along d σ s plane with contour plot at L o g ( F i t n e s s ) = 2 . (a,c,e) Varying different porosity to investigate inverse characterisation from the surface impedance Z s . (b,d,f) Varying different porosity to investigate inverse characterisation from two porous material properties [ Z c , k c ] .
Mathematics 10 03264 g009
Figure 10. 3D profile of the logarithm of fitness function versus the porosity ϕ and standard deviation σ s with the fixed fibre diameter d, and 2D projection along ϕ σ s plane with contour plot at L o g ( F i t n e s s ) = 2 , (a) d t = 15.07 , inversion from Z s , (b) d t = 15.07 , inversion from ( Z c , k c ) .
Figure 10. 3D profile of the logarithm of fitness function versus the porosity ϕ and standard deviation σ s with the fixed fibre diameter d, and 2D projection along ϕ σ s plane with contour plot at L o g ( F i t n e s s ) = 2 , (a) d t = 15.07 , inversion from Z s , (b) d t = 15.07 , inversion from ( Z c , k c ) .
Mathematics 10 03264 g010
Figure 11. 3D profile of the logarithm of fitness function versus the porosity ϕ and fibre diameter d with the fixed standard deviation σ s , and 2D projection along ϕ d plane with contour plot at L o g ( F i t n e s s ) = 2 , (a) σ s t = 0.21 , inversion from Z s , (b) σ s t = 0.21 , inversion from ( Z c , k c ) .
Figure 11. 3D profile of the logarithm of fitness function versus the porosity ϕ and fibre diameter d with the fixed standard deviation σ s , and 2D projection along ϕ d plane with contour plot at L o g ( F i t n e s s ) = 2 , (a) σ s t = 0.21 , inversion from Z s , (b) σ s t = 0.21 , inversion from ( Z c , k c ) .
Mathematics 10 03264 g011
Table 1. Range values of the transport parameters set in GA.
Table 1. Range values of the transport parameters set in GA.
ϕ d ( μ m ) σ s
( 0.5 , 1 ) ( 5 , 105 ) ( 0.1 , 0.6 )
Table 2. Physical parameters of six sintered metal fibre materials.
Table 2. Physical parameters of six sintered metal fibre materials.
Sample Number ϕ (%)d ( μ m)Thickness (mm)
1952217.77
2902219.26
3852219.52
4904029.54
5904019.86
690409.89
Table 3. Optimal standard deviation obtained from the inverse characterisation on the six felt samples.
Table 3. Optimal standard deviation obtained from the inverse characterisation on the six felt samples.
Sample Number123456
Optimal σ s 0.45020.46970.49610.44860.44310.5112
Table 4. Results of the non-acoustical parameters in the JCAL model for six felt samples.
Table 4. Results of the non-acoustical parameters in the JCAL model for six felt samples.
Sample Number ϕ σ
( N  ·  s  ·  m −4)
α Λ ( μ m ) Λ ( μ m ) k 0 ( μ m 2)
10.9541,954.061.4877.13113.8613041.30
20.9096,559.471.5353.3781.5627060.14
30.85173,335.581.6042.2767.8341938.14
40.9026,613.631.4799.33146.248339.47
50.9025,989.771.4699.92145.728379.13
60.9035,521.691.6592.41152.717871.56
Table 5. Randomly generated physical parameters of fibrous metal felts for numerical investigation.
Table 5. Randomly generated physical parameters of fibrous metal felts for numerical investigation.
Generated Sample ϕ d ( μ m ) σ s Thickness (mm)
A0.907495.57920.163546.1019
B0.653557.48240.407634.5052
C0.697715.07240.207435.1952
Table 6. Transport parameters identified using GA from the surface impedance Z s and two acoustic properties [ Z c , k c ] with relative errors between the theoretical and optimal values. The mean relative errors for each input parameter are given in boldface.
Table 6. Transport parameters identified using GA from the surface impedance Z s and two acoustic properties [ Z c , k c ] with relative errors between the theoretical and optimal values. The mean relative errors for each input parameter are given in boldface.
Generated SampleInput ParametersNumber ϕ ± Error (%)d ( μ m ) ± Error (%) σ s ± Error (%)Fitness Value
A Z s #10.9156 (+0.91)90.9148 (−4.88)0.1654 (+1.14)0.0074
#20.9123 (+0.55)92.8670 (−2.84)0.1653 (+1.11)0.0045
#30.9241 (+1.84)86.2075 (−9.81)0.1692 (+3.50)0.0151
Mean0.9173 (+1.09)89.9964 (−5.84)0.1666 (+1.90)
[ Z c , k c ] #10.9074 (+0.001)95.6981 (+0.12)0.1640 (+0.28) 3.2359 × 10 4
#20.9113 (+0.43)94.3968 (−1.24)0.1761 (+7.69)0.0035
#30.9166 (+1.02)92.5786 (−3.14)0.1912 (+16.93)0.0083
Mean0.9118 (+0.48)94.2245 (−1.42)0.1771 (+8.32)
B Z s #10.6344 (−2.93)60.2355 (+4.79)0.4067 (−0.22)0.0149
#20.6960 (+6.50)54.7718 (−4.72)0.4350 (+6.73)0.0233
#30.6492 (−0.66)63.8648 (+11.10)0.4476 (+9.83)0.0260
Mean0.6599 (+0.97)59.6240 (−3.73)0.4298 (+5.44)
[ Z c , k c ] #10.6545 (+0.14)55.8480 (−2.84)0.3934 (−3.48)0.0120
#20.6508 (−0.41)55.5301 (−3.40)0.3875 (−4.93)0.0130
#30.6483 (−0.80)56.9729 (−0.89)0.3963 (−2.77)0.0065
Mean0.6512 (−0.35)56.1170 (−2.38)0.3924 (−3.73)
C Z s #10.8859 (+26.97)13.4143 (−11.00)0.5449 (+162.71)0.0117
#20.8081 (+15.83)15.8471 (+5.14)0.4688 (+126.00)0.0089
#30.7293 (+4.53)15.2246 (+1.01)0.2976 (+43.48)0.0020
Mean0.8078 (+15.78)14.8287 (−1.62)0.4371 (+110.75)
[ Z c , k c ] #10.7048 (+1.02)19.0428 (+26.34)0.3826 (+84.43)0.0162
#20.7031 (+0.77)16.1080 (+6.87)0.2747 (+32.41)0.0084
#30.7012 (+0.50)13.9150 (−7.68)0.1122 (−45.92)0.0064
Mean0.7030 (+0.76)16.3553 (+8.51)0.2565 (+23.67)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, T.; Pan, B.; Song, X.; Sui, D.; Xiao, H.; Zhou, J. Heuristic Approaches Based on Modified Three-Parameter Model for Inverse Acoustic Characterisation of Sintered Metal Fibre Materials. Mathematics 2022, 10, 3264. https://doi.org/10.3390/math10183264

AMA Style

Zhao T, Pan B, Song X, Sui D, Xiao H, Zhou J. Heuristic Approaches Based on Modified Three-Parameter Model for Inverse Acoustic Characterisation of Sintered Metal Fibre Materials. Mathematics. 2022; 10(18):3264. https://doi.org/10.3390/math10183264

Chicago/Turabian Style

Zhao, Tianfei, Baorui Pan, Xiang Song, Dan Sui, Heye Xiao, and Jie Zhou. 2022. "Heuristic Approaches Based on Modified Three-Parameter Model for Inverse Acoustic Characterisation of Sintered Metal Fibre Materials" Mathematics 10, no. 18: 3264. https://doi.org/10.3390/math10183264

APA Style

Zhao, T., Pan, B., Song, X., Sui, D., Xiao, H., & Zhou, J. (2022). Heuristic Approaches Based on Modified Three-Parameter Model for Inverse Acoustic Characterisation of Sintered Metal Fibre Materials. Mathematics, 10(18), 3264. https://doi.org/10.3390/math10183264

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop