Underwater Spiral Wave Sound Source Based on Phased Array with Three Transducers

This paper realizes an underwater spiral wave sound source by using three omni-directional spherical transducers with three different phases. The pressure distribution of the sound field for a phased array is derived using the superposition theory of sound field. The generation of spiral wave field is presented, the relationship between the performance of phased array sound field and the array parameters is analyzed, and also verified by the finite element method (FEM). A spiral wave sound source with three spherical piezoelectric ceramic transducers is then designed and fabricated based on FEM simulation, and the performance of the sound source is analyzed. Measurements are made in a reverberation pool, and the result shows that the fabricated spiral wave sound source is capable of producing a spiral sound wave. Under a frequency of 3.5 kHz, the phase directivity has a fluctuation of ±21°, and the amplitude directivity range is 4.3 dB, which verifies the realization of the spiral wave sound source.


Introduction
As radio-based navigation is not available in deep seas, underwater acoustic navigation becomes significantly important, especially for autonomous underwater vehicles (AUVs). Acoustic navigation can be performed by long baseline, ultra-short baseline, and recently, a spiral wave method. Among these methods, navigation with spiral wave beacons offers the unique advantage of utilizing limited resources: only one hydrophone is required in an AUV, unlike the long baseline method which requires emitting equipment and the ultra-short baseline method which requires multiple hydrophones to form an array [1,2].
There are two methods to produce a spiral sound wave. One way is to use in-phase radial vibration of a spiral radiation surface. The other way is to use a circular array, where oscillators vibrate with different phases. The first method includes, the "Physical Spiral" transducer proposed by Dzikowicz and Hefner [3][4][5]. This transducer uses 1-3 piezocomposite material formed around a spiral-shaped backing to emit the spiral wave. The second method includes, the "Phase Spiral" transducer, also proposed by Dzikowicz and Hefner [4,5], the Cylindrical Transducer with divided electrodes proposed by Brown [6][7][8], and the Spiral Sound Transducer Based on Longitudinal Vibration Structure proposed by Wei Lu [9]. The "Phased Spiral" transducer also uses 1-3 piezocomposite material and backing, but the electrode is divided into 16 elements. The Cylindrical Transducer with divided electrodes uses a radially polarized piezoceramic hollow cylinder and the electrode is divided into four parts. Figure 1 presents the phased arrays in a polar coordinate system where the three array elements S 1 , S 2 , S 3 are placed at (d, 0), d, 2π 3 , d, 4π 3 , respectively.

Sound Pressure Field Produced by Phased Array with Three Transducers
Sensors 2018, 18, x FOR PEER REVIEW 2 of 15 elements, the three transducers can emit spiral waves. Of the above spiral wave transducers, the Cylindrical Transducer with divided electrodes has the best performance on phase directivity and transmitting voltage response. In addition to large fluctuation in phase directivity, other transducers suffer from other disadvantages. Too many driving voltages needed is one of these disadvantages. At present, to emit spiral sound waves by an array, it requires a large number of oscillators that vibrate with different phases. The larger the number of oscillators, the more signals are required, the more complex is the driving circuit, and larger the volume and weight of the transducer. This paper, therefore, analyzes a spiral sound source based on phased array with only three omni-directional transducers. Compared with the similar spiral sound sources, the proposed sound source requires the minimum number of transducers. First of all, the sound pressure field produced by a phased array with three transducers is analyzed theoretically. By separating the component of a spiral wave and the component of an interference wave in a sound field, the generation of a spiral wave field is discussed, and the influence of array parameters on the performance of the spiral wave source is analyzed. Then, based on the finite element simulation, a spiral wave sound source based on phased array with three spherical piezoelectric ceramic transducers is designed. The test result indicates that the designed sound source can emit spiral sound waves. Figure 1 presents the phased arrays in a polar coordinate system where the three array elements S , S , S are placed at , 0 , , , , , respectively. Assume that the working frequencies and vibration amplitudes of the three array elements are the same, with an initial phase of 0, , respectively, and all omni-directional.

Sound Pressure Field Produced by Phased Array with Three Transducers
When the size of the array element is much smaller than the wavelength, it can be regarded as a point sound source [10,11].
The sound pressure field of the point sound source can be expressed as: where, is the distance from point P to the point sound source, is the vibration angular frequency of the point sound source, is the medium wave number, and is the initial phase. According to the superposition theory of sound field [12,13], the sound pressure field produced by the phased array is given as: Assume that the working frequencies and vibration amplitudes of the three array elements are the same, with an initial phase of 0, 2π 3 , 4π 3 respectively, and all omni-directional. When the size of the array element is much smaller than the wavelength, it can be regarded as a point sound source [10,11].
The sound pressure field of the point sound source can be expressed as: where, r is the distance from point P to the point sound source, ω is the vibration angular frequency of the point sound source, k is the medium wave number, and ϕ 0 is the initial phase.
According to the superposition theory of sound field [12,13], the sound pressure field produced by the phased array is given as: There are the following expressions for r 1 , r 2 , r 3 : At the far field when the distance of point P from the origin r is greater, such that r d 2 λ , the difference in the amplitude of sound pressure is negligible. Therefore, the distance r 1 , r 2 , r 3 from three array elements S 1 , S 2 , S 3 to point P respectively, can be replaced with the distance r. On the other hand, the phase difference caused by the different distances to the array elements cannot be ignored.
Further expanding and rearranging, we can obtain: Ψ = cos(dk cos(ϕ)) + i sin(dk cos(θ))− cos 1 2 dk cos(−ϕ) cos 1 After extracting the spiral phase term, we can finally obtain: − 1 cos 1 2 dk cos(ϕ) Thus, the spiral wave field component p 1 of the total sound field can be written as: A r e j(ωt−kr) idke iϕ , (9) and the interference wave component p 2 can be written as: It can be noted that the spiral sound field component p 1 contains the phase component e iϕ , that changes linearly with azimuth.

The Performance of the Spiral Sound Field
The sound field of phased array contains the spiral sound field component p 1 and the interference component p 2 , along with a big phase characteristic difference between p 2 and p 1 , leading to fluctuation in phase directivity. The relationship between the performance of the spiral sound field and the product of the array radius and wavenumber is discussed from the variation of the relative magnitudes of two components with dk (the product of array and radius wavenumber).
The proportion of the interference component in the total sound field in the far field can be written as: Figure 2 shows a variation diagram of Ratio 1 (ϕ) with dk. It can be seen that the proportion of the interference component in the total sound field increases with the increase in dk. When dk is very small, the proportion of the interference component in the total sound field is very small, whereas the proportion of the spiral sound field component in the total sound field is very large. Then, with the increase of dk, there is a gradual increase in the proportion of the interference component, and a gradual decrease in the proportion of the spiral sound field component in the total sound field.
It can be noted that the spiral sound field component contains the phase component e , that changes linearly with azimuth.

The Performance of the Spiral Sound Field
The sound field of phased array contains the spiral sound field component and the interference component , along with a big phase characteristic difference between and , leading to fluctuation in phase directivity. The relationship between the performance of the spiral sound field and the product of the array radius and wavenumber is discussed from the variation of the relative magnitudes of two components with (the product of array and radius wavenumber). The proportion of the interference component in the total sound field in the far field can be written as: Figure 2 shows a variation diagram of Ratio with . It can be seen that the proportion of the interference component in the total sound field increases with the increase in . When is very small, the proportion of the interference component in the total sound field is very small, whereas the proportion of the spiral sound field component in the total sound field is very large. Then, with the increase of , there is a gradual increase in the proportion of the interference component, and a gradual decrease in the proportion of the spiral sound field component in the total sound field. The phase directivity function of the sound field of phased array can be written as: The slope of the phase directivity function of the sound field of phased array can be written as: Figure 3 shows a variation diagram of the phase directivity, and the slope of the phase directivity curve with . It can be seen from Figure 3 that when is very small, the phase directivity of phased array has a good linearity, and the slope of the phase directivity has a small amplitude fluctuation around 1. Then, with the increase in , the nonlinear fluctuation of the phase directivity of phased array increases gradually, and the fluctuation of the slope of the phase directivity increases gradually. The phase directivity function of the sound field of phased array can be written as: The slope of the phase directivity function of the sound field of phased array can be written as: Sensors 2019, 19, 3192 5 of 15 Figure 3 shows a variation diagram of the phase directivity, and the slope of the phase directivity curve with dk. It can be seen from Figure 3 that when dk is very small, the phase directivity of phased array has a good linearity, and the slope of the phase directivity has a small amplitude fluctuation around 1. Then, with the increase in dk, the nonlinear fluctuation of the phase directivity of phased array increases gradually, and the fluctuation of the slope of the phase directivity increases gradually. The ratio between the sound pressure amplitude of a phased array sound field and the sound pressure amplitude of the sound field produced by a single transducer can be written as: Figure 4 shows a variation diagram of Ratio with , and it can be seen that, when is very small, the amplitude of the phased array sound field is very small, whereas the amplitude directivity of the sound field is not obvious. Then with the increase in , the amplitude of the phased array sound field increases, and the amplitude directivity of the sound field increases. Therefore, we can conclude that in the range of ∈ 0,2 , when increases, the amplitude of the phased array sound field increases, the nonlinear fluctuation of phase directivity of the sound field increases, the phase directivity increases with the fluctuation of the azimuth rate, and the amplitude directivity becomes higher directivity rather than omni-directivity. The ratio between the sound pressure amplitude of a phased array sound field and the sound pressure amplitude of the sound field produced by a single transducer can be written as: Figure 4 shows a variation diagram of Ratio 2 (ϕ) with dk, and it can be seen that, when dk is very small, the amplitude of the phased array sound field is very small, whereas the amplitude directivity of the sound field is not obvious. Then with the increase in dk, the amplitude of the phased array sound field increases, and the amplitude directivity of the sound field increases. The ratio between the sound pressure amplitude of a phased array sound field and the sound pressure amplitude of the sound field produced by a single transducer can be written as: Figure 4 shows a variation diagram of Ratio with , and it can be seen that, when is very small, the amplitude of the phased array sound field is very small, whereas the amplitude directivity of the sound field is not obvious. Then with the increase in , the amplitude of the phased array sound field increases, and the amplitude directivity of the sound field increases. Therefore, we can conclude that in the range of ∈ 0,2 , when increases, the amplitude of the phased array sound field increases, the nonlinear fluctuation of phase directivity of the sound field increases, the phase directivity increases with the fluctuation of the azimuth rate, and the amplitude directivity becomes higher directivity rather than omni-directivity. Therefore, we can conclude that in the range of dk ∈ (0, 2), when dk increases, the amplitude of the phased array sound field increases, the nonlinear fluctuation of phase directivity of the sound field increases, the phase directivity increases with the fluctuation of the azimuth rate, and the amplitude directivity becomes higher directivity rather than omni-directivity.

Finite Element Simulation of Phased Array
The model of phased array in water was created in the finite element software COMSOL, to analyze the sound field produced by phased array. The simulation runs on a Windows PC. Figure 5 shows the finite element model of phases array.

Finite Element Simulation of Phased Array
The model of phased array in water was created in the finite element software COMSOL, to analyze the sound field produced by phased array. The simulation runs on a Windows PC. Figure 5 shows the finite element model of phases array. The monopole domain sources parameters were Q are 1, e , e respectively, for the three monopole sources. By adjusting the frequency parameter in the Study node, the sound field of the phased array with different was analyzed. Figure 6 and Table 1 compare the simulation result of phase directivity and amplitude directivity with the theoretical calculation results. The Pressure Acoustic, Frequency Domain Interface of COMSOL was used in the simulation. The monopole domain sources parameters were Q m are 1, e 2 3 iπ , e 4 3 iπ respectively, for the three monopole sources. By adjusting the frequency parameter in the Study node, the sound field of the phased array with different dk was analyzed. Figure 6 and Table 1 compare the simulation result of phase directivity and amplitude directivity with the theoretical calculation results.

Finite Element Simulation of Phased Array
The model of phased array in water was created in the finite element software COMSOL, to analyze the sound field produced by phased array. The simulation runs on a Windows PC. Figure 5 shows the finite element model of phases array. respectively, for the three monopole sources. By adjusting the frequency parameter in the Study node, the sound field of the phased array with different was analyzed. Figure 6 and Table 1 compare the simulation result of phase directivity and amplitude directivity  It can be observed from Figure 6 that when = 0.73, the theoretical calculation of the amplitude range in each direction and the finite element simulation analysis value are 3.17 dB and 3.22 dB respectively, with an error of 0.05 dB, and the theoretical calculation and finite element simulation of the phase directivity fluctuation are 20.39° and 20.82°respectively, with an error of 0.1%.
This shows that the theoretically-derived results can directly be used to describe the amplitude directivity and phase directivity of the phased array in water.

Basic Structure of Spiral Sound Source
From the theoretical and finite element simulation results, we can see that the amplitude of sound pressure can be increased with the increase in the product of array radius and wave number, with a cost of increase in phase directivity fluctuation, and increase of the amplitude directivity of the sound field. Given that the phase directivity fluctuation of the sound field is not very large, to keep the response of the sound source at a decent level, we designed the spiral wave sound source based on phased array with three transducers.
The spiral sound source consists of three transducers and a support. The transducer is a spherical piezoelectric ceramic transducer used for emitting omni-directional spherical waves. The array element transducers are fixed at three vertexes of an equilateral triangle by the support. Three signals are connected to three array element transducers respectively, for driving the sound source with the same frequency and same amplitude but phase difference of 2 3 , respectively.

Finite Element Simulation of the Spiral Sound Source
The commercial finite element software COMSOL was used to establish the finite element model of the transducer [14][15][16]. COMSOL software was used to analyze the emission performance of the spiral sound source. Different from the simulation in 2.3, this simulation uses The Acoustic-Piezoelectric Interaction, Frequency Domain Interface, to simulate the spiral wave sound source. As   This shows that the theoretically-derived results can directly be used to describe the amplitude directivity and phase directivity of the phased array in water.

Basic Structure of Spiral Sound Source
From the theoretical and finite element simulation results, we can see that the amplitude of sound pressure can be increased with the increase in the product of array radius and wave number, with a cost of increase in phase directivity fluctuation, and increase of the amplitude directivity of the sound field. Given that the phase directivity fluctuation of the sound field is not very large, to keep the response of the sound source at a decent level, we designed the spiral wave sound source based on phased array with three transducers.
The spiral sound source consists of three transducers and a support. The transducer is a spherical piezoelectric ceramic transducer used for emitting omni-directional spherical waves. The array element transducers are fixed at three vertexes of an equilateral triangle by the support. Three signals are connected to three array element transducers respectively, for driving the sound source with the same frequency and same amplitude but phase difference of 2 3 π, respectively.

Finite Element Simulation of the Spiral Sound Source
The commercial finite element software COMSOL was used to establish the finite element model of the transducer [14][15][16]. COMSOL software was used to analyze the emission performance of the spiral sound source. Different from the simulation in 2.3, this simulation uses The Acoustic-Piezoelectric Interaction, Frequency Domain Interface, to simulate the spiral wave sound source. As the structure of the spiral sound source is symmetrical, a half model of the sound source was adopted to carry out the emitting simulation of a spiral sound wave in water. PZT-5J material was chosen for the piezoelectric ceramic spherical shell of the array element transducers with a diameter of 30 mm and thickness of 1mm. The influence of the support and other structures on the spiral sound source was ignored. Figure 7 shows the finite element model of the spiral sound source. Piezoelectric Interaction, Frequency Domain Interface, to simulate the spiral wave sound source. As the structure of the spiral sound source is symmetrical, a half model of the sound source was adopted to carry out the emitting simulation of a spiral sound wave in water. PZT-5J material was chosen for the piezoelectric ceramic spherical shell of the array element transducers with a diameter of 30 mm and thickness of 1mm. The influence of the support and other structures on the spiral sound source was ignored. Figure 7 shows the finite element model of the spiral sound source. The paper studies the low-frequency emission characteristics for the spiral sound source, because is very large at resonance, causing series fluctuation of the phase directivity. The transmitting voltage response curve for the spiral sound source is shown in Figure 8. When the frequency is 3.5 kHz, the transmitting voltage response is 108.9 dB, and the transmitting voltage response of a single transducer is 106.8 dB. The simulation curve of the rate of change of the phase directivity of the sound source with the azimuth is shown in Figure 9. It can be seen that the rate of change of the phase directivity with the azimuth varies from 0.53 to 1.70, and is always greater than 0, indicating a continuous monotone variation of the phase directivity with the azimuth. The paper studies the low-frequency emission characteristics for the spiral sound source, because dk is very large at resonance, causing series fluctuation of the phase directivity. The transmitting voltage response curve for the spiral sound source is shown in Figure 8. When the frequency is 3.5 kHz, the transmitting voltage response is 108.9 dB, and the transmitting voltage response of a single transducer is 106.8 dB. the structure of the spiral sound source is symmetrical, a half model of the sound source was adopted to carry out the emitting simulation of a spiral sound wave in water. PZT-5J material was chosen for the piezoelectric ceramic spherical shell of the array element transducers with a diameter of 30 mm and thickness of 1mm. The influence of the support and other structures on the spiral sound source was ignored. Figure 7 shows the finite element model of the spiral sound source. The paper studies the low-frequency emission characteristics for the spiral sound source, because is very large at resonance, causing series fluctuation of the phase directivity. The transmitting voltage response curve for the spiral sound source is shown in Figure 8. When the frequency is 3.5 kHz, the transmitting voltage response is 108.9 dB, and the transmitting voltage response of a single transducer is 106.8 dB. The simulation curve of the rate of change of the phase directivity of the sound source with the azimuth is shown in Figure 9. It can be seen that the rate of change of the phase directivity with the azimuth varies from 0.53 to 1.70, and is always greater than 0, indicating a continuous monotone variation of the phase directivity with the azimuth. The simulation curve of the rate of change of the phase directivity of the sound source with the azimuth is shown in Figure 9. It can be seen that the rate of change of the phase directivity with the azimuth varies from 0.53 to 1.70, and is always greater than 0, indicating a continuous monotone variation of the phase directivity with the azimuth. The finite element simulation results of the phase directivity for the sound source is shown in Figure 10. It can be seen that as the azimuth is increased by 360°, the phase of the sound field is also increased by 360°, indicating that the sound field of the spiral sound source is a spiral sound field, and the nonlinear fluctuation of the phase directivity is ±16°. The finite element simulation result of the amplitude directivity for the sound source is shown in Figure 11. It can be seen that the amplitude directivity of the sound source is periodic with a range of 3.4 dB. The finite element simulation results of the phase directivity for the sound source is shown in Figure 10. It can be seen that as the azimuth is increased by 360 • , the phase of the sound field is also increased by 360 • , indicating that the sound field of the spiral sound source is a spiral sound field, and the nonlinear fluctuation of the phase directivity is ±16 • . The finite element simulation results of the phase directivity for the sound source is shown in Figure 10. It can be seen that as the azimuth is increased by 360°, the phase of the sound field is also increased by 360°, indicating that the sound field of the spiral sound source is a spiral sound field, and the nonlinear fluctuation of the phase directivity is ±16°. The finite element simulation result of the amplitude directivity for the sound source is shown in Figure 11. It can be seen that the amplitude directivity of the sound source is periodic with a range of 3.4 dB. The finite element simulation result of the amplitude directivity for the sound source is shown in Figure 11. It can be seen that the amplitude directivity of the sound source is periodic with a range of 3.4 dB.

TVR(dB)
According to the finite element simulation analysis, we conclude that the spiral wave sound source based on phased array with three transducers can emit the spiral sound wave where the phase continuously varies with the azimuth, such that when the azimuth varies by 360 • , the phase of the sound field also varies by 360 • . Sensors 2018, 18, x FOR PEER REVIEW 10 of 15 Figure 11. Finite element simulation results of the amplitude directivity for the spiral sound source.
According to the finite element simulation analysis, we conclude that the spiral wave sound source based on phased array with three transducers can emit the spiral sound wave where the phase continuously varies with the azimuth, such that when the azimuth varies by 360°, the phase of the sound field also varies by 360°.

Measuring Principle of Phase Directivity
If the sound source is placed at the origin of coordinates, the two hydrophones are placed at , 0, and , 0, respectively, then the azimuth difference can be written as: If the two hydrophones collect sound pressure signals at the same time, the signals collected are given as: = . (17) The phases of two signals at = 0 are then: The phase difference of the output signal between the two hydrophones can be written as: The sound field can be regarded as a spiral wave field if the experimental results prove that the phase difference ∆ between the output signals of two hydrophones equals the azimuth difference ∆ of two hydrophones.
To calculate the curve of the ratio KL with the azimuth, we first rotated the sound source, and then measured the phase difference ∆ between the output signals of two hydrophones in different directions. By integrating KL to , the phase directivity curve of the measured sound field was obtained.

Measuring Principle of Phase Directivity
If the sound source is placed at the origin of coordinates, the two hydrophones are placed at x 1 (r 0 , 0, ϕ 1 ) and x 2 (r 0 , 0, ϕ 2 ) respectively, then the azimuth difference can be written as: If the two hydrophones collect sound pressure signals at the same time, the signals collected are given as: The phases of two signals at t = 0 are then: The phase difference of the output signal between the two hydrophones can be written as: The sound field can be regarded as a spiral wave field if the experimental results prove that the phase difference ∆ϕ 2 between the output signals of two hydrophones equals the azimuth difference ∆ϕ 1 of two hydrophones.
To calculate the curve of the ratio KL with the azimuth, we first rotated the sound source, and then measured the phase difference ∆ϕ 2 between the output signals of two hydrophones in different directions. By integrating KL to ϕ, the phase directivity curve of the measured sound field was obtained.

Test Results of the Spiral Sound Source
A prototype of the spiral wave sound source based on phased array with three transducers was designed and fabricated according to the finite element simulation results. Figure 12 is the picture of spherical transducers and spiral sound source prototype. The resonant frequency of the array element transducer of the spiral sound source is 43 kHz, and the diameter of the transducer is 37.8 mm. Three transducers are fixed by the support which enables the three transducers to be both isometrically and circumferentially distributed, and the distance between the center of the support chuck and the center of the support is 50 mm.

Test Results of the Spiral Sound Source
A prototype of the spiral wave sound source based on phased array with three transducers was designed and fabricated according to the finite element simulation results. Figure 12 is the picture of spherical transducers and spiral sound source prototype. The resonant frequency of the array element transducer of the spiral sound source is 43 kHz, and the diameter of the transducer is 37.8 mm. Three transducers are fixed by the support which enables the three transducers to be both isometrically and circumferentially distributed, and the distance between the center of the support chuck and the center of the support is 50 mm.  Figure 13 shows the placement of the spiral wave sound source and hydrophones in the pool and equipment connection. The fabricated spiral wave sound source was tested in a reverberation pool. By measuring the direct pulse that the hydrophones received, the free field performance of the spiral wave sound source was tested. Two signal sources outputs three pulsed signals, where the phase difference of each signal is realized by triggering time delay. All the three signals were amplified to 200 Vpp through 3 power amplifiers and connected to the transducers. The frequency of the signal is = 3.5 kHz, and time delays are 0, , , respectively. The models of equipment involved are shown in Table 2.   Figure 13 shows the placement of the spiral wave sound source and hydrophones in the pool and equipment connection. The fabricated spiral wave sound source was tested in a reverberation pool. By measuring the direct pulse that the hydrophones received, the free field performance of the spiral wave sound source was tested. Two signal sources outputs three pulsed signals, where the phase difference of each signal is realized by triggering time delay. All the three signals were amplified to 200 Vpp through 3 power amplifiers and connected to the transducers. The frequency of the signal is f = 3.5 kHz, and time delays are 0, 1 3 f , 2 3 f , respectively. The models of equipment involved are shown in Table 2.

Test Results of the Spiral Sound Source
A prototype of the spiral wave sound source based on phased array with three transducers was designed and fabricated according to the finite element simulation results. Figure 12 is the picture of spherical transducers and spiral sound source prototype. The resonant frequency of the array element transducer of the spiral sound source is 43 kHz, and the diameter of the transducer is 37.8 mm. Three transducers are fixed by the support which enables the three transducers to be both isometrically and circumferentially distributed, and the distance between the center of the support chuck and the center of the support is 50 mm.  Figure 13 shows the placement of the spiral wave sound source and hydrophones in the pool and equipment connection. The fabricated spiral wave sound source was tested in a reverberation pool. By measuring the direct pulse that the hydrophones received, the free field performance of the spiral wave sound source was tested. Two signal sources outputs three pulsed signals, where the phase difference of each signal is realized by triggering time delay. All the three signals were amplified to 200 Vpp through 3 power amplifiers and connected to the transducers. The frequency of the signal is = 3.5 kHz, and time delays are 0, , , respectively. The models of equipment involved are shown in Table 2.   The experimental result of the transmitting voltage response of array element transducer is shown in Figure 14. The transmitting voltage responses of three array element transducers at a frequency of 3.5 kHz are 106.5 dB, 105.9 dB, 105.3 dB respectively, with a range of 1.2 dB, indicating good uniformity and conforming the simulation result.  The experimental result of the transmitting voltage response of array element transducer is shown in Figure 14. The transmitting voltage responses of three array element transducers at a frequency of 3.5kHz are 106.5 dB, 105.9 dB, 105.3 dB respectively, with a range of 1.2 dB, indicating good uniformity and conforming the simulation result. The experimental results of the horizontal phase directivity slope KL-curve of the spiral sound source are shown in Figure 15. KL fluctuates between 0.35 and 2.11 and is constantly above 0, indicating a monotone increase of the sound field phase with the increase of the phase. Compared with the simulation, we can see that the change trend of the slope of the phase directivity with the azimuth is basically in agreement with the simulation results. The test results of the horizontal phase directivity curve of the spiral sound source are shown in Figure 16. It can be seen that when the azimuth is increased by 360°, the phase of the sound field is The experimental results of the horizontal phase directivity slope KL-curve of the spiral sound source are shown in Figure 15. KL fluctuates between 0.35 and 2.11 and is constantly above 0, indicating a monotone increase of the sound field phase with the increase of the phase. Compared with the simulation, we can see that the change trend of the slope of the phase directivity with the azimuth is basically in agreement with the simulation results.  The experimental result of the transmitting voltage response of array element transducer is shown in Figure 14. The transmitting voltage responses of three array element transducers at a frequency of 3.5kHz are 106.5 dB, 105.9 dB, 105.3 dB respectively, with a range of 1.2 dB, indicating good uniformity and conforming the simulation result. The experimental results of the horizontal phase directivity slope KL-curve of the spiral sound source are shown in Figure 15. KL fluctuates between 0.35 and 2.11 and is constantly above 0, indicating a monotone increase of the sound field phase with the increase of the phase. Compared with the simulation, we can see that the change trend of the slope of the phase directivity with the azimuth is basically in agreement with the simulation results. The test results of the horizontal phase directivity curve of the spiral sound source are shown in Figure 16. It can be seen that when the azimuth is increased by 360°, the phase of the sound field is The test results of the horizontal phase directivity curve of the spiral sound source are shown in Figure 16. It can be seen that when the azimuth is increased by 360 • , the phase of the sound field is also increased by 360 • , indicating that the sound field of the spiral sound source is a spiral sound field, and the nonlinear fluctuation measurement value of phase directivity is ±21 • . also increased by 360°, indicating that the sound field of the spiral sound source is a spiral sound field, and the nonlinear fluctuation measurement value of phase directivity is ±21°. The test results of the horizontal amplitude directivity diagram of the spiral sound source are shown in Figure 17, it can be seen that the range of amplitude directivity in all directions is 4.3 dB. The test results show that the spiral wave sound source based on phased array with three transducers can emit spiral sound waves with the phase continuously varying with the azimuth, and the measured results are in agreement with the simulation results.

Conclusions
In this paper, the phased array with three omni-directional transducers was researched to apply in the emission of underwater spiral sound waves as a simple realization form of underwater spiral sound source. According to the superposition theory of sound field, the mathematical derivation of The test results of the horizontal amplitude directivity diagram of the spiral sound source are shown in Figure 17, it can be seen that the range of amplitude directivity in all directions is 4.3 dB. also increased by 360°, indicating that the sound field of the spiral sound source is a spiral sound field, and the nonlinear fluctuation measurement value of phase directivity is ±21°. Figure 16. Test results of the horizontal phase directivity curve of the spiral sound source.
The test results of the horizontal amplitude directivity diagram of the spiral sound source are shown in Figure 17, it can be seen that the range of amplitude directivity in all directions is 4.3 dB. The test results show that the spiral wave sound source based on phased array with three transducers can emit spiral sound waves with the phase continuously varying with the azimuth, and the measured results are in agreement with the simulation results.

Conclusions
In this paper, the phased array with three omni-directional transducers was researched to apply in the emission of underwater spiral sound waves as a simple realization form of underwater spiral sound source. According to the superposition theory of sound field, the mathematical derivation of The test results show that the spiral wave sound source based on phased array with three transducers can emit spiral sound waves with the phase continuously varying with the azimuth, and the measured results are in agreement with the simulation results.

Conclusions
In this paper, the phased array with three omni-directional transducers was researched to apply in the emission of underwater spiral sound waves as a simple realization form of underwater spiral sound source. According to the superposition theory of sound field, the mathematical derivation of far-field sound pressure distribution of the sound field produced by phased array was given, the generation of spiral sound wave emitted by the phased array was obtained, and the influence of array parameters on phase directivity, amplitude directivity, and amplitude of the sound field were discussed. On this basis, the spiral wave sound source based on phased array which takes the spherical piezoelectric ceramic transducers as the array elements was analyzed with finite element software simulation. The sound source was made, and the sound field was measured.
The theoretical, simulation and experimental results show that the phased array with three transducers can emit spiral waves, the performance of the phased array was affected by the product of array radius and wave number. When the product of array radius and wave-number increases, the sound field amplitude reduces, and the phase directional fluctuation decreases. The phase of the sound field produced by the prototype of the spiral sound source increases continuously with the azimuth at 3.5 kHz. The phase directional fluctuation was ±21 • . The measurement results of amplitude directivity and phase directivity are in agreement with the theoretical and simulation results, verifying the correctness of the generation of spiral wave by the phased array, and the realization of the spiral wave sound source.