A parameter optimization method of multi-layer nested slit resonator for the converter transformer noise in high voltage direct current converter station

Abstract

We propose a parameter optimization method of multi-layer nested slit resonator (MNSR) for the converter transformer noise in a high voltage direct current (HVDC) converter station. With the optimization objective of achieving quasi-perfect sound absorption at the target frequencies and the structural parameters of MNSR as the optimization variables, the optimization model of the structural parameters of MNSR is established and then solved by the sequential quadratic programming algorithm (SQP). The theoretical and simulation results of the four validation cases show that the optimized MNSR can achieve quasi-perfect sound absorption at all target frequencies, which is consistent with the optimization objective. The parameter optimization design method provides an efficient method to design a dual-frequency quasi-perfect sound absorber for the converter transformer noise.

Keywords

Parameter optimization metasurface converter transformer noise multi-layer nested slit resonator sound absorption

Introduction

High voltage direct current (HVDC) transmission can achieve long-distance, high-capacity power transmission with low line losses,¹ and has been playing an important role in China’s west-east power transmission and power system networking projects. In particular, the converter transformer is a critical piece in the HVDC transmission, and its noise is a main noise source in the HVDC converter station, which has a great impact on the physical and mental health of substation employees and surrounding residents and becomes an urgent environmental problem to be solved.

The noise generated by converter transformers is dominated by low-frequency harmonics,² but traditional sound-absorbing materials, such as porous materials,^3–5 have limited efficacy in controlling such low-frequency noise. Recent acoustic metasurfaces provide a promising solution for low frequency noise control.^6–10 Broadband absorbers are an alternative, but multi-harmonic frequency absorbers could be a more tailored solution for the converter transformers. Duan et al.¹¹ designed a tunable multi-frequency absorber using the series connection of multilayer Helmholtz Resonators (HRs), achieving near-perfect acoustic absorption at 110 Hz and 392 Hz. Liu et al.¹² proposed a dual-band absorber with the series connection of HRs and then constructed a broadband absorber that achieved near-perfect continuous absorption within 380 Hz–3600 Hz. The aforementioned literature only utilizes the first-order acoustic cavity modes of the structures, while the multilayer nested slit resonant previously proposed by the authors utilized the first-two orders of acoustic cavity modes of the structure to achieve dual-frequency quasi-perfect acoustic absorption,¹³ showing a good application prospect. However, the two frequencies affect each other, developing an efficient design approach to the structural parameters according to the target frequencies is a worthy research topic.

Currently, a number of researchers optimize the structural parameters of metasurface with the developed optimization algorithms to achieve the desired absorption coefficients at the target frequencies. Sun et al.,¹⁴ Ryoo et al.,¹⁵ and Romas et al.¹⁶ solved the structural parameter optimization model with the sequential quadratic programming algorithm (SQP). Yan et al.¹⁷ employed the particle swarm optimization algorithm to optimize the structural parameters of single-layer honeycomb microperforated plate to construct a broadband sound absorber, and the sound absorption coefficient is all above 0.85 within the range of 1000–3500 Hz. Li et al.¹⁸ optimally designed a convoluted spatial structure based on the genetic algorithm, and can reach a normal incident average absorption coefficient of 0.80 at [342–2000] Hz. In this paper, the SQP algorithm with high computational efficiency and good convergence is chosen to optimize the structural parameters of multi-layer nested slit resonator (MNSR).

In this paper, the theoretical acoustic impedance model of MNSR is first established based on the transfer matrix method. Taking the quasi-perfect sound absorption at the target frequencies as objective and the structural parameters of MNSR as optimization variables, the structural parameter optimization model is established. The SQP algorithm is applied to solve the optimization model, and the accuracy of the optimization model is verified by the finite element method.

Theory and verification

Theory

Figure 1(a) illustrates the structure of the MNSR. Figure 1(b) shows the air domain of the MNSR, which is separated into multiple air passageways by the nested partitions and connected sequentially by the connecting slits. The transfer matrices of the connecting slits, air passageways, and inner cavity will be derived with the equivalent parameter of air. All connecting slits and air passageways are set to be of the same width to simplify the theoretical derivation.

Figure 1.

The schematic diagram of MNSR. (a) Structural unit, (b) air domain of MNSR.

The transfer matrix $T_{c}$ of the connecting slit is expressed as¹⁹

T_{c} = [\begin{array}{l} \cos (k_{eff, c} l_{c}) & j Z_{eff, c} \sin (k_{eff, c} l_{c}) \\ jsin (k_{eff, c} l_{c}) / Z_{eff, c} & \cos (k_{eff, c} l_{c}) \end{array}]

(1)

where,

j = \sqrt{- 1}

is the imaginary unit,

l_{c}

is the length of the connecting slit,

k_{eff, c} = ω \sqrt{ρ_{eff, c} C_{eff, c}}

and

Z_{eff, c} = \sqrt{ρ_{eff, c} / C_{eff, c}} / S_{c}

are the equivalent wave number and equivalent impedance,

ω = 2 π f

is the angular frequency of incident sound and

f

is frequency,

S_{c} = w_{c} H

is the cross sectional area,

w_{c}

and

H

are the width and height. The equivalent density

ρ_{eff, c}

and equivalent compressibility modulus

C_{eff, c}

are expressed as¹⁰

ρ_{eff, c} = ρ_{0} \frac{v w_{c}^{2} H^{2}}{j 64 ω} {(\sum_{m = 0}^{\infty} \sum_{n = 0}^{\infty} {(α_{m}^{2} β_{n}^{2} (α_{m}^{2} + β_{n}^{2} + \frac{j ω}{v}))}^{- 1})}^{- 1},

(2)

C_{eff, c} = \frac{1}{P_{0}} - \frac{j 64 ω (γ - 1)}{P_{0} ν^{'} w_{c}^{2} H^{2}} \sum_{m = 0}^{\infty} \sum_{n = 0}^{\infty} {(α_{m}^{2} β_{n}^{2} (α_{m}^{2} + β_{n}^{2} + \frac{j ω γ}{ν^{'}}))}^{- 1},

(3)

where,

α_{m} = (2 m + 1) / w_{c}

β_{n} = (2 n + 1) / H

, m and n are nonnegative integer,

ν = μ / ρ_{0}

and

ν^{'} = κ / ρ_{0} C_{v}

are constant. The air physical parameters of specific heat ratio

γ (γ = 1.4)

, air density

ρ_{0} (ρ_{0} = 1.21 kg / m^{3})

, air pressure

P_{0} (P_{0} = 101.325 kPa)

, dynamic viscosity

μ (μ = 1.8 \times 10^{- 5} N \cdot s / m^{2})

, thermal conductivity

κ (κ = 0.02624 W / (m \cdot K))

, specific heat at constant volume

C_{v} (C_{v} = 7178.8 kJ / (kg \cdot K))

characterize the thermoviscous loss of MNSR.

The transmission matrix of sound pressure and volume velocity of the incident sound from the left side ( $p_{in, i}$ and $u_{in, i}$ ) of the ith air channel to the right side ( $p_{out, i}$ and $u_{out, i}$ ) are as follow,¹³

[\begin{array}{l} p_{in, i} \\ u_{in, i} \end{array}] = [\begin{array}{l} \cos (k_{eff 2, i} l_{t, i}) & j Z_{eff, i} \sin (k_{eff 2, i} l_{t, i}) / 2 \\ j 2 \sin (k_{eff 2, i} l_{t, i}) / Z_{eff 2, i} & \cos (k_{eff 2, i} l_{t, i}) \end{array}] \cdot [\begin{array}{l} p_{out, i} \\ u_{out, i} \end{array}] = T_{pa, i} \cdot [\begin{array}{l} p_{out, i} \\ u_{out, i} \end{array}],

(4)

where,

l_{t, i}

is the equivalent propagation distance of sound in the ith air passageway,

k_{eff 2, i}

and

Z_{eff 2, i}

are the equivalent wave number and equivalent impedance of ith air passageway. Since the air passageway has a rectangular cross section of width w and height

H

, the equivalent wave number

k_{eff 2, i}

and equivalent impedance

Z_{eff 2, i}

can be calculated by analogy with the connecting slit.

The transfer matrix $T_{ca}$ of the inner cavity with rectangular cross-section can be expressed as

T_{ca} = [\begin{array}{l} \cos (k_{eff 3} l) & j Z_{eff 3} \sin (k_{eff 3} l) \\ jsin (k_{eff 3} l) / Z_{eff 3} & \cos (k_{eff 3} l) \end{array}],

(5)

where, the equivalent wave number

k_{eff 3}

and equivalent impedance

Z_{eff 3}

can also be calculated by analogy with the connection slit, and

l

is the propagation distance in the inner cavity.

In addition, due to the section mutation between the outermost connecting slit and the air outside the metasurface, it should also include the correction transfer matrix $T_{Δ L_{slit}}$ , which can be expressed as²⁰

T_{Δ L_{slit}} = [\begin{array}{c} 1 & j Z_{eff, c} k w_{c} ψ \sum_{s = 1}^{\infty} \frac{\sin^{2} (s π ψ)}{{(s π ψ)}^{3}} \\ 0 & 1 \end{array}],

(6)

where,

s

is positive integer,

ψ = w_{c} / L

is slit rate,

k = ω / c_{0}

is wave number, and

c_{0} = 343 m / s

is sound speed.

The transfer matrix $T$ from the surface of MNSR to the inner cavity can be expressed as

T = T_{Δ L_{slit}} T_{c} T_{pa, 1} \cdot \cdot \cdot T_{c} T_{pa, (N - 1)} T_{c} T_{ca}

where,

N

is the number of nested air passageway. The acoustic impedance of MNSR is

Z = T (1,1) / T (2,1)

, and the sound absorption coefficient

A

are written as

A = 1 - {| (Z - Z_{0}) / (Z + Z_{0}) |}^{2}

where

Z_{0} = ρ_{0} c_{0} / S

is the characteristic acoustic impedance and

S

is sound absorbing area.

Dual-frequency quasi-perfect sound-absorbing metasurface with MNSR

According to the conclusion of the parametric analysis in the previous study,¹³ the length and width of connecting slits and the width of the nested air passageways are taken as the optimization parameters (where $H$ is 20 mm and $L$ is 60 mm), and the optimization model is established with the objective of maximizing the sound absorption coefficient at the target frequencies, as shown below,

\begin{array}{l} (w_{1}^{'}, . . ., w_{4}^{'}, l_{c 1}^{'}, . . ., l_{c 4}^{'}, w_{c 1}^{'}, . . ., w_{c 4}^{'}) = \arg \max_{w_{1}, . . . w_{4}, l_{c 1}, . . ., l_{c 4}, w_{c 1}, . . ., w_{c 4}} \sum_{i = 1}^{2} α (w_{1}, . . . w_{4}, l_{c 1}, . . ., l_{c 4}, w_{c 1}, . . ., w_{c 4}, f_{i}) \\ subject to f_{i} = f_{target, i} (i = 1,2) \end{array}

(8)

In this paper, the SQP algorithm is chosen to solve this optimization model for its good convergence, high computational efficiency, and powerful boundary search capability. Sequential quadratic programming algorithm converts the original problem into a series of quadratic programming problems, takes its optimal solution as the next search direction of the original problem, and then performs a constrained one-dimensional search for the objective function of the original constrained problem in that direction so that the approximate solution of the original constrained problem can be obtained. Repeating the above process, the optimal solution of the original problem can be obtained, and the specific algorithm flow is shown in the following Figure 2.

Figure 2.

The schematic diagram of SQP algorithm flow.

Since the noise energy of the converter transformer is concentrated at 50 Hz and its multiples, the four cases, (1) 200 and 350 Hz, (2) 200 and 400 Hz, (3) 200 and 450 Hz, and (4) 200 and 500 Hz, are chosen for verification. Taking the perfect sound absorption coefficient at the target frequencies as the optimization objective, the optimized model of MNSR is solved with SQP, and the optimized structural parameters are shown in Table 1 below.

Table 1.

Optimal structural parameters of the cases (in mm).

Cases	w ₁	w ₂	w ₃	w ₄	l _c1	l _c2	l _c3	l _c4	l _c5	w _c1	w _c2	w _c3	w _c4	w _c5
Case 1	5	5	1	1	5	1	1.34	5	5	1	5	5	1	1
Case 2	4.98	5	1.32	1	5	1	1	1	3.00	1	5	5	1	1
Case 3	4.61	5	1.04	3.81	3.64	1	1.29	1.46	1.37	1.07	1.91	4.01	3.98	3.18
Case 4	3.58	3.04	1.37	2.98	4.82	1.02	1.04	3.79	3.99	2.18	4.99	4.99	1.27	1.15

The simulation validation of the optimized MNSRs is performed using COMSOL Multiphysics 6.1, which is based on the finite element method. The simulation model is modeled using the “Pressure Acoustics, Frequency Domain Interface”, which is used to compute the pressure variations for the propagation of acoustic waves in fluids at quiescent background conditions. The normal plane wave radiation is used to simulate sound sources in the range of 50–700 Hz with intervals of 1 Hz. Acoustic hard boundary conditions are applied at the walls of the structure. Narrow region acoustics are applied to the internal air domain of the MNSR to account for thermal viscosity losses. The free tetrahedral is used to mesh the simulation model. The maximum mesh size of the air domain inside the MNSR and the waveguide are 3 mm, and 5 mm respectively, which meets the accuracy requirement. The absorption coefficient is calculated according to the incident sound pressure and the reflected sound pressure.

The corresponding finite element simulation model is constructed based on the above structural parameters, and the simulation results and the theoretical results are shown in Figure 3 below.

Figure 3.

Sound absorption coefficient versus frequency for the following cases, (a) target frequencies of 200 and 350 Hz, (b) target frequencies of 200 and 400 Hz, (c) target frequencies of 200 and 450 Hz, (d) target frequencies of 200 and 500 Hz.

The simulation results of the MNSR with optimized structural parameters are shown in Figure 3. As can be seen from the above figures, the simulation results are well consistent with the theoretical calculation results, indicating that the developed theoretical model can accurately characterize the acoustic response of the MNSR. However, there is still a slight error between the theoretical results and the simulation results due to the fact that the theoretical model simplifies the propagation pattern of sound in the air domain inside the MNSR. Furthermore, the theoretical and simulated absorption coefficients at all preset frequencies are greater than 0.95, which is consistent with the expected results, indicating that the sound absorption performance optimization model of MNSR and its solution method are an effective method to design dual-frequency sound absorber for the converter transformers noise.

Figure 4 depicts the intensity map of the sound pressure of the four optimized MNSRs at resonance frequencies. The sound pressure distribution is denoted by the relative sound pressure $| P / P_{0} |$ , while $P_{0}$ represents the incident sound pressure. From Figure 4(a), (c), (e) and (g), it can be seen that the first-order peak frequency of each optimized MNSR reaches its sound pressure maximum in the inner cavity, and has only one extreme value. The sound pressure at the second-order peak frequency of all the optimized MNSRs contains two sound pressure maxima, and the maxima are located in the inner cavity and the second air passageway according to Figure 4(b), (d), (f) and (h). The sound pressure maxima at the second-order peak frequency for Cases 2, 3, and 4 are located in the second air passageway, while those for Case 1 are located in the inner cavity. In addition, the sound pressure at the first-order peak frequency is greater than that at the second-order frequency in all cases.

Figure 4.

The intensity map of the sound pressure of the four optimized MNSRs at their peak frequencies. (a) 200 Hz and (b) 350 Hz in Case1, (c) 200 Hz and (d) 400 Hz in Case2, (e) 200 Hz and (f) 450 Hz in Case3, (g) 200 Hz and (h) 500 Hz in Case4.

Conclusion

In this paper, a structural parameter optimization method of a dual-frequency quasi-perfect absorber based on MNSR is proposed for the low-frequency harmonic noise of converter transformers in the HVDC converter station. With the optimization objective of achieving quasi-perfect sound absorption at the target frequency and the structural parameters of MNSR as the optimization variables, the optimization model of the structural parameters of MNSR is established and then solved by the SQP algorithm. According to the noise characteristics of the converter transformer, four different cases are set up in this paper to verify the established parameter optimization model. The theoretical and simulation results of all four cases match well, and all achieve quasi-perfect sound absorption at the target frequencies. We hope that the structural parameter optimization method of MNSR can provide guidance for the noise control of converter transformers in HVDC converter stations.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This project is supported by the Science and Technology Project of China Southern Power Grid (Grant No. GDKJXM20201968) and Graduate Research and Innovation Foundation of Chongqing, China (Grant No. CYB23019).

ORCID iD

Zhigang Chu

References

Long

Nilsson

. HVDC transmission: yesterday and today. IEEE Power Energy Mag 2007; 5(2): 22–31.

Liu

Luo

. Vibration and noise characteristics of the inductive filtering converter transformer. Electron Lett 2017; 53(10): 678–679.

Zhu

Sun

Tang

, et al. Gradient-structural optimization of metal fiber porous materials for sound absorption. Powder Technol 2016; 301: 1235–1241.

Chevillotte

. Controlling sound absorption by an upstream resistive layer. Appl Acoust 2012; 73(1): 56–60.

Lagarrigue

Groby

Tournat

, et al. Absorption of sound by porous layers with embedded periodic arrays of resonant inclusions. J Acoust Soc Am 2013; 134(6): 4670–4680.

Huber

. Recent progress in acoustic metamaterials and active piezoelectric acoustic metamaterials - a review. Appl Mater Today 2022; 26: 101260.

Gao

Zhang

Deng

, et al. Acoustic metamaterials for noise reduction: a review. Adv Mater Technol 2022; 7(6): 2100698.

Liang

Tao

, et al. Acoustic focusing by coiling up space. Appl Phys Lett 2012; 101(23): 233508.

Liu

Long-Xiang

Huang

, et al. Broadband acoustic absorbing metamaterial via deep learning approach. Appl Phys Lett 2022; 120(25): 251701.

10.

Long

Shao

Chen

, et al. Broadband near-perfect absorption of low-frequency sound by subwavelength metasurface. Appl Phys Lett 2019; 115(10): 103503.

11.

Duan

Shen

Wang

, et al. Acoustic multi-layer Helmholtz resonance metamaterials with multiple adjustable absorption peaks. Appl Phys Lett 2021; 118(24): 241904.

12.

Liu

Yang

, et al. Ultra-broadband acoustic absorption of a thin microperforated panel metamaterial with multi-order resonance. Compos Struct 2020; 246: 112366.

13.

Sun

Chu

. Dual-frequency anti-nonlinear sound-absorbing metasurface via multilayer nested microslit resonators. Phys Scripta 2023; 98(4): 045702.

14.

Sun

Chu

, et al. A parameter design method for multifrequency perfect sound-absorbing metasurface with critical coupled Helmholtz resonator. J Low Freq Noise Vib Act Control 2021; 40(4): 2054–2063.

15.

Ryoo

Jeon

. Dual-frequency sound-absorbing metasurface based on visco-thermal effects with frequency dependence. J Appl Phys 2018; 123(11): 115110.

16.

Ramos

Godinho

Amado-Mendes

, et al. Broadband low-frequency bidimensional honeycomb lattice metastructure based on the coupling of subwavelength resonators. Appl Acoust 2022; 199: 109038.

17.

Yan

Chen

, et al. Optimization design and analysis of honeycomb micro-perforated plate broadband sound absorber. Appl Acoust 2022; 186: 108487.

18.

Yan

. Acoustic transmission characteristics based on coiled-up space metamaterials. Appl Acoust 2023; 203: 109199.

19.

Lee

Kwon

. Estimation of the absorption performance of multiple layer perforated panel systems by transfer matrix method. J Sound Vib 2004; 278(4–5): 847–860.

20.

Shen

Yang

Guo

, et al. Low-frequency anechoic metasurface based on coiled channel of gradient cross-section. Appl Phys Lett 2019; 114(8): 083501.