METHOD AND DEVICE FOR GENERATING SYNTHETIC VORTEX SOUND FIELD WITH MORE MODE NUMBER

Abstract

A method and device for generating synthetic vortex sound field (SVSF) with more mode number includes the following steps: (1) a transducer array composed of N transducer units is constructed, and each transducer unit emits a sound field to generate an initial sound field; (2) at the same time, the position of the transducer unit and the phase of the sound field emitted by each transducer unit are changed, and each change produces a sound field, and thus changings times produces of sound fields, wherein the way to change the position of the transducer unit is to rotate the transducer array as a whole; (3) the initial sound field is superimposed with s of sound fields generated in step (2), to obtain SVSF with more mode number. The method and device for generating vortex sound field (VSF) can be used for underwater communication or acoustic imaging.

Claims

1. A method for generating synthetic vortex sound field (SVSF) with more mode number, characterized in that it includes the following steps: (1) A transducer array composed of N transducer units is constructed, and each transducer unit emits a sound field to generate an initial sound field; (2) At the same time, the position of the transducer unit and the phase of the sound field emitted by each transducer unit are changed, and each change produces a sound field, and thus changing s times produces s of sound fields, wherein the way to change the position of the transducer unit is to rotate the transducer array as a whole; (3) The initial sound field generated in step (1) is superimposed with s of sound fields generated in step (2), to obtain SVSF with more mode number; wherein N and s are integers of >0, and N*s is not less than 4.

2. The method according to claim 1, characterized in that the transducer array forms a virtual synthetic transducer array before and after rotation, and the number of array elements in the synthetic transducer array is N.sub.s, N.sub.s=(s+1)×N.

3. The method according to claim 2, characterized in that the array elements of said synthetic transducer array are arranged on one ring or a concentric ring consisted of at least two rings.

4. The method according to claim 3, characterized in that the array elements of said synthetic transducer array are arranged on one ring, and the phase of the sound field generated by the m-th array element in the synthetic transducer array is: ${\alpha }^{\prime }*\frac{2\pi \left(m-1\right)}{N_s},$ wherein 1≤m≤N.sub.s, m is an integer, α′is the mode number of said SVSF, $-\frac{N_s}{2}<{\alpha }^{\prime }<\frac{N_s}{2}.$

5. The method according to claim 1, characterized in that in the transducer array, the transducer units are arranged on a ring, and the rotation axis of the transducer array is the axis of the ring.

6. The method according to claim 5, characterized in that the phase of the sound field generated by the nth transducer unit at the initial position is: ${\alpha }^{\prime }*\frac{2\pi \left(n-1\right)}{N},$ wherein 1≤n≤N, n is an integer, and α′is the mode number of said SVSF, $-\frac{N_s}{2}<{\alpha }^{\prime }<\frac{N_s}{2};$ and/or the angle for each rotation of the transducer array is $\frac{2\pi }{N_s},$ and thus after the nth transducer unit rotates for the ith time, the phase of the generated sound field is: ${\alpha }^{\prime }*\frac{2\pi \left(n-1\right)}{N}+{\alpha }^{\prime }*\frac{2\pi }{N_s}*i,$ wherein 1≤i≤s, 1≤n≤N, i and n are integers, and α′ is the mode number of SVSF, $-\frac{N_s}{2}<{\alpha }^{\prime }<\frac{N_s}{2}.$

7. (canceled)

8. The method according to claim 1 which is used in underwater communication or acoustic imaging.

9. A device generating a SVFD with more mode number, which implements the method according to claim 1, and characterized in that it includes a rotating device and at least one transducer array composed of transducer units, and the rotating device is used to drive the transducer array to rotate

10. The device according to claim 8, which is used in the equipment of underwater communication or acoustic imaging.

11. The method according to claim 3, characterized in that the array elements on each ring are evenly arranged.

12. The generation device according to claim 9, characterized in that in the transducer array, the transducer units are evenly arranged on a ring; the rotating shaft by which the rotating device drives the transducer array to rotate passes through the center of the ring formed by the arrangement of transducer units; the rotating device is a precision rotating table for accurately controlling each rotation angle of the transducer array.

Description

DESCRIPTION OF FIGURES

[0027] FIG. 1. The schematic diagram of a uniform circular transducer array;

[0028] FIG. 2. The vortex ultrasonic field with a mode number of 1 (a), 2 (b) and 3 (c) obtained by 8 uniform circular transducer arrays in the prior art; and a schematic diagram of 8 transducer units unable to generate a VSF with a mode number of 4 (E);

[0029] FIG. 3. A schematic diagram of the state of the basic array formed by the transducer unit before rotation and after two rotations, and the synthetic transducer array formed by them;

[0030] FIG. 4. A schematic diagram of a VSF with a mode number of 4 generated by rotating the transducer basic array with N=8 once in the example and simulating N.sub.s=16 transducer arrays;

[0031] FIG. 5. A schematic diagram of a VSF with a mode number of 8 generated by rotating the transducer basic array with N=8 twice in the example and simulating N.sub.s=24 transducer arrays;

[0032] FIG. 6. The directivity of a VSF with a mode number of 3 directly generated by the transducer basic array with N=8 in the example;

[0033] FIG. 7. The directivity of a VSF with a mode number of 3 synthesized by this method using the transducer basic array with N=8 in the example;

[0034] FIG. 8. The directivity of a VSF with a mode number of 4 generated by the transducer basic array with N=8 in the example;

[0035] FIG. 9. The directivity of a VSF with a mode number of 4 synthesized using this method by rotating the transducer basic array with N=8 once in the example.

EXAMPLES

[0036] The technical solution of the present application was further illustrated in the following by specific examples.

[0037] In the prior art, the method generating a VSF by a uniform circular transducer array was as follows: Assuming that the uniform circular array composed of N circular transducers is located in the XOY plane, As shown in FIG. 1, the array radius is R, and the modulation azimuth of the nth transducer (i.e. the phase of the generated sound field) was ϕ.sub.n=2π(n-1)α/N,α is the number of topological charges (i.e. the number of modes). An excitation signal was applied to each transducer:

s.sub.n=A*exp(j2πft+jϕ.sub.n)   (1-1)

[0038] where A is the amplitude of the sound wave, f is the signal frequency, t is the time, and j is the imaginary unit.

[0039] Supposing that the coordinate of observation point T in a rectangular coordinate system was (x, y, z), and its coordinate in a spherical polar coordinate system is(r, φ, θ) (wherein r is the distance between the observation point and the coordinate origin, φ is the angle of the line between the observation point and the origin of coordinate axis from X-axis, and θ is the angle of the line between the observation point and the origin of coordinate axis from Z-axis), the sound pressure detected at the observation point is:

[00008]$\begin{array}{cc}{p}_T=\frac{A}{R_n}\exp \left(2\pi \mathrm{ft}+j{\varphi }_n+{\mathrm{jkR}}_n\right)& \left(1-2\right)\end{array}$

[0040] wherein k is the wave number, φ.sub.n is the spatial azimuth of the transducer in the spherical coordinates, φ.sub.n=ϕ.sub.n/α, R.sub.n is the distance from any transducer to the observation point T, R.sub.n could be expressed as:

[00009]$\begin{array}{cc}{R}_n=\sqrt{{\left(r\sin \mathrm{\theta cos}\phi -R\cos {\phi }_n\right)}^2+{\left(r\sin \mathrm{\theta sin}\phi -R\sin {\phi }_n\right)}^2+r\cos {\theta }^2}& \left(1-3\right)\end{array}$

[0041] When N of transducers are superimposed, the sound pressure at the detection point T(r, φ, θ) could be expressed as:

[00010]$\begin{array}{cc}{p}_n=\exp \left(-j\omega t\right)\underset{n=1}{\overset{N}{.Math.}}\frac{A}{R_n}\exp \left(j{\varphi }_n+{\mathrm{jkR}}_n\right)& \left(1-4\right)\end{array}$

[0042] The complex exponential form of formula (1-4) was expanded into a trigonometric function form:

[00011]$\begin{array}{cc}\stackrel{`}{p_n}=\frac{A}{R_n}\cos \left[{\mathrm{kR}}_n+2\alpha \pi \frac{\left(n-1\right)}{N}\right]+j\frac{A}{R_n}\sin \left[{\mathrm{kR}}_n+2\alpha \pi \frac{\left(n-1\right)}{N}\right]& \left(1-5\right)\end{array}$

[0043] After multiple transducers are superimposed, the amplitude expression of the sound field was as follows:

[00012]$\begin{array}{cc}P=\sqrt{\underset{n=1}{\overset{N}{.Math.}}{\left\{A\cos \left[{\mathrm{kR}}_n+2\alpha \pi \frac{\left(n-1\right)}{N}\right]/R_n\right\}}^2+\underset{n=1}{\overset{N}{.Math.}}{\left\{A\sin \left[{\mathrm{kR}}_n+2\alpha \pi \frac{\left(n-1\right)}{N}\right]/R_n\right\}}^2}& \left(1-6\right)\end{array}$

[0044] The phase expression of the formed sound field was as follows:

[00013]$\begin{array}{cc}\varphi =\arctan \frac{{.Math.}_{n=1}^N\left\{A\sin \left[{\mathrm{kR}}_n+2\alpha \pi \frac{\left(n-1\right)}{N}\right]/R_n\right\}}{{.Math.}_{n=1}^N\left\{A\cos \left[{\mathrm{kR}}_n+2\alpha \pi \frac{\left(n-1\right)}{N}\right]/R_n\right\}}& \left(1-7\right)\end{array}$

[0045] The experimental parameters used were the frequency f=1000 Hz, the sound velocity c=340 m/s, the sound amplitude A=1, the number of array elements N=8, the number of modes α=1, 2, 3, 4, and the array radius R=0.2 m. The VSF obtained by formula (1-6) and formula (1-7) is shown in FIG. 2. According to FIG. 2, when α=4, the VSF could not be formed. The characteristic of VSF was that the central sound intensity was 0, and the wave front in the propagation direction was spiral. Its characteristics were originated from the phase distribution of wavefront linear variation.

Example

VSF for the SOAM with More Mode Number According to the Present Invention

[0046] In this example, the parameters were defined as follows:

[0047] The number of original transducer units was N;

[0048] The number of transducer elements in the synthetic transducer array was N.sub.s; N.sub.s=(s+1)×N;

[0049] The synthetic mode number was α′,α′ was an integar, and met the following requirements:

[00014]$-\frac{N_s}{2}<{\alpha }^{\prime }<\frac{N_s}{2};$

[0050] If there are N of original transducer units, the number of vortex field modes that could be formed was α, which was an integer, and met the following requirements:

[00015]$-\frac{N}{2}<\alpha <\frac{N}{2}.$

[0051] If the synthetic mode number was α′, the modulation phase difference between two adjacent transducer units in the synthetic transducer array was:

[00016]$\Delta {\varphi }_s=2\pi \frac{{\alpha }^{\prime }}{N_s}.$

[0052] The rotation number of the transducer array was recorded as s; Synthetic transducer array denotes the array formed by taking the position of each transducer unit as an array, when each transducer unit used to synthesize a VSF generated a sound field. For example, in the prior art, when the transducer array was not rotated, the synthetic transducer array was the original transducer array. If the transducer array was rotated once (as shown in FIG. 3), the synthetic transducer array was a combination of the original transducer array and the transducer array obtained after rotation (as shown in the right figure of FIG. 3).

[0053] Therefore, in order to obtain a greater α′, it is necessary to increase the number N.sub.s of transducer units in the synthetic transducer array. In the traditional method, the number N of original transducer units must be increased. For the present method, it was only needed to increase the rotation times s of the original transducer array.

[0054] In particular, the operation method of this example was:

[0055] (1) N transducer units were evenly distributed on the ring with radius of R, and the obtained annular transducer array was controlled by the precision rotary table, which could drive the annular transducer array to rotate in the set direction (clockwise or counterclockwise).

[0056] (2) If a virtual vortex sound field with a mode number of α′was synthesized, at the initial position, the phase of the sound field generated by the nth transducer unit was

[00017]${\alpha }^{\prime }*\frac{2\pi \left(n-1\right)}{N},$

in which

[00018]$-\frac{N_s}{2}<{\alpha }^{\prime }<\frac{N_s}{2}.$

[0057] (3) If a synthetic transducer array with N.sub.s of synthetic array elements was required, the annular transducer array needed to be rotated k-1 times, to allow N.sub.s=kN. The ring was controlled by the precision rotary table, which would drive the transducer array to rotate in the set direction (clockwise or counterclockwise). The angle of the ring transducer array rotated each time is:

[00019]$\frac{2\pi }{N_s}.$

After the transducer array was rotated i times (1≤i≤s), the phase of the sound field successively transmitted by the nth transducer was:

[00020]${\alpha }^{\prime }*\frac{2\pi \left(n-1\right)}{N}+{\alpha }^{\prime }*\frac{2\pi }{N_s}*i,$

in which

[00021]$-\frac{N_s}{2}<{\alpha }^{\prime }<\frac{N_s}{2}.$

[0058] (4) By superimposing the original sound fields with different mode numbers formed at different positions of the array, the VSF with SOAM mode numbers could be synthesized.

[0059] The method of “sound field superposition” was: the vectors of the expressions (or measured values) for the initial sound field generated in step (1) and s of sound fields generated in step (2) were added, to obtain a new expression (a measured value), and the sound field represented by the new expression (the measured value) was the superimposed sound field. The expression denoted the sound pressure expression of the detection point T(r, φ, θ).

[0060] As the operation of the above method, the result was shown in FIG. 4. The sound field with a mode number of 4 (equal to N/2) was generated by the transducer basic array with N=8, as shown in FIG. 4 panel a, and then the sound field with a mode number of 4 was generated by rotating the basic array once, as shown in FIG. 4 panel b. Both of generated sound fields were superimposed to simulate the VSF with a mode number of 4 generated by N.sub.s=16 transducer arrays, as shown in FIG. 4 panel c. As shown in FIG. 5, the sound field with a mode number of 8 was generated directly with the transducer basic array of N=24, as shown in FIG. 5 panel a. If the basic transducer array with N=8 was used to generate the synthetic vortex sound field with a mode number of 8, the array needed to be rotated twice. After the spatial position of the array was rotated, the phase needed to be changed. The initial sound field generated by the transducer basic array with N=8 was shown in FIG. 5 panel b; after rotating the transducer base array with N=8 once, and the phase of the sound field emitted by each transducer unit was correspondingly changed, the generated sound field was shown in FIG. 5 panel c; after rotating the transducer base array with N=8 again, and the phase of the sound field emitted by each transducer unit was correspondingly changed, the generated sound field was shown in FIG. 5 panel d; after superimposing the sound fields emitted by the array at three different spatial positions, the synthetic vortex sound field with a mode number of 8 was obtained as shown in FIG. 5 panel e. Using the above method, the synthetic vortex sound field with a mode number of 8 was generated by 8 transducer units. As shown, this method could generate VSF having more mode number with less transducer units. Other parameters of this example were consistent with those used in the above method of VSF generated by uniform circular transducer array.

[0061] In order to demonstrate the advantages of the present application, the directivity of the sound field generated by this example was described below. The directivity function of the circular transducer array used in this example is:

[00022]$\begin{array}{cc}{D}_s=\frac{1}{N}\underset{n=1}{\overset{N}{.Math.}}{e}^{-j\frac{2\pi f*R}{c}\sin \mathrm{\theta cos}\frac{2\pi *n}{N}}.Math.\frac{2J_{\alpha }\left(\frac{2\pi f}{c}a\sin \theta \right)}{\frac{2\pi f}{c}a\sin \theta }& \left(1-7\right)\end{array}$

[0062] wherein R is the array radius, c is the sound velocity, j is the imaginary unit, and a is the radius of the transducer unit.

[0063] FIG. 6 is the directivity of a VSF with a mode number of 3 directly generated by the transducer basic array with N=8; FIG. 7 is the directivity of a VSF with a mode number of 3 synthesized by this method using the transducer basic array with N=8 in the example;

[0064] FIG. 8 is the directivity of a VSF with a mode number of 4 generated by the transducer basic array with N=8 in the example; FIG. 9 is the directivity of a VSF with a mode number of 4 synthesized using this method by rotating the transducer basic array with N=8 once in the example.

[0065] Obviously, by comparison of FIG. 6 and FIG. 7 as well as by directivity. Therefore, it had better imaging resolution and better transmission performance in the process of imaging and data transmission.

[0066] From the above example, it could be shown that the present application could synthesize a VSF with more mode number by rotating the transducer array composed of less transducer units, adjusting the phase of each transducer unit, and superimposing the VSF generated after each rotation with that generated before rotation. Compared with the prior art, the synthetic vortex ultrasonic field generated by the method of the present invention had better directivity. Applying this method to underwater communication, biomedical imaging and other equipment could reduce the number of transducer units and thus simplify the equipment. The information carrying capacity and imaging resolution could be increased as the increase of the mode number of VSF; the enhancement of directivity also made it have better imaging resolution and better transmission performance in the process of imaging and data transmission. Therefore, the application potential of the technology of the present invention was significant.