Environmental dust repelling from hydrophilic/hydrophobic surfaces under sonic excitations

Dust repelling from transparent polyvinyl chloride film surface via sonic excitation is examined and dynamics of repelled (inflight) dust particles are analyzed. An experimental rig is designed and built to assess the vibrational characteristics of the polyvinyl chloride film at different frequencies of sonic excitation. A high speed recording system and tracking program are utilized monitoring and evaluating the dynamics of the inflight particles. The dynamics of inflight particles are also simulated numerically and the predictions are compared with those of the experimental data. In order to examine the influence of dust particle adhesion on the dynamics of the inflight particles, the polyvinyl chloride film surface is hydrophobized through dip coating by functionalized nano-silica particles. Improvement of the optical transmittance of the dust mitigated film is determined via outdoor tests. The findings demonstrate that sonic excitation repels the particles from the film surface and it is more pronounced at 64 Hz excitation frequency while demonstrating that sonic excitation can be used for dust removal from transparent surfaces. The mitigation via sonic excitation improves the optical transmittance of the dusty surface by 77%, which becomes more apparent for hydrophobic surfaces.

method provided the effective cleaning of small particles, the scarcity of clean water in rural areas could limit the practical application of such arrangements in large scale solar energy harvesting farms. The use of sonic radiation could be extended to include fouling in tubing systems 13 and residuals removal in manufacturing, particularly additive manufacturing system 9,12 . In any case, the use of clean water becomes a necessity for such cleaning applications. Moreover, introducing mechanical excitations on surfaces via piezoelectric actuators and mechanical accelerators could clean the surfaces in ambiances where airborne particles are high present at high concentrations in air. However, it became necessary to use the additional mechanical system(s) to remove the repelled particles from surfaces such mechanical brushes 8 or to utilize gravitational potential, via surface inclination, to displace the repelled particles from surfaces 14 . The vibrational excitation of the antistatic coated surface could ease the mitigation of the dust particles from surfaces; however, the vibrational characteristics of the surface remained important for the efficient cleaning process. The combination of strong acoustic wave (sound wave) and mechanical vibration could be used effectively to clean the porous surfaces, which became particularly important for producing the biomedical parts 15 . However, practical applicability of such arrangements for dust mitigation is questionable because of relatively smooth surfaces involved in energy harvesting equipment. Moreover, the acoustic excitation towards removal of deposits from surfaces could be favorable for heat transfer enhancement 16 ; however, the sound waves generated could disperse and propagate in all directions in the liquid while resulting in loss of wave intensity, which lowered the mitigation of deposits from surfaces. Nevertheless, sound wave excitation in air ambient for the mitigation of environmental dust needs to be explored for practical applications.
On the other hand, environmental dust contains various elements with salt and oxide compounds. Some of the slat compounds in the dust particles do not conform stoichiometric ratio while creating additional ionic charges on the particle surfaces 17 . This adds to the interfacial adhesion, due to the van der Walls forces, between the particles and the settled surfaces. In addition, the charged particles form clusters-like edifices and they can attach onto the dust particle surfaces with large sizes 18 . The combination of particle clustering and ionic forces considerably enhances the work of adhesion required for dust removal from surfaces, i.e. increased contact area at the interface on the settled surface and strong interfacial forces enhance particle adhesion. One of the methods reducing the individual and clustered particle adhesion on the settled surfaces is to reduce the contact area between the particles and the settled surfaces. This can be achieved, possible, through proper texturing of the surface 19 . The surface texture with hierarchically distributed mico/nanopillars are favorable because such texture topology gives rise to hydrophobic state on the surfaces 20 . Hydrophobizing the surface lowers the transmittance of the optically transparent samples; however, overall reduction in optical transmittance of hydrophobized surface over the visible spectrum is almost 10% of the none-hydrophobized surface 21 . In environments subjected to frequent heavy dust settlements, hydrophobizing the surfaces eases dust mitigation from surfaces; in which case, overall optical transmittance can improve. Moreover, reducing the surface free energy of the textures causes a further reduction in the adhesion of the particles on surfaces. Hydrophizing the surface having low surface free energy, via texturing and chemical modification, can ease the mechanical repelling of dust particles from the surface due to low adhesion 18 . Hence, lowering the free energy of textured surfaces provides better opportunities to mitigate the environmental dust particles surfaces via acoustic excitations in terms of sonic waves. Acoustic excitations for dust mitigation have successfully applied in the previous studies to remove soot and ashes from tune banks, particularly at high-temperature applications 22 . The low frequency and high-intensity sound waves were effectively used removing the charged particles (powders) from the collection plate of an electrostatic precipitator 22 . Although acoustic levitation of dust particles enables are swept by the airflow, the flow forming a vortex structure is required for mitigating the levitated particles from the surfaces 23,24 . However, generating such flow structures sustaining the cleaning of large areas for prolonged durations remains difficult to achieve. In addition, the particle removal under sonic excitation was studied earlier and the focus was the removal of tooth 22 or charged particles 23 . However, dust particles have larger densities than sooth and low static charges unlike those of the cases reported in the early work 22 . Hence, the removal of dust particles from surfaces in dry environments under the sonic influence becomes interestingly new research and accomplishments of dust mitigation from surfaces by sonic excitation becomes fruitful despite the fact that the practical applications of sonic excitation for dust mitigation in: i) cleaning of sensor surfaces for autonomous systems, and ii) efficient operation of solar energy harvesting devices require further investigations. Nevertheless, the present work is proposed to investigate the dust removal from hydrophilic/hydrophobic and optically transparent polyvinyl chloride film surfaces using the sonic excitations pertinent to cleaning applications of protective layer for photovoltaic applications. The sample (polyvinyl chloride film) surfaces are hydrophobized by depositing functionalized nano-silica particles via dip coating. An experimental rig is designed and built for sonic excitations and dust removal. A high speed recording system is utilized monitoring and tracking the repelling dust particles from surfaces during the sonic excitations. An analytical approach is introduced to formulate the repelled dust particle dynamics and resulting predictions are compared with those of the experimental findings. The outdoor tests are carried out to evaluate the optical transmittance of dust repelled surfaces.

Experimental
Polyvinyl chloride film with 0.14 mm thickness and 120 mm diameter were used as samples. The sample surfaces were hydrophobized via a dip-coating technique through depositing the functionalized nano-silica-particles. The nanoparticles were synthesized in accordance with the early work 25 . The wetting of the coated surface was determined via goniometer (Kyowa, model DM 501) as similar to the previous work 26 . A high resolution camera and software were incorporated in the goniometer to capture and analyze the contact angle. De-ionized water was used in the measurements and droplet volume was controlled with an automatic dispensing system having a volume step resolution of 0.1 µL. The sessile droplet contact angle was evaluated using the high-precision  26 . Hence, the contact angle of the coating surface was 152° ± 3° with hysteresis 5° ± 2°. In addition, the contact angle of as received polyvinyl chloride film was measured as 82° ± 2° with contact angle hysteresis 41° ± 3°. The coating texture topology and dust particle adhesion on hydrophobic/hydrophilic surfaces are obtained using an atomic force microscope (AFM/SPM) probe with the friction mode. A fixture with two-axes freedom was built via using a 3D printer. The fixture accommodates both the circular sample holder and the loudspeaker with foam insulation in between them. Figure 1 shows the optical image. The loudspeaker (Edifier Inc.) operating at 9 V (DC) and 0.44 A was used to generate the sound waves at various frequencies. To minimize the mechanical disturbance of the polyvinyl chloride film by the loudspeaker mechanical vibration, vibration-reducing foams were used to isolate the loudspeaker from the polyvinyl chloride film holder. An accelerometer was incorporated ensuring the frequency measurement of polyvinyl chloride film during the sonic excitation by the loudspeaker at various frequencies. Initially, many tests were conducted securing the measurement repeatability and fixture stability at various frequencies of the sonic excitations. A highspeed camera (Speed Sense 9040) was utilized monitoring film oscillations when subjected to the sonic excitations and the movement of the repelled dust particles from the film surface. In the dust particles repelling experiments, www.nature.com/scientificreports/ the dust layer of almost 150 µm was deposited onto the polyvinyl chloride film surface. Initially, many tests were carried out calibrating the film response to the sonic excitations and dust particles' motion in terms of vertical and lateral motions. The tracker program was used to evaluate the dust particle dynamics from the recorded data. The high-speed records were obtained at 5,000 frames-per-second (fps) at a resolution of 1280 × 800 pixels. The pixel size of the images was 14 µm × 14 µm. The repeatability of the recording tests was assessed and the standard error estimated was in the order of 3%. The uncertainty (± u) of the measurements was evaluated incorporating the data measured (vertical and horizontal repelling heights of the particles). The confidence level of 95% was ensured basing the data repeatability. The uncertainty (σ u ) is 27 : σ u = x n x o (x − µ e ) 2 p(x)dx , here, µ e represents the mean/expected value of variable x, n corresponds to the number of points in the recorded data set, and p(x) resembles the probability distribution function. The probability distribution function was fitted in a Gaussian-function, which enabled to estimate the probability distribution function diameter. The standard uncertainty was obtained using the Gaussian fitting function. The number of pixels contributing to the crosscorrelation-peak was normalized for consistency of the recorded data in terms of pixels. The bias error was estimated to be about 0.5 pixels because of the complexity of the evaluation of small peaks in terms of sizing in the probability distribution function. The standard uncertainty was evaluated as 3%.
Dust was gathered from photovoltaic surfaces using soft brushes and kept in sampling bottles. They were characterized by utilizing scanning electron microscopy and energy dispersive spectroscopy (JEOL 6460), and X-ray diffraction (Bruker D8). The dust particles size distribution was evaluated using the particle size analyzer (Malvern Panalytical, Mastersizer 3000), which enabled to classify the size of the dust particles within 10-3.5 mm range through incorporating both red and blue light wavelengths.

Mathematical analysis
The mathematical analysis covers the film vibration under the sonic excitation and the repelled dust particle dynamics resulted from the film vibration.
Film vibration under sonic excitations. The circular polyvinyl chloride film is located above the sonic loudspeaker with a standoff distance (spacing between the loudspeaker and the film, Fig. 1) and sonic excitation of the film gives rise to vibrational motion of the film depending on the sonic excitation frequency and the amplitude. The film has a considerably low thickness (14 µm) and it possesses negligible flexural stiffness. The transverse deflection of the film because of the natural frequency of vibration can be formulated from the wave equation. It is worth to mention that the film has a circular shape and it is subjected to a nearly uniform radial force per unit length of in-plane (radial) tensile force ( T 0 ) at the film edges. Hence using the cylindrical coordinate system in line with Fig. 2, the equation governing the transverse deflection yields: Here, c = T 0 ρh is the speed of propagation of the transverse wave on the film surface, ρ is the film density, t is time, h is film thickness, w n (r, θ, t) is the transverse deflection. The boundary conditions should satisfy that, at the edges ( r = R ), w n (R, θ , t) = 0 , i.e. zero transverse displacements at the film edges. In addition, initially, the film is considered to be at rest, i.e. w n (r, θ, 0) = 0 and as the time approaches infinity time derivative of transverse displacement approaches to zero, i.e. ∂w n (r,θ,t) ∂t = 0 . The mathematical arrangement of the solution of Eq. (1) is provided in the Supplement (S1). The solution yields: Here, A 1 , B 1 , A 2 and B 2 are constants, ω = c is the angular frequency, m = 0, 1, 2, . . . . are constants that represent the number of diametral lines with zero deflection. J m is the Bessel functions of zero order. The natural frequency of the vibrating film is: Here: k mn can be obtained from the Bessel function of the first kind, m is an integer that represents the number of circumferential lines with zero deflection and n is an integer that represents the number of diametral lines with zero deflection. Hence, m and n can be used to characterize the modes' shape. It is worth to mention that: R(r) = J m ( mn r) for m = 0, 1, . . . , n = 1, 2, . . . and mn = k mn R and k mn is the n-th positive root of J m . The mode shapes of the vibrating film are also modeled numerically via solving Eq. (1) with the boundary conditions defined at r = R then w n (R, θ, t) = 0 . The Ordinary Differential Equation solver of the COMSOL Multiphysics finite element code is used to obtain the numerical solution. Since the vibrating film has a small thickness (low dimension as compared to its radius), it is represented as a 2D elastic film having properties given in Table 1. In the numerical solution, a uniform radial tensile force per unit length of 1.12 is applied at the edges of the film and the transverse deflection is constrained at the edges. The film is meshed with 5418 triangular elements and the grid independence tests demonstrate that the number of elements leads to the converged solution of Eigenvalues.
Since the film is excited by the sound waves via a loudspeaker, the forced vibrational analysis of the film needs to be considered. Hence, after obtaining the mode shapes of the film, a pulsating (sinusoidal) pressure is applied at the film bottom surface resembling the sound waves emanating from the loudspeaker. Generally, sound waves www.nature.com/scientificreports/ are considered to be pressure waves resulting from the change in pressure from that of ambient. The propagating sine wave can be represented as: Here, P = P − P atm is the change in pressure, P max is the amplitude (or maximum) pressure change, k = 2π is the wavenumber, ω = 2π T = 2πf is the angular frequency, x is spatial coordinate, t is time. Hence, the transverse deflection of the film under the forced vibration, w f (r, θ, t) , can be expressed as: Since the experiment is carried out at low-frequency ranges, the vibrational mode (0,1) dominates and the solution to w f (r, θ, t) can be considered to be axisymmetric. Hence, Eq. (5) becomes a function of r and θ , i.e.: T o is in-plane (radial) tensile force, P eff is the effective pressure created by sonic excitation, and transverse deflection of film.   6) can be expressed as 28 : Because of the consideration of consistency of numerical vibrational mode assessment of the film, the film displacement is also predicted numerically by solving Eq. (6) incorporating the conditions. The COMSOL Multiphysics Differential Equation Solver is used to solve Eq. (6). The vibrating film has a low thickness and it is represented as a 2D elastic film with properties given in Table 1. In addition, a sinusoidally varying force (source term), �P max sin(2πft) is introduced from the film bottom surface. The transverse deflection is constrained at the edges and the film is meshed with 5418 triangular elements. The grid independence tests are carried out securing the grid-independent results. Dynamic analysis of repelling dust particles. The dust particles are repelled from the film surface under the influence of film acceleration because of the sonic excitations. The inflight particle (repelled dust particles from the film surface) motion can be presented in the spherical coordinate system ( r, θ , φ, t ). Figure 3 shows schematically the film and inflight particle in the coordinate system. The formulation of the forces acting on the inflight particle is given in the Supplement (S2). Hence, the forces acting on the inflight dust particle at the onset of repelling from the film surface are: Here, F i is the inertia force, m is the particle mass, D is the drag force term according to Stokes' hypothesis, d p is the particle diameter, and F ac is the acoustic force term. Moreover, the airflow around the dust particle occurs at low Reynold number (Re) and the Stokes' hypothesis can be adopted formulating the drag force, i.e.: D = 3πµd p . In addition, the acoustic radiation force on a particle moving in a viscous fluid is formulated previously 29 and this formulation is used. The acoustic force yields F ac = π d 3 p 8 kfE ac , here k is the wave number, rms is acoustic energy density and k is the dipole scattering coefficient. The dipole scattering coefficient is taken as: f = , where, ρ = ρ p ρ a is the particle-to-air density ratio and δ = δ b r p www.nature.com/scientificreports/ is the ratio of boundary layer and particle radius. The boundary layer developed on a spherical particle can be obtained from δ b = 4.53d p √ Re 29 . The force term in Eqs. (8)(9)(10) can be formulated in terms of particle mass and acceleration, i.e.: m p where m p is the particle mass. In consistency with the solutions for the film vibrational motion, the numerical solution of Eqs. (8)(9)(10) can obtained using the Differential Equation solver of COMSOL Multiphysics finite element code. It is worth to mention that the second order Euler backward difference scheme is used to discretize the equations and the nonlinear solution is obtained with the aid of the Newton-Rapson method. Since the accuracy of the predictions is limited by the selection of the time increment, the time step is set at 10 −8 s in the numerical simulations. Table 2 gives the parameters used in the simulations.

Results and discussion
Environmental dust characteristics and dust particles mitigation from transparent polyvinyl chloride film are examined. The mechanisms of the dust particles repelling via sonic excitation are explored and the dynamics of the repelled particles are formulated. The findings of the repelled particle heights are compared with their counterparts obtained from the experiments. The optical transmittance of the dust mitigated surfaces is also evaluated.
Dust and surface characteristics. Figure 4a shows SEM microimages and Fig. 4b depicts the AFM line scan of the hydrophobized film surface. The silica particles have sizes almost 30 nm (Fig. 4a) and form a clustered layer on the sample surface. Small texture height is noted along the line scan (Fig. 4b) and the average roughness is about 40 nm. The peaks and valleys in the texture of the surface (Fig. 4b) demonstrates that silica particles clustered while forming small peaks on the surface, which creates a Lotus effect on the liquid droplets while reducing droplet contact angle hysteresis. The coated surface wetting is measured and the contact angle of the coated surface is about 150° ± 2° and hysteresis is 4° ± 2°, i.e. coating surface demonstrates superhydrophobic identity. On the other hand, dust is collected from PV panel surfaces in Dhahran, Saudi Arabia, and characterization tools are used to evaluate particle size, shapes, and elemental constitutes. Figure 5a,b show SEM microimages of dust particles while Fig. 5c shows dust particle distribution. Dust have different shapes and sizes (Fig. 5a). In addition, particularly, small particles agglomerate forming the clusters (Fig. 5b). The agglomeration and adherence of small particles are because of ionic forces formed on these particles. Elemental constitutes obtained from EDS analysis demonstrate that dust retains several elements, such as Si, Na, Ca, K, S, Fe, Cl, and O ( Table 3). The quantification of elemental constitutes does satisfy the stoichiometric ratio for small particles (2 µm ≤) particularly alkaline salt components (NaCl and KCl) as observed in Table 3, which is also reported in the earlier work 17 . This contributes to dust agglomeration (Fig. 5b). Dust shapes can be evaluated based on the shape factor ( A sh = P 2 4π A , where P is perimeter of the dust particle and A is area of dust cross-section) and aspect ratio ( A s = π(L l ) 2 4A , here L l is the largest projection length) 30 . The microscopic methods are utilized to determine the feature of dust geometry. The shape factor, mostly, ranges from 0.5 to 3. and one corresponds to dust with the almost circular feature. There is no clear distribution of shape factor with a range of dust sizes of 2-5 µm. As dust becomes greater than 5 µm, the shape factor becomes almost 3. However, it becomes almost one for small size dust (≤ 0.8 µm). The aspect ratio varies considerably with dust sizes and aspect ratio approaches almost one for small size dusts ((≤ 0.8 µm). In addition, a simple mathematical expression correlating the shape factor and the aspect ratio in terms of dust sizes could not obtained. The particle size of the dust varies in micrometer to nanometer ranges with an average of 1.2 µm (Fig. 5c). Figure 6 demonstrates the X-ray diffraction data for dust. The salt compounds (NaCl and KCl), calcite, silica peaks are apparent. The peaks of iron and silicon overlap and sulfur peak could be anhydrite or gypsum (CaSO 4 ) while iron (Fe) peak corresponds to clay-aggregated hematite (Fe 2 O 3 ). Moreover, dust adhesion on coated sample surfaces is calculated adopting the technique developed earlier 31 . The AFM probe deflection in friction mode is utilized determining the dust particle adhesion on the sample surface. Hence, F = kσ d �V , here k is spring constant (N/m) of probe tip, σ d is slope of probe deflection (Δz/ΔV, m/V), and ΔV is probe output recorded (mV) onset of deflection. The AFM probe used has the characteristics of kσ d = 5.80275 × 10 −15 N/mV. The particle of about 0.9 µm size on the hydrophilic sample surface (as received), the output probe voltage is recorded as 380 mV and equation ( F = kσ d �V ) yields the adhesion force of 2.21 × 10 −12 N. Similarly, for about 8 µm size particle, the probe output is 320 mV, which yields the adhesion force of 1.86 × 10 −12 N. Hence, the adhesion force increases almost 16% as the particle size reduces from 8 to 0.9 µm, which demonstrates the increased adhesion of small particles on the sample surface. The adhesion force calculations are repeated for the hydrophobized sample surfaces. The probe output reading for 0.9 µm size particle on the hydrophobized (coated) sample surface is 205 mV, which gives the adhesion force of 1.45 × 10 −12 N, Table 2. Parameters used in the particle dynamic analysis.

Parameter Value
Density of air, ρ kg m 3 1. www.nature.com/scientificreports/ and for 10 µm size particle on coated sample surface yields the probe output of 170 mV, i.e. the adhesion force for 10 µm size particle is 0.98 × 10 −12 N. Therefore, hydrophobizing the sample surface reduces the particle adhesion force by almost 35% for 0.9 µm particle and about 47% for 8 µm particle. Consequently, hydrophobizing the surface becomes more effective towards reducing the particle adhesion on the surfaces, particularly for large size particles. Adhesion force experiments are repeated 10 times to secure accurate data. The standard deviation where N is the number of repeats, x i is the output recorded (mV), − x is the mean value of probe output) is estimated at about 7%.
Dynamics of repelling dust particles. The polyvinyl chloride (PVC) circular thin film with 14 µm thickness is excited by the sonic waves (sonic radiation) within 35-300 Hz and 80-100 dB amplitude from the bottom of the film while the particles are located on the top surface of the film (Fig. 1). The film possesses very negligible flexural stiffness and the mode of deflection of the film can be obtained from Eqs. (2) and (3). Moreover, initially, the assessment of the mode of film vibration is carried out both numerically (Eqs. 2 and 3) and experimentally. Figure 7a shows the mode of vibration of the film at two different sonic excitation frequencies (64 Hz and 104 Hz) while Fig. 7b depicts the corresponding modes obtained from the experiments. Tables 4  and 5 give the data for the mode of shapes and natural frequencies of film vibration. At low excitation frequencies (45 Hz-85 Hz), single-mode results while at high excitation frequency four-mode results. Since the dust particles have small sizes, the mode of film vibration results in the clustering of the dust particles on the film surface (Fig. 7b). The transverse displacement (along the y-axis) of the film surface is evaluated experimentally using high speed camera data incorporating the tracker program at low frequencies one mode shape (1,0) of the film vibration (Fig. 7b). Figure 8 shows the transverse displacement of the film surface with different excitation sound frequencies and resulting sound power. The displacement remains maximum for the frequency of 64 Hz because of the resonant frequency. To evaluate the effect of the location of the sound excitation reference to the film location (standoff distance), various tests are conducted. Figure 9 shows the sound power variation with the frequency. The local peak sound power occurs at a frequency of 64 Hz. Further tests are carried out to assess the influence of the input voltage of the sound wave generator on the dynamic characteristics of the film. Table 6 gives the data for film dynamic characteristics. In addition, the amplitude of film vibration at various frequencies is also obtained from the analytical (Eq. 7) and numerical (Eq. 6, COMSOL simulations) approach and the values are included in Table 6. The findings reveal that the input voltage of 8.4 V for loudspeakers with d = 30 mm standoff distance (Fig. 1) results in the maximum sound power, which is used for dust mitigation study, i.e. the sound power 104.9 dB is resulted (Table 6). Moreover, the velocity of the film in the transverse direction is determined from experiments and the findings are shown in Fig. 10. To evaluate the influence of sonic excitation frequency on the film transverse velocity, three frequencies are included in Fig. 10 for comparison. The velocity of the film in the transverse direction follows the oscillatory motion of the film with the same frequency. The maximum amplitude of the transverse velocity is about 1 m/s, which occurs at 64 Hz of the external excitation frequency. It is worth mentioning that the sound power is the maximum for 64 Hz excitation frequency with a standoff distance of 30 mm (Table 6).
Under the sonic excitation, rapid displacement of the film causes dust particle acceleration. The trajectory of the repelled particles are formulated through Eqs. (8)(9)(10) and predicted using Differential Equation solver (DE) of COMSOL multi-physics code. Figures 11 and 12 show the temporal behavior of the inflight particle vertical (Fig. 11) and horizontal heights (Fig. 12) obtained from the experimental data using the tracker program (incorporating the high-speed recorded data for hydrophilic surface) and predicted from the numerical simulations, at which the particle is located at (200, 5.74°, 5°) on the film surface for four sonic excitation frequencies (45 Hz, 64 Hz, 75 Hz, and 104 Hz), respectively. Here, the coordinate system ( r, θ , φ ) is according to Fig. 3. It is worth to mention that as the particle acceleration reaches a value that takes over the adhesion and weight of the particle, the particle can depart from the film surface. The vertical height of the inflight dust particle follows a parabolic rise and fall with time. This is because of the force balance between the particle inertia, gravity, and drag forces. www.nature.com/scientificreports/ The gravitational force over the drag force of the dust particle is in the order of 10 −3 which demonstrates that the gravitational force has a major influence on the inflight particle deceleration with time. The maximum vertical height of the inflight dust particle occurs at about 0.02 s from its departure for 45 Hz excitation; however, the  www.nature.com/scientificreports/ horizontal location continues to increase with time for a longer duration (0.1 s). Moreover, as the excitation frequency changes, the time occurrence of the maximum vertical displacement differs, i.e. 0.1 s for 64 Hz and 0.08 s for 104 Hz of excitations. This indicates that the inflight time of the repelled particle remains longer in the air for 64 Hz excitation. Hence, the repelled particle trajectory has two-stages such that in the first stage, the particle rises both vertically and horizontally under the influence of the sonic excitation while in the second stage the particle follows a falling trend from its maximum vertical height. This allows the repelled particles landing on the film surface significantly away from their repelling locations on the film, which becomes more apparent for 64 Hz of sonic excitation. This indicates that the particle can be removed from surface of the film while creating the multiple sonic excitations. For inflight particle velocity, the horizontal component of the particle increases while the vertical velocity component becomes zero at the maximum peak location, particularly at 64 Hz excitation frequency. As comparing the experimental findings with the predictions of the inflight particle heights in vertical and horizontal directions, both results are in good agreement. The small differences are because of the experimental errors (7%) and the consideration of uniform density particles, which may change slightly from particle to particle because of varying elemental composition (Table 3). Nevertheless, both findings are  www.nature.com/scientificreports/ in agreement. Figure 13a,b show an optical image of the inflight particles at various times and different sonic excitation frequencies for hydrophilic and hydrophobic film surfaces, respectively. The tracking of a particle at different time steps is marked in Fig. 13a,b. It can be observed that at low (45 Hz) and high (105 Hz) sonic excitation frequencies, the heights of the inflight particles are lower than that of the 64 Hz. In addition, at these frequencies, some particles remain on the film surface, i.e. particle inertia force remains less than the adhesion and weight forces of the particles. This situation can also be seen from Fig. 14a,b, in which the temporal variation of horizontal and vertical heights of the inflight dust particle is shown for various sonic excitation frequencies.
In the case of the hydrophobic surfaces, vertical and horizontal heights of the dust particles change. Figure 15 shows the temporal variation of horizontal and vertical heights of the inflight particle at the initial location (200, 5.74°, 5°) on the hydrophilic and hydrophobic film surfaces for different sonic excitation frequencies (45 Hz, 64 Hz, and 75 Hz). It is worth to note that inflight particle heights are obtained experimentally. The inflight particle height on the hydrophobic film becomes larger than its counterpart corresponding to the hydrophilic surface at low excitation frequency 45 Hz. As the excitation frequency increases, the inflight particle height difference becomes small along the vertical line (normal to the film surface); however, some small decrease in the horizontal height of the particle is observed for the hydrophobic surface. Nevertheless, the inflight particle heights in both vertical and horizontal directions become similar to whether the particle is initially located on the hydrophobic or hydrophilic film surfaces as excitation frequency increases to 64 Hz. This is attributed to the inertial force generated on the particle located on the film surface. As the frequency reduces (45 Hz), the film displacement and velocity reduce significantly in vertical and horizontal directions. This lowers the inertial force created on the particle to be repelled from the film surface. Hence, at low excitation frequency, the influence of the particle adhesion force on the particle acceleration becomes critically important on the particle inflight dynamics for the case of the hydrophobic surface. As the excitation frequency increase 64 Hz, the inertial force generated on the particle becomes much larger than the adhesion force. Therefore, the inflight dynamics of the repelled particle becomes independent of the surface wetting characteristics, i.e. hydrophobic or hydrophilic. The optical transmittance of the hydrophobic and hydrophilic film surfaces is carried out after the dust repelling experiments incorporating the sonic excitations. It is worth to mention that the duration of the experiment is kept 20 s for hydrophobic and hydrophilic surfaces. The relative optical transmittance is evaluated adopting the transmittance ratio, i.e. the relative transmittance (T rel ) is: , here T film is the optical transmittance of the film, T cleaned is optical transmittance of the film after the dust is repelled. Table 7 gives the transmittance ratio at 64 Hz of sonic excitation for hydrophobic and hydrophilic film surfaces. The relative transmittance is almost in the order of 57% for hydrophilic surface and it becomes higher (87%) for the hydrophobic surface. The relative transmittance is higher for hydrophobic surface than hydrophilic surface after dust mitigation by sonic excitation. It is worth to mention that the optical transmittance of hydrophobic sample is about 90% of clean (none-dust settled) sample for the wavelength range 600 nm 800 nm 21 . Hence, the use of hydrophobic surface improves the overall optical transmittance of the samples after dust mitigation via sonic excitation. www.nature.com/scientificreports/ However, a care must be taken to use the hydrophobic coating on the optically transmitting surfaces because of almost 10% reduction in the optical transmittance after hydrophibizing the surface occurs 21 . Hence, hydrophobizing of surfaces lowers the optical transmittances in none-dusty regions; however, it improves overall optical transmittance as the samples surfaces being subjected to regular heavy dust depositions. Moreover, the attainment of a low ratio of optical transmittance is because of the mode of film vibration at 64 Hz, which generates a singlemode and repelled dust particles in the outer region of the film moves towards the film center with multiple excitations during 20 s. Nevertheless, as the film surface tilted more than 30° about the vertical axis, the transmittance ratio is expected to increase significantly.

Conclusion
Dust repelling from transparent polyvinyl chloride film under sonic excitation is examined in relation to dust mitigation from transparent surfaces. An experimental rig is designed and built evaluating the dynamics of the repelling dust particles from the film surface. Environmental dust particles adhesion on hydrophilic and hydrophobic surfaces are evaluated through atomic force microscopy. The film surface is hydrophobized via dip-coating by the functionalized nano-silica particles. High-speed monitoring system and the tracker program are used to measure the vertical and horizontal heights of the infight dust particles for various frequencies of the sonic excitation. The flexural characteristics of the film due to applied sonic power are formulated and the mode of film vibration is determined for different sonic power frequencies. The trajectory of dust particles emanating from the film surface is formulated in a spherical coordinate system and the resulting equations are solved using COMSOL Differential Equation solver. The findings of the dynamic characteristics of the inflight dust particles Figure 11. Vertical displacement (height) of inflight dust particle (repelled) predicted from numerical simulations and obtained from experiment. The dust particle initial location on the film surface is (200, 5.74°, 5°).
Scientific Reports | (2020) 10:19348 | https://doi.org/10.1038/s41598-020-76418-2 www.nature.com/scientificreports/ are compared with those of the experimental counterparts. It is demonstrated that the inflight dust particles heights along the vertical and horizontal directions agree with those of the experimental findings. The film displacement and transverse velocity become the maximum for the excitation frequency of 64 Hz with the standoff distance between the film and the sonic excitation source of 30 mm. In this case, the film vibration results in a single principle vibration mode (0,1), which occurs at the film center. Hydrophobizing the film surface lowers the dust particle adhesion considerably; hence, adhesion force reduces by almost 35% for 0.9 µm particle and about 47% for 8 µm particle. The vertical height of the dust particle is higher for excitation frequency of 64 and reducing the excitation frequency to 45 Hz or increasing to 104 Hz, vertical height reduces significantly. The inflight particle displacement is large in both vertical and horizontal axes, which implies that the particles from the film surface can be repelled away from the film. Hence, sonic excitation can mitigate the dust particles from the transparent polyvinyl chloride surface under a proper selection of the standoff distance and the frequency of the excitation. The hydrophobic surface provides longer repelling distances for dust particles under low (55 Hz) Figure 12. Horizontal displacement (height) of inflight dust particle (repelled) predicted from numerical simulations and obtained from experiment. The dust particle initial location on the film surface is (200, 5.74°, 5°).
Scientific Reports | (2020) 10:19348 | https://doi.org/10.1038/s41598-020-76418-2 www.nature.com/scientificreports/ www.nature.com/scientificreports/ and high (75 Hz) frequency sonic excitations. This is because of low inertia force generated on the dust particle, which is comparable to the force of adhesion of the particle. However, as the sonic excitation frequency is set at 64 Hz, the inertia force generated on the dust particle becomes considerably larger than the adhesion force. Hence, the height of the inflight dust particle does not influence notably with the adhesion force. Dust mitigation