# 1 HIGH-AMPLITUDE ACOUSTIC FIELD IN A DISC-SHAPED RESONATOR

## Texto completo

(1)

1

### RESONATOR

M. Cervenka, R. Balek. M. Bednarik

CTU in Prague, FEE, Technicka 2, 166 27 Prague 6, Czech Republic e-mail: cervenm3@feld.cvut.cz

Abstract

This paper is concerned with theoretical study of high-amplitude shock-free acoustic field in a gas-filled piston-driven cylindrical resonator where the diameter exceeds its height in a substantial way. Due to this geometry, only radial modes are driven. Since these modes are non-equidistant, resonance conditions are not fulfilled for higher harmonics of the acoustic field resulting from nonlinear effects, they are suppressed, the acoustic saturation effect is avoided and Q-factor of the resonator increases. Such a resonant cavity behaves similarly as a

variable-cross-section resonator, but the geometry is much simpler which gives rise to its greater utilization e.g. in ultrasonic processing applications. For the purpose of numerical simulation of the problem, set of Navier-Stokes equations formulated in cylindrical coordinates is used. The equations are integrated by using of a central semi-discrete difference scheme, characteristic boundary conditions are implemented. Numerical results confirm excitation of shock-free high-amplitude waveform in the resonant cavity.

Keywords: acoustic resonator, Navier-Stokes equations, numerical integration

### Introduction

Nowadays, we can observe rapid development of nonlinear acoustics that is the branch of acoustics concerned with acoustic waves, whose amplitudes cannot be treated as infinitesimally small ones. This situation is caused by wide application area of these waves and also by fast development of powerful computer technology allowing their theoretical study.

One of the application areas is connected with utilization of high-amplitude acoustic standing waves generated in resonators driven by means of vibrating piston or inertial force. When a standing wave is driven into high amplitude, nonlinear effects (convection, nonlinear relation between density and pressure) cause distortion of originally harmonic wave. As the thermo-viscous attenuation is proportional to the square of frequency, dissipation of acoustic energy increases in a substantial way due to excitation of higher harmonics and acoustic pressure amplitude is no longer proportional to amplitude of driving.

(2)

2

Philosophy of this work issues from the firstly mentioned approach (shaped resonators), whereas it overcomes one its disadvantage - complicatedness of the acoustic resonator cavity. Here a thin, piston-driven cylindrical resonator is used, nevertheless, transversal mode is driven instead of a mode longitudinal. As higher eigenfrequencies are not integer multiples of the eigenfrequency fundamental for cylindrical waves, resonance conditions are not fulfilled for higher harmonics and thus appearance of shock wave and intense nonlinear dissipation can be avoided. Geometry of the problem is shown in Figure 1.

Figure 1 – Geometry of the problem formulated using dimensionless coordinates.

### Model equations

Model equations of nonlinear acoustics are represented by nonlinear partial differential equations. For the sake of simulation of the above mentioned problem, the appropriate model equation must be a multidimensional one, some of them have recently been proposed in literature, see e.g. [6-8]. Here, set of dimensionless equations issuing from Navier-Stokes equations is used. With the dimensionless quantities

2 2

0 0 0 0

where

### ,

are axial and radial coordinates, l is length of the resonant cavity,

is time,

### ω

is driving frequency,

### ,

0 are the total and ambient fluid density,

### ,

are the axial and radial components of acoustic velocity,

### c

0 is the small-signal speed of sound and

### ,

are the acoustic pressure and energy respectively, the set model equations reads

2 2 2

2 2 2

2 2 2 2

2 2 2 2

,

( ) 1

,

( ) 1

,

U V V

T Z R R

U P U UV UV G U U U

T Z R R Z R R R

V UV P V V G V V V V

T Z Z R Z R R R R

E T

### π

∂Λ+ ∂ Λ⎛ ⎞+ ∂ Λ⎛ ⎞= −Λ

⎜ ⎟ ⎜ ⎟

∂ ∂ ⎝ Ω ⎠ ∂ ⎝ Ω ⎠ Ω

⎞ ⎞

⎛ ⎛

∂ Λ + ∂ + Λ + ∂ Λ⎛ ⎞= −Λ +++

⎟ ⎜ ⎟ ⎟

⎜ ⎜

∂ ∂ ⎝ Ω ⎠ ∂ ⎝ Ω ⎠ Ω Ω ∂⎝ ∂ ∂ ⎠

⎞ ⎞

⎛ ⎛

∂ Λ + ∂ Λ⎛ ⎞+ ∂ + Λ = −Λ +++

⎟ ⎟

⎜ ⎟ ⎜ ⎜

∂ ∂ ⎝ Ω ⎠ ∂ ⎝ Ω ⎠ Ω Ω ∂⎝ ∂ ∂ ⎠

+

∂ ∂

( ) ( ) ( )

.

E P U E P V E P V

Z R R

+ ∂ + +

⎡ ⎤+ ⎡ ⎤= −

ΩΩΩ

⎣ ⎦ ⎣ ⎦

(3)

3

The equations are supplemented by the pressure-energy relation

2

2

where

### γ

is the ratio of specific heats. Here

### /

0 is ratio of the driving frequency and frequency of the fundamental longitudinal mode

0

0

and

0

0 0

### c

2 is the attenuation parameter where b is the diffusity of sound.

Left sides of the equations (1) are formulated in the conservative form, the right sides represent geometrical sources (due to formulation of the problem with use of curvilinear coordinates) and dissipative terms. Model of the acoustic energy dissipation process is in the equations (1) somewhat simplified, it is assumed that the acoustic field is non-vortical (the terms containing rotation of velocity are omitted) and the small third-order terms in the energy transfer equation are dropped. All the dissipation mechanisms are included in heightened value of the parameter G.

For the purpose of numerical integration of the set (1), high-resolution second-order central semi-discrete difference scheme  was used. The characteristic boundary conditions were implemented as the adiabatic free-slip ones  (hard resonator's walls at

0,

### 1

), time-varying position of vibrating piston (boundary at Z=0) was implemented to avoid nonphysical matter flux through the boundary. For R=0, symmetry of the problem was utilized.

### Eigenfrequencies of a cylindrical resonator

In order to obtain the first insight into the problem, the equations (1) can be linearized and solved analytically. However, it is simpler and more straightforward to issue from linearized wave equation that can be written in dimensionless form

2 2

2 2

2 2

(2)

where

### c l

0 is dimensionless velocity potential and i= −1.

Solution of the Eq. (2) can be easily found using the Fourier series. Using homogeneous Neumann boundary conditions at the walls and

p

### ∂

at the piston wall (Z=0), where

p

### R

is distribution of the vibration velocity, the solution can be written in the form

### ]

0 0 0

0

( , ) J / cos (1 ) ; 0,1 , 0, ,

sin

k

k k

k k k

R Z

R R

Z Z R R

=

Φ =

## ∑

− ∈ ∈ (3)

where

### J ()

0 is the first-kind zero-th order Bessel function,

### ν

k represents the k-th root of equation 1

2 2 2 2 2

0

k

k

and

### α

k are coefficients of the Bessel-Fourier decomposition of

### U

P, Individual eigenfrequencies can be found from the condition

2 2 2 2 2 2 0

k

### π

, where n is an integer, in form 2 2 2 2

, / 0.

n k n

k R

Ω = + (4)

(4)

4

shows that

0,k

k

0

### 0,

which means that the higher eigenfrequencies are not integer

multiples of the eigenfrequency fundamental. For example, if

0

### 6

(radius of the resonator is six times greater than its length), we obtain

0,k

### Ω =

For driving of the k-th transversal mode, it must be fulfilled that

k

### 0.

It should also be noted that for reaching high-amplitude acoustic field in the k-th transversal mode, value of the coefficient

### α

k must be maximized. For example, if the velocity of the left resonator's wall (piston) is given as

0

( ) H( ); 0,1 ,

p p

U R =A pRR p∈ (5)

where Ap, p are constants, H is the Heaviside function and the 1st transversal mode is driven, the maximum value for

### α

1 is obtained for

0

!

1 0 1 0

1 0

x

x pR

=

## ∫

where

### μ

1 is the 1st root of equation

0

.

### Numerical results

Some numerical results obtained by means of numerical integration of Eqs. (1) are presented in this section. For all the cases, hard-walled piston-driven cylindrical resonator filled with air at ambient room conditions is assumed. Parameter G accounting for all the attenuation processes was set to

3

### =

− . Two geometrical configurations are compared here. Firstly, it is a disc-shaped resonator with

0

### 6,

where the 1st transversal mode is driven using piston velocity distribution

0

p p

### T

which is close to "the optimal driving". And secondly, it is an elongated thin resonator with

0

### 1/ 6,

where the 1st longitudinal mode is driven using piston velocity distribution

p

p

### T

Acoustic pressures presented in the following figures are "measured" at the center of vibrating piston (

### 0

).

0 0.5 1 1.5 2 2.5 3

0.8 1 1.2 1.4 1.6

x 105

p [Pa]

Period [−]

Longitudinal mode Transversal mode

Figure 2 – Acoustic pressure in the case of longitudinal and transversal mode driving.

Figure 2 shows comparison of time courses of longitudinal (

### 1

) and transversal

(Ω =0.2064) mode pressures for the same driving amplitude 5 104

p

(5)

5

conditions. It can be seen that the transversal mode attains substantially higher amplitude compared to the highly distorted longitudinal mode, where shock wave appears.

Figure 3 shows comparison of amplitude characteristics (the 1st harmonics of acoustic pressure) in the case of longitudinal and transversal mode driving at resonance conditions. As it is common in nonlinear systems, the pressure amplitudes are not proportional to amplitude of driving. This situation is caused by excitation of higher harmonics that undergo stronger the nonlinear attenuation that is proportional to the square of frequency.

3 4 5 6 7 8 9 10

x 10−4 0

1 2 3 4

x 104

p 1

[Pa]

Ap [−]

Transversal mode

Longitudinal mode

Figure 3 – Amplitude of the 1st harmonics of pressure with respect to driving velocity in the case of longitudinal and transversal mode driving.

3 4 5 6 7 8 9 10

x 10−4 1

1.01 1.02 1.03

Ω

/

Ω LIN

[

]

Ap [−]

Figure 4 – Shift of resonance frequency as a function of driving velocity amplitude (transversal mode).

Figure 4 shows dependence of transverse-mode resonance frequency on amplitude of driving. The linear theory predicts this resonance frequency as

LIN

### 0.2033.

It can be seen from the

figure that the nonlinear resonance frequency increases with amplitude of driving. Such a behavior is due to the mechanical analogy called the hardening-spring resonance frequency shift. This resonance frequency shift is not observed for the case of longitudinal mode driving.

### Conclusions

(6)

6 Acknowledgements

This work was supported by grant MSM 6840770015 and GACR grant project No. 202/06/0576.

### References

 Atchley, A. A.; Gaitan, D. F. Finite amplitude standing waves in harmonic and anharmonic tubes. J. Acoust. Soc. Am. 93, 1993, pp. 2489-2495

 Illinski, Yu. A.; Lipkens, B.; Lucas, T. S.; Van Doren, T. W.; Zabolotskaya, E.. Nonlinear standing waves in an acoustical resonator. J. Acoust. Soc. Am. 104, 1998, pp. 2664--2674

 Luo, C.; Huang, X. Y.; Nguyen, N. T. Generation of shock-free pressure waves in shaped resonators by boundary driving, J. Acoust. Soc. Am. 121, 2007, pp. 2515

 Huang, P. T.; Brisson, J. G.. Active control of finite amplitude acoustic waves in a confined geometry. J. Acoust. Soc. Am 102, 1997, pp. 3256-3268

 Gusev, V. E.; Bailliet, H.; Lotton, P.; Job, S.; Bruneau, M. Enhancement of the Q of a nonlinear acoustic resonator by active suppression of harmonics. J. Acoust. Soc. Am. 103, 1998, pp. 3717-3720

 Vanhille, C.; Campos-Pozuelo, C. Numerical and experimental analysis of strongly nonlinear standing acoustic waves in axisymetric cavities, Ultrasonics43, 2005, pp. 652

 Cervenka, M.; Bednařík, M. Nonlinear standing waves in 2-D acoustic resonators. Ultrasonics. 2006, vol. 44, Suppl. 1, pp. 773-776

 Cervenka, M.: Numerical Simulation of Two-Dimensional High-Amplitude Acoustic Field in Cylindrical Resonators Using Kuznetsov's Equation. 19th International Congress on Acoustics [CD-ROM]. Madrid: Sociedad Espaňola de Acústica, 2007

 Kurganov, A.; Tadmor, E. High-Resolution Central Schemes for Nonlinear Conservation Laws and Convection-Diffusion Equations, J. Comput. Phys.160, 2000, pp. 241

## Figure Figure 1 – Geometry of the problem formulated using dimensionless coordinates. p.2 Figure 2 – Acoustic pressure in the case of longitudinal and transversal mode driving p.4 Figure 3 – Amplitude of the 1st harmonics of pressure with respect to driving velocity in the case of longitudinal and transversal mode driving p.5 Figure 4 – Shift of resonance frequency as a function of driving velocity amplitude (transversal mode) p.5

## Referencias

Actualización...

Related subjects : Nonlinear effects
{"ft_lang_label":"__label__en","ft_lang_prob":0.8258629,"math_prob":0.956351,"size":12305,"snap":"2022-05-2022-21","text_gpt3_token_len":3823,"char_repetition_ratio":0.13543615,"word_repetition_ratio":0.06555944,"special_character_ratio":0.2763917,"punctuation_ratio":0.14029364,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.98324776,"pos_list":[0,1,2,3,4,5,6,7,8],"im_url_duplicate_count":[null,null,null,null,null,null,null,null,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2022-05-19T05:17:23Z\",\"WARC-Record-ID\":\"<urn:uuid:6d0bd20c-ea73-4736-a337-d7eb5344223f>\",\"Content-Length\":\"263427\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:a644ea17-5c9e-4bd2-8cf0-dab41b83d62d>\",\"WARC-Concurrent-To\":\"<urn:uuid:991dfef3-310f-43ac-934a-e1ac5118508f>\",\"WARC-IP-Address\":\"164.90.139.124\",\"WARC-Target-URI\":\"https://1library.co/document/eqook55q-high-amplitude-acoustic-field-in-disc-shaped-resonator.html\",\"WARC-Payload-Digest\":\"sha1:33XFNP4KLJKAXPWT5JKA6IDJP4MM6HTT\",\"WARC-Block-Digest\":\"sha1:GLDO43Y4SEG53LJN2CWH4UKSVMRSLFCG\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2022/CC-MAIN-2022-21/CC-MAIN-2022-21_segments_1652662525507.54_warc_CC-MAIN-20220519042059-20220519072059-00096.warc.gz\"}"}