Optimization analysis of non-contact ultrasonic degradation system based on synchronous resonance of multiple modes

2.1. System constituents

Non-contact ultrasonic degradation system is mainly composed of an ultrasonic excitation module, a degradation module and an auxiliary module, the schematic illustration of the system is shown in Fig. 1.

Fig. 1Schematic illustration of the non-contact ultrasonic degradation system: 1 – Ultrasonic generator, 2 – ultrasonic transducer, 3 – ultrasonic horn, 4 – diaphragm, 5 – sample solution, 6 – cell, 7 – degradation vessel, 8 – altitude adjuster, 9 – base

The ultrasonic excitation module is composed of an ultrasonic generator, an ultrasonic transducer and an ultrasonic horn. The ultrasonic generator converts electric supply into high-frequency AC electrical signals, which are subsequently transformed into high-frequency mechanical vibration via a transducer [11, 12].

The degradation module is a specially processed reaction container, at the bottom of which, a diaphragm made of elastic material is placed. The diaphragm is capable of vibrating at a high frequency along the longitudinal direction. A preload force should be applied between the diaphragm and the ultrasonic horn to press the ultrasonic horn against the external surface of the diaphragm tightly, as a result of which, the ultrasonic energy is able to transmit into the sample solution, causing the degradation of the cell and the release of the target material inside the cell.

Auxiliary module mainly includes a height adjuster and a base, which not only supports the whole system but also adjusts the preload force between the ultrasonic horn and the diaphragm.

2.2. Key issue

The precondition that the non-contact ultrasonic degradation system effectively disrupts the cell is that the ultrasonic wave can smoothly pass through the diaphragm and transmit into the sample solution to be lysed as much as possible. The flexural vibration performance of the diaphragm is of great importance because it directly determines the efficiency of the non-contact ultrasonic degradation system.

Ultrasonic wave is a type of longitudinal mechanical vibration whose frequency is greater than 20 kHz, and it can transmit through the medium [13, 14]. To transmit as much ultrasonic energy into the sample solution as possible, the diaphragm should resonate with the ultrasonic horn, which means the natural frequency of the diaphragm should be approximately equal to the operating frequency of the ultrasonic horn. Hence, to ensure the natural frequency of diaphragm matches the operating frequency of the ultrasonic horn, both the structure of the diaphragm and the preload force between the diaphragm and the ultrasonic horn should be optimized; this optimization is the key issue to be solved in this paper.

3.1. Simplified model of the non-contact ultrasonic degradation system

Compared with a flat diaphragm, the diaphragm with dome-shaped structure has a greater natural frequency and is apt to vibrate. The advantage of the dome-shaped design of the diaphragm is that the dome shape increases the natural frequency of the diaphragm without causing the diaphragm to be so stiff that it cannot deflect to transmit the vibratory movements and the ultrasonic energy into the sample solution. Fig. 2 is a simplified model of the non-contact ultrasonic degradation system, whose curvature radius, width and thickness are represented by $R$, $W$ and $T$, respectively. In addition, the preload force between the diaphragm and the ultrasonic horn is represented by $F$.

Fig. 2Simplified model of the non-contact ultrasonic degradation system

3.2. Theoretical principle of flexural vibration

For the dome-shaped diaphragm, the differential equation of flexural vibration is given as [15]:

where $M$, $C$ and $K$ are the mass matrix, damping matrix and stiffness matrix of the diaphragm, respectively. $\ddot{X}$, $\dot{X}$ and $X$ are the acceleration vector, velocity vector and displacement vector, respectively, and $F$ is the excitation force vector.

It is the physical coordinates that the equation set of motion is described by in Eq. (1). The physical coordinates of all of the mass points of the diaphragm are included in each equation, as a result of which, there is much difficulty in solving the equation set directly. Therefore, modal analysis based on the finite element method (FEM) is used to calculate the natural frequency of the diaphragm.

3.3. Finite element model of the diaphragm

The natural frequency of each mode of the diaphragm can be calculated by finite element method. Take a dome-shaped diaphragm made of acetal for example, whose structural parameters and preload force are as follows: $R=$22 mm, $W=$12 mm, $T=$ 0.4 mm, $F=$3 N. And the solid model of the diaphragm is designed by SolidWorks 2013, as show in Fig. 3.

Then the solid model of the diaphragm is imported into ANSYS 11.0. The element type of the diaphragm is defined as the four-noded shell element ‘SHELL63’. ‘SHELL63’ has both bending and membrane capabilities. Both in-plane and normal loads are permitted. The element has six degrees of freedom at each node, translation is allowed in all coordinate directions as well as rotation around the coordinate axes. Stress stiffening and large deflection capabilities are included.

Fig. 3Solid model of the diaphragm

Fig. 4Finite element model of the diaphragm

Fig. 5Eigenshapes of the first ten modes of the diaphragm

a) The first mode

b) The second mode

c) The third mode

d) The fourth mode

e) The fifth mode

f) The sixth mode

g) The seventh mode

h) The eighth mode

i) The ninth mode

j) The tenth mode

In the step of meshing the model, the free meshing is used, and the size level is set to 6 by default. The finite element model, with a total of 1104 elements and 1106 nodes, is illustrated in Fig. 4. In the non-contact ultrasonic degradation system, the annular boundary part of the diaphragm is fixed on the degradation vessel, so all the DOFs of the annular boundary part’s external areas are set to be constrained.

The dome-shaped diaphragm is made of acetal, the material properties of which are shown in Table 1.

3.4. Natural frequencies and eigenshapes of the dome-shaped diaphragm

Based on ANSYS 11.0, static analysis of the dome-shaped diaphragm is performed, after which modal analysis is performed. Therefore, the natural frequencies and the eigenshapes of the first ten modes are extracted, as shown in Table 2 and Fig. 5.

3.5. Synchronous resonance of multiple modes

According to the data in Table 2, taking the order of mode as the abscissa and the natural frequency as the ordinate, the distribution of the first ten natural frequencies of the diaphragm is plotted, as shown in Fig. 6. The first ten natural frequencies of the dome-shaped diaphragm are found to have the characteristics of a ladder-like distribution. The frequencies can be divided into three gradients: the first gradient (represented by G1) is composed of M1 and M2; the second gradient (represented by G2) is composed of M3, M4, M5 and M6; and the third gradient (represented by G3) is composed of M7, M8, M9 and M10.

Fig. 6The distribution of the first ten natural frequencies of the diaphragm

For the application of the resonance, the characteristics of the ladder-like distribution of the first ten natural frequencies of the diaphragm has the benefit to enhance the vibration intensity. When the operating frequency of the ultrasonic horn is equal to the frequency of a certain gradient (here, the frequency of gradient is defined as the mean value of the natural frequencies of the modes that are included in the gradient), the synchronous resonance of multiple modes will occur, as a result of which, the vibration intensity of the diaphragm will be enhanced.

If the operating frequency of the ultrasonic horn is equal to the frequency of the first gradient, then the diaphragm will operate in the first gradient, which means there will be only two modes (M1 and M2) violently vibrating. If the operating frequency of the ultrasonic horn is equal to the frequency of the third gradient, then the diaphragm will operate in the third gradient, and there will be four modes (M7, M8, M9 and M10). However, the frequency of the third gradient is so high that the effect of vibration damping is great. As a consequence, the second gradient should be applied to transmit as much ultrasonic energy into the sample solution as possible. Here, the operating frequency is equal to the frequency of the second gradient, and the synchronous resonance of M3, M4, M5 and M6 will occur.

4.1. Optimization analysis of the curvature radius

The thickness $T$ of the diaphragm and the preload force $F$ between the ultrasonic horn and the diaphragm should be kept constant when the effect of the curvature radius $R$ on the natural frequency of each mode and their distribution characteristics is analyzed. Here, the thickness $T$ is set as 0.5 mm, and the preload force $F$ is set as 0 N.

The finite element modal analysis is conducted on the diaphragm having values of the curvature radius $R$ of 15 mm, 20 mm, 25 mm, 30 mm, 35 mm and 40 mm. As a result, the natural frequencies of the first ten mode of the diaphragm with different curvature radius is obtained, as listed in Table 3. In addition, the distribution of the first ten natural frequencies of the diaphragm is plotted, as shown in Fig. 7.

Fig. 7 shows that the first ten natural frequencies of the dome-shaped diaphragm with different curvature radius have characteristics of a ladder-like distribution, which is divided into three gradients: the first gradient (represented by G1) is composed of 1 and M2; the second gradient (represented by G2) is composed of M3, M4, M5 and M6; and the third gradient (represented by G3) is composed of M7, M8, M9 and M10.

According to the data in Table 3, the frequencies of the three gradients are calculated and then listed in Table 4. To evaluate the distribution characteristics of the first ten natural frequencies of the dome-shaped diaphragm with different curvature radius, an evaluation index-fluctuating frequency of gradient, is introduced. The fluctuating frequency of gradient is defined as the standard deviation of the natural frequencies of the modes that are included in the gradient. The fluctuating frequency of the gradient is capable of reflecting the deviation of the natural frequencies of the modes that are included in the gradient. Obviously, the smaller the fluctuating frequency of gradient is, the smaller the deviation of the natural frequencies is, which is more beneficial for synchronous resonance of multiple modes that are included in the gradient. According to the data in Table 3, the fluctuating frequencies of the three gradients are calculated, as listed in Table 4.

Fig. 7The distribution of the first ten natural frequencies of the diaphragm with different curvature radiuses

According to Table 4, the effects of the curvature radius on the frequencies and the fluctuating frequencies of the three gradients are illustrated in Figs. 8-9, respectively.

Fig. 8Effect of the curvature radius on the frequencies of the three gradients

Fig. 9Effect of the curvature radius on the frequencies of the three gradients

As stated in the preceding part of this paper, the diaphragm should operate in the second gradient (G2) to transmit as much ultrasonic energy into the sample solution as possible. From Figs. 8-9, it can be concluded that the frequency of the second gradient (G2) decrease gradually with the increasing curvature radius. Meanwhile, the fluctuating frequency of the second gradient (G2) reaches the minimum when the dome-shaped diaphragm has a curvature radius $R$ of 35 mm. Therefore, the curvature radius $R$ of the dome-shaped diaphragm is optimized to be 35 mm.

4.2. Optimization analysis of the thickness

The curvature radius $R$ of the diaphragm and the preload force $F$ should be kept constant when the effect of the thickness $T$ on the natural frequency of each mode and their distribution characteristics is analyzed. Here, the curvature radius $R$ is set as 35 mm, and the preload force $F$ is set as 0 N.

The finite element modal analysis is conducted on the diaphragm having values of the thickness $T$ of 0.3 mm, 0.4 mm, 0.5 mm, 0.6 mm, 0.7 mm and 0.8 mm. As a result, the natural frequencies of the first ten mode of the diaphragm with different thickness is obtained, as listed in Table 5. In addition, the distribution of the first ten natural frequencies of the diaphragm is plotted, as shown in Fig. 10.

Fig. 10The distribution of the first ten natural frequencies of the diaphragm with different thicknesses

Fig. 10 shows that the first ten natural frequencies of the dome-shaped diaphragm with different thickness have characteristics of a ladder-like distribution, which is divided into three gradients: the first gradient (represented by G1) is composed of M1 and M2; the second gradient (represented by G2) is composed of M3, M4, M5 and M6; and the third gradient (represented by G3) is composed of M7, M8, M9 and M10.

According to the data in Table 5, the frequencies and the fluctuating frequencies of the three gradients are calculated, as listed in Table 6.

According to Table 6, the effects of the thickness on the frequencies and the fluctuating frequencies of the three gradients are illustrated in Figs. 11-12, respectively.

Fig. 11Effect of the thickness on the frequencies of the three gradients

Fig. 12Effect of the thickness on the fluctuating frequencies of the three gradients

For the non-contact ultrasonic degradation system, the diaphragm should operate in the second gradient (G2) to transmit as much ultrasonic energy into the sample solution as possible, as stated above. Fig. 11 indicates that the thickness $T$ has a great effect on the frequency of the second gradient (G2). Further studies find that the frequency of the second gradient has a linear relationship with the thickness of the diaphragm. The linear-regression analysis is conducted using the least square method, and the equation of linear regression is obtained as follows:

where $f_{G2}\left(T\right)$ is the frequency of the second gradient (G2), $T$ is the thickness of the diaphragm. The square of correlation coefficient is $r^2=\text{0.99946}>\text{0.8}$, which proves $f_{G2}\left(T\right)$ is linearly proportional to $T$.

The equation of linear regression has a significant reference value in engineering application. On one hand, if the thickness of the diaphragm is known, then the frequency of the second gradient can be estimated. In this case, we can choose the ultrasonic horn with a certain operating frequency to match the diaphragm. On the other hand, if the operating frequency of the ultrasonic horn has been set, then the thickness of the diaphragm can be calculated to have the frequency of the second gradient match the operating frequency of the ultrasonic horn.

It can be concluded from Fig. 12 that, when the thickness $T$ of the diaphragm is in the range of 0.4 to 0.6 mm, the fluctuating frequency of the second gradient is relatively smaller. Hence, the thickness $T$ of the diaphragm should be designed to be in the range of 0.4 to 0.6 mm.

4.3. Optimization analysis of the preload force

The curvature radius $R$ and the thickness $T$ of the diaphragm should be kept constant when the effect of the preload force $F$ on the natural frequency of each mode and their distribution characteristics is analyzed. According to the results optimized above, the curvature radius $R$ and the thickness $T$ of the diaphragm are set as 35 mm and 0.5 mm, respectively.

The finite element modal analysis is conducted on the diaphragm under values of the preload force $F$ of 0 N, 2 N, 4 N, 6 N, 8 N and 10 N. As a result, the natural frequencies of the first ten mode of the diaphragm under different preload force is obtained, as listed in Table 7. Besides, the distribution of the first ten natural frequencies of the diaphragm is plotted, as shown in Fig. 13.

Fig. 13The distribution of the first ten natural frequencies of the diaphragm under different preload forces

From Table 7, when the preload force ranges between 0 N and 10 N, the frequency of each mode changes little, as a result of which, the six distribution curves in Fig. 13 are interlaced with each other. However, it is clear that the first ten natural frequencies of the dome-shaped diaphragm under different preload forces have characteristics of a ladder-like distribution, which can be divided into three gradients: the first gradient (represented by G1) is composed of M1 and M2; the second gradient (represented by G2) is composed of M3, M4, M5 and M6; and the third gradient (represented by G3) is composed of M7, M8, M9 and M10.

According to Table 7, the frequencies and the fluctuating frequencies of the three gradients are calculated, as listed in Table 8.

According to Table 8, the effects of the preload force on the frequencies and the fluctuating frequencies of the three gradients are illustrated in Figs. 14-15, respectively.

Fig. 14Effect of the preload force on the frequencies of the three gradients

Fig. 15Effect of the preload force on the fluctuating frequencies of the three gradients

The diaphragm should operate in the second gradient (G2), as stated above. It can be concluded from Fig. 14 that the preload force $F$ has a slight effect on the frequency of the second gradient (G2). Further studies indicate that the frequency of the second gradient has a linear relationship with the preload force between the diaphragm and the ultrasonic horn. The linear-regression analysis is conducted using the least square method, and the equation of linear regression is obtained as follows:

where $f_{G2}\left(F\right)$ is the frequency of the second gradient (G2), and $F$ is the preload force between the diaphragm and the ultrasonic horn. The square of the correlation coefficient is $r^2=\text{0.99949}>\text{0.8}$, which proves $f_{G2}\left(F\right)$ is linearly dependent on $F$.

The equation of linear regression also has a significant reference value in engineering application. If the structure of the diaphragm has been set, that is, the curvature radius $R$ and the thickness $T$ of the diaphragm have both been set, then the frequency of the second gradient can be fine-tuned to match the operating frequency of the ultrasonic horn through changing the value of the preload force $F$.

5.1. Experimental methods

The concentration of the raw sample solution of bacillus atrophaeus is 10⁶ cfu/mL. Then take 1500 µL and divide them into three parts on average, which are marked as Group A, B and C. That is to say the volume of each group is 500 µL. The main experimental procedure is illustrated in Fig. 17.

Fig. 17The main experimental procedure

1) Group A is degraded by the non-contact ultrasonic degradation system.

Group A is degraded by the experimental platform of the non-contact ultrasonic degradation system for 60 seconds. It’s impossible to count the amount of the bacillus atrophaeus that have been degraded by the experimental platform, while the amount of the bacillus atrophaeus that remain alive can be counted through plate cultivation. After degradation, the amount of bacillus atrophaeus that remain alive is represented by $N_{\mathrm{A}\left(\mathrm{1,500}\right)}$, here ‘1’ means the sample solution after degradation is not diluted, ‘500’ means the volume of the sample solution is 500 µL.

$N_{\mathrm{A}\left(\mathrm{1,500}\right)}$ cannot be obtained directly since there will be loss in volume after degradation. So 100 µL of the sample solution after degradation is taken to do the following experiments. And the amount of acillus atrophaeus that remain alive is represented by $N_{\mathrm{A}\left(\mathrm{1,100}\right)}$. Therefore, the relation between $N_{\mathrm{A}\left(\mathrm{1,500}\right)}$ and $N_{\mathrm{A}\left(\mathrm{1,100}\right)}$ is as follows:

However, $N_{\mathrm{A}\left(\mathrm{1,100}\right)}$ cannot be obtained directly through plate cultivation since the value of $N_{\mathrm{A}\left(\mathrm{1,100}\right)}$ is too huge to count after plate cultivation, which means dilution is a necessary step before plate cultivation. So the 100 µL of the sample solution after degradation is diluted, the dilution factor of which is 1:100. After dilution, take 150 µL and divided them into three parts on average, namely the volume of each part is 50 µL. Then three parts are inoculated on three agar plates independently, which are placed into the incubator at a temperature of 36 °C and then inverted. 24 hours later, the three agar plates are taken from the incubator, and then the colonies of bacillus atrophaeus are counted. The amount of bacillus atrophaeus in each agar plate is represented by $N_{\mathrm{A}\left(\mathrm{100,50}\right)}^{\text{'}}$, $N_{\mathrm{A}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}}$ and $N_{\mathrm{A}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}\text{'}}$, and the mean value of $N_{\mathrm{A}\left(\mathrm{100,50}\right)}^{\text{'}}$, $N_{\mathrm{A}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}}$ and $N_{\mathrm{A}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}\text{'}}$ is represented by ${\stackrel{-}{N}}_{\mathrm{A}\left(\mathrm{100,50}\right)}$. Here ‘100’ means the dilution factor is 1:100, and ‘50’ means the volume of bacillus atrophaeus solution inoculated on agar plate is 50 µL. Therefore, the relation between $N_{\mathrm{A}\left(\mathrm{1,100}\right)}$ and ${\stackrel{-}{N}}_{\mathrm{A}\left(\mathrm{100,50}\right)}$ is as follows:

2) Group B is degraded by the conventional ultrasonic processor.

Group B is degraded by Branson Sonifier^® S-250D digital ultrasonic processor for 60 seconds. After degradation, the amount of bacillus atrophaeus that remain alive is represented by $N_{\mathrm{B}\left(\mathrm{1,500}\right)}$, here ‘1’ means the sample solution after degradation is not diluted, ‘500’ means the volume of the sample solution is 500 µL.

$N_{\mathrm{B}\left(\mathrm{1,500}\right)}$ cannot be obtained directly. So 100 µL of the sample solution after degradation is taken to do the following experiments. And the amount of acillus atrophaeus that remain alive is represented by $N_{\mathrm{B}\left(\mathrm{1,100}\right)}$. Therefore, the relation between $N_{\mathrm{B}\left(\mathrm{1,500}\right)}$ and $N_{\mathrm{B}\left(\mathrm{1,100}\right)}$ is as follows:

$N_{\mathrm{B}\left(\mathrm{1,100}\right)}$ cannot be obtained directly through plate cultivation, and dilution is a necessary step before plate cultivation. So the 100 µL of the sample solution after degradation is diluted, the dilution factor of which is 1:100. After dilution, take 150 µL and divided them into three parts on average, namely the volume of each part is 50 µL. The three parts are inoculated on three agar plates independently, which are placed into the incubator at a temperature of 36 °C and then inverted. 24 hours later, the three agar plates are taken from the incubator, and then the colonies of bacillus atrophaeus are counted. The amount of bacillus atrophaeus in each agar plate is represented by $N_{\mathrm{B}\left(\mathrm{100,50}\right)}^{\text{'}}$, $N_{\mathrm{B}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}}$ and $N_{\mathrm{B}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}\text{'}}$, and the mean value of $N_{\mathrm{B}\left(\mathrm{100,50}\right)}^{\text{'}}$, $N_{\mathrm{B}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}}$ and $N_{\mathrm{B}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}\text{'}}$ is represented by ${\stackrel{-}{N}}_{\mathrm{B}\left(\mathrm{100,50}\right)}$. Here ‘100’ means the dilution factor is 1:100, and ‘50’ means the volume of bacillus atrophaeus solution inoculated on agar plate is 50 µL. Therefore, the relation between $N_{\mathrm{B}\left(\mathrm{1,100}\right)}$ and ${\stackrel{-}{N}}_{\mathrm{B}\left(\mathrm{100,50}\right)}$ is as follows:

3) Group C is not degraded.

Group C is not degraded, and the amount of bacillus atrophaeus in group C is represented by $N_{\mathrm{C}\left(\mathrm{1,500}\right)}$. Here ‘1’ means group C is not diluted, ‘500’ means the volume of group C is 500 µL. Obviously, $N_{\mathrm{C}\left(\mathrm{1,500}\right)}$ can be regarded as the amount of bacillus atrophaeus in group A or B before degradation.

Take 100 µL from group C, the amount of bacillus atrophaeus in it is represented by $N_{\mathrm{C}\left(\mathrm{1,100}\right)}$. The relation between $N_{\mathrm{C}\left(\mathrm{1,500}\right)}$ and $N_{\mathrm{C}\left(\mathrm{1,100}\right)}$ is as follows:

However, $N_{\mathrm{C}\left(\mathrm{1,100}\right)}$ cannot be obtained directly through plate cultivation, and dilution is a necessary step before plate cultivation. This time the dilution factor is 1:1000 because the value of $N_{\mathrm{C}\left(\mathrm{1,100}\right)}$ is huger. After dilution, take 150 µL and divided them into three parts on average, namely the volume of each part is 50 µL. The three parts are inoculated on three agar plates independently, which are placed into the incubator at a temperature of 36 °C and then inverted. 24 hours later, the three agar plates are taken from the incubator, and then the colonies of bacillus atrophaeus are counted. The amount of bacillus atrophaeus in each agar plate is represented by $N_{\mathrm{C}\left(\mathrm{100,50}\right)}^{\text{'}}$, $N_{\mathrm{C}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}}$ and $N_{\mathrm{C}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}\text{'}}$, and the mean value of $N_{\mathrm{C}\left(\mathrm{100,50}\right)}^{\text{'}}$, $N_{\mathrm{C}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}}$ and $N_{\mathrm{C}\left(\mathrm{100,50}\right)}^{\text{'}\text{'}\text{'}}$ is represented by ${\stackrel{-}{N}}_{\mathrm{C}\left(\mathrm{1000,50}\right)}$. Here ‘1000’ means the dilution factor is 1:1000, and ‘50’ means the volume of bacillus atrophaeus solution inoculated on agar plate is 50 µL. Therefore, the relation between $N_{\mathrm{C}\left(\mathrm{1,100}\right)}$ and ${\stackrel{-}{N}}_{\mathrm{C}\left(\mathrm{1000,50}\right)}$ is as follows:

4) Calculation of the degradation efficiency.

The degradation efficiency is defined as the ratio of the amount of the bacillus atrophaeus that has been degraded to the total amount of the bacillus atrophaeus without any treatment and is represented by $\gamma$. Therefore, the degradation efficiencies of the experimental platform and the Branson Sonifier® S-250D digital ultrasonic processor are as follows:

According to Eqs. (4)-(12), the degradation efficiencies can be given as:

5.2. Experimental results

The results are listed in Tables 9-11, from which the degradation efficiency can be calculated.

${\gamma }_{\mathrm{A}}=$86.7 %, ${\gamma }_{\mathrm{B}}=$82.5 %.

The results prove that the non-contact ultrasonic degradation system developed in this paper has the ability to degrade bacillus atrophaeus.