Numerical simulation of the dimensional transformation of atomization in a supersonic aerodynamic atomization dust-removing nozzle based on transonic speed compressible ﬂow

To simulate the transonic atomization jet process in Laval nozzles, to test the law of droplet atomization and distribution, to ﬁnd a method of supersonic atomization for dust-removing nozzles, and to improve nozzle efﬁciency, the ﬁnite element method has been used in this study based on the COMSOL computational ﬂuid dynamics module. The study results showed that the process cannot be realized alone under the two-dimensional axisymmetric, three-dimensional and three-dimensional symmetric models, but it can be calculated with the transformation dimension method, which uses the parameter equations generated from the two-dimensional axisymmetric ﬂow ﬁeld data of the three-dimensional model. The visualization of this complex process, which is difﬁcult to measure and analyze experimentally, was realized in this study. The physical process, macro phenomena and particle distribution of supersonic atomization are analyzed in combination with this simulation. The rationality of the simulation was veriﬁed by experiments. A new method for the study of the atomization process and the exploration of its mechanism in a compressible transonic speed ﬂow ﬁeld based on the Laval nozzle has been provided, and a numerical platform for the study of supersonic atomization dust removal has been established.


Introduction
Based on newly modeled data from the World Health Organization, 92% of the world population lives in places where air quality levels exceed the WHO ambient air quality guidelines for the annual mean of particulate matter with a diameter of less than 2.5 micrometers (PM2.5) (World Health Organization 2016). This kind of dust with a very small particle size is called respiratory dust, and it comes from industries such as mining and machining and civilian sources (Betts 2002;Chen et al. 2012;Shi et al. 2009). In addition to endangering human health through inhalation, it can also attach to many toxic and harmful substances, with which it can undergo some complex chemical reactions (Bern et al. 2009;Maertens et al. 2008;Jang and Kamens 1999). Aerodynamic spray is one of the main methods in the field of dust removal, and its efficiency mainly depends on the atomization fineness and particle speed. In the process of collecting fog droplets, the droplet size needs to be close to that of the target dust particles, which results in an extremely high requirement for the atomization efficiency of the nozzle (Dust control theory; Huang et al. 2015;Park et al. 2017;Yang et al. 2018;Jin et al. 2018).
The Laval nozzle is a kind of device used for fluid control and is widely used in air acceleration. Under certain pressures, the velocity of air in the expanding section breaks the speed of sound (Cengel and Cimbala 2013;Boyanapalli et al. 2013;Kuan and Witt 2013). Because of this characteristic, the structure is often used in jet atomization nozzles (Yang et al. 2003). In this process, the jet fluid instantly breaks down into tiny droplets by overcoming the surface tension in the supersonic flow field due to the relative velocity difference between the flow and the liquid jet (Ö zkar and Finke 2017;Ghosh and Hunt 2000;Yang et al. 2012). The ultrasonic dusting atomizer is a type of Laval nozzle that is used to produce ultrasonic speeds and break droplets, but this method cannot be considered a good application of the abovementioned method. We must question whether there is a way to avoid the waste resulting from this supersonic process and to find a supersonic atomization dusting nozzle that uses a different process than ultrasonic atomization (Soo 1973). Therefore, to improve the atomization performance of dusting nozzles and explore a kind of ultrafine atomizing nozzle with an ultrahigh efficiency, it is necessary to study the atomization process in Laval nozzles, which is a process in which supersonic gas breaks up droplets; however, no scholars have done this, and the process is difficult to study through experiments. The following analogous studies have been performed.
Several scholars used the gas-liquid two-phase flow method and simulated the atomization process of kerosene jet liquid injection in supersonic transverse laminar flow, but the atomization particle size distribution and gas flow acceleration process were not reflected in the results (Jing andXu 2006, Jing 2010;Shunhua and Jialing 2008;Li et al. 2015). Wu (2016) studied the jet atomization process of liquid hydrocarbon fuel in a scramjet combustor by combining experimental research and theoretical analysis. However, no achievements have been made in crushing atomization in the complex transonic process of compressible fluids. Deshpande and Narappanawar (2016) studied the structure and internal flow field characteristics of the Laval nozzle and its related effects; good or bad results were determined only on the basis of air flow parameters such as the velocity of flow, but droplet atomization and the particle size distribution were not considered. Sinha et al. (2015) and Alhussan and Teterev (2017) studied the spray jet process by numerical methods and considered that the nature of instability depends on the density of the jet fluid and the ambient fluid and the velocity of the jet. Inside the Laval nozzle, the air density and temperature are different at different locations, and atomization is a complicated process (Introduction to aerodynamics 2007). However, none of these studies combines the Laval air compression acceleration process with the gas-liquid atomization process.
Supersonic atomization in Laval nozzles is unique, and the density, temperature and speed of the air flow should not be constant. Thus, although it has already performed very well, atomization simulation with transverse laminar flow cannot approximate the real situation inside the atomizer in existing simulations.
In addition, most of the existing Laval studies simulate incompressible airflow, and the dimensions of the simulation model are all two-dimensional. To approximate the real jet atomization process of the transonic compressible flow field, it is necessary to calculate in the 3D model space. This has not yet been achieved by any researchers, mainly because of the following difficulties.
(1) The flow of supersonic gas in the Laval nozzle is a kind of compressible flow, and the density and temperature of the gas are different in different positions.
(2) The flow process in the Laval nozzle is 2D axisymmetric. In CFD calculations, the 2D axisymmetric model should be applied in the case of COMSOL, which is essentially a 2D calculation process. In data processing, the dataset of results is rotated with the symmetric axis. This is the correct way to obtain 3D data. However, the problem is that the results of the droplet atomization model cannot be treated this way because after the rotation, the original data points will develop into rings. Because of this, we cannot observe the distribution of particles in 3D space nor can we obtain the particle size distribution and concentration distribution of a certain interface. Detailed instructions are mentioned later in Chapter 4. (3) In the 3D model of COMSOL, it is impossible to calculate the compressible airflow acceleration process in the Laval nozzle because it is unable to be defined axisymmetrically, and the grid cannot be divided axisymmetrically. In addition, it is not accurate to apply symmetric boundaries to the calculation. Detailed instructions are mentioned later.
Therefore, the calculation of the flow field and atomization cannot be realized simultaneously in both 2D axisymmetric and 3D models. In this way, the characterization of jet atomization in transonic compressible flows in Laval nozzles becomes a seemingly impossible task. However, an approximate method has been found in our study. The rationality of the method is verified by experiments. A new idea for related research and software development will be developed. Although this combination has achieved good results, due to professional limitations, the two-way coupling between droplets and compressible flow fields cannot be realized by this simulation, which will need to be solved by software and theoretical researchers. The chapters of the paper are summarized below.
After the introduction, the article contains four chapters. In chapter two, the basic theoretical equations and derivations used in the simulation process are introduced. In the chapter three, the model details and grid division of the numerical simulation, including the grid division and optimization of the 2D axisymmetric model and the grid and grid division of the 3D symmetric model, are introduced. In chapter four, the results and a discussion of the simulation are highlighted. By combining the simulation results of several methods to discuss their infeasibility and reasons, the paper introduces the simulation research method, which realizes the numerical simulation process of droplet atomization in the transonic flow of a compressible gas in 3D space by transforming the dimension. Physical quantitative analysis and atomization experiments are carried out to verify the reliability of the simulation. Finally, the dust reduction test is used to prove the guiding significance of the method for practical application by analyzing the dust removal efficiency. The last section presents the conclusions of the article.

Theory of transonic compressible flow
In this work, the simulation is realized by using the finite element method as implemented in CFD of COMSOL Multiphysics software. COMSOL Multiphysics is a special multifield-coupled numerical simulation software, which plays an important role in theoretical research and practical simulation in various fields and can realize the coupling between any multiphysical fields. It is able to establish onedimensional to three-dimensional geometric models and interface with other software. In the mesh division, the direct mesh generator can be used to divide the triangle and tetrahedron units, and some detailed parts can be refined. Local encryption can achieve better results than ANSYS software. In the definition of the variables and equations, embedding can be realized quickly, and the interface is more friendly, for the drawing and other subsequent processing, COMSOL Multiphysics also has a powerful postprocessing function, which can carry out point calculations and sections, draw the parameter evolution graph of the body, and realize the integral operation of the trajectory and variables.
According to the characteristics of the Laval nozzle, some assumptions have been made through the Spalart-Allmaras module.
The conservation of mass is a fundamental concept of physics. Within some problem domains, the amount of mass remains constant; mass is neither created nor destroyed. The mass of any object is simply the volume that the object occupies times the density of the object. For a fluid (a liquid or a gas), the density, volume, and shape of the object can all change within the domain with time, and mass can move through the domain (Deshpande and Nitin 2015). The conservation of mass (continuity) tells us that the mass flow rate m through a tube is a constant and equal to the product of the density q, velocity ũ, and flow area A: Considering the mass flow rate equation, we could increase the mass flow rate indefinitely by simply increasing the velocity for a given area and a fixed density. In real fluids, however, because of compressible effects, the density does not remain fixed with increasing velocity. We must account for the change in density to determine the mass flow rate at higher velocities. When the specific heat ratio is c, the gas constant is R s , and the temperature is T. If we start with using the isentropic effect relations, the equation of state and the mass flow rate equation given above, we can derive a compressible form of the mass flow rate equation based on the definition of the Mach number M and the speed of sound c: here M is the Mach number. Substituting Eq.
(2) into Eq. (1), when q represents the mass density of the fluid, its mass flow rate can be defined by: when P is the pressure, the equation of state is: Substituting Eq. (4) into Eq. (3), the mass flow rate is defined by: when the total pressure is P tot and the total temperature is T tot , from the Isentropic flow equations: Substituting Eq. (6) into Eq. (5), the mass flow rate can be indicated as: Numerical simulation of the dimensional transformation of atomization in a supersonic… Another isentropic relation gives: Substituting Eq. (8) into Eq. (7): This equation as shown relates the mass flow rate to the flow area A, total pressure P tot and temperature T tot of the flow, the Mach number M, the ratio of specific heats of the gas c (value of 1.4 in the study), and the gas constant R s (Mass Flow Rate Equation). Appropriate initial variables are also defined by the equations above in the simulations in our work.
The fluid is assumed to be a Newtonian fluid, which satisfies the Navier-Stokes (N-S) and the continuity equations.
The compressible form of the Navier-Stokes (N-S): where F is the external force per unit volume, and l is the coefficient of viscosity. For compressible fluids, generally, there are two coefficients l and l 0 , similar to Hook's elastomer. Both measure the magnitude of the internal friction work caused by the expansion or contraction of a fluid.
Except for extreme cases such as with high temperatures and high frequency sound waves, the motion of gases can be approximately assumed to have a l' = 0 because l 0 r Á ũ is much smaller than P. This fact has been indicated in the theory of molecular motion but was originally a hypothesis put forward by Stokes.
l varies with the temperature, which conforms to Sutherland's law.
In this work, l ref is determined by T l,ref (substituted with the static temperature inlet T stat ), and S l usually valued at 111 K.
The density and the velocity of the fluid further satisfy the continuity equation: In addition, the fluid should be a steady, axisymmetric, isotropic, and adiabatic flow. In other words, ou ot ¼ 0, so Eq. (10) can be simplified as follows: To close the equation, some equations need to be introduced as follows: Equation (14) indicates energy conservation. The thermal conductivity of the fluid satisfies Fourier's Law: The Prandtl number P r is a dimensionless parameter in fluid mechanics to characterize the relative importance of the momentum exchange and the heat exchange in the fluid flow. Except for the critical point, the P r of the gas is almost independent of the temperature and pressure, and it takes the value of 0.72 in this study.
The boundary condition for the wall in the simulation can be restricted as follows: Equation (17) assumes that there is no slip and no heat exchange on the outside wall of the nozzle (Patel et al. 2016).
Based on the above assumptions, the physical process of the supersonic flow field in a Laval nozzle can be simulated by solving the above equations when the total pressure P tot , the total temperature T tot and the Mach number M in at the inlet are given. In this case, the inlet has a total pressure of 6 atm, a total temperature of 293 K, and a Mach number of 0.01. All parameters are listed in Table 1.

Theory of atomization simulation (K-H)
First, to represent the continuity of the fluid, the jet was regarded as a continuous droplet column of particles released by several calculation time steps in the transient calculation, as first proposed by Reitz. The particles' attributes are as follows: the dynamic viscosity is l p , the surface tension is r p , and the density is q p . Because the cross-sectional area of the liquid column is the cross- Inlet T tot (K) 293 Inlet M in 0.01 Outlet P out (atm) 1 R s (J/(kgÁK)) 287 sectional area of the particle, when the particle has a mass of m p , the initial radius R is controlled by the m p when the particles' material attribute is fixed: when the rate of the particle velocity change over time is dv dt , according to Newton's law, the force on the particle in the present is: The particle receives a drag force of air in the flow field. When the velocity of the air flow is u, the dynamic viscosity is l, the density is q, and the diameter of the particle is d p , the velocity is v, as described by Stokes, and the Shear stress is s p ¼ q p d 2 p 18l ; then, the force of the particle can be expressed as: when the velocity difference between the air flow and the droplet at their contact surface is Du = |u -v|, and because q/q p ( 1 by the Kelvin-Helmholtz instability model (K-H), the following is known.
When the gas phase Weber number is We g ¼ The radius of the sub-droplets produced when the mother liquid droplets with a radius of r p are broken is: when r KH is the K-H broken particle radius, the shear force is s KH ¼ 3:788B KH K KH X KH : B KH is an empirical constant \ 10, which represents the atomization time of particles; the smaller the value is, the earlier the breakup time is, and the faster the atomization rate is (Patel et al. 2015). Then, the particle atomization rate is: Gas-phase properties were derived from compressible flow simulations, and other parameters in the droplet atomization model were from experimental measurement and set in Table 2.
By determining the rate of mass flow m r and the initial radius of particle d 0 , the initial velocity of particle V 0 and the transient step size were determined, which ensured the continuity of the jet and tracking requirements for particles.

Introduction of the new method: dimensional transformation
The velocity U of any point in a 3D space (x, y, z) can be decomposed into V z and V r in a 2D axisymmetric coordinate system. In addition, in a 3D space it can be decomposed into U x , U y and U z . Their relationship is shown in Fig. 1.

Fig. 1 Vector U of dimension transformation
Numerical simulation of the dimensional transformation of atomization in a supersonic… : The actual variable setting in the software can be replaced by an interpolation function. The velocity field is set as follows: 3 Model details and simulation mesh 3.1 Model details Figure 2 shows a sketch of the simulation setup composed of the Laval Nozzle and the calculation area of the external atmosphere. In this setup, the length of the Laval nozzle and that of the external calculation area were chosen as 30 and 1000 mm. The distances from the throat to the nozzle inlet and outlet are 10 and 20 mm, respectively. The nozzle has an inlet diameter and an outlet diameter of 10.2 and 10 mm (inlet [ outlet), respectively, because of energy conservation and the interior resistance. The atmosphere has a diameter of 1000 mm. The throat has a diameter of 2 mm. All the above parameters were set for the sake of closely approximating reality.

Two-dimensional mesh
The geometry of the nozzle was reproduced with COM-SOL Multiphysics software. As described in detail by Moses and Stein (1978), only half of the computational domain was required due to symmetry, as shown in Fig. 3. The mesh was constructed to follow the flow, facilitating convergence of the solution (Singhal and Mallika 2013;Venkatesh and Jayapal 2015;Patel et al. 2016).
In the software, the sidewall curve of the nozzle was drawn in the two-dimensional working plane by the parametric equation. A closed nozzle geometry was generated by drawing perpendicular lines from the two ends of the curves to the center axis, the axis between the two ends and the sidewall curve. Then, a rectangle with a side length of 500 9 1000 mm was drawn to indicate the atmospheric area outside the nozzle.
The mesh was divided into triangular meshes because these can capture the variation of the turbulent flow field. The boundary layer and angular refinement used stripshaped quadrilateral meshes, which facilitate spreading along the boundary. The following diagram is a scale graph of the grid; the chromatographic legend on the right represents the mesh cell scale relationship in each subgraph.
An extreme refinement fluid dynamics mesh was used in the nozzle area, a refinement fluid dynamics mesh was used near the nozzle outlet, and a fluid dynamics coarsening mesh was used in the other areas. An additional refinement method based on boundary control was used to refine the throat. All the grid refinement details in the diagram, such as the angle refinement, maximum cell, minimum cell and curvature number, are marked in Fig. 3. As shown in Fig. 4, the flow field has clear velocity boundaries and was no longer changed by mesh refinement, which ensured that the speed of the microchange could be captured.

Three-dimensional mesh
As shown in Fig. 5, the 3D mesh for the 1/4 symmetric model was divided by a tetrahedral mesh. The entire mesh was divided into four parts to match the calculations. Similarly, the boundary-controlled refinement method was applied, and the specific parameters are shown.
However, in later research, it was proven that the above method was not appropriate. As shown in Fig. 6, a further optimization was made, with the mesh cells presenting an approximate axisymmetric distribution.
4 Results and discussion 4.1 Simulation under the 2D axisymmetric module space

Simulation of the 2D flow field
As shown in Fig. 7, the velocity, pressure, density and dynamic viscosity of flow in the Laval nozzle are different at each position. Accordingly, based on verification using experiments and formulas, the error of the simulation was found to be 3.2%.

2D simulation of droplet atomization
Under the same model space, the previously calculated flow field results were called. The results of the atomization simulation are shown in Fig. 8.  In Fig. 8, the fog shows a state of a ''breathing wave''. The particle size near the axis of symmetry is approximately 2 microns, and they accelerate at supersonic speeds in the nozzle. The particle size away from the axis of symmetry is approximately 10 microns, and the particle breaking position away from the axis of symmetry is the subsonic band. Therefore, the spatial distribution of smallsized particles is inside, while that of large-sized particles is outside.
Thus, particles close to the axis of symmetry have higher velocities. In addition, their particle radius is five times that of the particle near the axis of symmetry, and the resistance during the movement is much greater than that of the particle with a smaller radius. This allows particles close to the axis of symmetry to catch up in a very short period with the outer particles that were sprayed earlier, so the wave-like spray distribution is observed macroscopically.

Problems of the simulation in the 2D axisymmetric module space
When we want to further discuss the 3D distribution law of particles, we find that the biaxial symmetry model cannot be realized. Just as everyone knows, for a 2D axisymmetric model, if you want to get a 3D result you need to rotate along the axis, but in Fig. 9, by postprocessing the rotated data set, the results are terrible. Rings were formed by particle rotation, which not only increased the difficulty of analyzing and determining the statistics of the particle distribution but also does not conform to the actual situation. In addition, in this model space, the particles released at the sidewall of nozzle correspond to the water injection mode of the annular discrete slot. Thus, the 3D visualization of particles cannot be realized in a 2D axisymmetric model. Therefore, we want the model space to be threedimensional in the first place.

Simulation of the flow field
In Fig. 10, the 1/4 of the flow field was calculated by applying symmetric boundaries in three dimensions, and the other 3/4 was achieved by postprocessing; its mesh is shown in Fig. 5. Although this method produced a correct calculation result on the symmetric boundary, the results between the boundary were wrong. This is because the mesh of the boundary at symmetry cannot approximate axisymmetric conditions. Thus, the next step is to further approximate the mesh to axisymmetric conditions.     In Fig. 11 with the mesh optimized along the axis in the 3D model (mesh in Fig. 6), although the results are better than those obtained by the previous method, some tinpot edges and corners still exist in the flow field, which is not satisfactory.
Because the flow field of the Laval nozzle is totally axisymmetric, the mesh setting cannot be perfect in any case, and it is not appropriate to calculate the 2D axisymmetric problem in a 3D space.
In a word, it is not accurate to calculate the compressible flow field in the Laval nozzle under the 3D model only, regardless of whether or not it is approximated as symmetric or axisymmetric.

Dimension transform numerical simulation of atomization in supersonic atomizing dustremoving nozzle based on transonic speed compressible flow
Generate the flow distribution data with a 2D axisymmetric spatial basement into 3D basement data. Import the above data with a 3D basement, including the velocity, density, dynamic viscosity and temperature of the flow field, into the 3D model by means of an interpolating function. In the model settings, some details were set according to the rules of spatial transformation by Eq. (24).
In the process of calculating droplet atomization, the above functions were called, and the particle states were calculated by a transient. The simulation results of particle atomization in the nozzle are shown in Fig. 12.
Calculate the atomization angle and compare it with the experiment result in the next chapter. The position of the released particles is z = 7.5 mm from the outlet, and the diameter of the particle concentration section is d & 2.5 mm, so tan a = d/z = 0.356, and a = 20°.
As shown in Fig. 12, particles with a large size adhere to the wall of the nozzle, and the micron-sized particles are instantly produced and blown out of the nozzle wall at a very high speed, according to the initial velocity direction when the particle is released. The results are basically consistent with the two-dimensional results.
To enhance the effect of the simulation analysis, increase the number of particles released. As shown in Fig. 13, a large number of micron-sized particles present delamination and spray into the outer space of the Laval nozzle.
The phenomenon of a breath wave is obvious. In threedimensional space, the particle coverage cannot reach full Fig. 12 The simulation results of atomization a in the nozzle, b out of the nozzle Fig. 13 Atomization results out of the nozzle 360°coverage, which is due to the position of particle release on the sidewall surface that has a certain angle and initial velocity. This phenomenon is called an irregular jet, which is consistent with the results of the subsequent atomization experiment. As shown in Fig. 14, the particles are released from the sidewall at 0 s and move slowly. When the flow of particles is in contact with the supersonic band inside the nozzle at 4 9 10 -6 s, the jet droplets are instantly broken into microns and accelerated to subsonic speed (approximately 230-300 m/s) dispersions outside the nozzle.
The closer the axis of symmetry is to the nozzle, the greater the velocity difference between the particles and the flow field is, the more intense the acceleration and breakage are, and the faster the velocity is; however, the latter is ejected out of the nozzle. The axisymmetric stratified distribution of the airflow field causes this phenomenon.
According to the steady state calculation report, the error of the separation group is estimated to be 0.00087,0.00061; the residual error of the separation group is estimated to be 50, 4.9 9 10 3 ; the maximum number of iterations of separation 1(se1) is 1000; the damping factor of the separation step (ss1) is 0.5; the damping factor of the turbulent variables is 0.35; the tolerance factor is 1.0; the solver is PARDISO; the relative tolerance is 0.0012; and the transient calculation absolute tolerance is 1.0 9 10 -6 .

Rationality analysis by experiment
Aerodynamic atomization will produce an ultrafine fog curtain and a high-speed wind curtain, fog droplets in the fog curtain will capture the free dust, and the turbulent air layer formed by the wind curtain will restrain the dust diffusion Liu et al. 2019;Cai et al. 2018). Therefore, the accuracy of the simulation is verified by the atomization and dust removal experiments and the guiding significance of the research on the design of dustremoving nozzles.
The experimental parameters were set according to the boundary conditions of the simulation, and the experimental results were compared. The nozzle used in the experiment was designed and manufactured according to the simulation parameters.
In Fig. 15, based on the comparison, the simulated fog droplets were mainly concentrated along the axial direction and had a very high velocity, which was similar to the actual spray situation. The flame shape of the fog and the characteristic of the breath-waves in the process of droplet migration were in good agreement with the experimental results. The main distribution of the particles at the exit position was between 2 and 100 microns, and the spray angle was basically consistent with the experimental results. In both the simulation and the experiment, the water droplets with larger particle sizes were ''blown'' out along the sidewall of the nozzle outlet.  In Fig. 16, based on a particle size test by laser diffraction, it is basically proven that the D50 is approximately 8 microns at 250 mm from the nozzle outlet. This is because the data in the nozzle cannot be measured accurately under the conditions of the existing test instrument, and the accuracy of the simulation can only be verified by the results of the longer distance data outside the nozzle.
The results agree well, which proves that the simulation method is true and reliable.
According to the existing dust removal theory, it is considered that the particle size of droplets used in dust harvesting should be close to that of dust particles (coal dust particles), so for conventional respiratory dust (particle size less than 10 microns) in the coal mining industry, such atomization methods should have a high dust removal efficiency. Therefore, a dust control experiment was carried out to investigate the dust control efficiency of the nozzle on the respiratory coal dust.
In Fig. 17, the experimental process included simulating the air pressure, setting the water flow dust removal experimental conditions, taking dust samples through the dust sampler, weighing the sample before and after drying the mass, calculating the dust concentration, and carrying out the dispersion test. Finally, according to the measured dust concentration results and dispersion test results, the Fig. 16 Atomization experiment to verify a the particle size distribution of the simulation results at 250 mm from the outlet of the nozzle, b the particle size distribution of the particle size test by laser diffraction at 250 mm from the outlet of the nozzle, and c the particle size test by laser diffraction transonic atomization dust removal was compared with the conventional respiratory dust removal efficiency.
After the dust removal experiment, this kind of atomization method was shown to achieve a respiratory coal dust removal efficiency of 78.9% in 1 min, and the capture effect of respiratory dust was much higher than that of the conventional method (52.03%), which proves that the simulation method has a basic and guiding significance for the development and research of new technology.

Conclusions
(1) It is proved that the 2D axisymmetric, 3D symmetry and 3D approximate axisymmetric models cannot accurately simulate the process of droplet fragmentation and atomization in the transonic flow field in three-dimensional space.
(2) Through transforming the dimension, based on the COMSOL software, the droplet crushing and atomization process in the transonic flow field in 3D space is simulated in the 3D model, and the atomization experiment is carried out to verify the analysis. The particle distributions of the mesoscopic and macroscopic phenomena are all consistent with a high degree of coincidence. (3) The physical process, macro phenomena and particle distribution of supersonic atomization are analyzed in combination with this simulation. (4) Through the atomization and dust removal experiments, the significance of the simulation was proven.
Funding Supported by the National Natural Science Foundation of China (NO: 51704146, 51274116, 51704145).

Compliance with ethical standards
Conflict of interest We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work, there is no professional or other personal interest of any nature or kind in any product, service and/or company that could be construed as influencing the position presented in, or the review of, the manuscript entitled author.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.