Harmonic Response Analysis of Tank Design Effect on Ultrasonic Cleaning Process

: Several ultrasonic cleaning tanks (UCTs) had a problem: a manufacturer complained that there were damages to cleaning objects, they were unclariﬁed, and it needed to be abruptly solved. To investigate and solve the problem, a small UCT ﬁlled with 3.92 L of water, with a frequency of 28 kHz, two horn transducers, and a total power of 100 W was built for simulation and experiment. A built tank body material of UCT can be adjustable to acrylic, glass, and stainless steel. Since the cavitation causing the cleaning relates to acoustic pressure, harmonic response analysis (HRA) in ANSYS software was employed to calculate the acoustic pressure inside the UCT for different designs such as mentioned materials, power, thickness, volume, and frequency. The HRA results revealed uneven acoustic pressure depending on the tank designs, consistent with foil corrosion and power concentration experiments. Furthermore, using the tank body material with acrylic, glass, and stainless steel provided the highest, moderate, and lowest acoustic pressure levels, respectively. The uneven acoustic pressure resulted from the differences in material transmission coefﬁcients. In addition, the damage occurred because of improper tank design, resulting in excessive acoustic pressure. Therefore, the tank design is indispensable in designing high-efﬁciency UCTs to reduce damage and meet customer requirements.


Introduction
Thailand is one of the technology leaders in manufacturing ultrasonic cleaners in Southeast Asia. Many manufacturers are capable of producing ultrasonic cleaners ranging from simple functioning devices for household usage to large machines with complex systems for industrial usage. The advantages of domestic manufacturers are that they can design and produce ultrasonic cleaners and components that match customers' requirements at a low price. In addition, they can quickly solve any usage problems complained by the customers and suggest ways to improve the cleaning performance. Behind the success of leadership is using finite element analysis (FEA). Research on the development of ultrasonic cleaning based on FEA was divided into two groups: cleaners and processes.
First, since the cleaners consist of three essential components, i.e., transducer, tank body, and controller (generator), most research has focused on developing the first two rather than the last because transducer and tank body developments can lead to high-efficiency, rapid cleaning, meet the requirements of a wide range of users, and are worth the investment. Meanwhile, the controller development is relatively slow because the existing technology can be practical on an industrial scale. FEA to develop the transducer and tank body impacting the industries can be exemplified as follows. The structural analysis [1,2], modal analysis [3][4][5][6], and harmonic response analysis (HRA) [6] were employed to investigate stress, deformation, and related results for optimizing the design of front masses for transducers to suit usages. Furthermore, since the transducer consists of front mass, piezoelectric (PZT) material, and back mass, the PZT is also studied and developed for many purposes [7]. As a result, HRA confirmed that PZT4 provided more intense acoustic pressure than PZT8 [8,9].
Second, for the cleaning process development, the structural analysis helped to find a suitable fouling removal method [10,11]. In addition, the HRA determined suitable positions of cleaning objects [12] and optimized transducer placements [13][14][15] to achieve higher ultrasonic cleaning efficiency.
Because ultrasonic cleaning efficiency depends on many factors such as cleaning fluid, temperature, standing waves, power, and frequency [16,17], optimized design concerning these factors ensures optimum efficiency. Therefore, from the literature review mentioned above and the authors' experience, it can be concluded that the HRA is more suitable than the structural analysis and modal analysis to achieve higher ultrasonic cleaning efficiency. For example, the HRA can display more accurate acoustic pressure, leading to fluid cavitation, when adjusting the related factors than the other analyses. In addition, the structural and modal analyses require boundary conditions of external sources with specific frequencies. Still, the HRA can report the results at all frequencies, does not require boundary conditions from external sources, and is easy to use. Therefore, most of the recent research prefers to use the HRA to help design and develop ultrasonic cleaning technology [8,9,[12][13][14][15].
The origin of this research came from customers complaining to manufacturers. They informed us that the cleaning objects were damaged and not clean in some ultrasonic cleaning tanks (UCTs). In fact, repositioning the transducers to the proper position, as suggested in previous research [13,14], is already a potential solution. However, the mentioned method cannot identify the cause of this problem; therefore, finding the root cause, alternative solution, and ways to improve efficiency should be urgently addressed and were aims of this research. Preliminary investigation revealed that a material typically used to invent the tank body could cause the problem.
An actual UCT in a laboratory was built to help solve the problem. The tank body material could be adjusted to typical material used in the cleaning process, such as acrylic, glass, and stainless steel. Using water as a cleaning fluid and three mentioned materials as variables, the HRA was employed to determine the acoustic pressure inside the UCT. The HRA simulation results were compared with the experimental results obtained from foil corrosion and power concentration experiments to confirm reliability and methodology. Next, the HRA was applied to calculate acoustic pressure by varying tank designs such as tank body material, thickness, power, volume, and frequency. All simulation results were analyzed to find the root cause and a solution to achieve high-efficiency UCT. The highlight is that boundary conditions and settings in this research are based on actual cleaning processes in the factory. In addition, the built UCT was exclusively designed for this research. Finally, the research methodology was explained step by step.

Harmonic Response Analysis (HRA) and Ultrasonic Cleaning
The HRA is a linear dynamic analysis employed to determine the reaction of materials responding to sinusoidal wave loads in the steady state. Recently, it has been widely applied in developing the UCTs [8,9,[12][13][14][15] and is briefly described in this section.
The HRA can be calculated in ANSYS software by dividing domains of the UCT materials into interfaces, water, PZT4, aluminum alloy, stainless steel, and tank body [9,14,15], as shown in Figure 1. However, the division is that equations governing the domains are different. When alternating current is applied to PZT4 inside the transducer, the PZT4 vibrates with a natural frequency. Then, the vibration is transferred through different components of UCT until entering the water and resulting in acoustic pressure.
In the PZT4 domains and PZT4/solid interface, {u} can be calculated from a coupling between structural and electrical responses as  p} are the first and second derivatives of nodal acoustic pressure, respectively.
In the water domain, {p} is an acoustic pressure of wave propagating through a medium and has a time-harmonic relation as p = p 0 e iωt , which can be calculated in the HRA of ANSYS software based on the Navier-Stokes equations by where ω is the angular frequency (driving frequency), [M] is the acoustic mass matrix, [C] is the acoustic damping matrix, [K] is the acoustic stiffness matrix, and {F} is the acoustic load vector. The subscript F denotes the fluid, which was water in this research. The {p} is a crucial result representing the cavitation effect of the cleaning process. The larger the value of {p}, the greater the intensity of cavitation but does not imply the bubble size. The bubble size can be estimated from the ω. The higher the value of ω, the smaller the cavitating bubble size [9,18]. For the most efficient cleaning, the bubble size should be close to the size of the contaminant and dirt. Generally, factories will start the cleaning process at 28 kHz. This frequency gives a bubble size of approximately 14 µm, which can be used to clean many types of dirt, microorganisms, grease, and dust, etc. The mixing frequency produces bubbles of various sizes, which is a benefit and helps to clean more kinds of dirt [15]. Therefore, if accurate HRA simulation results are obtained, it will lead to high-efficiency UCTs matching customers' requirements. The description above is a brief example of using the HRA to develop the cleaning process.

Reflection, Transmission, and Vibration, in the UCT
As mentioned in Section 1, the cause of the problem may be the material to make the tank body. In physics, cavitation caused by vibration, reflection, and transmission recurring inside the UCT can be explained in Figure 2.
where ρ is the material density and c is the sound speed. The Z relates reflection (α R ) and transmission (α T ) coefficients as in Equations (6) and (7) In order to link the physics in Figure 1 to the cavitation, the transducer in this figure consists of materials: PZT4 as a vibrator, aluminum alloy as a front mass, and stainless steel as a back mass. The tank body contains water and can be adjusted to materials such as stainless steel, acrylic, and glass, depending on usage purposes. Therefore, when the PZT4 vibrates, it causes ultrasonic waves to pass through the water, producing reflected and transmitted waves and cavitation. Generally, the front mass is designed with a less dense material to allow easy wave transmission to the front, while the back mass is designed with a dense material to reflect the waves from the back to the front, and the PZT is designed to match the vibration of the front mass shape. As a result, changing the tank designs such as tank body material, power, volume, and frequency alters the variables in Equations (1)- (7), resulting in different {p}. Therefore, the {p} is a key result to be calculated.

Methodology
The research methodology is shown in Figure 3. The yellow and blue colors are for the HRA simulation and the experiment, respectively. It can be described as follows.

An Experimental UCT
Since there are no UCTs from the manufacturers whose properties correspond to the aims of this research, a prototype UCT in a laboratory was designed and built. To understand the functioning of built UCT, two parts needed to be mentioned: an electric circuit and a tank body. Figure 4 reveals (a) diagrams of an electric circuit and (b) a CAD model of UCT, specifically modified and built, respectively, supporting this research. The functions can be summarized as follows. Firstly, an electric circuit for 40-60 kHz designed for an ultrasonic generator presented in [20] was employed as a prototype and adjusted to suit 28 kHz, as shown in Figure 4a. The 28 kHz is a natural frequency of the horn transducer generating a vibration with longitudinal amplitude shown in the direction of Figure 1, enhancing the cleaning efficiency. On the contrary, other natural frequencies produce a vibration with less amplitude in the direction obstructing the cleaning efficiency.
Secondly, the two transducers were attached to a rectangular stainless steel base placed a supporting structure. The transducer details are in Figure 1, connecting to the electric circuit of Figure 4a. For the supporting structures, stands and beams were aluminum alloys welded to the base and transducers to allow vibration as much as possible, as shown in Figure 4b. Moreover, a tank body with a thickness of 2.0 mm did not weld to the supporting structures but could be adjusted to stainless steel, acrylic, and glass, materials generally used by UCT manufacturers.

Experiment
This paragraph explains the equipment and its setup. For convenience, we did not rebuild a circuit board entirely but modified it from a commercial product of the manufacturer using our design presented in Figure 4a. The modified circuit board was connected to the supporting structures and two transducers functioning as the ultrasonic generator applying ultrasonic waves of 28 kHz, which was tested for safety and precision before being employed in the actual experiment. For the tank body, three rectangular tanks made of glass, acrylic, and stainless steel with a dimension of 14 × 20 × 18 cm 3 (W × L × H) 5.04 L, were built. Figure 5 shows the equipment, which are outcomes of Figure 4: (a) a modified circuit board and (b) a built UCT with an adjustable tank body to glass, acrylic, and stainless steel. In the experiment, water as a cleaning fluid of 14 × 20 × 14 cm 3 , 3.92 L was filled into the tank body placed on the supporting structures shown in Figure 5b. When the power of the ultrasonic generator was turned on, cavitation occurred in the water, according to the physics principles with the theoretical background described in Section 2. The experiment changed the tank body material to determine the effect on the ultrasonic cleaning. The experiment was divided into foil corrosion and power concentration experiments.

Foil Corrosion Experiment
This experiment was conducted to test the suitable location of the ultrasonic cleaning in the fluid. The same technique used previously in [12][13][14] was applied again. The foil sheet was stretched with wires; it was then submerged to the water inside the tank body with a sonication time of 3 min. Cavitation caused corrosion of the foil sheet, which can be employed for acoustic pressure and tank design analyses. Figure 6a reveals an example of foil corrosion after the experiment using the tank body made of acrylic, which can be changed to glass and stainless steel later.

Power Concentration Experiment
This experiment was conducted to determine the cavitation intensity. Again, the same technique employed previously in [9,13] was used. The NGL measurer (Model UPC 30000) measured the power concentration in three sampling points: P 1 , P 2 , and P 3 . P 1 and P 2 were in a middle plane of the UCT located near and far away from transducers, respectively. P 3 is near the tank body at an aside. Each sampling point was repeated five times to calculate an average value and a standard deviation. The NGL measurer with a resolution of ±0.1% W/L is a popular tool among manufacturers to determine cavitation intensity of the UCTs of frequency below 68 kHz. At a frequency higher than 68 kHz, the NGL measurer is not appropriate because the measurement sensitivity does not match the cavitation intensity. The higher the measured power concentration, the greater the cavitation intensity. In academics, the cavitation intensity can be significantly presented in terms of acoustic pressure [9,13]. Figure 6b shows an example of the power concentration experiment using the tank body made of acrylic and sampling points.

Simplified and Mesh Models
The designed CAD model in Figure 4b was simplified by ignoring components that did not affect the simulation results, leading to a simplified CAD model shown in Figure 7a. It consists of water, the tank body, PZT4, stainless steel, and aluminum alloy domains. The simplified model in Figure 7a was employed as a prototype to create a mesh model, as shown in Figure 7b. Because of symmetry, a half-mesh model was therefore created to calculate faster. According to the mesh criteria, at least six mesh elements are required in one wavelength [9,15,21]. Since the transducer frequency was 28 kHz and the sound velocity in water was 1507 m/s, the wavelength was 53.8 mm. Therefore, the largest mesh size that produced accurate results enough was approximately 9 mm. In this research, the mesh size was 3.5 mm, smaller than the mesh criteria. As a result, this mesh size can be applied to simulate up to 68 kHz with high accuracy. In this research, the maximum frequency to be investigated was 40 kHz, so the created mesh size was smaller than the criteria, ensuring accuracy. In the mesh model, elements over 96% were hexahedrons, while 4% were wedges and pyramids. A picture on the right shows an enlarged of the elements in fine detail. After using the mesh independent analysis to find suitable mesh models, models which provided accurate results proper to computational resource usage, we found the mesh models with 88,531-150,523 elements, 368,529-47,3736 nodes, 0.58-0.70 skewness was the suitable mesh model. The number of elements and nodes varies slightly depending on the changing variables. Figure 8 shows boundary conditions defined in the HRA labeled as A-H. A small picture reveals boundary conditions for transducers in detail. Since the UCT is approximately 3.92 L, which is a small volume compared with the UCTs in [13][14][15] with a capacity of 10-20 L, the critical settings, therefore different from [13][14][15], are an acoustic radiation boundary (B) assigned to the topwater surface and an acoustic impedance boundary (D) assigned to the tank body. The assignments of B and D were made to mimic the reflection, transmission, and radiation of the ultrasonic waves and achieve high accurate simulation results consistent within the aims of this research. The UCTs in [13][14][15] had higher water volume, so without setting of both B and D did not affect the results and conclusion. In addition, the symmetry condition enabled rapid calculation. As a result, the simulation results across the symmetrical plane were identical. Mechanical properties of materials such as acrylic, glass, stainless steel, aluminum alloy, PZT4, and water were defined according to Table 1. When the boundary conditions and material properties were defined entirely, the calculation displayed the acoustic pressure in shade color for easy analysis.

Results and Discussion
This section divides the HRA simulation results of acoustic pressure into two parts: validation and tank design effect. First, the acoustic pressure distribution patterns were compared with the foil corrosion and the power concentration experiments to validate the results and methodology of this research. Second, the acoustic pressures at the sampling points obtained by varying variables such as tank body material, thickness, volume, and frequency were reported and analyzed to investigate the tank design effect on the ultrasonic cleaning process. Figure 9 shows a comparison between the acoustic pressure distribution patterns on the left and foil corrosion on the right for (a) the xz and (b) the yz planes using acrylic as a tank body. Both planes were located in the middle position of the UCT, as shown in small pictures on the right. The red and blue colors present the maximum and minimum acoustic pressures, respectively. The comparison results revealed the consistency of both results in marked areas and indicated that the center position of UCT was suitable for the ultrasonic cleaning process, as expected. In addition, the acoustic pressure distribution depended on transducer positioning, consistent with the reports in [8,9,[12][13][14][15]. Moreover, the comparison results of stainless steel and glass tank bodies remained the same but differed only in acoustic pressure and corrosion intensities.  Figure 10 presents a comparison between the acoustic pressure and power concentration for tank bodies made of acrylic, glass, and stainless steel at the sampling points (a) P 1 , (b) P 2 , and (c) P 3 , in the same plane as Figure 9a. The blue is for the acoustic pressure, and the red is for the power concentration. Again, consistency of both results was found for all sampling points. Both results depended on the measured points. The error bar was calculated from the standard deviation of measurements. The difference in acoustic pressure at each sampling point, the effect of the tank body on acoustic pressure, and an explanation in physics are discussed at the end of this section. Results and discussion reported above of Figures 9 and 10, including the methodology, confirmed the reliability of this research. Therefore, the HRA simulation by varying variables also yielded accurate results.

Tank Design Effect
This section reports the effect of the tank design on the acoustic pressure by varying variables such as material, power, thickness, volume, and frequency. Similar boundary conditions and settings were assigned to all UCT models but different only mentioned variables. All results were collected from the plane in Figure 9a

Materials
The acoustic pressures were plotted against the tank body materials at sampling points P 1 , P 2 , and P 3 , as shown in Figure 11. As expected, the acoustic pressure depended on the measured position with the same distribution pattern as displayed in Figure 9. By comparing the measured points, the acoustic pressures were the highest at P 1 , moderate at P 2 , and the lowest at P 3 , for all sampling points. In-depth, acoustic pressure inequality can be explained as P 1 closest to transducers, P 2 inferior, and P 3 farthest. P 1 and P 2 provided greater acoustic pressure than P 3 , so they had more cavitation; therefore, they were proper for the ultrasonic cleaning process. In addition, P 3 was not a proper position for the cleaning since it had less cavitation.  Figure 12 shows the maximum acoustic pressures: (a) comparison results for varying power of transducers at 60-140 W and (b) a sample of the distribution pattern of acrylic with the power of 100 W. As expected, the maximum acoustic power increased with increasing the power, consistent with [12,15]. However, changing the power and tank body material did not affect the acoustic pressure distribution pattern. Instead, it remained the same as in Figure 12b but only differed in the maximum acoustic pressure. Likewise, acrylic provided the highest acoustic pressure, followed by glass and stainless steel accordingly.

Thickness
Similarly, Figure 13 presents the maximum acoustic pressure: (a) comparison results for varying materials of the tank body and thicknesses at 1.0-2.5 mm and (b) a sample of the distribution pattern of stainless steel with a thickness of 2.5 mm. The original tank body thickness is 2.0 mm. The maximum acoustic pressure increased with increasing the thickness. As expected, the increasing thickness did not affect the acoustic pressure distribution pattern since it depends on the frequency and transducer placement [7][8][9][12][13][14][15]. Again, the highest, moderate, and lowest acoustic pressures were found in acrylic, glass, and stainless steel tank bodies, respectively. It is worth noting that the maximum acoustic pressure was almost unchanged for stainless steel. According to the density of materials in Table 1, it might be because stainless steel has a density of 7550 kg/m 3 , 7 times higher than acrylic, and 3.5 times greater than glass, approximately. Therefore, it can be said that the UCT made with the tank body's lower density provides higher acoustic pressure than the one with higher density. More in-depth supporting reasons are discussed at the end of this section. However, the tank body with lower density may have a lower durability not suitable for such applications requiring frequent relocating or long-continued usage.

Volume
The water in the original UCT has a dimension of 14 × 20 × 14 cm 3 (W × L × H) with a volume of 3.920 L; therefore, the HRA simulation was performed by extending the length from 20 mm to 22 mm and 24 mm. The extending changed the water volume to 4.312 L and 4.704 L, respectively. Figure 14 presents the maximum acoustic pressure: (a) comparison results for varying volume and tank body materials and (b) a sample of the distribution pattern of glass with a water volume of 3.920 L. The maximum acoustic pressure decreased with increasing volume. This result is consistent with the manufacturer's report that for the same UCT model and specification, the tank body with a larger volume had lower cleaning efficiency than the smaller one. The larger the volume, the lower the efficiency of cleaning. Therefore, if the customers require a larger UCT, the manufacturers must increase the number of transducers to achieve a proper acoustic pressure level with the highest cleaning efficiency and meet customers' requirements. Moreover, comparing Figure 14b to Figures 13b and 12b, the distribution patterns still remained the same, so extending the length of the tank body in this range did not affect the distribution pattern. Additionally, acrylic, glass, and stainless steel provided superb, medium, and the least acoustic pressures, respectively.

Frequency
The transducers were changed to 15 kHz and 40 kHz using exact specifications provided by the manufacturer to investigate the frequency effect on the acoustic pressure. Figure 15 compares the maximum acoustic pressure at 15 kHz, 28 kHz, and 40 kHz by varying tank body materials. As expected, the maximum acoustic pressure increased with increasing frequency, consistent with [12,15]. The higher the frequency, the greater the maximum acoustic pressure. Moreover, varying tank body materials also affected the acoustic pressure. A similar trend remained the same, from the highest to the lowest: acrylic, glass, and stainless steel, accordingly.   The top, middle, and bottom rows display the HRA simulation results of acrylic, glass, and stainless steel, respectively. Again, increasing the frequency affected the distribution pattern. Higher acoustic pressure with smaller areas and more distribution throughout the UCT were observed. When increasing the frequency, the HRA simulations revealed higher acoustic pressure with the smaller area and more dispersed throughout the UCT, caused by increasing the number of wave superpositions. The higher the frequency, the greater the number of wave superpositions. The distribution patterns in this figure were consistent with the reported in [12,15]. Therefore, increasing the frequency increases the cavitation intensity and reduces the size of bubbles, which are dispersed more throughout the UCT. The section after this part discusses the maximum acoustic pressure inequality provided by acrylic, glass, and stainless steel, including the cause of the problem. By comparing all materials in Sections 4.2.1-4.2.5, acrylic was shown to yield the highest acoustic pressure, followed by glass and stainless steel, respectively. Therefore, the acoustic pressure depends on the tank body material. In factories, each tank body material is implemented in a different situation. For example, the tank body made of stainless steel is suitable for general cleaning objects. However, it is not suitable for acidic or base cleaning fluids. In contrast, glass is specially employed for acidic or base cleaning fluid. Furthermore, acrylic is often used in processes requiring timely inspection of cleanliness in a factory or an academic laboratory.
We determine the reason in physics why acrylic produced the highest acoustic pressure, and glass and stainless steel decreased, respectively. The reflection (α R ) and transmission (α T ) coefficients were calculated using Equations (6) and (7) and acoustic impedances in Table 1. The calculation results are reported in Table 2. Clearly, the acoustic pressure directly relates to the transmission coefficient. The tank body material with high α T yielded more significant acoustic pressure than the one with less. The higher the α T , the greater the acoustic pressure. Therefore, if the cleaning process requires intense cavitation, the UCT should have a tank body material with high α T . The HRA simulation results also indicated that the tank designs, such as power, body material, thickness, volume, and frequency, affected the acoustic pressure. Therefore, the complaint problem would be concluded that came from the tank designs that were not suitable for washing the products. As a result, the cavitation effect generated by the problematic UCTs might be too high to cause damage or too low to make unclean on the products. Indeed, choosing suitable tank designs and the HRA simulation not only solves the complaint problem but also ensures high cleaning efficiency. All research findings have already been applied as essential information for developing a new generation of high-efficiency UCTs by the manufacturers.

Conclusions and Further Study
This article aimed to investigate the cause of ultrasonic washing damage problems reported by the customers whose factories applied the UCTs in the cleaning process and find ways to achieve high cleaning efficiency. First, an electric circuit and a CAD model of the UCT were designed and then actually built in a laboratory for experiment and a computer for the simulation. For the experiment, the built UCT had a tank body of UCT that could be adjusted to acrylic, glass, and stainless steel materials. For the HRA simulation, the tank designs such as body material, power, thickness, volume, and frequency could be modified. Next, the foil corrosion and power concentration were applied to the experiments, and the HRA was employed in the simulation. Since acoustic pressure relates to cavitation and ultrasonic cleaning, it was analyzed to reach the aims. The HRA results revealed that the acoustic pressure was comparable to the foil corrosion and power concentration experiments, confirming the validity of simulation results and the research methodology. Then, by varying the tank designs, the HRA results revealed that the acoustic pressure increased with increasing the power, thickness, and frequency, but decreased with increasing volume. The tank body made of acrylic provided the highest acoustic pressure, moderate for glass, and the lowest for stainless steel, accordingly. In addition, changing frequency affected the acoustic pressure distribution pattern while varying the material, power, thickness, and volume did not alter the one, as expected. The acoustic pressure difference relates to the transmission coefficient of the material. The higher the transmission coefficient, the greater the acoustic pressure. Finally, all analysis results indicated that the complaint problem was caused by the tank design that was unsuitable for usage. Therefore, the HRA simulation, together with a proper tank design, will lead to a high-efficiency ultrasonic cleaning process in accordance with the customers' requirements. The findings from this research are how to apply the HRA in developing the UCT and tank design effect on cavitation for interested persons to apply in their work.
Water is not a usage cleaning fluid in an actual ultrasonic cleaning process. Instead, most factories prefer to use deionized (DI) water. The DI water is mixed with surfactants and cleaning chemicals to activate the cavitation effect quickly. The mixing ratio depends on the type of dirt, sonication time, and cleaning objects to be washed, which is a secret of each factory and chemical vendor. With a proper mixing ratio, cavitation occurs rapidly and violently compared with regular water to avoid recontamination [14] and reduce sonication time. However, the proper mixing ratio comes from the test repeatedly, which consumes workload and budget. Unfortunately, the HRA simulation cannot accurately predict the acoustic pressure after the mixing. From the authors' experience, adding the chemicals to aid in the ultrasonic cleaning process relates to fluid dynamics coupled with chemical reactions, which the HRA does not cover. Therefore, in further study, finding computer simulations and mathematical models that eliminate weakness and help support the HRA will surely enhance the cleaning efficiency and promote the domestic manufacturers to always be the leader in this field.