Effect of slit length on linear and non-linear acoustic transfer impedance of a micro-slit plate

– The effect of the slit length on the acoustic transfer impedance of micro-slit plates (MSPs) is investigated in the linear and non-linear regime for a speci ﬁ c slit geometry. This geometry is inspired by slits obtained by cutting and bending the plate. MSPs are plates with arrays of slit-shaped perforations, with the width of the order of the acoustic viscous boundary layer thickness. Impedance tube measurements on two accurately manufactured plates are compared to numerical solution of the Linearized Navier-Stokes equations and to analytical limits. The impedance of the plate is obtained by the impedance of a single slit divided by the plate porosity. The resistance of a slit is independent on the slit length and on the plate porosity. In the linear regime the resistance is accurately predicted by a two-dimensional numerical model. In the non-linear regime, the resistance is strongly dependent on the amplitude of the acoustic waves. The inertance of the slit is weakly dependent on the slit length and on the plate porosity, for low and moderate amplitudes. For high amplitudes, a complicated amplitude dependency of the inertia of short slits is found. One expects that most of the conclusions obtained can be generalised to other slit geometries.


Introduction
Micro-slit absorbers and plates (MSAs, MSPs) have been proposed as sound absorbers at low frequencies, providing lightweight and compact solutions to substitute conventional materials [1].MSPs are plates with slit widths in the sub-millimeter range and low porosity (order of 1%).In conventional designs, micro-slit plates are backed by a cavity forming micro-slit absorbers (MSAs).One of the advantages of slits with respect to circular perforations is that, for equal porosity, a single slit replaces a large number of circular perforations.Furthermore, a slit can be used to delimit flexible structures whose vibration can contribute to the sound absorption [2,3].However, the manufacturing of slits is difficult and can be an obstacle in industrial applications.A possible manufacturing process is to cut the plate, bending the two portions close to the cut, as displayed in Figure 1.A slit is created without removing material from the plate and can lead to new designs.One of the advantages of this geometry is that the edges in contact with the slits are protected from external agents in harsh environments.Another advantage is the possibility to reach sub-millimeter slit widths.This manufacturing technique is used to produce Ò Acustimet plates by Sontech [4,5].In this work, a geometry inspired by the geometry of Figure 1 is studied.Impedance tube measurements are used to investigate the effect of the slit length in two accurately manufactured micro-slit plates.The edges of the slits are kept as sharp as possible.Both plates have the same porosity and total slit perforation length.In the linear regime, experimental results are compared to numerical solutions of the Linearized incompressible Navier-Stokes equations.Microperforated plates (MPPs) and micro-slit plates can be designed to obtain excellent linear acoustic properties but, at high amplitudes, the non-linear effects deteriorate the performance of the absorbers [6][7][8].In practical applications, the acoustic particle velocity in the slits can reach high amplitudes.For this reason, the change of resistance and inertance of the slits due to non-linear effects for long and short slits has been studied in this this work.In literature, several manufacturing techniques are employed to create slits.In classical applications, a slit can be created by removing material from the plate [2,3,9].Slits can also be generated by mating two slotted layers [6,10].Alternative designs of MSAs have been reported in the literature [2,3,6,11,12].Several publications concern the acoustic impedance of MPPs at high amplitudes in presence of a bias flow or for sound-excited flows [13][14][15][16][17].In Section 2, the theoretical background is presented for the linear and nonlinear regimes.In Section 3, the experimental setup and the geometry of the plates are discussed.In Section 4, the twodimensional numerical model is described.In Section 5, the comparison between numerical and experimental results is presented in the linear and non-linear regimes.As explained by Cummings and Eversman [18], non-linear losses at very high amplitudes are due to the formation of a quasi-steady jet flow.This jet displays a contraction after flow separation from the edges of the slit: the so-called vena contracta.In Appendix A, the vena contracta factor for a simplified model of the geometry is calculated.In Appendix B, the quasi-steady incompressible model is presented.In Appendix C, a correction of the quasi-steady model for the viscous boundary layer thickness is discussed.

Theoretical background
In the linear regime, dissipation of acoustic energy takes place in the oscillating boundary layer of thickness , where x = 2pf is the angular frequency, q is the air density (q = 1.18 kg/m 3 at 25 °C and atmospheric pressure) and l is the dynamic viscosity of air (l = 1.85 Â 10 À5 kg/ms at 25 °C).The ratio between the slit width b and the thickness of the viscous boundary layer d v is the Shear number In a micro-slit plate, typical Shear numbers, in the frequency range of interest, are of order unity.In the nonlinear regime, for moderate excitation amplitudes, vortices form locally at the edges of the slits.At very high amplitudes, this leads eventually to the formation of jets.Additional dissipation of acoustic energy is involved.The behaviour of the plate can be studied as a function of the Strouhal number (St b ), defined as the ratio between the slit width b and the amplitude of the oscillating particle displacement at the slits.In formulas, where U p is the cross-sectional surface averaged acoustic velocity amplitude at the slit (u p ðtÞ ¼ Re½û p expðixtÞ ¼ U p cosðxtÞ, for harmonic oscillations).For St b ) 1 (linear regime) the particle displacement is smaller than the slit width and vortices are not formed.For St b ( 1 (strongly non-linear), vortices are formed and they move away from the slit forming a free jet.For St b % 1 (moderate nonlinear) vortices form at the edges of the slits and they remain local.Alternatively, one can define a Strouhal number St t based on the plate thickness at the slit (t = t p À t d in Fig. 3) in order to compare the vortices displacement with the thickness of the plate at the slit [19].

Transfer impedance in linear regime
In the linear regime, the concept of transfer impedance is introduced in the frequency domain (for purely harmonic oscillations) of frequency f.At a distance large compared to the slit width b but small compared to the acoustic wavelength k = c/f, the flow can be described in terms of plane acoustic waves.As this region is compact, the corresponding complex amplitude û of the acoustic velocity uðtÞ ¼ Re½û expðixtÞ ¼ U cosðxtÞ is the same on both sides of the plate.The transfer impedance of the plate is defined as the ratio between the complex acoustic pressure difference Áp and the amplitude of the acoustical velocity û in a cross-section upstream of the plate.The pressure is found by the extrapolation of the plane wave solutions (on both sides) to the sample surface (a formal discussion can be found in Sect.4).The dimensionless transfer impedance of the plate is where q is the density of air and c is the speed of sound in air.Note that the plate transfer impedance is a complex quantity.The dimensionless transfer impedance of the plate z plate is with i 2 = À1, Re[z plate ] the resistive part of the transfer impedance of the plate (or resistance of the plate) and Im[z plate ] the reactive part of the transfer impedance of the plate (or inertance of the plate).The transfer impedance of a slit is defined as the ratio between the complex pressure difference p and the amplitude of the crosssectional acoustical velocity ûp in the slit.In formulas, where ûp ¼ û=U and U the porosity of the plate.Therefore, in the first order of approximation the transfer impedance of the plate can be obtained by the impedance of a single slit.In practical applications, the micro-slit plate is backed by a cavity.As the slit forms the neck of a Helmholtz resonator with a portion of the back cavity as volume, the flow within the slit will be considered as locally incompressible up to the first resonance frequency of the resonator.In the audio range, the square of the Helmholtz number is small, i.e.He 2 ¼ xb c À Á 2 < 10 À1 .Therefore, one can assume a frictionless flow of the incompressible flow.
Thermal effects in the slit are neglected [20].

Non-linear regime
In the non-linear regime, the resistance due to vortex shedding dominates the absorption mechanism, as shown by Ingard and Ising [21].Cummings and Eversmann [18] Figure 1.Representation of a typical geometry obtained with the process of cutting the plate and bending the extremities.assume a quasi-steady flow behaviour to describe the behaviour of perforations at high Shear numbers and very high amplitudes of acoustic particle velocity.The model assumes that the acoustic flow separates at the edges of the slits and forms a free jet with a cross-section smaller than the perforation area.The ratio between the cross-section of the jet and the cross-section of the perforation is called the vena contracta factor C v .In their model, using the Bernoulli equation one can derive the relationship between pressure change across the plate and particle velocity Áp % where Real[z plate,L ] is the plate resistance in the linear regime.Derivation of the theoretical limit can be found in Appendix B. The vena contracta coefficient for the geometry discussed in this study is assumed to be C v = 0.82.This value is found in the potential flow limit for small porosity using the hodographic method [22,23].Derivation of the vena contracta factor C v can be found in Appendix A.
The theoretical limit can be corrected for the effect of the quasi-steady viscous boundary layer.This causes a reduction of the effective porosity that leads to an increase of the resistance.The correction for the thickness of the viscous boundary layer is in Appendix C. For the inertial part of the transfer impedance Im[z plate ], Ingard and Ising assume that in the upstream of the flow separation the flow remains identical to the potential flow prevailing in the linear case [21].In the downstream free jet, the inertia is negligible.Hence, the inertia should be reduced by a factor 2. In other words, DIm½z plate;NL ¼ ÀIm½z plate;L =2.The factor onehalf is explicitly discussed in Morse and Ingard [24].This simple limit will be compared to experimental results in Section 5.2.To study non-linear effects, the linear contribution is subtracted from the non-linear resistance and inertance.In order to analyze the effect on a single slit, the porosity U is introduced.

Impedance tube setup
The experimental setup used in this study is an impedance tube with 6 pre-polarized 1/4 inch microphones (type BWSA, sensitivity 50 mV/Pa).The tube is made of aluminium with an inner diameter D i = 50 mm, a wall thickness t w = 10 mm and length l t = 1000 mm.The excitation system is a 25 W loudspeaker.The six microphones are equally placed at a distance of 175 mm.A relative calibration is performed on the microphones using the microphone closest to the end of the tube (sample side) as the reference microphone.The position of this microphone with respect to the end of the impedance tube is x ref = 47.7 mm.Details on the setup and the calibration system can be found in the works of Temiz et al. [25] and Kojourimanesh et al. [26].The micro-slit plates are positioned at the end of the impedance tube through a sample holder.For this study two sample holders are used to compare the effects of three-dimensional effects for a plate confined by the impedance tube from two sides and from one side (L h 1 ¼ 50 mm and L h 2 ¼ 9 mm).The impedance tube termination with the sample holders is shown in Figure 2.Both the sample holders have a groove for an o-ring to guarantee air-tightness from both the sides of the sample.A script built-in NILabView software controls the signal processing and data acquisition during the measurements.For this study, the sampling rate is 20 kHz for the excitation signal and 10 kHz for recording the input signal.The amplitude of the excitation signal is adjusted automatically until it satisfies the pre-determined pressure value for the reference microphone pðx ref Þ within an accuracy of 2%.This amplitude is also used to derive the acoustic velocity at the sample.The calculation of the reflection coefficient at the sample is based on the plane wave assumption.For the evaluation of the reflection coefficient the method from Jang and Hi [27] is used.For each frequency, every microphone records the complex pressure amplitude pðxÞ at position x, with pþ and pÀ respectively the amplitudes of the wave travelling in the positive and in the negative directions x = 0 corresponds to the end of the impedance tube (sample side), k is the complex wavenumber.Taking viscothermal effects into account as proposed in Peters et al. [28], the complex wavenumber is where i is the imaginary unit, Pr is the Prandtl number (Pr = 0.72), c is the heat capacity ratio (c = 1.4) and Sh D the Shear number based on the impedance tube diameter D i defined as the ratio of the tube diameter Experimentally, the closed pipe termination at x = 0 is used as a reference for the accuracy of the measurements.For an amplitude pðx ref Þ ¼ 2 Pa, the maximum deviation from the theoretical value R = 1.000 is less than 0.3%.The closedend measurements are performed for several excitation amplitudes.It appears that for amplitude p x ref ð Þ > 23 Pa the accuracy of the measurements reduces to 1% up to f = 400 Hz.For f > 400 Hz this accuracy reduces to around 3%.Hence, the study of the amplitude dependence of the measurements is restricted to 0.4 Pa < p x ref ð Þ < 23 Pa.For low frequencies (f < 200 Hz), no significant effects are found increasing the measurement time (number of samples).For f > 700 Hz the results appear to be less reliable.Therefore, the frequency range of the measurements is restricted to 20 Hz < f < 700 Hz corresponding to 1 < Sh < 6.

Measuring the transfer impedance of the sample and acoustic velocity
To measure the transfer impedance of the plate the following procedure is followed.At the impedance tube termination an additional pipe segment of the same length of the sample holder is added.The open pipe termination is located at x open = x 0 +L h + t p , with x 0 = 0 is the right side of the plate in Figure 2, L h the length of the holder and t p the thickness of the sample.The reflection coefficient R o for the open pipe termination is measured.The dimensionless radiation impedance is calculated using The tube is then loaded with the sample by the use of a sample holder.The reflection coefficient of the sample-loaded termination R s is measured.The sample-loaded impedance can be calculated in the same way as the radiation impedance as z s = (1 + R s )/(1 À R s ).The samples have relatively low porosity.Therefore, the radiation impedance is expected to be much lower than the impedance of the plate.Nevertheless, the radiation effects are taken into account and the dimensionless transfer impedance of the plate z plate is calculated as the difference between the sample-loaded impedance and the radiation impedance, in formula z plate = z s À z rad .The radiation of the room is close to that of a free field.Deviations due to unwanted changes in room acoustics are taken into account repeating the open pipe termination experiments before each set of measurements.One does observe some systematic deviation from free-space radiation as a result of room resonances.Assuming the radiation impedance is in series with the plate impedance, the room effect is corrected for by measuring the radiation impedance.In this open pipe radiation impedance measurement, the sample is replaced by a ring in the sample holder so that the geometry (pipe length, position in the room) is exactly the same as when measuring with a sample.The analysis for non-linear studies is performed for f 200 Hz.In this analysis of non-linear effects the amplitude jûj of the flow velocity is a key parameter.The magnitude of the acoustic velocity at the plate is calculated as jûj ¼ pðx ref Þ=jz plate j.Hence, it is assumed that pðx ¼ 0Þ % pðx ref Þ.This is certainly accurate at low frequencies (f < 340 Hz) given x ref % 50 mm.

Specifications of the samples
The acoustical behaviour of short and long slits in micro-slit plates is compared by considering two samples.The geometry of the plates is inspired to the plate in Figure 1.The plates are realized in brass with milling process [29].Sketches of the plates showing the main parameters are shown in Figure 3.The plates are shown in Figure 4.The external diameter of the plates is D p = 70 mm, allowing hosting the plates in the holder.
The effective diameter of the portion of the plate where the slits are located is D i = 50 mm, with D i the internal diameter of the impedance tube.The total plate thickness is t p = 5 mm and the nominal slit width is b = 0.5 mm.The external width of the ditch is w d,e = 5 mm, the internal width of the ditch is w d,i = 2.25 mm.The ditch thickness is t d = 2.75 mm.The thickness of the plate at the slit is t = t p À t d = 2.25 mm.The slit length is l s,short = 15 mm for short slits and l s,long = 35 mm for long slits.The angle between the internal ditch and the outside ditch is 45°.The plates are realized in such a way that the total length of the slits P is the same.The total slit length is P = 7 Â l s, short = 3 Â l s,long % 105 mm.The porosity (of the portion of the plate in the impedance tube) is U = 2.7%.In order to provide access to the rotary cutter, the ditch length is longer than the slit length, l d = l s + 2 mm.The actual widths of the single slits are measured experimentally by means of a digimatic indicator.Maximum deviations from the prescribed dimensions are of the order of 2% (± 1 lm) for the plate with short slits and 4% (± 1 lm) for the plate with long slits.The measured average of b over the slit length has been used b = 0.505 mm for short slits and b = 0.520 mm for long slits.Therefore, the porosity is U short = 2.7%, for the plate with short slits and and U long = 2.78%, for the plate with long slits.The edges in contact with the slits are kept as sharp as possible to

Numerical model
In this section, a two-dimensional numerical model is proposed.Solution of the linearized Navier-Stokes equations for an incompressible flow is considered.Consider the cross-section of the plates shown in Figure 3.A single slit of width b in a confinement channel of width a is considered.The confinement channel width is chosen such as the porosity of the two-dimensional model U 2D is equal to the actual porosity of the plate U.The plate of thickness t p is enclosed between an upstream and a downstream channel of lengths L u and L d .The computational domain is shown in Figure 5.The domain is divided into three sub-domains for meshing purposes.Continuity of pressure and velocity is assumed at the boundaries between the calculation domains.All the lengths are made dimensionless with the slit width b.The dimensionless plate thickness is t Ã p ¼ t p = b ¼ 5 mm=0:5 mm ¼ 10 (Distance BC in Fig. 5).The confinement channel width is a* = a/b = 18 mm/ 0.5 mm = 36 (Segment AH in Fig. 5).The lengths of the channels upstream (Segment AB in Fig. 5) and downstream the channel (Segment CD in Fig. 5) are chosen to have The upstream face of the plate is located at x Ã BG ¼ À4:5 and the slit opening at x* = 0.The downstream face of the plate is located at x Ã CF ¼ 5:5 and the end of the slit at x* = 1.The upstream side of domain AH is located at The downstream side of domain DE is located at x Ã DE ¼ ½À76:5; 77:5 and y Ã ¼ y=b ¼ ½0; 36.
The low He number approximation is made.The incompressible Linear Navier-Stokes equations for a 2D domain in a dimensionless form in the frequency domain are hereby presented for a radial frequency x: with x* = x/b and y* = y/b.The dimensionless velocity ðû Ã ; vÃ Þ is ðû=bx; v=bxÞ and the dimensionless pressure is pÃ ¼ p=ðq bx ð Þ 2 ).These equations are implemented in COMSOL Multiphysics as user defined equations (PDE) and solved.For the inlet segment AH and outlet segment DE uniform pressures are imposed, p Ã AH ¼ 1 and p Ã DE ¼ 0. At the segments BC and FG (walls) no-slip boundary conditions are applied, i.e. (u*, v*) = (0, 0).At the segments AB, CD, FE, GH slip boundary conditions are implemented to simulate the hydro-dynamical interaction between neighbouring slits in micro-slit plate with multiple slits.In the present case, the slits are very thin (b/a ( 1) therefore a slip boundary condition is used instead of a periodic boundary condition.A slip condition implies an equality of the velocities at the corresponding boundaries.Due to the dimensions of the slits, the deviation from a symmetric case is small.An unstructured mesh of quadratic triangular elements is used, with a finer mesh at the walls with no-slip conditions.The mesh inside the plate domain (Domain 2) is finer with the maximum element size is M el /b = 5 Â 10 À2 and the minimum is m el /b = 1 Â 10 À4 .A mesh convergence study shows a quadratic convergence of the results using the average velocity in a cross-section (line) calculated at a location x* = À1.7a*.Far from the slit, the acoustic pressure is uniform in the cross-section and the amplitude of the pressure depends linearly on the position along the duct  (parallel flow behaviour described in Aulitto et al. [20]).For À1.7a* < x* < À1.1a*, one has The complex constants can be determined by a linear fit of the pressure data obtained by numerical simulations for these regions far from the discontinuity.For a Sh b = 2.5 (corresponding to f = 120 Hz for a slit width b = 0.5 mm), the linear fit gives a coefficient of determination 1 À R 2 = 10 À6 [30].The transfer impedance of the slit Z Ã slit is determined by Z Ã slit ¼ bÀ d Û Ã with Û Ã being the flux calculated in a generic section of the channel upstream the slit, defined as slit one can derive the dimensionless transfer impedance of the slit z slit as with Sh b the Shear number based on the slit width.Tests for several channel widths a are performed.It appears that the resistance of the slit is independent of the porosity of the plate.Reducing of a factor 7 the porosity, the resistance increases by less than 0.7%.The inertia of the plate changes with changing the porosity.For a drastic reduction of the porosity, by a factor 7, the inertia increases by 30%.For small changes of the porosity around the nominal value, as the difference found between the two samples, the change is negligible.One can conclude that the transfer impedance of the slit appears to be only weakly dependent on the porosity U.

Results in linear regime
In Figure 6, the resistive and the reactive part of the impedance are shown for long and short slits in the range 2.5 < Sh b < 6 corresponding to 120 Hz < f < 700 Hz.The slit impedance of a single slit is displayed, z slit = z plate U, calculated using U short and U long for the short slits and long slits, respectively.The amplitude of the acoustic waves is pðx ref Þ ¼ 2 Pa.Frequencies below 120 Hz are excluded because of the presence of non-linear effects that will be discussed in Section 5.2.Frequencies above 700 Hz are ignored due to uncertainties in the measurements.The results are presented for a sample holder L h = 9 mm.
It appears that the maximum deviation between the resistance of long and short slits is of the order of 4%.This deviation is most probably due to the difference in slit width between the two plates and could be due to differences in edge sharpness.At low Sh b the resistance scales with 1/b 2 .An uncertainty of 2% in b explains a difference of the order of 4%.The influence of the length of the slits is negligible.The difference between the inertance of short and long slits is of the order of 10%.The inertance of short slits is lower than the inertance of long slits.This is due to three-dimensional effects due to the geometry of the slits.Three-dimensional effects do not change significantly for the plate confined from both sides (holder with L h = 50 mm) or for a free plate (holder with L h = 9 mm).The effect of the sample holder length is negligible for f < 500 Hz.At higher frequencies the data with the long holder show oscillating frequency dependency.This is to expected to be connected to the presence of a table in the measurement room.Therefore, the small holder is chosen to display the results.In Figure 6, the two-dimensional numerical model for a single slit is compared to the experimental results.It appears that the 2D numerical model of a single slit predicts (within few percentage of accuracy) the impedance (both resistive and reactive part) of long slits and the resistance of short slits.In the same figure, the semi-analytical model for high Shear numbers of Aulitto et al. [20] for a plate with square sharp edges with thickness t p = 0.5 mm is shown.For the semianalytical model, the plate thickness is assumed to be the same as the slit width.It appears that the model predicts reasonably well the resistance of the plate with the geometry proposed in this work.This confirms that the resistance of micro-slit plate is a local effect, strongly affected by the geometry of the edge and less sensitive to the global geometry of the plates [20].The inertance obtained when assuming b = t p = 0.5 mm, on the other hand, is much lower than that of the plates used in this study.In conclusion, in the linear regime (Sh b > 2.5), the difference between the resistance (or resistive part of the impedance) for short and long slits is small.The inertance (or reactive part) of short slits is smaller (within 10%) than long slits due to three-dimensional effects.Experiments with two sample holders exclude a dependence on the confinement of the plate.It appears that the two-dimensional numerical model well predicts the resistance of long and short slits.

Results in non-linear regime
For Shear numbers (Sh b < 2.5) deviations from the numerical model for both short and long slits appear.These deviations depend on the amplitude of the acoustic waves.This is shown for long slits for several amplitudes in Figure 7.These deviations at low frequencies are due to non-linear effects and decrease for increasing Shear numbers.At Sh b < 2.5 for p x ref ð Þ ¼ 2 Pa one has St b < 0.4.Non-linear effects impact both the inertance and the resistance of the slit.In particular, at the lowest frequency, the resistance for the highest amplitude (p x ref ð Þ ¼ 23 Pa) is 7 times higher than the resistance at the lowest amplitude (p x ref ð Þ ¼ 0:4 Pa).The inertance is reduced by almost a factor 2, as expected from the model of Ingard and Ising [21].
In this study, sound generation as higher harmonics due to non-linearity is not considered.For long slits, the plots of DR NL and DI NL are provided as function of the Strouhal number for several amplitudes of the acoustic waves.On the lower horizontal scale, the Strouhal number based on the slit width St b is used.On the upper horizontal scale, the Strouhal number based on the plate thickness at the slit St t appears.In Figure 8, results are shown for high amplitudes of the acoustic waves (pðx ref Þ ! 10 Pa).Both the change in resistance and inertance due to non-linearity increase (in absolute value), as expected, for decreasing Strouhal number (St b ¼ xb=U p ).For 1/St b ) 1 the change in resistance DR NL is approaching the theoretical quasisteady potential flow limit with correction for the viscous boundary layer for St b ?0. The change in inertance approaches the value DI NL = À0.5.The correction in the inertance due to non-linear effects is almost half the inertance in the linear case.Figure 9 compares the non-linear resistance and inertance for long slits to that of short slits for high amplitudes (p x ref ð Þ ! 10 Pa).In Figure 10, the results are shown for low and moderate amplitudes for short and long slits (p x ref ð Þ 6 Pa).These results display some Shear number dependency, which is less pronounced at higher amplitudes.
It can be seen that the resistance changes due to nonlinearity are almost identical for long and short slits, both for high and moderate amplitudes.Non-linear effects on the inertance of short slits are, in absolute value, smaller than for long slits.At moderate amplitudes (1/St b < 20) one observes a weak non-linear behaviour reported by Ingard and Labate [19].The vortices are formed at the edges, but they remain close to the slit.For 1/St b > 20 and for higher amplitudes, the vortices start moving away from the slits.One observes in this region differences between long and short slits.For very high amplitudes 1/St b > 50 the behaviour of the inertia of short and long slits is completely different.The flow for the short slits becomes essentially three-dimensional while it remains approximately two-dimensional for the long slits because of the confinement in the impedance pipe and hydrodynamical interactions between slits.Differences for high amplitudes are most probably due to different behaviour of the synthetic jet (zero net mass flow) in 3D depending on the length of the aperture [31].Difference in the slits (such as small surface perturbations in the edges) can generate different behaviour of the jet and different spatial  evolution of three-dimensional vortices.Examples of complex three-dimensional behaviour of jet formed by a slit are the axis switching and forking [32,33].The behaviour of short and long slits for low St b is different.In axis switching, the lateral ends of the jet will curve towards the symmetry axis, so that, within a distance comparable to the jet width, a almost plane jet will be formed in a direction normal to the original jet.In forking, the planar jet breaks down into separate jets.These observations confirm that while the resistance is determined locally, the inertance is a more global flow effect.

Discussion and conclusions
In this work, the effect of slit length on the linear and the non-linear acoustic transfer impedance of two accurately manufactured micro-slit plates has been investigated.The study is limited to sharp edges.Experimental data are obtained by impedance tube measurements.The frequency range of interest is 20 Hz < f < 700 Hz corresponding to Shear numbers in the range 1 < Sh b < 6.As for circular micro-perforations, the dissipation of sound occurs mainly at the edges of the micro-slits [20,25].Both in the linear and non-linear regime, the resistance of a single slit is independent of the porosity of the plate and on the slit length.Therefore, by dividing the transfer impedance of a slit by the porosity of the plate, one obtains an accurate prediction of the plate resistance for micro-slit plates.In the linear regime, an accurate prediction of the resistance of the plate is obtained using a two-dimensional numerical solution of the linearized Navier-Stokes equations in a single slit.In the non-linear regime, the plate resistance is strongly dependent on the amplitude of the acoustic waves.A simple quasi-steady model provides an order of magnitude for the asymptotic value of the non-linear resistance.Combining the results obtained in the present work for a specific geometry with the earlier studies on linear acoustical properties of micro-slits (see Aulitto et al. [20] and Temiz et al. [25]) one can conclude that the independence of resistance on the slit length is due to the fact that the resistance is mainly determined by the local acoustic flow around the edges.Therefore, this observation is independent of the exact slit geometry.The inertance of slits is sensitive to the acoustic flow outside the perforations.This explains the strong difference between the behaviour of slits and circular micro-perforations [20].Also, one observes a small but yet significant reduction of the inertance of short slits with respect to long slits.The linear inertance of long slits is accurately predicted by the locally two-dimensional incompressible numerical model.For moderate amplitudes, the non-linear inertance is not strongly affected by the slit length.For high amplitudes, one observes a strong deviation between the inertance of short and long slits.For long slits, the inertance at very high amplitudes is reduced by a factor two with respect to the linear case, as predicted by the intuitive model of Ingard and Ising [21].For short slits, the amplitude dependency of the inertance is more complex due to three-dimensional effects.The weak dependency of the inertia on the slit length is also expected to be independent of the exact geometry of the slits.The non-linear behaviour depends on the vena contracta factor, which is dependent on the details of the slit geometry [8,34].However, the fact that the slit length influences more the nonlinear effects on the inertia than the resistance is expected to be independent of details of the slit geometry.Hence, most conclusions drawn from the present study are expected to be quite general.

Appendix A
A.1 Vena contracta factor for potential flow At high amplitudes the flow detaches at the edges and then forms a free jet that contracts.The computation of jet flow is difficult, but in plane potential flow the problem of free jet can be solved by conformal mapping.The contraction factor of the jet (vena contracta factor) can be calculated for several geometries following approaches found in literature [22,23].In Figure A.1 the geometry of a single slit is compared with a simplified model that will be used to calculate the vena contracta factor using the method of Spurk and Aksel [22].As defined in Section 4, the porosity of the two-dimensional model is U 2D = b/a.The simplified geometry is connected to the idea that, for low b/a, the main parameter is the 45°angle between the two sides of the slit, because the interaction of the jet with the walls can be neglected.
The emerging jet contracts from the cross-section B -B 0 of width b to the cross-section C -C 0 of width C v b, with C v the vena contracta factor.At section C -C 0 the pressure inside the jet is equal to the ambient pressure, since the curvature of the streamlines vanishes.From the Bernoulli's equation, the velocity on the boundary of the jet is where p 1 and p 0 are the pressure before the slit and after the slit, respectively.The duct can be associated to a region in the complex z-plane by z = x + iy, with i 2 = À1 and spatial coordinates (x, y).In order to determine the shape of the free jet, a mapping resulting from the definition of the complex conjugate velocity can be used Using conformal mapping, the flow region in the duct can be mapped into a velocity plane, the so-called hodograph plate.For small porosity (b/a ( 1) the flow can be generated from the superposition of a source of strength 4C v bu 0 at w = (u À iv) = 0 and sinks of 2C v bu 0 at w = À1, +1, Ài, +i, with u 0 the velocity at the edge of the jet.The complex potential is In the present work the limit for b/a ( 1 is considered.The mapping function z = z(f) has to be calculated to determine the free surface in the z-plane.From equation (A.2) follows that The solution to this integral is ðA:5Þ The integration constant can be found assuming w ¼ ð1 þ iÞu 0 = ffiffi ffi 2 p at z = 0 + ib and w = u 0 at z = 1 + iC v b/2.It follows that the vena contracta factor in the limit b/a ( 1 is C v = 0.82.

Appendix B B.1 Quasi-steady incompressible flow
Assuming a quasi-steady incompressible flow with a free jet of vena contracta factor C v , the pressure difference Dp(t) across the plate resulting from the oscillating flow velocity u p ðtÞ ¼ U p cosðxtÞ is: where q is the air density and u j is the free jet velocity.
For the continuity of the velocity one has that ûj ¼ u p =C v and u p ¼ u=U where u p ðtÞ ¼ U p cosðxtÞ is the cross-sectional averaged acoustical velocity in the perforation, uðtÞ ¼ U cosðxtÞ is the cross-sectional averaged acoustical velocity in the pipe upstream the plate and U is the porosity.Hence, the pressure difference across the plate is ðB:2Þ The instantaneous power dissipated P w is given by P w ¼ Áp uA i ; ðB:3Þ where where T is the period of the harmonic oscillation.Substituting equations (B.1) and (B.2), this becomes

A.
Aulitto et al.: Acta Acustica 2022, 6, 6 and the viscous boundary layer thickness in the tube (Sh d = D i /d v ).The reflection coefficient at the end of the tube (x = 0) is[27]

Figure 2 .
Figure 2. Impedance tube termination with a) short sample holder of length L h = 9 mm and b) long sample holder of length L h = 50 mm.

Figure 4 .
Figure 4. Picture of the samples: a) plate with short slits and b) plate with long slits.

Figure 3 .
Figure 3. Sketches of the plate and of the cross-section of a single slit with parameters.

Figure 5 .
Figure 5. Two-dimensional computational domain for the plate of thickness t p and slit width b enclosed between an upstream and a downstream channel.Three sub-domains (1-2-3) are defined for meshing purposes.

Figure 6 .
Figure 6.Comparison of the a) slit resistance and b) slit inertance for short (--) and long (---) slits as a function of Shear number in linear regime.The numerical model (&) and the semi-analytical model for a plate with square sharp edges (. . ..) are shown [20].

Figure 9 .
Figure 9. Change in a) slit resistance and b) slit inertance due to non-linearity for long ( ) and short ( ) slits as function of St b and St t for several amplitudes at the reference microphone (p x ref ð Þ ! 10 Pa).Lines are defined as in Figure 8.

Figure 10 .
Figure 10.Change in a) slit resistance and b) slit inertance due to non-linearity for long ( ) and short ( ) slits as function of St b and St t for several amplitudes at the reference microphone (p x ref ð Þ !6 Pa).Lines are defined as in Figure 8.

Figure A. 1 .
Figure A.1.Schematic representation of the slit and zoomed simplified model for the contraction of the free jet after the plate.

2 À Á 2
is the pipe cross-section.Assuming a harmonically oscillating velocity in the pipe u p , one finds for the time-averaged dissipated power "

1 1 2
Correction for boundary layer thickness For pðx ref Þ ¼ 40 Pa (really high amplitude) and C v = 0.82, we have jû p j ¼ C v ffiffiffiffiffiffi 2Áp q q ¼ 6:7 m=s.For b = 0.5 mm, the Reynolds number Re b ¼ qU p b=l ¼ 223.The steady viscous boundary layer has a thickness of the order of d v =b % ffiffiffiffiffi Re b q % 0:1.Therefore, one can expect the frictionless theory to have an accuracy of the order of 10%.For very high amplitudes, using the Twaites solution [35] of the integral boundary layer equation of Von Karman one can obtain an estimation of the viscous boundary layer thickness.One has h momentum thickness of the viscous boundary layer of thickness d v and x ¼ 0 is at the slit neck.Neglecting the boundary layer thickness and assuming a uniform velocity U ðxÞ the mass conservation law becomes U ðxÞðb À 2xÞ % U p b ðC:2Þ because the slit angle is p/2 (Fig. A.1) and for t ) b 12 qû p jû p j and ûp , with u p ðtÞ ¼ Re½û p expðixtÞ ¼ U p cosðxtÞ.For St b ( 1 one can assume a quasi-steady incompressible flow with a free jet of vena contracta factor 2 v U 2 : ðB:5Þ One can define the time-averaged non-linear plate resistance as R plate,NL,t using the expression Using as reference the expression in Temiz et al. [8] (R ref ¼ ðqU Þ=ð2C 2 v U 2 Þ) one finds the theoretical asymptote: The value of R plate,NL,t is used to derive the analytical asymptote in Figures 8-10.This value is fairly close to the asymptote 0.7 < R plate,NL,t /R ref < 0.8 found in Temiz et al. [6] and Auriemma [8].Temiz et al. [8] is considering circular perforations with sharp square edges.Auriemma is considering slit shaped perforations with right angled edges.In the geometry proposed in this work, the channel length is zero.
Assuming a linear velocity profile in the boundary layer of thickness d v , one has a displacement thickness d Ã v ¼ d v =2 and a momentum thickness h ¼ d v =6.This implies: Given Re b ¼ qbU p =l % 223 one has d Ã v % 0:05b.This implies a reduction of the power because of the reduction of the porosity U eff ¼ Uð1 À 2d Ã v =bÞ so that R plate;NL;t Cite this article as: Aulitto A. Hirschberg A. Arteaga I.L. & Buijssen E.L.R.H. 2022.Effect of slit length on linear and non-linear acoustic transfer impedance of a micro-slit plate.Acta Acustica, 6, 6.