Broadband design of acoustic metasurfaces for the stabilization of a Mach 4 boundary layer flow

A piecewise acoustic metasurface is designed to suppress the first mode while marginally amplifying the Mack second mode in a Mach 4 flat-plate boundary layer (BL) flow. The results of linear stability theory (LST) and the e N method demonstrate the stabilization effect and transition delay performance, respectively. However, the direct numerical simulation (DNS) results indicate that the designed broadband acoustic metasurface actually weakly excites the first mode with a slightly larger fluctuating pressure amplitude at the surface, which is in contrast to the analysis of LST. The discrepancies are found to lie in the ‘roughness’ effect caused by the recirculation zones inside the microslits and the alternating expansion and compression waves induced at the slit edges, which significantly amplifies the first mode. For further clarification of the competitive mechanism between the acoustic stabilization and ‘roughness’ destabilization effects of metasurfaces on the first mode, a carefully designed metasurface is installed at the maximum growth rate region, which excites the first mode on the metasurface but inhibits its development downstream.

Extensive studies have been performed to stabilize the Mack second mode and inhibit the BL transition. Kimmel [7] identified the following categories of transition control techniques that could damp the Mack second mode: (1) passives, such as global or local shaping [7,8], porous coatings [2]; and (2) actives, such as retuned blowing and suction [9], CO 2 injection [10], and plasma actuators [10]. Because of the severe heat fluxes and high temperatures around the high-speed vehicle surface, it is difficult to implement active control techniques, and therefore, passive techniques are of primary interest. Malmuth et al. [6] first proposed the porous coating concept that utilizes the viscosity in narrow pores to efficiently dissipate the acoustic Mack second mode. Fedorov et al. [2] deduced the surface impedance boundary condition to model the acoustic effect of realistic microstructures and applied linear stability theory (LST) to investigate the stabilization effect of various inflow conditions and pore shape parameters. Their theoretical finding was then verified in the GALCIT T-5 shock tunnel [11], with a 5° halfangle sharp cone at free-stream Mach numbers ranging from 4.59 to 6.4. The BL on the half cone covered by the porous surface was mostly laminar, whereas transition was observed halfway along the other half cone with an untreated solid surface. With the motivation that a practical porous coating should be symbiotic with the actual materials of TPS, Fedorov et al. [12] experimentally and theoretically studied the effect of a feltmetal coating with random microstructures on hypersonic BL instabilities and showed that this type of porous coating strongly stabilized the second mode and marginally destabilized the first mode. Maslov [13] experimentally studied the stabilization effects of both regular porous coatings and random felt-metal porous coatings and found that both types stabilized the second mode while the first mode was excited. Wagner et al. [14] pioneered the use of carbon-carbon (C/C), an intermediate state of C/C-SiC already employed in hypersonic vehicles for thermal protection, to control Mack second mode waves. Wind tunnel experiments have been conducted at Mach 7.5, and the stabilization of the second mode as well as an increase in the laminar portion of the BL over the porous surface has been observed [14]. Tu et al. [15] investigated the optimal porosity and pore radius for stabilizing the most unstable second mode waves. Zhu et al. [16,17] investigated the effect of a permeable wall on the second mode. They found that the suppression of the near-wall disturbances by the permeable wall can disrupt the phaselocked relationship and prevent the growth of secondary instabilities. Due to the coating's marginal effect on the mean flow and great suppression of the most unstable BL instability, this technique has been widely investigated by theoretical analysis [18][19][20], wind-tunnel experiments [11,12,[20][21][22], and direct numerical simulations (DNS) [23,24]. More updated progresses could be found in Refs. [25][26][27][28]. Recently, Zhao et al. [29] expanded the research subjects by introducing the acoustic metasurface concept. Acoustic metasurfaces are planar metamaterial structures constructed with monolayer or multilayer stacks of subwavelength building blocks, which have significantly broadened the horizon of acoustic wave manipulation from wave-front modulation to sound insulation and absorption. According to this definition, the aforementioned porous coating belongs to the absorptive acoustic metasurface. Contrary to the long-held belief that acoustic wave energy dissipation is indispensable in suppressing the Mack second mode, Zhao et al. [29] proposed a nearly nondissipative impedance-near-zero acoustic metasurface, which inhibits the growth of the Mack second mode by the out-of-phase behavior that minimizes the near-surface acoustic pressure. They further designed a reflection-controlled metasurface for suppressing the Mack second mode by controlling the reflection directions of disturbance waves [30].
Compared to the abundance of studies for stabilizing the Mack second mode, there seem to be limited options for controlling the first mode. Kimmel [7] noted that this mode could be stabilized by wall cooling, suction and a favorable pressure gradient. Another way to damp the first mode is to use a very thin perforated sheet stretched over a plenum chamber. When the TS waves propagate along the BL, the fluctuating pressure forces air in and out of the plenum chamber that produces a near-wall region of negative Reynolds stress. Carpenter & Porter [31] conducted LST analyses for the incompressible BL on a flat plate and showed that the TS waves are stabilized when the wall admittance phase is close to π/2. Whether the experiments will confirm the theoretical predictions remains to be seen. They indicated that this theory leads to severe constraints on the perforated sheet thickness, and the pressure variations may lead to areas of quasi-steady inflow and outflow, thereby modifying the mean flow and its stability characteristics. Wang & Zhong [32] analyzed the porous coating admittance and noticed that its phase angle was closely related to the destabilization of the first mode, i.e., the destabilization effect decreased as the phase angle of admittance decreased. However, they used an artificial porous coating boundary for numerical demonstration instead of realistic coating structures. Tian et al. [25] found that the first mode is amplified as the admittance phase θ varies from 0.5π to π, whereas the second mode is destabilized when θ tends to 0.5π and suppressed as θ ≥ 0.75π; furthermore, the stabilization or destabilization effects on Mack modes are intensified by the increase in the amplitude of admittance. Following this principle, they theoretically designed an acoustic metasurface that could damp the second mode in a wide frequency range with the first mode slightly destabilized. Zhao et al. [33] designed a metasurface that suppressed the first mode in a Mach 4 BL flow with the Mack second mode not considered. Furthermore, Wang et al. [34] proposed a broadband design strategy that damps the first mode and Mack second mode in the high-frequency region through LST analyses.
In the present study, a piecewise acoustic metasurface aimed at suppressing the first mode while marginally exciting the Mack second mode is designed by LST and subsequently validated by DNS in a Mach 4 BL flow. This paper is organized as follows. Section 2 introduces the numerical methods employed in the present research. After that, the high-speed BL stabilization problem and broadband design strategy are described in Section 3. Section 4 compares the results from LST and DNS. The discrepancies are further verified and discussed. The conclusions are summarized in Section 5.

Numerical methods
In the current research, the two-dimensional (2-D) DNS is employed to directly resolve the detailed BL flow fields among the metasurface microstructures, which are used to verify the effectiveness of the impedance boundary condition when employed in the LST and the stabilization performance of the designed broadband metasurface. The 2-D LST and e N methods are mainly used to analyze the increment ratio of the unstable modes and determine the transition delay performance on the metasurface. It is worth noting that the three-dimensional (3-D) first mode is the most unstable in a supersonic BL on a flat plate [35]. In the Mach 4 BL flow, the wave angles [Θ ≡ arctan(β/α r ), β and α r are respectively spanwise and streamwise wavenumbers] of these unstable 3-D waves calculated at x * = 0.1 m (Fig. 7a) fall in the range of 0 ~ 3.49, which is much smaller than M 2 ∞ − 1 . That implicates that the first modes considered here belong to the inviscid regime [35]. According to the previous studies [36][37][38], the 2-D analyses of the metasurface effect on the first mode could apply equally to the 3-D cases. Additionally, the ideal metasurface with 2-D microslits can be regarded as an equivalent 3-D one with the same acoustic characteristics, provided that these have the same hydraulic diameter d h , porosity n, and depth H [23].

Direct numerical simulation
The DNS governing equations are 2-D Navier-Stokes equations in dimensional conservation form: where Q * is a vector of conservative variables and E * and F * are flux vectors in the x and y directions, respectively. These vectors are expressed as Here, ρ * is the density; u * , v * are velocity components in Cartesian coordinates; e * = p * / (γ-1) + ρ * (u *2 + v *2 )/2 is the total energy; and p * is the pressure. The superscript * denotes dimensional variables. τ * and q * are the stress tensor and heat flux with the components where T * denotes the temperature. The dynamic viscosity μ * is calculated using Sutherlands law. The specific heat ratio and Prandtl number of perfect gas flows are assumed to be γ = 1.4 and Pr = 0.72, respectively. R is a constant with R = R 0 /M, where gas constant R 0 = 8.314 J/(mol·K); molar mass M = 0.029 kg/mol.

Linear stability theory
In this study, 2-D instabilities are considered, and the small nondimensionalized perturbation is expressed in harmonic wave form: The dimensionless coordinates x and y are normalized by the BL thickness scale l * = µ * e x * /ρ * e u * e . The flow velocity components u * , v * , temperature T * and pressure p * are normalized by u e * , T e * and ρ e * u e *2 , respectively. The subscript e denotes the mean flow , α is the dimensionless streamwise wavenumber normalized by 1/l * , and ω is the angular frequency normalized by u e * /l * . With an assumption of a local-parallel boundary layer, the small disturbances governing equations is derived by substituting disturbance quantities into compressible Navier-Stokes equations and equations of state and then subtracting the mean flow parts. The details can be found in Malik [39]. Substituting Eq. (5) into the disturbances governing equations and neglecting the nonlinear terms yields the linearized dispersion relation where nonzero elements of the matrix L are evaluated by the local mean flow velocity and temperature, which can be obtained from a compressible self-similar solution [39]. The boundary conditions are where A is admittance of acoustic metasurface. Notably, for a smooth wall, A is equal to zero.
The spatial theory is able to demonstrate the changes in amplitudes of eigenmodes with distance; therefore, the dispersion problem is converted into a spatial problem, in which ω is given and α is unknown. α is obtained by employing the single-domain spectral collocation method [39]. The growth rate of a wave perturbance is σ = − α i , where α i is the imaginary part of α. The disturbance wave is unstable with a negative α i . A regular acoustic metasurface structure consisting of microslits is shown in Fig. 1, with width 2b, spacing s, and depth H nondimensionalized by l * .
The acoustic metasurface admittance for normal incident waves was derived in Ref. [18,40] (6) L(α, ω)� = 0, where A is the metasurface admittance, Z is the metasurface impedance, n (≡2b/s) is the porosity of the acoustic metasurface, and the dimensionless viscous wavenumber k v and thermal wavenumber k t are given by Prk v , respectively [41]. The admittance of the acoustic metasurface is a complex number and can be expressed as A = |A|e iθ , where |A| denotes the admittance magnitude and θ denotes the admittance phase. Accordingly, the disturbance energy flux contributed by the acoustic metasurface per unit area is where c.c. denotes the corresponding complex conjugate. To suppress the disturbances, i.e., E w < 0, the admittance phase angle should be in the range of θ∈[0.5π, 1.5π].
The e N method is a semiempirical method based on LST, which predicts the transition location by accumulating the growth rate of unstable waves [42]. This method is applied to evaluate the transition delay performance by the designed acoustic metasurface. In the e N method, the logarithmic amplification ratio of the amplitude versus its initial value for each wave is: where x 0 is the location of onset of instability, and A 0 is the disturbance amplitude there. The envelope of these ratios, i.e., the N factor that represents the amplitude evolution of the most amplified disturbance, is:

Problem description
In previous studies, an acoustic metasurface was conventionally designed with relatively deep microslits or pores to absorb the acoustic energy of Mack second mode waves [24]. However, the first mode dominates the BL transition at a lower Mach number. In such cases, the acoustic metasurface should be designed to stabilize the first mode while marginally promoting the second mode. The stabilization problem is simplified as a highspeed flat-plate BL flow at zero angle of attack, as shown in Fig. 2a the leading edge and plate surface with approximately 50 grids within the BL to resolve the disturbance accurately. To achieve global optimum stabilization, the broadband acoustic metasurface is cut into 35 pieces with different sets of geometrical dimensions of microslits (Fig. 2b). According to the grid independence study in Ref. [43], the shallowest slit has a grid dimension of 10 × 15 (width × depth), while the deepest slit has a grid dimension of 15 × 20. The detailed dimensions of broadband acoustic metasurface will be introduced in Section 3.2 later. For comparison, a conventional metasurface with microslit depth H * = 1.0 × 10 − 3 m, width 2b * = 1.0 × 10 − 4 m (aspect ratio, Ar ≡ 2b * / H * = 0.1), and porosity n (≡ 2b * /s * ) = 0.5 is also meshed with a uniform grid dimension of 10 × 30 (Fig. 2c).
An unsteady disturbance source is introduced in the calculation of DNS by a slot of periodic blowing and suction between x 1 * = 0.02 m and x 2 * = 0.03 m with the expression where q w * is the mass flow flux in the normal wall direction. A small forcing amplitude ε = 1.0 × 10 − 3 is adopted, and the forcing frequency f * is fixed at 60 kHz. At this frequency, the first mode is gradually amplified along the surface.

Broadband design strategy
Tian et al. [25] and Wang et al. [34] theoretically investigated the influence of the admittance phase and amplitude on the first and Mack second modes under Mach 4 flow conditions. Tian et al. [25] designed a metasurface with gradient microslits to suppress the second mode across a wide frequency band without distinctly amplifying the first mode. However, the primary concern should focus on the stabilization of the first mode, and the gradient design is challenging for engineering manufacturing. Wang et al. [34] further designed a piecewise acoustic metasurface to obtain the local optimal admittance amplitude and phase to delay transition. The broadband design strategy outlined in Ref. [34] is reviewed below for the sake of completeness. Figure 3 shows the growth rates of unstable modes in the BL over the smooth wall analyzed by LST. Both the first mode (ω < 0.12) and second mode (ω > 0.23) coexist in the BL. Since the dimensionless angular frequencies ω of their maximum growth rates basically remain constant in the streamwise direction, ω can be chosen to characterize the BL instabilities. As shown in Fig. 4, the stabilization or destabilization effect on the unstable modes was closely dependent on the metasurface admittance phase θ, and the effect was intensified by an increase in the admittance magnitude |A|. The first mode is marginally stabilized if θ approaches 0.5π and is greatly destabilized near θ = π as |A| increases. In contrast, the second mode is restrained when θ ≥ 0.75π and amplified as θ approaches 0.5π. For the present study, the metasurface was cut into 35 pieces in the streamwise direction to enable engineering realization. The admittance A of each piece should simultaneously meet the following three principles through a careful selection of the geometrical parameters in Eq. (8); that is: 1). As the dimensionless angular frequency ω < 0.12, the admittance phase θ should be designed to approach 0.5π with a small |A| to ensure the stabilization effect on the first mode (Fig. 4a); 2). In the frequency region where the first and second modes overlap (0.12 < ω < 0.23), for all phases, amplitude |A| should be designed as small as possible to minimize the destabilization effect on the second mode; 3). When the frequency ω > 0.23, the admittance phase should satisfy θ > 0.75π, and the amplitude effect on the growth rate is negligible (Fig. 4b and c).
As shown in Fig. 5, the required geometry parameters are achieved by solving Eq. (8) with program loops in the following order: loop 1 is 0.1 × 10 − 3 m ≤ H * ≤ 1 × 10 − 3 m, loop 2 is 0.1 ≤ Ar ≤ 1.5, and loop 3 is 0.1 ≤ n ≤ 0.8. These parameters span the range of Fig. 3 Contour of the growth rates of the unstable modes as a function of the dimensionless angular frequencies ω and the coordinates along the streamwise x * for a smooth solid wall interest for practical applications. The selected parameters should first satisfy the admittance phase θ requirement (Step I in Fig. 5), and then the amplitude |A| should be small enough to marginally excite the second mode (Step II in Fig. 5). After optimization, the final microstructure parameters at different streamwise directions x * are shown in Table 1. The increase in H * along the streamwise ensures the stabilization effect on the first and the high-frequency second modes. Ar increases to 0.7 to detune the admittance phase θ to approach 0.5π for first mode, whereas a further increasement will lead to the increase in the amplitude |A|. The porosity n increases gradually along the streamwise direction, but restrains small values to obtain a small |A|. For the sake of simplification, the optimization processes are only performed at 7 locations in the streamwise direction, and the optimal parameters are then interpolated into the midpoints of 35 pieces of the designed metasurface.

Verification
The accuracy of DNS and the self-similarity compressible Blasius solution used in LST are first verified in the BL flow. As shown in Fig. 6, the profiles of velocity and temperature of DNS are consistent with the self-similar solutions at x * = 0.24 m. There are no  visible discrepancies between the velocity and temperature profiles of smooth walls and the two above-mentioned acoustic metasurfaces, which indicates that acoustic metasurfaces have minimal effects on the mean flow.

LST analysis
As shown in Fig. 7a, the broadband acoustic metasurface successfully damps wide frequencies of first mode disturbances compared with the smooth wall. For Mack second mode disturbances with low frequencies, as depicted in Fig. 7b, the designed metasurface weakly enhances the Mack second mode and promotes the position of the maximum growth rate upstream. The growth rates of the high-frequency second mode disturbances become less than zero when two metasurfaces are used, which indicates that the high-frequency disturbance is completely suppressed (Fig. 7c). For comparison, the conventional acoustic metasurface damps all the second mode perturbations but strongly excites the first mode. As shown in Fig. 8, the e N method is employed to verify the transition delay performance of the acoustic metasurface. In the frequency range of 40 ~ 80 kHz, the broadband acoustic metasurface damps the first mode disturbances, while the conventional metasurface amplifies them. It can also be observed that 60 kHz is the dominant frequency, of which the disturbance reaches the largest N at the end of the flat plate. This frequency corresponds to the dimensionless angular frequency ω with the range 0.03 ~ 0.9 in the streamwise direction, and that indicates the typical first mode as denoted in Fig. 3. With respect to the Mack second mode of low frequencies (160 ~ 240 kHz), the N factor curves of both acoustic metasurfaces move upstream, while the maximum N factor of the broadband metasurface is only slightly reduced compared to that of the smooth wall. The N factor in the high-frequency region (320 ~ 400 kHz) of both metasurfaces is less than zero, which indicates that both metasurfaces could completely suppress the high-frequency second mode. Figure 9 depicts the N-factor envelopments of different surfaces. Although the N factor of the broadband acoustic metasurface in the upstream region (x * = 0.05 ~ 0.1 m) is higher than that of the smooth wall, the transition would not be triggered because of the small N factor. Moving downstream, the N factor envelopment of the broadband acoustic metasurface is lower than that of the smooth wall, which indicates its potential to delay the transition. In contrast, the conventional one tends to prompt transition with a higher N factor envelopment under the Mach 4 inflow condition.

DNS verification
In the above LST analyses, the impedance boundary condition was applied to represent the acoustic effect of the real microstructures on the BL but ignored the recirculation flows inside the slits and the alternating expansion and compression waves induced at the slit edges [43] (Fig. 10). To take into account this simplification, we perform DNS of the broadband and conventional acoustic metasurfaces (see Fig. 2) to verify the stabilization effect at the dominant frequency of 60 kHz. As shown in Fig. 11, compared with the smooth wall, the designed broadband acoustic metasurface actually destabilizes the first modes with larger fluctuating pressure amplitudes at the surface, which is in contrast to the analysis of LST in Section 4.2 (Fig. 7a). The conventional metasurface largely excites the first mode, as expected. For clarification, the impedance boundary condition for the conventional one is also implemented in DNS for comparison. The impedance model  could generally represent the effect of realistic metasurface microstructures for the first mode, but the results are not consistent with the comparisons for the Mack second mode [43]. The fluctuating pressure amplitude is obviously larger on the realistic metasurfaces with microslits, which indicates that the local flow perturbations would represent a 'roughness' effect and excite the first mode. In the previous research on the stabilization of the Mach second mode, the roughness effect of metasurfaces is considered to be negligible. Compared to the roughness insensitivity of the second mode, the first mode is easier to be excited [12,44]. This work shows that the carefully designed metasurface acoustic characteristics by LST could stabilize the first mode, but the DNS results revealed that the local flow perturbance induced by the microstructures will amplify the first mode. It can be concluded that the acoustic metasurface has simultaneous opposite influences on the first mode, and quantitative analyses need further research.
To prove the stabilization effect of the acoustic metasurface on the first mode by DNS, the wall temperature increases to T w = 400 K to further amplify the first mode under the same freestream flow conditions as those in Section 3.1. The periodic suction-blowing disturbance is chosen to be 80 kHz in the unsteady calculation, which is the dominant frequency of the first mode for this wall temperature condition after LST analysis. The microslits are designed with uniform parameters: depth H * = 4.0 × 10 − 4 m, width 2b * = 6.0 × 10 − 4 m, and porosity n = 0.536. These parameters ensure that the metasurface admittance phase is equal to 0.52π to stabilize the first mode at 80 kHz. Figure 12 shows the disturbance growth rates predicted by LST. The hot wall condition approximately doubles the growth rate of the first mode compared to that in Fig. 7a. With the designed metasurface installed along the plate surface (0.05 m < x * < 0.4 m), the first mode is seen to be significantly suppressed. In the DNS study, the metasurface is only installed at the maximum growth rate region, i.e., 0.08 m < x * < 0.2 m to avoid the 'roughness' effect on the downstream flow. As shown in Fig. 13, the wall fluctuating pressure amplitudes at the metasurface region are larger than those on the smooth wall, which is contrary to the LST prediction again. The 'roughness' effect of the metasurface plays a dominant role in exciting the first mode. However, the amplitudes gradually decrease downstream of the metasurface, which indicates that the acoustic metasurface could substantially stabilize the amplification process of the first mode. On the other hand, 'roughness' is limited to the local influence. The results indicate that the acoustic metasurface should be installed in a confined region where the maximum growth ratio of the first mode occurs to minimize the 'roughness' effect.

Conclusion
In this paper, a piecewise acoustic metasurface is designed to suppress the first mode across a wide frequency band without amplifying the second mode distinctly in a Mach 4 flat-plate flow. The LST proves its stabilization effect, and the e N method demonstrates Fig. 12 Dimensional growth rates along the hot wall with a temperature of 400 K and disturbance frequency of 80 kHz Fig. 13 Comparisons of fluctuating pressure amplitude distributions along the hot wall with a temperature of 400 K and disturbance frequency of 80 kHz the trend of delay transition. However, the DNS simulations show contrary results. The discrepancies lie in the neglect of the local flow perturbation caused by the recirculation flows inside the microslits and the alternating expansion and compression waves induced at the slit edges, which behave as 'roughness' to excite the first mode. The competitive mechanism between the acoustic stabilization and 'roughness' destabilization effect of the metasurface on the first mode was also found, which was different from that for stabilizing Mack second mode. In the metasurface region, 'roughness' plays a leading role in exciting the first mode, while downstream, the disturbance amplification process is inhibited because of acoustic modulation by the metasurface. This not only means that the metasurface microstructures should be carefully designed but also that the installment location should be considered when applying the metasurface to delay the BL transition triggered by the first mode. However, the 'roughness' mechanism on the first mode should be further clarified. Additionally, the impedance model should be improved by including the 'roughness' effect to reach a consistent result between LST and DNS.