Which Is the Motion State of a Droplet on an Inclined Hydrophilic Rough Surface in Gravity: Pinned or Sliding?

The motion state of a droplet on an inclined, hydrophilic rough surface in gravity, pinned or sliding, is governed by the balance between the driving and the pinned forces. It can be judged by the droplet’s shape on the inclined hydrophilic rough surface and the droplet’s contact angle hysteresis. In this paper, we used the minimum energy theory, the minimum energy dissipation theory, and the nonlinear numerical optimization algorithm to establish Models 1–3 to calculate out the advancing/receding contact angles (θa/θr), the initial front/rear contact angles (θ1−0/θ2−0) and the dynamic front/rear contact angles (θ1−*/θ2−*) for a droplet on a rough surface. Also, we predicted the motion state of the droplet on an inclined hydrophilic rough surface in gravity by comparing θ1−0(θ2−0) and θ1−*(θ2−*) with θa(θr). Experiments were done to verify the predictions. They showed that the predictions were in good agreement with the experimental results. These models are promising as novel design approaches of hydrophilic functional rough surfaces, which are frequently applied to manipulate droplets in microfluidic chips.


Introduction
Droplets on surfaces are a phenomenon observed in everyday life, as well as in many environmental or industrial applications: coating processes [1][2][3], combustion processes [4][5][6], printing [7,8], self-cleaning surfaces [9][10][11], self-catchment surfaces [12][13][14], protein adsorption chips, etc. [15][16][17]. The study of the motion states of droplets on an inclined, hydrophilic rough surface in gravity is a fundamental problem in the mechanics of wetting and spreading [18][19][20][21][22][23][24][25], which facilitates a better understanding of how to manipulate a droplet on a rough surface. Obviously and simply, a small droplet on an inclined hydrophilic rough surface has two main motion states: pinned and inchworm sliding. However, predicting the motion state of a droplet on an inclined, hydrophilic rough surface is more complicated and difficult, because it concerns the surface inclination, the droplet shape, the droplet's contact angle hysteresis, and the dynamic behavior of the droplet's three-phase contact line.
There have been many studies to predict the pinned or sliding state of a Newtonian fluid (water, glycerol etc.) droplet on an inclined, smooth surface. The earliest theoretical work was attributed to Furmidge et al. [26]. In their theory, when a droplet stays on an inclined, smooth surface, the component of the gravitational force along the inclined surface F g can be expressed by where mg is the gravitational force. For a Newtonian fluid droplet, the pinned force F p is equal to the capillary force, and can be expressed by F p = γ lv w × [cos θ r − cos θ a ] (2) where γ lv the liquid-vapor surface tension, w the width of the drop perpendicular to the motion, θ a is the advancing contact angle, and θ r is the receding contact angle, respectively. If F g ≥ F p the droplet slips, whereas if F g < F p the droplet is pinned. Subsequently, Hashimoto et al. [27] studied the motion states of a droplet on an inclined, rough surface. They used almost the same method as Furmidge's method to predict the droplet's motion state, only replacing θ a and θ r on a smooth surface with that on a rough surface. In the Furmidge and Hashimoto methods, w is always replaced with the contact circle diameter of a droplet on a horizontal plane. Their predictions did not consider the shape change of a droplet's shape when the droplet stays on an inclined surface; therefore, their prediction results contain many errors. After these works, Masao Doi et al. [28][29][30] used the minimum energy dissipation principle to analyze the evolution of the droplet's shape when a droplet begins to stay on an inclined surface. They derived the droplet's motion state by solving a series of equations of contact line evolution. However, for a rough surface the equations of the droplet's three-phase contact line evolution are very complicated and not easily solved. Then, Legendre et al. [31][32][33] developed a volume of fluid (VOF) method and corresponding JADIM software solver, which could numerically simulate the changes of the droplet shapes and the dynamic front/rear contact angles (θ 1− * /θ 2− * ), when a droplet stays on an inclined, smooth surface. By comparing θ 1− * /θ 2− * with θ a /θ r , they judged the motion state of a droplet on an inclined surface in gravity. However, for a rough surface, the VOF method needs to finely mesh the bottom of a droplet, due to micro or nano structures on the hydrophilic rough surface. Maybe it is complicated and trivial. Frechette et al. [34] did many experiments of droplets on the inclined, smooth surfaces and showed the relationship between the droplets' motion states and the changes of contact angles. However, they did not give the theoretical models.
Our research is also for Newtonian fluid droplets. We used the minimum energy theory and the minimum energy dissipation theory to analyze the motion states of a droplet on an inclined, hydrophilic rough surface in gravity, and gave the corresponding prediction method. We did the following steps. First, we set up Model 1 to calculate θ a and θ r . Second, we set up Model 2 to calculate out the initial droplet profile Ω 0 , the initial droplet front contact angle θ 1−0 , and the initial droplet rear contact angle θ 2−0 , when a droplet begins to stay on an inclined hydrophilic rough surface. Third, we set up Model 3 to calculate out the dynamic droplet profile Ω * (* represents every position during droplet motion), the dynamic front contact angle θ 1− * , and the dynamic rear contact angle θ 2− * when the droplet stretches or contracts its three-phase contact line on the inclined, rough surface. Fourth, we gave out the prediction for the motion states of a droplet on an inclined, hydrophilic rough surface. Finally, we did many experiments to verify the predictions, and found that the predictions are in good agreement with the experimental results.
We gave a simple description for the methods of Models 1-3 and the prediction, which will be described in detail in the following sections. Model 1 was based on the minimum energy theory. As is shown in Figure 1a, a droplet was imaged to stay on a flat hydrophilic rough surface. When we imaged to continuously add the volume of the droplet, but fixed the three-phase contact line, we could gain the potential energy ∆E. When ∆E is equal to the energy barrier E barr , preventing the contact line from moving, the contact angle is θ a . In contrast, we could calculate θ r by imaging to continuously decrease the volume. Model 2 was based on the minimum energy theory. As is shown in Figure 1b, when a droplet was initially plated on an inclined surface, the sum of gravitational energy and the interface free energy is minimal. We gave out the integral expression of the sum energy, numerically dispersed the droplet, minimized, and got the initial shape Ω 0 , θ 1−0 /θ 2−0 . Model 3 was based on the minimum energy dissipation theory. As is shown in Figure 1c, when a droplet moves on an inclined surface, the droplet has the local minimum sum energy on each point of droplet inchworm motion. We dispersed the droplet, minimized the sum energy with the constraint of contact line length, and calculated out the dynamic shape Ω * and θ 1− * /θ 2− * . As is shown in Figure 1d,e, the prediction was based on comparing θ 1−0 /θ 2−0 and θ 1− * /θ 2− * with θ a /θ r . We thought that the front contact line of droplet moves if θ 1−0 (θ 1− * ) ≥ θ a and the rear contact line of droplet moves if θ 2−0 (θ 2− * ) ≤ θ r ; otherwise, they are pinned. Furthermore, the droplet can keep sliding if both the front and the rear contact lines move. Otherwise, the droplet will be pinned finally when both the front and the rear contact lines stop. barr tact angle is a  . In contrast, we could calculate r  by imaging to continuously decrease the volume. Model 2 was based on the minimum energy theory. As is shown in Figure 1b, when a droplet was initially plated on an inclined surface, the sum of gravitational energy and the interface free energy is minimal. We gave out the integral expression of the sum energy, numerically dispersed the droplet, minimized, and got the initial shape 0  , 10   / 20   . Model 3 was based on the minimum energy dissipation theory. As is shown in Figure 1c, when a droplet moves on an inclined surface, the droplet has the local minimum sum energy on each point of droplet inchworm motion. We dispersed the droplet, minimized the sum energy with the constraint of contact line length, and calculated out the dynamic shape *  and 1*   / 2*   . As is shown in Figure 1d     In the initial droplet state, θ 1−0 > θ a and θ 2−0 > θ r , we predicted the front contact line moving and the rear contact line pinned. (e) In one point of droplet motion, θ 1− * < θ a and θ 2− * > θ r , both the front and the rear contact lines are pinned, the droplet will finally be pinned, and the motion state is predicted as "stretching-to-pinned". We gave a simple example to predict the droplet motion state of "stretching-to-pinned" on an inclined rough SiO 2 surface, which will be described in detail in Section 4.3.1. A SiO 2 surface was patterned by circular microstructures (d = 6 µm, h = 12 µm, and a = 60 µm). The droplet had the volume of 40 µL, and the surface was 39 • inclined to the horizontal plane. As shown in Table 1, using Models 1-3 we got θ a = 75.61 • , θ r = 42.91 • , θ 1−0 = 77.71 • , θ 2−0 = 46.43 • , θ 1− * = 75.48 • , and θ 2− * = 44.05 • . Because θ 1−0 > θ a and θ 2−0 > θ r , the rear end of the droplet is pinned and the front end advances, initially leading to drop stretching. Stretching increases the three-phase contact line length and decreases θ 1− * and θ 2− * . In one point of motion, θ 1− * = 75.48 • < θ a and θ 2− * = 44.05 • > θ r , which lead to the droplet being pinned. The droplet motion state was regarded as "stretching-to-pinned".
In this work, the buoyancy force can be ignored due to the low air density; similarly, the fluid drag force can be ignored due to the near-zero slip velocity. In other practical scenarios, such as slurry Taylor droplets on inclined surfaces, both the buoyancy force and fluid drag forces should be considered [35,36]. Table 1. Numerical results and motion state predictions for the droplet on an inclined, rough SiO 2 surface. A droplet on the homogeneous hydrophilic rough surface is always in Wenzel state. We only calculated θ a and θ r of a droplet on a rough surface in Wenzel state. For an equilibrium droplet on a horizontal rough surface, when the droplet volume V1 decreases or increases to the interval V1 r < V1 < V1 a , but the three-phase contact line keeps immobile, the apparent contact angles (ACAs) corresponding to the critical V1 r and V1 a are the receding contact angle θ r and the advancing contact angle θ a , respectively.

Surface Tilt
As is shown in Figure 2, the droplet on a rough horizontal surface and the relative energy of the system E w−1 can be expressed by where r gh = 1 + πdh (a+d) 2 is the roughness factors, r b−1 the radius of contact circle on the surface, ρ the density of the liquid, g the gravitational acceleration, ϕ is the angle between the radius vector and the positive z-axis, r(ϕ) is the length of radius vector, γ lv is the interface tension coefficient of the liquid vapor, and θ e is the equilibrium contact angle on a smooth flat surface (for details, see Sections S3.1 and S3.2 in Supplementary Materials). As is shown in Figure 3a, we increased the droplet volume V  but kept the contact circle immobile   b 1 r n r ; the droplet has the local minimum relative energy v+Δv fix E   .
Also, we increased the droplet volume V  and let the contact circle be mobile to get the whole minimum relative energy v+Δv E . When the droplet keeps the contact circle immobile, the droplet has the planar potential energy v+Δv fix v+Δv The energy barrier Figure 2. The sketch of calculation for the relative energy of the system E wmin−1 . (a) A droplet stays on a horizontal rough surface decorated by circular pillars; h is the height of the pillar, d is the diameter of the pillar, and a is the periodic spacing of pillars. (b) A profile of the half-droplet; h a is the distance between the top of the droplet and the surface, h amax is the maximum h a , r b−1 is the radius of contact circle on the surface, r bmax is the maximum r b−1 , and ϕ 0 ≤ ϕ ≤ π 2 is the angle between the radius vector r(ϕ) and the positive z-axis.
The droplet staying on a rough surface reaches the minimum relative energy. Using the finite-difference method and the nonlinear optimization algorithm, we calculated the equilibrium of relative total energy E wmin−1 and the equilibrium contact circle radius r b (for details, see Sections S3. 3

and S3.4 in Supplementary Materials).
As is shown in Figure 3a, we increased the droplet volume ∆V but kept the contact circle immobile r[n + 1] = r b ; the droplet has the local minimum relative energy E v+∆v−fix . Also, we increased the droplet volume ∆V and let the contact circle be mobile to get the whole minimum relative energy E v+∆v . When the droplet keeps the contact circle immobile, the droplet has the planar potential energy ∆E = E v+∆v−fix − E v+∆v . The energy barrier E barr , which prevents the droplet's three-phase contact line from moving, can be calculated by [37][38][39] where U is the adhesive friction between liquid and solid, and l act the actual length of the contact line. The moment ∆V increases to make ∆E = E barr , the three-phase contact circle begins to advance, and the corresponding contact angle is θ a . On the other hand, as is shown in Figure 3b, the moment ∆V decreases to make ∆E = E barr , the three-phase contact circle begins to recede, and the corresponding contact angle is θ r . (for details, see Sections S3.5-S3.8 in Supplementary Materials).    As is shown in Figure 4, the droplet initially stays on the hydrophilic rough surface inclined to the horizontal plane with a tilt angle φ. We selected the rough surface as the XOY plane. One point on the symmetry axis of the droplet base was defined as the origin O, the direction going ascent along the inclined surface was defined as the positive Y direction, the direction vertical to the surface and going to droplet curvature was defined as the positive Z direction, and the inside direction vertical to the YOZ plane was defined as the positive X direction. Also, we defined the vector from O to one point on the surface of the droplet as → r (β, α), where the azimuth angle β(−π ≤ β ≤ π) was the angle from the positive X axis to the projection of → r (β, α) on the XOY plane, and the zenith angle α 0 ≤ α ≤ π 2 was the angle from the positive Z axis to As is shown in Figure 3a, we increased the droplet volume V    The relative total energy E w−2 of the system on the inclined rough surface can be expressed as the following (for the deduction, see Sections S5.1 and S5.2 in Supplementary Materials): where r(β, π 2 ) is the length variable of the radius vector, with the zenith angle → r (β, π 2 ). The droplet forms its initial shape at the minimum relative energy. Using the finitedifference method and the nonlinear optimization algorithm, we simulated the initial droplet profile Ω 0 , the initial front contact angle θ 1−0 , the initial rear contact angle θ 2−0 , and the initial droplet contact line length l 0 for a droplet on the inclined hydrophilic rough surface (for details, see Sections S5.3 and S5.4 in Supplementary Materials). The initial state of the droplet is the whole minimum relative energy state, but it is certainly unstable if θ 1−0 ≥ θ a or θ 2−0 ≤ θ r . As is shown in Figure 5a, when θ 1−0 ≥ θ a and θ 2−0 > θ r , the front end of the droplet advances while the rear end stays pinned, leading to drop stretching. Alternatively, as shown in Figure 5b, when θ 1−0 < θ a and θ 2−0 ≤ θ r , the rear end of the droplet retracts while the front end stays pinned, leading to drop contracting. When the droplet's contact line moves, the front contact angles and the rear contact angles will change with the change of the droplet contact line length. Every contact line length corresponds to the droplet's dynamic front contact angle θ 1− * and the droplet's dynamic rear contact angle θ 2− * . According to minimum energy dissipation principle, the droplet has local minimum total potential energy at every contact line length, whether stretching or contracting. In Model 3, we first set the length as l given . Second, with the constraint l given , we minimized the local total potential energy E w−2 and calculated the dynamic droplet shape Ω * , θ 1− * , and θ 2− * corresponding to l given . Using the algorithm for continuously changing the set l given , we can calculate θ 1− * and θ 2− * at every point during the droplet moving period. (for the algorithm and the flow chart, see Sections S6. 1

  When a Droplet Evolves Its Contact Line Length on an Inclined, Hydrophilic Rough Surface
The initial state of the droplet is the whole minimum relative energy state, but it is certainly unstable if

The Prediction Method of the Droplet: Pinned or Sliding
When the droplet stays on an inclined, hydrophilic rough surface, we can predict the front end of the three-phase contact line, moving or not, by the comparison between  . According to the motion state of the front and the rear end of the three-phase contact line, we finally predict the droplet being pinned or sliding. The prediction procedure is divided into five steps, as follows (for the flow chart, see Section S7 in Supplementary Materials):


Step 1: We used Model 1 to calculate the advancing angle a  and the receding angle

The Prediction Method of the Droplet: Pinned or Sliding
When the droplet stays on an inclined, hydrophilic rough surface, we can predict the front end of the three-phase contact line, moving or not, by the comparison between θ 1−0 /θ 1− * and θ a , and predict the rear end of the three-phase contact line moving or not by the comparison between θ 2−0 /θ 2− * and θ r . According to the motion state of the front and the rear end of the three-phase contact line, we finally predict the droplet being pinned or sliding. The prediction procedure is divided into five steps, as follows (for the flow chart, see Section S7 in Supplementary Materials):

SiO 2 Rough Surface Fabrication
The fabrication process for SiO 2 rough surfaces started from 4 inch, n-type (100) silicon wafers. Firstly, the AZ4620 photoresist was patterned for the diameter d and the periodic spacing a of the periodic circular microstructures. Secondly, deep reactive ion etching (DRIE) was processed for the height h of circular microstructures. Thirdly, the wet thermal oxidation was processed to grow 500 nm thick silicon dioxide covered on the surfaces.

SiO 2 Rough Surface Measurement
The morphologies of SiO 2 rough surfaces were measured by field emission scanning electron microscopy (FE-SEM, S4700, Hitachi, Japan). As is shown in Figure 6, the parameters of microstructures (d,a and h) were defined in the SEM images. The fabricated SiO 2 rough surfaces were decorated by the microstructures with parameters d = 6 µm, h = 12 µm, and a = 60 µm.

Characterization of Droplet Equilibrium Contact Angles
We used an SDC-80 (Sindin, China) profile and contact angle measurement to characterize the droplet equilibrium contact angles on smooth surfaces. Droplets were set on the horizontal table. We took pictures for the equilibrium droplets and gained the equi-

Characterization of Droplet Equilibrium Contact Angles
We used an SDC-80 (Sindin, China) profile and contact angle measurement to characterize the droplet equilibrium contact angles on smooth surfaces. Droplets were set on the horizontal table. We took pictures for the equilibrium droplets and gained the equilibrium contact angles. As is shown in Figure 7, the equilibrium contact angle on the smooth PMMA surface was 74.73 ± 0.74 • , and the angle for the smooth SiO 2 surface was 65.57 ± 1.27 • .

Characterization of Droplet Equilibrium Contact Angles
We used an SDC-80 (Sindin, China) profile and contact angle measurement to characterize the droplet equilibrium contact angles on smooth surfaces. Droplets were set on the horizontal table. We took pictures for the equilibrium droplets and gained the equilibrium contact angles. As is shown in Figure 7, the equilibrium contact angle on the smooth PMMA surface was 74.73 ± 0.74°, and the angle for the smooth SiO2 surface was 65.57 ± 1.27°.

Characterization of Droplet Motion on Inclined Surfaces
As is shown in Figure 8, also on a SDC-80, a hydrophilic rough surface was laid on an inclined table. The droplets were emitted from a needle to the rough surface. The droplet shapes and motion states were recorded by a side video camera, whose optical axis was perpendicular to the trajectory of the droplets. By different droplets' volumes, different tilt angles, and different hydrophilic rough surfaces, we got different experimental data about droplets' profiles ( 0  , *  ), droplets' front contact angles (

Characterization of Droplet Motion on Inclined Surfaces
As is shown in Figure 8, also on a SDC-80, a hydrophilic rough surface was laid on an inclined table. The droplets were emitted from a needle to the rough surface. The droplet shapes and motion states were recorded by a side video camera, whose optical axis was perpendicular to the trajectory of the droplets. By different droplets' volumes, different tilt angles, and different hydrophilic rough surfaces, we got different experimental data about droplets' profiles (Ω 0 ,Ω * ), droplets' front contact angles (θ 1−0 ,θ 1− * ), and droplets' rear contact angles (θ 2−0 ,θ 2− * ).

The Advancing and Receding Contact Angles(θ a /θ r ) Change with the Droplet Volume
We introduced the droplet bond number (for water, Bo = (2.73 × 10 −3 ) 2 ( 3V 4π ) 2 3 ) to study the relationship between the droplet volume and the advancing/receding angle (θ a /θ r ). The bigger the droplet Bond number is, the larger the droplet volume is. As shown in Figure 9, when the bond number increases, θ a goes down and θ r goes up. The so-called contact angle hysteresis is expressed by CAH = θ a − θ r . The larger the volume of droplet is, the smaller CAH is, and the easier it is for the droplet to move.

The Advancing and Receding Contact Angles( a  / r  ) Change with the Droplet Volume
We introduced the droplet bond number (for water, ) to study the relationship between the droplet volume and the advancing/receding angle ( a  / r  ). The bigger the droplet Bond number is, the larger the droplet volume is. As shown in Figure 9, when the bond number increases, a  goes down and r  goes up. The socalled contact angle hysteresis is expressed by ar CAH  . The larger the volume of droplet is, the smaller CAH is, and the easier it is for the droplet to move.

Droplets on an inclined Rough SiO2 Surface
As is shown in Table 2 and Figure 11, when the tilt angles of the rough SiO2 hydro-

Droplets on an inclined Rough SiO 2 Surface
As is shown in Table 2 and Figure 11, when the tilt angles of the rough SiO 2 hydrophilic surface change from 12 • to 52 • , the final motion states of droplets on the inclined surface change from "pinned" to "sliding". The results from Model 1 show that the θ a and θ r of the SiO 2 rough surface, decorated by microstructures (d = 6 µm, h = 12 µm, and a = 60 µm), are 75.61 • and 42.91 • , respectively. Table 2. Numerical results and motion state predictions of droplets on the inclined rough SiO 2 surface a . For the surface with tilt angle 12 • , results from Model 2 and the prediction were θ 1−0 = 64.85 • < θ a and θ 2−0 = 58.74 • > θ r , respectively. We predicted that the motion state of the droplet would be "pinned" on the 12 • inclined SiO 2 surface.
For the surface with tilt angle 52 • , results from Model 2 and the prediction in the droplet's initial state were θ 1−0 = 81.21 • > θ a and θ 2−0 = 40.19 • < θ r , respectively. We predicted the motion state of the droplet is "sliding" on the 52 • inclined SiO 2 surface.
We did experiments to verify predictions. Experimental results showed that when the surfaces tilt angles were in the range 40 • < φ < 42 • , the motion state of the droplet is "contracting-to-pinned", which has 1~2 • errors with the prediction results 39 • ≤ φ < 44 • . It was showed that the experimental data agreed well with the prediction results (for video, see Video S1-S4 in the Supplementary Materials). For the surface with tilt angle 12°, results from Model 2 and the prediction were 10   = 64.85° < a  and 20   = 58.74° > r  , respectively. We predicted that the motion state of the droplet would be "pinned" on the 12° inclined SiO2 surface.
For the surface with tilt angle 39°, the results from Models 2 and 3 and the prediction, in the droplet's initial state, were that

Droplets on an Inclined, Smooth PMMA Surface
As is shown in Table 3 and Figure 12, when the tilt angles of the smooth PMMA hydrophilic surface change from 10 • to 28 • , the final motion states of droplets on the inclined surface change from "pinned" to "sliding". Results from Model 1 indicate that the θ a and θ r of the smooth PMMA surface are 81.62 • and 65.10 • , respectively.  We did experiments to verify our predictions. Experimental results showed that when the surface tilt angles are in the range 19° <  < 23° , the motion state of droplet is "contracting-to-pinned"; this has 4~5° errors, with prediction results at 14° <  < 19°. It was shown that the experimental data agreed well with the prediction results (for video, For the surface with tilt angle 10 • , the results from Model 2 and the prediction were that θ 1−0 = 76.20 • < θ a and θ 2−0 = 67.56 • > θ r , respectively. We predicted that the motion state of droplet is "pinned" on the 10 • inclined PMMA surface. For the surface with tilt angle 14 • , the results from Models 2 and 3 and the prediction, in the droplet's initial state, were θ 1−0 = 77.48 • < θ a and θ 2−0 = 63.89 • < θ r , leading to drop contracting, respectively. At one point (*), θ 1− * = 80.71 • < θ a and θ 2− * = 65.13 • ≈θ r lead to drop pinning. We predicted the motion state of droplet would be "contracting-to-pinned" on the 14 • inclined PMMA surface.
For the surface with a tilt angle of 28 • , the results from Model 2 and the prediction, in the droplet's initial state, were that θ 1−0 = 86.85 • > θ a and θ 2−0 = 55.71 • < θ r , respectively. We predicted the motion state of droplet would be "sliding" on the 28 • inclined PMMA surface.
We did experiments to verify our predictions. Experimental results showed that when the surface tilt angles are in the range 19 • < φ < 23 • , the motion state of droplet is "contracting-to-pinned"; this has 4~5 • errors, with prediction results at 14 • < φ < 19 • . It was shown that the experimental data agreed well with the prediction results (for video, see Video S5-S8 in the Supplementary Materials).

Conclusions
All this work contributes to understanding the wetting and the spreading properties of a droplet on an inclined, hydrophilic rough surface, both theoretically and practically. In this paper, we used the minimum free energy theory, the minimum energy dissipation theory, and nonlinear optimization algorithms to model and calculate the advancing/receding contact angles (θ a /θ r ), the initial front/rear contact angles (θ 1−0 /θ 2−0 ), and the dynamic front/rear contact angles (θ 1− * /θ 2− * ) of a droplet on an inclined hydrophilic rough surface. Also, we predicted the droplet motion state by comparing θ 1−0 (θ 2−0 ) and θ 1− * , (θ 2− * ) with θ a (θ r ). Additionally, experiments were done to verify the predictions. The experimental data were found to agree with the predictions. Our method can be used to optimize the hydrophilic rough surface, which can be used to exploit devices like variable-focus lances, electronic displays, and micro-fluidic systems.