The static shape of drop under electrowetting actuation is well studied and recent electrowetting theory and experiments confirm that the local contact angle (microscopic angle) is unaffected while the apparent contact angle (macroscopic angle) is characterized by the Lippmann-Young equation. On the other hand, the evolution of the drop motion under electrowetting actuation has received less attention. In this paper, we investigate the motion of a conducting water drop on an electrowetting device (EWD) using the level set method. We derive a contact line two-phase flow model under electrowetting actuation using energy dissipation by generalizing an existing contact line model without the electric field. Our model is consistent with the static electrowetting theory as the dynamic contact angle satisfies the static Young's equation under equilibrium conditions. Our steady state results show that the apparent contact angle predicted by our model satisfies the Lippmann-Young's relation. Our numerical results based on the drop maximum deformation agree with experimental observations and static electrowetting theory. Finally, we show that for drop motion our results are not as good due to the difficulty of computing singular electric field accurately. Nonetheless, they provide useful insights and ameaningful first step towards the understanding of the drop dynamics under electrowetting actuation.

The paper proposes a general framework allowing the analysis of wetting problems in thesituation when interfacial tensions depend on external fields. An equation predictingapparent contact angles of sessile droplets deposited on rough surfaces in the presence ofexternal fields is derived. The problem of wetting is discussed in the framework of thevariational approach. Derivation of a general equation generalizing the Cassie and Wenzelapproaches is presented. The effects related to the line tension which are important fornano-structured surfaces are considered.

The electrowetting process is commonly used to handle very small amounts of liquid on a solid surface. This process can be modelled mathematically with the help of the shape optimization theory. However, solving numerically the resulting shape optimization problem is a very complex issue, even for reduced models that occur in simplified geometries. Recently, the second author obtained convincing results in the 2D axisymmetric case. In this paper, we propose and analyze a method that is suitable for the full 3D case.