Preprint
Article

Influence of Acoustic Streams on the Efficiency of Ultrasonic Particle Agglomeration

Altmetrics

Downloads

74

Views

26

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

13 November 2023

Posted:

13 November 2023

You are already at the latest version

Alerts
Abstract
The article is devoted to the study of ultrasonic agglomeration of PM2.5 in uniform and non-uniform ultrasonic fields. The possibility of increasing the efficiency of ultrasonic agglomeration by initiating acoustic streams in a resonant non-uniform ultrasonic field is shown. A non-uniform ultrasonic field with zones of high and low sound pressure levels formed using a bending-oscillating disk transmitter made it possible to initiate acoustic vortex-type streaming that promote the movement of particles into the nodal areas of a standing wave and between them. The efficiency of ultrasonic treatment was determined indirectly by the degree of inertial trapping of agglomerates. It has been established that the impact of an non-uniform ultrasonic field in comparison with the impact of a uniform ultrasonic field makes it possible to increase the efficiency of inertial capture of agglomerates PM2.5 from 89% to 95%, PM1.5 from 85% to 92%, PM0.5 from 76% to 85%. The results obtained made it possible to propose a design for an agglomeration chamber, which allows for a multiple increase in the productivity of ultrasonic processing of a gas-dispersed flow.
Keywords: 
Subject: Physical Sciences  -   Applied Physics

1. Introduction

The main efforts of environmental services are aimed at ensuring high-quality purification of ambient air polluted by industrial production and man-made natural phenomena. This is due to the fact that the presence and uncontrolled distribution of aerosols of various substances in the air has a negative impact on humans, flora and fauna [1,2]. It is generally accepted that the most dangerous particles are those measuring 2.5 microns or less, since the respiratory system introduces such particles into the bloodstream and accumulates toxic material in human organs [3,4]. The extremely dangerous nature of fine aerosols has been confirmed by numerous studies, in particular those conducted by L. Calderon-Garciduenas, A. Solt, C. Enriquez-Roldan and others. As a result of these studies, a connection has been established between long-term exposure to fine aerosols and sudden deaths of healthy children and young people [5]. The proof of such a correlation explains the exceptional importance of finding ways to protect people and developing systems for purifying air from fine aerosols.
All existing cleaning systems (separators) based on the inertial collection method practically do not capture particles smaller than 2.5 microns [6,7,8,9], and systems capable of capturing such particles (filters and electric precipitators) are characterized by low dust capacity or operate with the formation of very harmful particles. environment and human health nitrogen oxides and ozone.
The most effective way to solve the problem is the use of high-intensity ultrasonic vibrations, the effect of which on aerosols causes processes of convergence and agglomeration of small particles into large ones. This process is called ultrasonic agglomeration [10,11]. Fine particles combined due to ultrasonic action can subsequently be easily captured by conventional gas cleaning devices [12,13].
However, as follows from results of studies by various authors (J.A. Gallego-Juarez, R.R. Andres, C. Sheng), ultrasonic agglomeration loses its advantages when aggregating particles smaller than 2.5 μm for their subsequent removal, while providing a significant increase in efficiency (up to 99%) in the removal of particles larger than 2.5 μm, as per [14,15,16].
The reason for this lies in the very physical principles underlying the sonic agglomeration process. It is known that the impact of ultrasonic vibrations on gas media generates two types of nonlinear effects [17] affecting dispersed particles in gases:
1) Effects arising due to the phenomena of momentum transfer of the gas phase between particles;
2. Effects arising due to oscillatory motions of solid particles of different sizes relative to each other.
The first group of effects contributes to the hydrodynamic mechanism (first described by W. Koenig in 1891) of particle agglomeration, when equivalent forces of hydrodynamic interaction (arising as a result of mutual distortion of flow fields around particles) bring particles closer to each other and promote their agglomeration. In this case, the probability of collision of particles is proportional not only to the interaction force of particles, but to their concentration, primarily (and, in fact, to the distance between particles).
The second group of effects is realized in orthokinetic mechanism (first described by O. Brandt, H. Freund, E. Hiedemann), which consists in the fact that the particle doesn’t oscillate with the gas phase synchronously, while it slides relative to the gas phase and oscillates with a smaller amplitude of displacement than the gas itself. The sweep effect is related to the inertia of the particle [18]. The ratio of the particle oscillation amplitude to the gas oscillation amplitude is called the drag coefficient. The drag coefficient depends on the particle size, and particle collision occurs due to the fact that particles of different sizes oscillate with different amplitude (with non-zero amplitude relative to each other) and randomly collide and aggregate.
For fine particles smaller than 2.5 μm, the drag coefficient modulus is close to 1, i.e., all particles oscillate almost in-phase with approximately equal amplitude, regardless of the sound pressure value. Accordingly, the orthokinetic interaction is weakly realized particle [18]. In addition, for particles to aggregate, the distance between them should not exceed the displacement amplitude of gas molecules in the sound field.
For these reasons, the low probability of particle collision, especially at low concentrations (large distances between particles) leads to the fact that, even at the maximum permissible sound pressure level (before the destruction of the formed agglomerates), increasing the efficiency of agglomeration of particles smaller than 2.5 μm due to ultrasonic vibrations becomes virtually impossible [19].
For example, in [19] it was shown that studies of ultrasonic agglomeration were carried out under conditions of detection of high particles in concentrations from 0.2 to 5 g/m3. Such conditions are typical during the production of any product and its capture in the form of highly dispersed particles. At the same time, due to ultrasonic exposure, the cleaning efficiency increased from 30-50% to only 40-85%.
However, research results indicate that at relatively low concentrations (less than 0.2 g/m3, which is typical for industrial emissions, for example at thermal power plants), the efficiency of agglomeration will be very low, the degree of agglomeration will not exceed 2-3.
Various agglomeration systems utilizing standing wave aerosol impacts are proposed to address this problem.
It is known that Gallego-Juárez et al. developed a semi-industrial pilot-scale system combining an electrostatic precipitator and an acoustic agglomeration chamber to remove submicron particles (average particle size 0.8 µm) from a flue gas flow [20]. In this system, four ultrasonic transducers were arranged at the bottom of the agglomeration chamber, and the standing wave sound field was generated perpendicular to the gas flow to ensure agglomeration. The disadvantage of this design is low efficiency at low concentration. This is due to the fact that during generation of the standing wave the particles almost do not interact with each other even within the node areas.
Miura et al. developed several types of airborne ultrasonic transducers that generate powerful sound waves in the air [21,22]. Among them is an ultrasonic radiator in which a rigid wall was attached around the circumference of a circular transversely vibrating plate. This ultrasonic radiator shapes an intense standing wave sound field into a cylindrical sealed space perpendicular to the plane of the circular transversely vibrating plate. Therefore, the fine particles can be retained in enclosed space in which the powerful sound field is formed. Using this sound wave agglomeration equipment, the effect of fine particle concentration on agglomeration was studied, and the agglomeration effect was higher at higher particle concentrations. However, ultrasonic wave agglomeration for PM2.5 achieved an agglomeration efficiency of no more than 63%. This is due to large distances between particles, at which the ultrasonic interaction processes do not take place.
This determines the urgent need to find new ways to further increase the efficiency of ultrasonic agglomeration of particles smaller than 2.5 μm.
Understanding the actual mechanisms of ultrasonic agglomeration suggests that, if it is impossible to further increase the sound pressure level (the value of direct impact on the particles), the probability of collision of small particles can be increased by giving them additional movements that contribute to their collision and aggregation. Such an additional impact on small particles in the ultrasonic field can be provided only through utilization of secondary effects arising under certain conditions of generation of high-intensity ultrasound impact.
One of such secondary intensifying effects is the so-called "ultrasonic wind". This is a second order effect due to momentum transfer, represented by the Reynolds voltage. Theoretical studies of acoustic flows are well developed [23,24] and carefully modeled [25,26].
However, when ultrasonic (US) impact is implemented in an open space or a long gas duct, the radiation pressure generated enables the particles to be moved only in one direction without significant increase in efficiency of interaction between particles [27].
For this reason, in order to increase the efficiency of small particles interaction in the maximum ultrasonic field in terms of sound pressure level, it is essential to generate the motion of particles in different directions, so that their mutual movement and interaction is realized.
In addition, it is required to increase the concentration of particles in certain areas by moving them to significantly increase the efficiency of particle agglomeration, especially at their low concentration [8,11,20,28]. This can be provided by generation of vortex flows, because only vortex flows will increase the efficiency of interaction (aggregation) of particles due to the following factors:
1. Increasing the speed of mutual movement of particles of different sizes.
2. Increasing duration of interaction of particles due to their confinement in vortex areas;
3. Generating zones with increased concentration of particles.
Therefore, the research results presented below are focused on identification of conditions for generation of vortex flows in ultrasonic fields with the maximum sound pressure level of ultrasonic impact. Comparative studies of the process of agglomeration of particles with a size of 2.5 μm in the presence of vortex flows and without them allows to establish the actual values of increasing the efficiency of ultrasonic agglomeration due to vortex formation of gas-dispersion flow by acoustic streaming.

2. Materials and Methods

2.1. Test bench

To study the agglomeration process of 2.5 μm particles and to conduct comparative studies in the presence of vortex flows and without them, a bench implementing ultrasonic agglomeration and subsequent capture of aggregated particles was used. The test bench consists of agglomeration chamber (Figure 1) and sequentially installed cyclone based on the design of VNIIOGAZ CN-15 with a capacity of 10 m3/hour. The cyclone outlet tube is connected to the air handling unit.
The developed and manufactured agglomeration chamber of the bench (RF Patent No. 2759506) is suitable for implementation of various conditions and modes of exposure to high-intensity ultrasonic vibrations due to application of different types of ultrasonic transducers, and owing to possible resonance amplification of vibrations by changing the distance between radiator and reflector [28].
Two ultrasonic oscillating systems (No. 1, Figure 1) with disk radiators (No. 4, Figure 1) of two different types creating uneven or uniform distribution of vibrations for ultrasonic impact with and without flow generation were used as a source of ultrasonic impact in the agglomeration chamber. A reflector (No. 5, Figure 1) with centrally placed inlet tube (No. 6, Figure 1) is installed across the disk radiator.
The process of particle aggregation in the agglomeration chamber takes the following course. The gas flow enters the tube (No. 6, Figure 1) installed in the center of the reflector (No. 5, Figure 1). After that, the gas-dispersion flow spreads uniformly from the central area to the periphery. In the process of propagation, the flow is affected by ultrasonic vibrations generated between the surfaces of radiator (No. 4, Figure 1) and the reflector (No. 5, Figure 1). The gas flow (No. 8, Figure 1) and aggregated particles (No. 9, Figure 1) are discharged through outlet tubes (No. 10, Figure 1).
The design of the agglomeration chamber allows adjusting height of the reflector position to achieve fine tuning of the resonance mode (standing wave mode in a gas medium). When the distance between the radiator and the reflector ensures generation of standing wave, the effectiveness of the treatment is increased by means of increasing the sound pressure level.
The test bench was equipped with measuring and auxiliary equipment. A noise meter Ecofizika-110A was used to measure the main parameter (i.e. sound pressure level). The TIPAS-1 aerosol dispersion characteristics measuring instrument, which enables monitoring the aerosol parameters, was used to measure the dispersion properties of the aerosol. TIPAS-1 consists of two units:
1. The first unit implements the optical method of small-angle scattering (radiator, reflector) and is designed to measure micron aerosol particle sizes and concentrations in the range of 1-100 µm and 0.009-8.0 g/m3;
2. The second unit implements the optical spectral transparency method (radiator and reflector) and is designed to measure submicron aerosol particle sizes and concentrations in the range of 0.1 – 2 µm and 0.009-8.0 g/m3.
Sampling to determine the mass concentration of particles was carried out at the inlet of the agglomeration chamber and at the outlet of the cyclone. The dispersed particles were supplied by pneumatic atomizer of ejection type installed in front of the inlet of the agglomeration chamber.
A laser plane plotter is installed in the chamber to track particle trajectories. The laser plane passes through the axis of symmetry of the radiator. The camera is focused on the laser plane and enables video recording of the particles travelling between the radiator and the reflector. The directions of gas flows are determined indirectly by trajectories of light-scattering (tracer) particles.

2.2. Ultrasonic transducers

Two ultrasonic radiators consisting of a piezoelectric Langevin transducer and a disk-shaped bending and oscillating radiator were used for the research. At the same time, due to excitation of oscillations on the non-basic mode, oscillation maxima and minima are generated on the frontal surface of the radiator, interchanging at a distance from each other corresponding to half the wavelength of bending oscillations in the radiator material [29,30]. To generate an inhomogeneous ultrasonic field and its impact on the gas-dispersion flow, a bending-oscillating disk radiator with a flat frontal surface (hereinafter referred to as the flat radiator) was used. In turn, to generate a homogeneous ultrasonic field, a radiator with phase-equalizing collars with a height equal to half of the wavelength in the gas (hereinafter referred to as a step-function radiator) was used.
Figure 2 shows diagrams explaining the principle of operation of the two types of disk radiators used for comparative studies. The radiators were made with equal-sized radiation surfaces and with the test frequencies as close to each other as possible.
As follow from Figure 2(a), the neighboring annular domains of the disk generate equal-amplitude oscillations in antiphase. Such a radiator enables generation of inhomogeneous ultrasonic field, in which vortex flows between the regions oscillating in counter-phase can be formed. In this version of the radiator, the characteristic dimensions of the vortices generated will be determined by the wavelength of ultrasonic vibrations in the gas medium, and their number – by the number of neighboring areas of the ultrasonic field oscillating in counterphase and the distance between the radiator and the reflector [28,29].
To ensure phase synchronism (generation of homogeneous ultrasonic field) of vibrations emitted into the gas medium from each annular domain, the front surface of this type of radiators has a stepwise profile (Figure 2(b)). The arrangement and size of the projections on the front surface determine the key parameters of the ultrasonic field generated by the radiator. Therefore, to generate a uniform ultrasonic field, step-type transitions are made in the node circle zones with a height equal to half of the wavelength of ultrasonic vibrations in the gas medium. Thus, generation of a uniform ultrasonic field is ensured both when the reflector is installed at a resonant distance from the oscillating surface of the radiator and when it’s not available.
Figure 3 shows the vibration distributions and pictures of the fabricated radiators of both types.
Each piezoelectric transducer connected to the radiator is equipped with an electronic generator for its power supply, which maintains constant amplitude of oscillations. Technical specifications of ultrasonic transducers are presented in Table 1.

3. Results

3.1. Calculation of sound pressure distribution in coagulation chamber

The ultrasonic field parameters were estimated in order to establish conditions of the most effective impact and comparative analysis of sound pressure level distributions generated by means of different radiators, provided the distance from the radiator to the reflector (for both flat surface and step-function type radiators) at which the standing wave occurs is ensured.
Modeling of vibration distribution was carried out by means of finite element analysis using ‘Harmonic Acoustics’ harmonic acoustic analysis module. During modeling, the boundary conditions were set on the basis of obtained experimental data regarding parameters of developed and manufactured ultrasonic radiators (frequency of vibrations and vibration amplitudes distribution of the radiating surface). The surfaces of the internal walls of the chamber are specified as reflecting boundaries with a certain absorption level. By comparing the results of modeling with the values of sound pressure level measured experimentally, the level of absorption of ultrasonic vibrations by the walls of agglomeration chamber was determined, which amounted to 11%. This value is accepted in further calculations.

3.2. Analysis of sound pressure distribution

The sound pressure level distributions for both radiators are shown in Figure 4. It follows from presented results that there is a sound pressure level between the radiator and the reflector (160-172 dB) sufficient for particle agglomeration. The distance (wavelength λG = 15 mm) at which the maximum sound pressure level is reached is also determined.
It established that the sound pressure level between the reflector and the radiator when using the step function radiator is 168-171 dB, a uniform ultrasonic field is generated therewith, and when a flat radiator is used it is 168-172 dB, the generated field is inhomogeneous, the areas of the minimum sound pressure level, which are located near the zeros of vibrations of the disk radiator are clearly pronounced.
As follows from the presented results, for a step radiator, the sound pressure amplitude gradient has a non-zero component only along the acoustic axis. This means that the equivalent force that can initiate the “sonic wind” is also directed along the acoustic axis. However, since the reflector located opposite the radiator serves as an obstacle to the “sound wind,” no acoustic flows are actually formed in the air gap with a step radiator. In this case, the equivalent force leads to the formation of a static pressure gradient in the gas, counteracting the force (a single vortex flow may occur due to the inhomogeneity of the field along the perimeter of the radiator).
It follows that the step radiator is not capable of creating conditions that further increase the efficiency of particle agglomeration. Choosing a certain distance between the radiator and the reflector can lead to an increase in the sound pressure level due to the phenomenon of resonance. However, even with a resonant layer thickness (15 mm), the created sound field will not ensure effective agglomeration of particles due to their short residence time in the gas gap.
For a flat radiator, the sound pressure amplitude gradient (and since there is absorption of vibrations, there will also be a phase gradient) has non-zero components along all coordinate axes. In this case, depending on the position of the observed point in space, the amplitude gradient can be directed at different angles (from 0 to 2π) to the surface of the radiator. This means that the equivalent volume force capable of initiating the “sonic wind” has similar directions, creating a vortex motion. In this case, the surface of the radiator and reflectors do not interfere with the movement of the gas flow and, therefore, the occurrence of acoustic flows with a flat radiator is possible.

3.3. Particle path analysis

At the first stage, the paths of particles in the agglomeration chamber were analyzed using the developed test bench. Particle paths were tracked on the basis of frame-by-frame tracking of tracer particles filmed with a camera. The output power of the generator was set so that the sound pressure level for both types of disks used was at least 165 dB. The resulting frame images are shown in Figure 5.
Based on the analysis of the obtained particle path data, diagrams of paths travelled by particles in the agglomeration chamber using both types of radiators were plotted (Figure 6 and Figure 7).
When using a flat radiator (at distances of 7.5 and 15 mm), a standing wave is generated with two and three nodal planes of particle velocity, respectively, into which solid particles move, agglomerate and where they are contained. Herewith, intense vortex flows generated due to adjacent regions of the ultrasonic field oscillating in counter-phase originate. As a result, the particles are captured by the vortex motion both within one node areas and between them through the low sound pressure level zones located in the vicinity of (opposite to) the zeros of oscillations of the disk radiator. In this case, in addition to the agglomeration of particles in the node areas of vibrational velocity their additional reciprocal vortex movement and interaction with each other take place, which further results in increase in the efficiency of agglomeration
Changing the distance to a greater or lesser extent, relative to the resonance distance, results in origination of a traveling wave and disappearance of node areas. At the same time, vortex flows are practically not generated, reducing the efficiency of agglomeration. It is also established that at a distance equal to λ, the intensity of vortex flows is higher than at a distance equal to λG/2. At the distances equal to ¼ λG and ¾ λG, conditions at which the sound pressure level reaches its minimum value become effective.
For comparison, Figure 7 shows a flow diagram for a radiator generating a uniform ultrasonic field.
As in the previous case, when a step-function radiator is used, a standing wave with two and three nodal planes of vibrational velocity into which the solid particles travel, agglomerate and are retained, respectively, is generated when the distance between the radiator and the reflector is 7.5 mm, as well as when the distance is 15 mm.
Further particle agglomeration does not occur due to containment of already formed particles in the standing wave nodes. At the same time, there is no mutual movement of particles within the nodal zone without any additional forced gas flow. Continuous supply of aerosol results in additional increase in the size of particles retained by the ultrasonic field in the node areas. At the same time the vortex flows are not observed.

3.4. Estimation of particle capture efficiency based on gas-dispersion flow and ultrasonic impact parameters

To identify the possibility of increasing the capture efficiency using the developed stand, experimental studies were carried out. Research has made it possible to determine the influence of various factors (flow rate of gas-dispersed flow, sound pressure level, mass concentration) on the efficiency of agglomeration of dispersed particles. The efficiency of agglomeration was assessed by comparing the sizes of the initial particles entering the inlet of the agglomeration chamber with the sizes of particles entering the inlet of the cyclone (sampling was carried out in the pipe that connects the cyclone and the agglomeration chamber).
To assess the contribution of ultrasonic agglomeration to the efficiency of the process of collecting dispersed particles, during research, the efficiency of collecting particles entering the inlet of the agglomeration chamber was determined.
The capture efficiency was calculated using the following equation: η
η = 1 N o u t l e t N i n l e t · 100 %
where Noutlet – is a particle concentration at the cyclone outlet, g/m3, Ninlet is a particle concentration at the agglomeration chamber inlet, g/m3.
In order to determine the contribution of effects arising from ultrasonic exposure to the process of particle agglomeration, a control experiment was conducted without ultrasonic exposure. In all experiments, the hydrodynamic parameters of the gas-dispersed flow and the characteristics of the particles were the same. This made it possible to exclude the influence of factors related to the properties of particles (adhesion, extrusion, etc.) and flow parameters (hydrodynamic turbolization, etc.).
To exclude the influence of other effects arising in the sound field on the agglomeration efficiency (for example, the appearance of higher and lower harmonics, distortion of the waveform at high sound pressure levels), the parameters of ultrasonic exposure for both types of radiators within the same experiment were the same.

3.5. Determining the optimum gas flow rate

At the next stage, the gas-dispersion flow rate, at which the maximum efficiency of particle capture is ensured, both with and without acoustic impact was determined. The studies were conducted under the following conditions: sound pressure level was 165 dB; initial mass concentration Ninlet = 2∙10-2 g/m3. The gas flow rate varied from 2.5 to 10 m3/hour. Microtalc Jetfine T1 CA with particle diameter d50 =1 μm, d95 =3 μm was used as a dispersed material in the research.
The bulk density of microtalc is 0.9 g/cm³, the true density is 2.8 g/cm³. The studies were carried out at an air temperature of 22 degrees and a relative humidity of 35%. Figure 8 shows the dependence of particle capture efficiency on gas flow rate and the dependence of the particle size at the outlet of the agglomeration chamber on the flow rate.
The dependence of the efficiency of particle collection on the velocity of the gas-dispersed flow is extreme. The presence of the extremum is due to the mutually opposite dependence of the size of the formed particle agglomerates and the efficiency of particle collection in the cyclone on the air flow speed. As the flow speed increases, the efficiency of the cyclone increases due to an increase in the centrifugal force acting on the particles in the cyclone (Figure 8A). However, at the same time, the size of the aggregates formed in the agglomeration chamber decreases due to a decrease in the residence time of the particles in the ultrasonic field (Figure 8B). This leads to a limitation in the growth of the efficiency of the agglomeration chamber – cyclone system.
It follow form the obtained curves that the best gas flow rate, at which the maximum capture efficiency is provided, is 6.2 m3/hour both when using flat and step disk radiator.
However, without ultrasonic impact, the optimum flow rate was 8-10 m3/hour, which corresponds to the optimum flow rate for a cyclone.
With a further increase in gas flow (above 6.2 m3/hour), in the presence of ultrasonic influence, the capture efficiency begins to decrease. This is due to the fact that reducing the size of the formed agglomerates (due to a further reduction in the residence time of particles in the agglomeration chamber) becomes a determining factor affecting efficiency.
When the flow rate exceeds 8-10 m3/hour, the decrease in efficiency increases. This can be explained by the additional influence of the cyclone itself, the effectiveness of which is reduced due to increased turbulization of the flow, which leads to mixing of the flow in the cyclone, reducing stratification
This indicates that the efficiency of particle agglomeration has a greater influence on the capture efficiency than establishing the optimum gas flow rate.
Thus, the efficiency of particle agglomeration under the impact of ultrasonic vibrations depends on the gas flow rate (on duration of ultrasonic exposure). It is obvious that the longer the period of ultrasonic exposure, the higher the efficiency of agglomeration.
In addition, analysis of the relationships presented in Figure 8B showed that there is a nonlinearity in agglomeration efficiency versus flow rate. In the case of using a stepped radiator, nonlinearity is associated with an increase in the probability of collision of large particles held at the nodes of a standing wave with newly arrived particles.
An increase in flow speed leads to the separation of large particles from the nodal areas and their entrainment into the cyclone.
When using a flat radiator, nonlinearity is associated both with the factor described above and with a decrease in the intensification of agglomeration due to acoustic flows. In this case, the movement of particles between nodal regions by acoustic flows is reduced. As a result, with increasing gas flow, the agglomeration efficiency curves for both types of radiators approach each other.
Further studies were conducted with an optimum flow rate of 6.2 m3/hour.

3.6. Effect of sound pressure level on particle capture efficiency

To determine the influence of the sound pressure level on the capture efficiency and particle size at the outlet of the agglomeration chamber, dedicated experimental studies were carried out at a concentration of Ninlet = 2∙10-2 g/m3. Figure 9 shows the obtained dependencies.
It follows form analysis of the obtained curves that at the initial sound pressure level L = 130 dB the capture efficiency practically does not depend on the used radiator type, there is almost no efficiency gain due to vortex flows, which indicates their insufficient development and influence.
As the sound pressure level increases, an increase in particle capture efficiency is observed. At the same time, the increase in efficiency is more significant when exposed to inhomogeneous field. This is due to the intensification of vortex flows with increasing sound pressure level.
The growth of the capture efficiency continues until the sound pressure level reaches the values of 165 dB, herewith for the step radiator the efficiency is less than 90%. In turn, the vortex acoustic flows generated by the flat radiator at the sound pressure level of 165 dB allow increasing the capture efficiency by 6%.
It is known that when exposed to gas flows or in a non-uniform ultrasonic field, large agglomerates begin to disintegrate due to the influence of acoustic flows [31]. It can be seen from Figure 9B that particle size does not increase significantly above 160 dB SPL. Therefore, we can conclude that no further increase in particle size is observed due to the destruction of agglomerates.

3.7. Effect of mass concentration on particle capture efficiency

To determine the influence of particle mass concentration on the capture efficiency and particle size at the outlet of the agglomeration chamber, experimental studies were carried out at a sound pressure level of 165 dB. The curves below (Figure 10) were obtained as a result of these studies.
As follow from the presented curves, at concentrations up to Ninlet = 2∙10-2 g/m3 the particle capture efficiency increases significantly. At the same time, the capture efficiency is higher when using a flat radiator, due to increased efficiency of particle agglomeration. At concentrations of Ninlet = 3∙10-2 g/m3 and above, the increase in efficiency slows down, and at concentrations above Ninlet = 4∙10-2 g/m3 the capture efficiency becomes almost the same with both types of radiators. This happens because at the specified concentrations the distances between particles appear sufficient for realization of orthokinetic and hydrodynamic mechanisms of particle agglomeration, and the influence of vortex flows has little effect on increasing the particle collision probability.
Analysis of the dependencies presented in Figure10B shows that the impact of a flat radiator on particles leads to a significant increase in particles already at a concentration of Ninlet = 0.2 10-2 g/m3, while particles with a size of 2.5 μm increase to a size of 7.4 μm. At the same time, due to the impact of a stepped radiator on the particles, the agglomerates at the exit from the agglomeration chamber increase to only 4.2 microns. Such a significant difference in the sizes of agglomerates is due to an additional intensifying factor: acoustic flows arising in a non-uniform ultrasonic field generated by a flat radiator. This allows particles to interact not only within the nodal regions, as happens when exposed to a uniform field generated by a stepped radiator, but also to interact with particles moving from one nodal region to another.
As the concentration increases, the difference between the sizes of the formed agglomerates decreases and at a concentration of Ninlet = 5∙10-2 g/m3 it is less than 10%.

3.8. Determination of fractional efficiency

To confirm the efficiency of capturing particles of different sizes, experimental studies were carried out using a mixture of microtalc of different grades (Jetfine T1 CA d50 =1 μm, Luzenac St30 d50 =6 μm, Luzenac H50 d50 =9 μm) which had polydisperse composition. Comparative experiments were carried out at initial particle concentration Ninlet = 2∙10-2 g/m3, at three different modes of agglomeration chamber operation: exposure with a flat radiator; exposure with a step-function radiator; without ultrasonic exposure. As a result, the fractional efficiency data shown in Figure 11A were obtained. Figure 11 B shows the dependence of the particle size at the outlet of the agglomeration chamber on the particle size at the inlet of the coagulation chamber.
The analysis of the obtained curves shows that:
1. Particle capture efficiency without ultrasonic impact (curve 3, Figure 11A) is significantly lower than with exposure to it.
2. When using a step-function radiator (curve 2, Figure 11A), when a standing wave is generated, the particles hardly interact with each other within the node areas, and at very low concentration the agglomeration efficiency is reduced to zero.
3. The use of a flat radiator (curve 1, Figure 11A) enables generation of vortex flows between adjacent regions oscillating in counter-phase, which significantly increases the interaction of particles located in nodal zones.
4. The dependences presented in Figure11B allow us to conclude that when exposed to a flat radiator, the maximum increase in agglomeration efficiency is achieved on the smallest particles
Thus, the results of experimental studies have confirmed that increasing the efficiency of ultrasonic agglomeration can be reached owing to generation of acoustic streaming in the form of vortex flows. In this case, the efficiency of capturing particles of 2.5 µm in size increases by 6% - from 89% to 95%, particles of 1.5 µm in size by 7% - from 85% to 92%, and particles of 0.5 µm in size by 9% - from 76% to 85 %.

3.9. Increasing the productivity of gas-dispersed flow processing

The considered design of the agglomeration chamber can significantly increase the efficiency of coagulation, which subsequently increases the efficiency of capturing dispersed particles.
However, the processing capacity of the gas-dispersed flow of such a chamber (with a disk diameter of 320 mm) is no more than 10 m3/hour. To increase productivity, the authors proposed the design of an agglomeration chamber with several parallel oscillating plates, which are installed at a resonant distance from each other and implement the identified optimal conditions and modes of ultrasonic exposure. A sketch of the proposed sintering chamber is shown in Figure 12.
The module is based on an ultrasonic oscillatory system (No 5, Figure 12) with a plate radiator (No 4, Figure 12), which are installed in the housing (No 1, Figure 12). A reflector is installed opposite the emitter (No 6, Figure 12). The sketch and vibration shape of the plate radiator are shown in Figure 13.
The radiator of the oscillatory system is made in the form of parallel plates (No 1, Figure 13). The plates undergo bending vibrations (Figure 13 B). When installing an oscillating system and a reflector in the agglomeration chamber, the distance between the plates of the radiator and reflector is equal to the wavelength of ultrasonic vibrations in the air.
With a resonant distance between the plates of the emitter and reflector in the air gap, conditions are created for the emergence of intense acoustic flows and vortices that promote the coagulation of dispersed particles and the deposition of their aggregates.
In general, the module works as follows. A gas flow containing highly dispersed particles enters the inlet pipe (No 2, Figure 12). Further, along the direction of gas movement, flow distributors are installed (No 7, Figure 12), which evenly distribute the gas-dispersed flow throughout the chamber volume.
Next, highly dispersed particles enter the area of ultrasonic influence - the gaps between the plates of the ultrasonic emitter (No 4, Figure 12) and reflector (No 6, Figure 12). In this area, high-intensity ultrasonic action occurs on particles and liquid droplets in a thin gas gap.
This design of the agglomeration chamber makes it possible to increase the productivity of processing a gas-dispersed flow in proportion to the area of the radiator.

4. Conclusions

The studies of agglomeration process of particles of 2.5 µm and less in ultrasonic field generated by oscillating disk radiators in resonant gaps were carried out in terms of this paper. The conducted studies and presented results of agglomeration of particles of 2.5 µm and less in a uniform field of ultrasonic vibrations confirmed that there is no technical possibility to implement the process with an efficiency of more than 89% for particles smaller than 2.5 µm.
To solve this problem, a new concept of increasing the efficiency of agglomeration through generation of alternating zones in the resonance gaps, in which particles oscillate in counter-phase through the use of bending oscillating ultrasonic disk radiators was proposed. Fabrication and application of such radiators ensured generation of vortex type acoustic streaming capable of moving particles within the node areas of the standing wave and between them, which increased the probability of their collision.
The conducted studies and obtained experimental results allowed us to identify conditions and modes of increasing the efficiency of ultrasonic agglomeration owing to generation of acoustic streaming capable of moving particles within and between the node areas of the standing wave.
It is shown that the new approach to ultrasonic agglomeration and its practical implementation provide an increase in the rate of inertial capture for particles of 2.5 μm size by 6% - from 89% to 95%, for particles of 1.5 μm size by 7% - from 85% to 92%, and for particles of 0.5 μm size by 9% - from 76% to 85%.
The obtained research results made it possible to propose a design for an agglomeration chamber, which allows for a multiple increase in the productivity of ultrasonic processing of a gas-dispersed flow.

Author Contributions

Vladimir N. Khmelev: Supervision, Conceptualization, Reviewing and Editing; Andrey V. Shalunov: Methodology, Validation, Investigation, Reviewing and Editing; Viktor A. Nesterov: Formal analysis, Methodology, Software, Investigation, Validation, Reviewing and Editing; Sergey A. Terentiev: Investigation, Reviewing and Editing.All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Russian Science Foundation (project No. 19-19-00121).

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, X.; Zhang, S.; Xing, J.; Wang, Y.; Chen, W.; Ding, D.; Hao, J. Progress of air pollution control in China and its challenges and opportunities in the ecological civilization era. Engineering 2020, 6, 1423–1431. [Google Scholar] [CrossRef]
  2. Miller, L.; Xu, X. Ambient PM2. 5 human health effects—Findings in China and research directions. Atmosphere 2018, 9, 424. [Google Scholar] [CrossRef]
  3. Pope, C.A.; Ezzati, M.; Cannon, J.B.; Allen, R.T.; Jerrett, M.; Burnett, R.T. Mortality risk and PM 2.5 air pollution in the USA: An analysis of a national prospective cohort. Air Quality. Atmosphere & Health 2017, 11, 245–252. [Google Scholar] [CrossRef]
  4. Halonen, J.; Lanki, T.; Yli-Tuomi, T.; Tiittanen, P.; Kulmala, V.; Pekkanen, J. Particulate air pollution acute cardio respiratory hospital admissions and mortality among the elderly. Amrican Journal of Epidemilogy 2009, 20, 143–153. [Google Scholar] [CrossRef]
  5. Calderón-Garcidueñas, L.; Solt, A.C.; Henríquez-Roldán, C.; Torres-Jardón, R.; Nuse, B.; Herritt, L.; Reed, W. Long-term Air Pollution Exposure Is Associated with Neuroinflammation, an Altered Innate Immune Response, Disruption of the Blood-Brain Barrier, Ultrafine Particulate Deposition, and Accumulation of Amyloid β-42 and α-Synuclein in Children and Young Adults. Toxicologic Pathology 2008, 36, 289–310. [Google Scholar] [CrossRef] [PubMed]
  6. Fu, P.B.; Wang, F.; Yang, X.J.; Ma, L.; Cui, X.; Wang, H.L. Inlet particle-sorting cyclone for the enhancement of PM 2.5 separation. Am J Environ Sci. 2017, 51, 1587–1594. [Google Scholar] [CrossRef]
  7. Hoekstra, A.J. Gas flow field and collection efficiency of cyclone separators. Ph.D. Thesis, Technical University Delft. 2000.
  8. Leith, D.; Licht, W. The Collection Efficiency of Cyclone Type Particle Collectors-A New Theoretical Approach. AIChE Symp. 1972, 68, 196–206. [Google Scholar]
  9. Qian, F.; Wu, Y. Effects of the inlet section angle on the separation performance of a cyclone. Chem. Eng. Res. Des. 2009, 87, 1567–1572. [Google Scholar] [CrossRef]
  10. Riera, E.; González-Gomez, I.; Rodriguez, G.; Gallego-Juarez, J.A. Ultrasonic agglomeration and preconditioning of aerosol particles for environmental and other applications. Power Ultrasonics 2015, 1023–1058. [Google Scholar] [CrossRef]
  11. Ng, B.F.; Xiong, J.W.; Wan, M.P. Application of acoustic agglomeration to enhance air filtration efficiency in airconditioning and mechanical ventilation (ACMV) systems. PLoS ONE 2017, 12, 1–26. [Google Scholar] [CrossRef]
  12. Chen, H.; Liu, W.; Li, J.; Xun, X.; Shen, X. Experimental Study on Acoustic Agglomeration of Fine Particles from Coal Combustion. In Proceedings. In Proceedings of the 2010 International Conference on Digital Manufacturing and Automation, ICDMA. 1, Changcha, China, 18-20 December 2010; pp. 702–705. [Google Scholar] [CrossRef]
  13. Fan, F.; Zhang, S.; Wang, W.; Yan, J.; Su, M. Numerical investigation of PM2.5 size enlargement by heterogeneous condensation for particulate abatement. Process Safety and Environmental Protection 2019, 125, 197–206. [Google Scholar] [CrossRef]
  14. Gallego-Juarez, J.A.; Rodriguez, G.; Acosta, V.; Riera, E. Power ultrasonic transducer with extensive radiator for industrial processing. Ultrasonic Sonochemistry 2010, 17, 954–964. [Google Scholar] [CrossRef]
  15. Andres, R.R.; Acosta, V.M.; Lucas, M.; Riera, E. Modal analysis and nonlinear characterization of an airborne power ultrasonic transducer with rectangular plate radiator. Ultrasonic 2018, 82, 345–356. [Google Scholar] [CrossRef] [PubMed]
  16. Sheng, C.; Shen, X. Simulation of acoustic agglomeration processes of poly-disperse solid particles. Aerosol Science and Technology 2007, 41, 1–13. [Google Scholar] [CrossRef]
  17. König, W. Hydrodynamisch-akustische Untersuchungen. Über die Kräfte zwischen zwei Kugeln in einer schwingenden Flüssigkeit und über die Entstehung der Kundtschen Staubfiguren. Ann. Phys. und Chem 1891, 42, 549. [Google Scholar] [CrossRef]
  18. Brandt, O. , Freund, H. , Hiedemann, E. Zur Theorie der akustischen Koagulation. Kolloid-Zeitschrift. 1936, 77, 103–115. [Google Scholar]
  19. Khmelev, V.N.; Shalunov, A.V.; Nesterov, V.A. Improving the separation efficient of particles smaller than 2.5 micrometer by combining ultrasonic agglomeration and swirling flow techniques. PLoS ONE 2020, 15, 1–16. [Google Scholar] [CrossRef] [PubMed]
  20. Gallego-Juarez, J.A.; Riera-Franco De Sarabia, E.; Rodriguez-Corral, G.; Hoffmann, T.L.; Gálvez-Moraleda, J.C.; RodrÍguez-Maroto, J.J.; Acha, M. Application of Acoustic Agglomeration to Reduce Fine Particle Emissions from Coal Combustion Plants. Environmental Science & Technology 1999, 33, 3843–3849. [Google Scholar] [CrossRef]
  21. Uchiyama, Shinnosuke & Asami, Takuya & Miura, Hikaru. Development of an aerial ultrasonic sound source with a truncated cone-shaped reflective plate on a circular transverse vibrating plate. Acoustical Science and Technology 2022, 43, 43–49. [Google Scholar] [CrossRef]
  22. Motoi, R.; Asami, T.; Miura, H. Agglomeration of aerosol using intense standing wave field of cylindrical shape. In Proceedings of the 23rd International Congress on Acoustics : Integrating 4th EAA Euroregio 2019, Aachen, Germany, 9-13 September 2019. [Google Scholar] [CrossRef]
  23. Nyborg, W.L. Acoustic Streaming. Phys. Acoustics 1965, 2, 265–331. [Google Scholar]
  24. Lighthill, Sir J. Acoustic Streaming. J. Sound Vib 1978, 61, 391–418. [Google Scholar] [CrossRef]
  25. Hamilton, M.F.; AIlinskii, Y.; Zabolotskaya, E.A. Acoustic Streaming Generated by Standing Waves in Two-Dimensional Channels of Arbitrary Width. J. Acoust. Soc. Am. 2003, 113, 153–160. [Google Scholar] [CrossRef] [PubMed]
  26. Aktas, M.K.; Farouk, B. Numerical Simulation of Acoustic Streaming Generated by Finite-Amplitude Resonant Oscillations in an Enclosure. J. Acoust. Soc. Am. 2004, 116, 2822–2831. [Google Scholar] [CrossRef]
  27. Yong, S.; Sangmo, K. Acoustic Streaming. Encyclopedia of Microfluidics and Nanofluidics 2014, 25–33. [Google Scholar] [CrossRef]
  28. Bochenkov, A.S.; Golykh, R.N.; Nesterov, V.A.; Tertishnikov, P.P.; Khmelev, V.N.; Tsyganok, S.N.; Shalunov, A.V. Method of ultrasonic. coagulation. Patent RF2759506, B01D 51/08, 15 November 2021. [Google Scholar]
  29. Khmelev, V.N.; Shalunov, A.V.; Nesterov, V.A. Improving the performance of air purification efficiency from fine-dispersed particles by ultrasonic exposure in swirling flow. Int. J. Environ. Sci. Technol 2020, 17, 3927–3934. [Google Scholar] [CrossRef]
  30. Khmelev, V.N.; Shalunov, A.V.; Nesterov, V.A. Summation of high-frequency Langevin transducers vibrations for increasing of ultrasonic radiator power. Ultrasonics 2021, 114, 1–12. [Google Scholar] [CrossRef]
  31. Baillot, F.; Blaisot, J.; Boisdron, G.; Dumouchel, C. Behaviour of an air-assisted jet submitted to a transverse high-frequency acoustic field. Journal of Fluid Mechanics 2009, 640, 305–342. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the ultrasonic agglomeration chamber: 1 – ultrasonic vibrating system; 2 – piezoelectric transducer; 3 - concentrator; 4 – disk radiator; 5 - reflector;6 – inlet tube; 7 – polluted gas; 8 – outgoing gas-dispersion flow (to cyclone); 9 – aggregated particles; 10 – outlet tubes; 11 – amplitude distribution; 12 – vortex flows; 13 – housing; 14 – radiator flange; λG – distance between radiator and reflector.
Figure 1. Schematic diagram of the ultrasonic agglomeration chamber: 1 – ultrasonic vibrating system; 2 – piezoelectric transducer; 3 - concentrator; 4 – disk radiator; 5 - reflector;6 – inlet tube; 7 – polluted gas; 8 – outgoing gas-dispersion flow (to cyclone); 9 – aggregated particles; 10 – outlet tubes; 11 – amplitude distribution; 12 – vortex flows; 13 – housing; 14 – radiator flange; λG – distance between radiator and reflector.
Preprints 90332 g001
Figure 2. Diagrams of generation of ultrasonic radiation using bending and oscillating ultrasonic disk radiators: (a) radiator with a flat frontal surface (flat radiator). (b) radiator with phase-aligning collars (step-function radiator).
Figure 2. Diagrams of generation of ultrasonic radiation using bending and oscillating ultrasonic disk radiators: (a) radiator with a flat frontal surface (flat radiator). (b) radiator with phase-aligning collars (step-function radiator).
Preprints 90332 g002
Figure 3. Shape of vibrations and photo of ultrasonic disk radiators: (a) flat radiator. (b) step-function radiator.
Figure 3. Shape of vibrations and photo of ultrasonic disk radiators: (a) flat radiator. (b) step-function radiator.
Preprints 90332 g003
Figure 4. Sound pressure level distribution inside the agglomeration chamber: (a) flat radiator; (b) step-function radiator.
Figure 4. Sound pressure level distribution inside the agglomeration chamber: (a) flat radiator; (b) step-function radiator.
Preprints 90332 g004
Figure 5. Particle distribution in the agglomeration chamber: (a) flat radiator; (b) step-function radiator.
Figure 5. Particle distribution in the agglomeration chamber: (a) flat radiator; (b) step-function radiator.
Preprints 90332 g005
Figure 6. Diagram of particle paths in the agglomeration chamber when using a flat radiator: 1 – chamber body; 2 – radiator; 3 – reflector; 4 – inlet tube; 5 – annular domain for emission of the sounded gas-disperse flow; 6 – areas of agglomerates formation; 7 – gas-dispersion flow; 8 – distribution of amplitudes of vibrations of the disk radiator; 9 – zeros of vibrations of the radiator; 10 – areas of maximum vibration amplitudes; 11 – vortex areas.
Figure 6. Diagram of particle paths in the agglomeration chamber when using a flat radiator: 1 – chamber body; 2 – radiator; 3 – reflector; 4 – inlet tube; 5 – annular domain for emission of the sounded gas-disperse flow; 6 – areas of agglomerates formation; 7 – gas-dispersion flow; 8 – distribution of amplitudes of vibrations of the disk radiator; 9 – zeros of vibrations of the radiator; 10 – areas of maximum vibration amplitudes; 11 – vortex areas.
Preprints 90332 g006
Figure 7. Diagram of particle paths in the agglomeration chamber when using a step-function radiator: 1 – chamber body; 2 – radiator; 3 – reflector; 4 – inlet tube; 5 – annular domain for emission of the sounded gas-disperse flow; 6 – areas of agglomerates formation; 7 – gas-dispersion flow; 8 – distribution of amplitudes of vibrations of the disk radiator; 9 – zeros of vibrations of the radiator; 10 – areas of maximum vibration amplitudes.
Figure 7. Diagram of particle paths in the agglomeration chamber when using a step-function radiator: 1 – chamber body; 2 – radiator; 3 – reflector; 4 – inlet tube; 5 – annular domain for emission of the sounded gas-disperse flow; 6 – areas of agglomerates formation; 7 – gas-dispersion flow; 8 – distribution of amplitudes of vibrations of the disk radiator; 9 – zeros of vibrations of the radiator; 10 – areas of maximum vibration amplitudes.
Preprints 90332 g007
Figure 8. Dependences of particle capture efficiency on gas flow rate (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the flow rate (B): 1 – flat radiator; 2 – step radiator; 3 – without ultrasonic influence.
Figure 8. Dependences of particle capture efficiency on gas flow rate (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the flow rate (B): 1 – flat radiator; 2 – step radiator; 3 – without ultrasonic influence.
Preprints 90332 g008aPreprints 90332 g008b
Figure 9. Dependence of the capture efficiency on the sound pressure level (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the sound pressure level (B): 1 – flat radiator; 2 – step radiator.
Figure 9. Dependence of the capture efficiency on the sound pressure level (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the sound pressure level (B): 1 – flat radiator; 2 – step radiator.
Preprints 90332 g009aPreprints 90332 g009b
Figure 10. Dependence of capture efficiency on particle mass concentration (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the mass concentration of particles (B): 1 – flat radiator; 2 – step radiator.
Figure 10. Dependence of capture efficiency on particle mass concentration (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the mass concentration of particles (B): 1 – flat radiator; 2 – step radiator.
Preprints 90332 g010
Figure 11. Fractional efficiency (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the particle size at the inlet of coagulation chamber (B): 1 – flat radiator; 2 – step radiator; 3 – without ultrasonic influence.
Figure 11. Fractional efficiency (A) and the dependence of the particle size at the outlet of the agglomeration chamber on the particle size at the inlet of coagulation chamber (B): 1 – flat radiator; 2 – step radiator; 3 – without ultrasonic influence.
Preprints 90332 g011aPreprints 90332 g011b
Figure 12. – Design of the ultrasonic agglomeration chamber: 1 – housing; 2 – inlet pipe; 3 – outlet pipe; 4 – radiator; 5 – ultrasonic piezoelectric transducer; 6 – reflector; 7 – flow distributors; N = 1…3; λ – wavelength.
Figure 12. – Design of the ultrasonic agglomeration chamber: 1 – housing; 2 – inlet pipe; 3 – outlet pipe; 4 – radiator; 5 – ultrasonic piezoelectric transducer; 6 – reflector; 7 – flow distributors; N = 1…3; λ – wavelength.
Preprints 90332 g012
Figure 13. – Ultrasonic oscillatory system: A – Sketch of an ultrasonic oscillatory system; B – vibration mode of the radiating element; 1 – radiating element; 2 – piezoelectric transducer; 3 – housing; 4 – fan; 5 – flange; 6 – power cable.
Figure 13. – Ultrasonic oscillatory system: A – Sketch of an ultrasonic oscillatory system; B – vibration mode of the radiating element; 1 – radiating element; 2 – piezoelectric transducer; 3 – housing; 4 – fan; 5 – flange; 6 – power cable.
Preprints 90332 g013
Table 1. Specifications of ultrasonic transducers.
Table 1. Specifications of ultrasonic transducers.
Disk transducer type Flat radiator Step-function radiator
Diameter of the radiator, mm 320 320
Power consumption, W 245 240
Oscillation frequency, kHz 22.1 22.2
Surface oscillation amplitude max/min 51/45 51/47
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated