1. Introduction
It is important to understand the statics and dynamics of bodies under the effect of gas–liquid interface in nature, ranging from standing and locomotion of insects on water (Bush & Hu Reference Bush and Hu2006) through water entry of bodies (Truscott, Epps & Belden Reference Truscott, Epps and Belden2014) to mineral flotation (Ueda et al. Reference Ueda, Tanaka, Uemura and Iguchi2010). The design of micromachines or microrobots on the gas–liquid interface is becoming increasingly popular in research, for example, processing platforms floating on a liquid surface (Bowden et al. Reference Bowden, Terfort, Carbeck and Whitesides1997; Barbot et al. Reference Barbot, Tan, Power, Seichepine and Yang2019; Basualdo et al. Reference Basualdo, Bolopion, Gauthier and Lambert2021; Jung et al. Reference Jung2021), water walking microrobots (Hu, Chan & Bush Reference Hu, Chan and Bush2003; Hu et al. Reference Hu, Lum, Mastrangeli and Sitti2018) and strider-like water jumping microrobots (Koh et al. Reference Koh, Yang, Jung, Jung, Son, Lee, Jablonski, Wood, Kim and Cho2015).
Floating and sinking are the two main states of a small object on the liquid surface in the absence of an initial velocity of the body, considering the effect of surface tension and gravity. Vella (Reference Vella2015) presented a review of research on floating and sinking states of a small object. A small object with larger density than water can float at the water surface under the effect of surface tension or sink in the water. Whether the state is floating or sinking is determined by some parameters, e.g. the density, size, shape and surface properties of the object, and the distance away from nearby other objects. Pierson & Magnaudet (Reference Pierson and Magnaudet2018) conducted experimental research on the inertial settling of a sphere through an interface between two immiscible liquids with moderate density difference, and investigated the floating and sinking regimes. In addition to different sinking modes, due to the absence of initial sphere velocity, only the static state of flotation of sphere was illuminated but the transient state of initially sinking and then floating has not been mentioned.
The impact of bodies onto liquid can lead to rich physical phenomena, which involve the inertial force, the surface tension force, the hydrodynamic force and the gravitational force. Research on this topic focuses on the impact of spherical particles onto liquid. For a relatively low impact speed, the particle may either sink into the liquid with the pulled cavity pinching off or rebound with the pulled cavity not pinching off (Chen et al. Reference Chen, Liu, Lu and Ding2018), depending on the properties of the particle and the liquid. For a high impact speed, a region of high pressure in liquid can appear and form a splash in the early stage of the impact (Thoroddsen et al. Reference Thoroddsen, Etoh, Takehara and Takano2004; Marston, Li & Thoroddsen Reference Marston, Li and Thoroddsen2012). Wagner's theory (Wagner Reference Wagner1932) describes this important early stage, and is often applied in many studies concerning the high-speed impact between two phases (solid–liquid or liquid–liquid), e.g. Cointe & Armand (Reference Cointe and Armand1987), Howison, Ockendon & Wilson (Reference Howison, Ockendon and Wilson1991) and Cimpeanu & Moore (Reference Cimpeanu and Moore2018).
Recently, Galeano-Rios et al.'s (Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021) experiments on impact of a superhydrophobic sphere onto water for three densities of particle (1.2, 2.2 and 3.2 g cm−3), three radii of particle (0.83, 1.24 and 1.64 mm) and different impact velocities (30–110 cm s−1) found a new dynamical phenomenon that the pulled cavity pinches off and eventually the particle rebounds. This phenomenon was called resurrection, which occurs in a narrow range of impact velocities and for a sufficiently small density difference between the particle and the water. The observed resurrection was reproduced by direct numerical simulations (DNS) conducted for the particle density of 1.2 g cm−3. The observation of the interesting resurrection phenomenon may be crucial to the design of strider-like water jumping microrobots (Koh et al. Reference Koh, Yang, Jung, Jung, Son, Lee, Jablonski, Wood, Kim and Cho2015). Up to date, the ranges of values in main parameters for the occurrence of resurrection have not been determined due to the limited cases of experiments and more limited cases of DNS.
A strider has long cylinder-like legs, and the cylinder-like shape can increase the surface tension force in order to bear the weight of the strider and the impact especially as the strider jumps. Vella, Lee & Kim (Reference Vella, Lee and Kim2006) performed an experiment to investigate the descending process of a horizontal cylinder on the air–water interface, and developed a simple hydrodynamic model based on the experimental data. Using the developed model, they predicted that the sinking time taken for the cylinder to reach the fully immersed state due to gravity is directly proportional to $\sqrt {{{l}_{c\textrm{a}}}/g} $ where lca is the capillary length and g is gravitational acceleration. Vella & Metcalfe (Reference Vella and Metcalfe2007) developed a nonlinear two-dimensional (2-D) mathematical model of equilibrium of a floating infinite cylinder and studied the criterion for the cylinder sinking. According to Bhatnagar & Finn (Reference Bhatnagar and Finn2006) and Chen & Siegel (Reference Chen and Siegel2018), a floating cylinder under the surface tension effect may reach equilibrium states, but some of these states are unstable. Janssens, Chaurasia & Fried (Reference Janssens, Chaurasia and Fried2017) studied a partly submerged cylinder under the effect of the surface tension imbalance and discussed the Marangoni propulsion for the floating cylinder, which possibly inspires a new way to manipulate long cylinders at the liquid–gas interface. By comparing the stabilities of confined and unconfined floating cylinders, Zhang & Zhou (Reference Zhang and Zhou2023) found the confinement by two hydrophobic plates with a small spacing can assist with the interfacial floatation of a cylinder with a larger weight.
Different from the floatation behaviour of a 2-D body, the impacting dynamical behaviour of a 2-D body or a long cylinder onto a liquid bath is much richer, dependent on the impacting velocity and the particle-to-liquid density ratio. Some experimental research about a long horizontal cylinder falling or impacting onto water was conducted. Greenhow & Lin (Reference Greenhow and Lin1983) and Greenhow (Reference Greenhow1988) experimentally observed the evolution of air–water interface as a horizontal long circular cylinder descends through the interface. Colicchio et al. (Reference Colicchio, Greco, Miozzi and Lugni2009) conducted experimental and numerical investigations of the water-entry and water-exit of a light (lighter than water) circular cylinder either freely falling at a height above the interface or exiting the water at a depth below the interface. Goharzadeh & Molki (Reference Goharzadeh and Molki2012) designed a release unit of horizontally releasing a horizontal cylinder to make the horizontal cylinder vertically descend through the interface and experimentally observed the impact dynamics of a cylinder and the surface wave. Wei & Hu (Reference Wei and Hu2014) experimentally investigated the three-dimensional (3-D) effect on water entry with an impacting horizontal circular cylinder. Lyu et al. (Reference Lyu, Wei, Tang, New and Li2015) experimentally studied the water entry of a horizontal circular cylinder into flowing water with an impact velocity. The horizontal displacement was found to increase as the water velocity increases, and the trajectories of the hollow cylinder collapsed at the initial stage of submerging in water when the ratio of impact velocity to water velocity was the same as the solid cylinder. When the impact speed is relatively high, Wagner's (Reference Wagner1932) theoretical framework can be used to study the impacting dynamical behaviour. For the liquid entry of an arbitrary 2-D body, Zhao & Faltinsen (Reference Zhao and Faltinsen1993) presented a simple asymptotic solution based on the Wagner theory, without considering the surface tension effect. Korobkin (Reference Korobkin1997) used the method of matched asymptotic expansions to analyse the high-speed liquid–solid impact problem while the role of the surface tension can be estimated by this method.
Some numerical research about a 2-D body or horizontal long cylinder falling or impacting onto water was also conducted. Zhu, Faltinsen & Hu (Reference Zhu, Faltinsen and Hu2007) numerically simulated the deformations of water surface induced by the water entry and exit of a horizontal 2-D circular cylinder with both forced and free vertical motions, and the simulated results were found to be in favourable agreement with the experimental results by Greenhow & Lin (Reference Greenhow and Lin1983). Hafsia et al. (Reference Hafsia, Mnasri, Mohamed and Maalel2009) conducted a 2-D numerical simulation of a horizontal cylinder entering into and pulled out from water and well predicted the deformation of the interface as simulated by Lin (Reference Lin2007). Liu, Gao & Ding (Reference Liu, Gao and Ding2017) conducted 2-D modelling and simulations of fluid–structure interaction involving dynamic wetting as a horizontal circular cylinder sinks through the interface. Kiara, Paredes & Yue (Reference Kiara, Paredes and Yue2017) carried out a 2-D numerical investigation of the water entry of horizontal cylinders without and with spin. Iranmanesh & Passandideh-Fard (Reference Iranmanesh and Passandideh-Fard2017) performed 3-D numerical simulations on the water entry of a horizontal circular cylinder for low Froude numbers, and the numerical results were found to be in a good agreement with the experimental data by Wei & Hu (Reference Wei and Hu2014).
Up to date, the work about the occurrence or non-occurrence of the resurrection (observed for an impacting superhydrophobic sphere; see Galeano-Rios et al. Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021) for a small cylinder impacting onto a liquid bath has not been reported. In this situation, some intriguing questions have been raised. Will the interesting resurrection occur for a small 2-D cylinder impacting onto a liquid bath? What are the ranges of the density of the impacting cylinder, the impact velocity, and the surface tension coefficient for the occurrence of the resurrection? What are the occurring mechanism and condition of the resurrection? Without doubt, these questions are crucial to the research about the impacting of small particles onto a liquid bath and the related industry.
In this paper, in order to answer the questions well, the impact of a small superhydrophobic cylinder onto a liquid bath is investigated using DNS based on the pseudo-solid method. The dynamic modes of the cylinder are studied and classified into four types. Various states of the impacting cylinder are comprehensively identified in a space of two parameters (the density ratio of the cylinder to the liquid and the Weber number). The occurring mechanism and condition of the resurrection are investigated.
2. Problem formation and mathematical methods
Figure 1 shows a schematic of a 2-D superhydrophobic circular cylinder of diameter D impacting downwards onto a liquid bath from air in a computational domain (spanning 10D in width and 20D in length) in a downward gravity field in Cartesian coordinates $(\bar{x},\bar{y})$. The initial cylinder velocity of non-zero magnitude is V (when the cylinder contacts the undisturbed liquid–gas surface), which is chosen as the reference velocity for our system. The related physical parameters used in our 2-D simulations can be found in table 1. To simplify the simulation of the fluid–solid interaction, the pseudo-solid approach is adopted (see the details in Appendix A), which has been shown as a satisfactory model for the impact of a solid sphere with a very large contact angle (Galeano-Rios et al. Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021). Within the pseudo-solid approach, a special stiff 2-D droplet with high viscosity (restricting the liquid motion inside) and high surface tension coefficient (restricting the relatively large deformation during impact) is used to simulate the solid cylinder. In our 2-D calculations, the densities of the pseudo-solid, the liquid and the gas are denoted by ρs, ρl and ρg, respectively, while the dynamic viscosities of the pseudo-solid, the liquid and the gas are denoted by $\mu_{s}$, $\mu_{l}$ and $\mu_{g}$, respectively. The surface tension coefficients of the pseudo-solid surface and the liquid surface are denoted by σs and σl, respectively. The dynamic viscosity of the pseudo-solid (i.e. the cylinder) is 250 times as large as that of the liquid in the bath (see $\mu_{s}/\mu_{l}$ in table 1), while the surface tension coefficient of the pseudo-solid surface is 150 times as large as that between the liquid and the gas (see σs/σl in table 1).
Within the pseudo-solid method, it is unable to impose a value for the contact angle of the pseudo-solid. Instead, the contact angle of the pseudo-solid remains 180° naturally, due to the specific relation among the surface tensions of the solid–liquid, solid–gas and liquid–gas interfaces (see the details in Appendix A). Applying the pseudo-solid method, we find that there is a very thin film of gas with the thickness of approximately three meshes (less than D/100) between the pseudo-solid and the liquid and there is no triple contact line (CL; contact point in two dimensions) on the pseudo-solid. However, the liquid-entry behaviour with the superhydrophobic surface appears to be little influenced by the microscopic details of the contact with this surface, according to Galeano-Rios et al. (Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021). We have explained this from the perspective of energy in Appendix B.
The physical parameters are non-dimensionalised as follows:
where u and v are the dimensionless horizontal and vertical velocities of the fluid, respectively, $v_{s}$ is the dimensionless vertical velocity of the cylinder mass centre, and τ is the dimensionless time. The governing Navier–Stokes equations in dimensionless form are given by
where the definitions of the dimensionless numbers Re, We and Fr are listed in table 1, the Nabla operator $\boldsymbol{\nabla }\textrm{ = }\textrm{(}\partial /\partial x\textrm{,}\partial /\partial y\textrm{)}$, the velocity vector u = (u, v), κ denotes the curvature of the liquid–gas interface (the liquid surface) or the pseudo-solid surface, n denotes the unit normal of the interface and δ denotes the Dirac delta function at the interface. The Bond number Bo (with Bo = ρgD 2/σ) is commonly employed for interfacial floatation (Vella Reference Vella2015). To employ the Bond number in this paper is also feasible with the relationship Bo = We/Fr 2. To simplify our analysis, the Froude number Fr is set as a moderate constant (table 1), under which the highlighted resurrection is likely to occur. By changing We, the relative importance of the gravity and the surface tension (measured by Bo) varies. Thus, for the liquid-entry problem in this paper (with Fr fixed), there is no essential difference between using Bo and using We.
To mark the interface location, the volume fraction functions cl(x, y, t) and cs(x, y, t) are used, with which the value of cl, cs and 1 − cl − cs denote the volume fraction of the liquid, the pseudo-solid and the gas, respectively, in a computational cell. The domains of cl = 1, cl = 0 and 0 < cl < 1 represent the liquid, the gas and the liquid surface, respectively, while the domains of cs = 1, cs = 0 and 0 < cs < 1 represent the pseudo-solid, the liquid or gas and the pseudo solid surface, respectively. The volume fractions cl and cs are governed by the advection equations:
For the whole computational domain, the density ρ (dynamic viscosity μ) can be interpolated with the pseudo-solid density ρs (dynamic viscosity $\mu_{s}$), the liquid density ρl (dynamic viscosity $\mu_{l}$) and the gas density ρg (dynamic viscosity $\mu_{g}$):
To reduce computational cost, the symmetry boundary condition is set at the left side of the computational domain. Only the right half of the physical domain will be computed (see figure 1). In addition, the outflow boundary condition is employed at the top of computational domain as
while no-slip boundary conditions are employed at the right side and the bottom of the computational domain as
The DNS are implemented using the well-known, open-source, volume-of-fluid (VoF) package Gerris (Popinet Reference Popinet2003, Reference Popinet2009). The adaptive mesh refinement is adopted and the local mesh sizes are adjusted depending on the interfacial location, the velocity gradient and the vorticity. An independence analysis is conducted to determine the appropriate domain size and grid refinement (see Appendix C).
3. Results and discussion
In order to show the validity of our DNS results, we simulate the impact of a superhydrophobic sphere onto a liquid bath using the 2-D axisymmetric model and compare the calculated results at different times with the experimental data by Galeano-Rios et al. (Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021), as shown in figure 2. The parameters of liquid and solid in Galeano-Rios et al. (Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021) are employed in these simulations while the parameters of gas in theirs are not given. We employ ${\rho _g} = 10\ \textrm{kg}\ {\textrm{m}^{ - 3}}$ (with ${\rho _g}/{\rho _l} = {10^{ - 2}}$) and ${\mu _g} = {10^{ - 5}}\ \textrm{Pa}\ \textrm{s}$ for the gas in these simulations. The gas–liquid density ratio of ${\rho _g}/{\rho _l} = {10^{ - 2}}$ is larger than that in reality (approximately 10−3) while the influence of gas is still negligible. In addition, we find the density ratio of ${\rho _g}/{\rho _l} = {10^{ - 2}}$ (also in table 1) can possibly tune the Gerris solver for better convergence. The numerical results are in a good agreement with the published experimental data. Notably, the equilibrium contact angle of the sphere is 160° for the experiments in Galeano-Rios et al. (Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021) while the contact angle in our simulations is naturally 180°, which may account for the error of our simulations in figure 2(a). In addition, the deformation of the pseudo-solid used to simulate the sphere is generally less than 5 % of the original sphere diameter, as shown in figure 2(b). This proves validity and applicability of the pseudo-solid approach used here, as it was in Galeano-Rios et al. (Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021). It should be noted that surface tension and viscosity of the pseudo-solid for the 2-D axisymmetric model are adopted from Galeano-Rios et al. (Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021). For the 2-D cases focused on in this paper, the pseudo-solid parameters (i.e. σs/σl and $\mu_{s}/\mu_{l}$ in table 1) can also achieve the maximum deformation less than 5 % of the original cylinder diameter.
3.1. Four flow regimes
For the simplicity of our investigation, the Reynolds number for the liquid Rel is fixed in the following (as in table 1). Four typical flow regimes are identified in our simulations. Flow characteristics (including the interface evolution and the liquid velocity vector) of the four are shown in figure 3. For a given density ratio λρ, the four regimes (floating, bouncing, resurrecting and sinking) appear in sequence as the liquid Weber number Wel (which is more representative than the Weber number for the gas or pseudo-solid) increases. For the different regimes in figure 3, the corresponding mass centre location and velocity of the cylinder are shown in figure 4(a,b), respectively. An interesting observation can be made for the highlighted case of resurrecting (see figure 3(c) and the green dot-dashed line in figure 4a). When the cylinder reaches the lowest position (with the minimum value of y), the cylinder is not yet fully immersed and the surface tension still plays a role. The cylinder begins to ascend after passing the lowest position and then becomes fully immersed. It is founded that this interesting process typically accounts for approximately 10–15 % of the entire impact (from initial touch of τ = 0 to the moment of resurfacing, see the cross in figure 4).
For the floating regime (in figure 3a), the cylinder remains descending after coming into contact with the liquid, and then starts to ascend at the lowest position (see the image of τ = 5.6 in figure 3a) due to the recovery of the liquid–gas interface. Meanwhile, the cylinder remains attached to the liquid–gas interface. For the bouncing regime (in figure 3b), the cylinder movement is similar to that of the floating regime, except that the cylinder separates from the liquid during the ascent motion (see the image of τ = 19.9 in figure 3b). For the resurrecting regime (in figure 3c), the cylinder also undergoes a descent motion and then an ascent motion. During the ascent motion, the air cavity pinches off and the cylinder is fully submerged (see the image of τ = 12.8 in figure 3c). However, the cylinder will finally break through the interface (see the image of τ = 37.8 in figure 3c). For the sinking regime (in figure 3d), the cylinder will never break through the interface after fully submerged.
Figure 5 shows the phase diagram of the flow regimes for a cylinder impacting on the liquid surface in a parameter space (λρ, Wel). Only the cases of λρ ≥ 1 are considered, because the cylinder with its density smaller than the liquid density (λρ < 1) will manifest itself eventually even if the impacting velocity is extremely high. When the density of the cylinder is slightly smaller than that of the liquid, the impacting cylinder can become fully immersed at first while it may take a significant time for the cylinder to resurface, depending on the impact velocity V. For a specific density ratio λρ, as Wel increases from 0 to a large enough value, the four regimes, i.e. floating, bouncing, resurrecting and sinking, appear sequentially. It can be seen from figure 5 that the resurrection occurs in a very narrow region (between the solid curve and the dot-dashed curve). For a large enough density ratio λρ (larger than 1.3 or so), the resurrection will never occur.
In most cases of the resurrecting regime, we find that the cylinder will remain afloat after resurfacing (see figure 3(c) and black crosses in figure 5). Meanwhile, in very rare cases of the resurrecting regime, the cylinder will bounce and get detached from the liquid surface after resurfacing (denoted by red crosses in figure 5). However, for the superhydrophobic sphere impacting onto the liquid, the regime of resurrecting with bouncing has been found (see figure 11 in Galeano-Rios et al. Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021) but the regime of resurrecting with floating has not been found so far. Though the cases of resurrecting (with bouncing or floating) are rare for both the impacting cylinder and the impacting sphere, the resurrection with floating seems more possible to occur for the cylinder while the resurrection with bouncing seems more possible to occur for the sphere.
An interesting result can be found in figure 5 that the large density ratio λρ is helpful to the occurrence of bouncing rather than floating for the impacting cylinder (at a given value of Wel, e.g. Wel = 20). A similar result can also be seen in figure 7 of Lee & Kim (Reference Lee and Kim2008). Only the cases of bouncing and floating are considered here. When the superhydrophobic cylinder descends at the beginning, its mechanical energy is absorbed by the liquid. After the cylinder reaches the lowest point, the surface tension force lifts the superhydrophobic cylinder up and the energy returns to the cylinder (with a dissipation of the energy). Similar to the cases of the impacting sphere (Galeano-Rios et al. Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021), the cylinder with higher density (larger value of λρ) recovers relatively more energy during impact than the cylinder with lower density. Thus, the large density ratio λρ promotes the occurrence of bouncing.
3.2. Force analysis
For the force analysis of the cylinder, we use the dimensional forms for all physical quantities to aid in the comprehension. During the impact, the forces acting on the cylinder are the capillary force, the gravitational force and the hydrodynamic force. The hydrodynamic force can be regarded as the sum of the form drag, the hydrostatic pressure and the added inertia force.
Though there is no microscopic three-phase CL (point) in our simulations with the pseudo-solid approach (due to the very thin gas film between the solid and the liquid), the interface configurations in our simulations are relatively accurate because of the same essence in surface energies (see Appendix B). To analyse the forces on the cylinder theoretically, the position of the CL is estimated based on the configurations from the simulations. The position on the cylinder where the thickness of the thin gas film rapidly increases (i.e. the end of the thin gas film) is deemed as the contact point.
From the Newton's second law, the motion of the cylinder can be described as
where the mass of the cylinder is given by ${m_s} = {\rm \pi}{D^2}{\rho _s}/4$ and ma denotes the added mass. Here ${\bar{f}_p},{\bar{f}_d},{\bar{f}_\sigma }$ and $\textrm{ }{\bar{f}_G}$ are the vertical components of the hydrostatic pressure force, the form drag force, the surface tension force and the gravitational force, respectively. The vertical component of the hydrostatic pressure force ${\bar{f}_p}$ is equal to the total weight of liquid displaced by the wetted portion of the cylinder and the cavity above the wetted portion of the cylinder. Following Vella & Li (Reference Vella and Li2010) and Bush & Hu (Reference Bush and Hu2006), the vertical component of the hydrostatic pressure force can be expressed as
where ${S_w} = (\alpha - \sin \alpha \cos \alpha ){\rho _l}{D^2}/4$ is the volume of liquid displaced by the wetted portion of the cylinder, ${S_c} = {\rho _l}D\sin \alpha ( - h + (D/2)\cos \alpha )$ is volume of liquid displaced by the cavity above the wetted portion of the cylinder (see figure 6) and h is the vertical coordinate of the cylinder centre. The other three force components can be expressed as
where Cd denotes the drag coefficient. The added mass ${m_a}$ denotes the mass of liquid displaced by the wetted portion of the cylinder (where ${m_a} = {\rho _l}{S_w}$). With (3.2) and (3.3a–c), we can rewrite (3.1) as
For the cases of medium Wel in this paper, the impact process in each flow regime can be divided into three stages: slamming, pinning and sliding. Figure 7(a) shows the time evolution of the vertical components of the hydrodynamic force, the capillary force and the gravitational force exerted on the cylinder for the three stages. In figure 7(b), sinα can indicate the position of the CL and the magnitude of the vertical capillary effect. The slamming stage is the early period of impact process, at which the cylinder velocity decreases rapidly due to the hydrodynamic force. The CL rapidly climbs to the cylinder surface of α ≈ 90°, where the vertical component of the capillary force reaches almost its maximum value. In this stage, the movement of the cylinder is dominated by the hydrodynamic force. After climbing to the cylinder surface of α ≈ 90°, the CL is pinned at this position, and the pinning stage starts. The capillary force remains nearly constant due to pinning (see figure 7a) and the hydrodynamic force decreases slowly from a relatively low value. At this pinning stage, the most dominant force is the capillary force.
At the end of the pinning stage, the downward velocity of the cylinder further decreases to approximately zero. At the sliding stage, the cylinder begins to ascend under the effect of the capillary force. For a relatively small Wel (in the floating regime), the CL begins to slide downward from the pinned position (i.e. α decreases from 90°, see figure 3a). Meanwhile, for a relatively large Wel (in the resurrecting or sinking regime), the CL begins to slide upward from the pinned position (i.e. α increases from 90°, see figure 3c,d). In the bouncing regime, the CL may slide downwards continuously from the pinned position before the cylinder becomes detached from the liquid (e.g. Wel = 18 in figure 7b), or the CL will slide upwards from the pinned position and then slide downwards continuously before the detachment (e.g. Wel = 36 in figure 7(b), see the corresponding insets).
3.3. Early stage motion analysis
In this section, we present a simplified model to predict the early-stage motion of an impacting cylinder for the cases in the regime of floating and the regime of bouncing according to the above force analysis.
At the slamming stage, the hydrodynamic force is dominant. From Ji, Song & Yao (Reference Ji, Song and Yao2017), the magnitude of the hydrostatic pressure force is much smaller than the form drag and the added inertia force. The hydrostatic pressure force can be neglected throughout this stage. Due to the relatively high velocity of the cylinder during this stage (Wel > 5 in this paper), the form drag is significantly greater than other forces. Therefore, we only consider the form drag at this stage, and (3.4) can be rewritten as
Integrating (3.5) with the initial condition ${\bar{v}_s} =- V$, we can obtain the cylinder velocity as a function of the time at the slamming stage:
Previous studies on the drag coefficient of a partially immersed superhydrophobic cylinder remain inadequate. The experiments in Hunt et al. (Reference Hunt, Zhao, Silver, Yan, Bazilevs and Harris2023) indicate that the drag coefficient is 0–1.5 for a partially immersed superhydrophobic sphere, depending on the impacting velocity, the diameter and the vertical position. When fully immersed, the drag coefficient of the cylinder is approximately double that of the sphere with the same diameter for Re ~ 103 (White Reference White2016). Taking these into account, we find that the value of 1.5 can be appropriate for Cd at the slamming stage in the cases of this paper.
At the pinning stage of the cylinder, the capillary force starts to play a role and its vertical component can be seen as a constant (see figure 7a). In the meantime, the form drag and the added inertia force also play roles, because the cylinder velocity is still large enough at the whole stage. The volume of liquid displaced by cylinder is equal to half the volume of the cylinder, leading to α ≈ 90° or sin α ≈ 1 (see figure 7b). The gravity force is also taken into consideration. At this stage, (3.4) is therefore reduced to
Integrating (3.7), we can obtain the time evolution of the cylinder velocity at the pinning stage:
where C denotes a parameter to be determined.
Figure 8 shows the DNS result of the cylinder velocity for the case of figure 3(a), which is also compared with the theoretical prediction of (3.6) and (3.8). For this case, when $\tau = tV/D \approx 0.7$, which corresponds to the start of the pinning stage, the cylinder velocity is given by ${v_s} = {\bar{v}_s}/V =- 0.68$. Accordingly, we can obtain the parameter C = 0.846 in (3.8). For the drag coefficient of the partial immerged body at the pinning stage, we use Cd = 0.6 (in Vella & Li (Reference Vella and Li2010), Cd = 0.6 was also used in a similar case) in theoretical calculations.
3.4. Critical condition for resurrection
As shown in figure 5, the transition between the resurrecting regime and the sinking regime only exists for λρ ≤ 1.3. From figure 3(c,d), the cylinder will be completely submerged after the cavity above the cylinder pinches off. Then the cylinder continues to accelerate upward under the inertial propulsion of the surrounding fluid for a short time. Subsequently, the cylinder starts to decelerate due to the drag of the surrounding fluid. During the decelerating process, the forces acting on the cylinder are the buoyancy force, the form drag, the gravitational force and the added inertia force. Therefore, following Kim et al. (Reference Kim, Hasanyan, Gemmell, Lee and Jung2015), the motion of the fully immersed cylinder can be expressed as
where ${\bar{f}_b} = {\rm \pi}{D^2}{\rho _l}g/4$. For the fully submerged cylinder moving perpendicular to the axis of itself, by substituting $\alpha = {\rm \pi}$ into ${m_a} = (\alpha - \sin \alpha \cos \alpha ){\rho _l}{D^2}/4$, the added mass becomes ${m_a} = {\rm \pi}{D^2}{\rho _l}/4$ (see also White Reference White2016). The drag coefficient is given by Cd ≈ 1.5, which is the appropriate value for a fully submerged cylinder (see Vella & Li Reference Vella and Li2010). Equation (3.9) can therefore be rewritten as
We consider the critical cases that the cylinder velocity will precisely decelerate to zero at the moment when it breaks through the interface (see the cross in figure 4). At t = ts, the cylinder just reaches its maximum upwards velocity ${\bar{v}_{cr}}$ (see figure 9a). This velocity is critical for the regime of sinking and the regime of resurrecting. Solving the ordinary differential equation (3.10), it can be obtained that
where
For the case that is precisely at the transition between the resurrecting regime and the sinking regime, the cylinder velocity will precisely decelerate to zero at the moment when it breaks through the interface. By letting ${\bar{v}_s} = 0$ in (3.11), we can get the moment te when the cylinder just breaks through the interface. Integrating the velocity ${\bar{v}_s}$ in (3.11) from t = ts to t = te, we can obtain
where $\Delta {\bar{H}_{cr}}$ denotes the critical height difference between the point where the cylinder is just reaches its maximum upwards velocity vmax and the point where the cylinder just breaks through the interface (see figure 9a).
We find that the critical height difference $\Delta {\bar{H}_{cr}}$ almost remains constant for different density ratios, as is shown in figure 9(b). Notably, for a very large density ratio λρ (larger than 1.3 in figure 5), the regime of resurrecting does not exist, leading no critical case. As depicted in figure 5, the allowable parameter region for resurrecting is relatively narrow for both λρ and Wel, which may account for the constancy of the value of $\Delta {\bar{H}_{cr}}$. Thus, the critical velocity ${\bar{v}_{cr}}$ of the fully immersed cylinder can be obtained from (3.13) as
The critical velocity ${\bar{v}_{cr}}$ for different λρ calculated by (3.14) is compared with the corresponding DNS results, as shown in figure 10. During the whole entry, if the maximal velocity of the cylinder after fully immersed is larger than ${\bar{v}_{cr}}$, the resurrection will occur for the cylinder.
4. Conclusions
We have numerically investigated the impact of a superhydrophobic cylinder onto a liquid bath. The results demonstrate that with an increase in the Weber number, the system exhibits four distinct flow regimes, i.e. floating, bouncing, resurrecting and sinking regimes. The highlighted resurrecting regime pictures that, even if the cylinder is fully submerged in the liquid during the impact, it can still emerge from the liquid. With a further analysis, we find that when the cylinder reaches the lowest position, the cylinder is not fully immersed yet and the surface tension still plays a role. After reaching this position, the cylinder begins to move up when influenced by the surface tension force. Subsequently, the cylinder is fully immersed in the liquid and the surface tension effect on the cylinder disappears. For the resurrecting regime, the duration for the cylinder ascending from the lowest position to the full immersion position is approximately 10–15 % of the entire process (from initial touching the liquid surface to resurfacing). The fully immersed cylinder then continues moving up and resurfaces eventually, even if its speed is slowed down by the viscous drag.
From further numerical results, the phase diagram for the four regimes is presented in a parameter space of the Weber number and the particle-to-liquid density ratio. The resurrection regime for the impacting cylinder can only exist at a relatively narrow parameter space (where the Weber number is medium and the cylinder density is not very large). A force analysis is conducted for the impacting cylinder, and three stages (i.e. slamming, pinning and sliding) can be divided according to the effect of the surface tension (or the CL position). We find that the form drag dominates in the slamming stage while the surface tension force and the gravity force also play roles in the pinning stage, which can be applied to analyse the early stage motion of an impacting cylinder.
In order to predict whether the interesting resurrection phenomenon will occur for the impacting cylinder, the dynamics is investigated for the cylinder when it is fully immersed. We find that the critical height for the cylinder to resurface is almost constant for different cylinder densities. With this finding, the critical velocity of the fully immersed cylinder is derived, which can be used for the prediction of resurrection.
Funding
This research was supported in part by the National Natural Science Foundation of China (no. 11972170).
Declaration of interests
The authors report no conflict of interest.
Appendix A
To aid in the comprehension of the pseudo-solid approach, the volume fractions of liquid and solid (cl and cs) for the three phases are shown in figure 11. At the triple-phase region, the (surface tension) force balance can be expressed as
where σlg, σsg and σsl denote the surface tension coefficients for the liquid–gas interface, the solid–gas interface and the solid–liquid interface, respectively.
For the VoF method, the theoretical surface tension coefficient σl (σs), as presented in table 1, works in the domain where 0 < cl < 1 (0 < cs < 1). Thus, we have σlg = σl for the liquid–gas interface, while having σsg = σs for the solid–gas interface. The two points meet our expectations. However, for the solid–liquid interface, both 0 < cl < 1 and 0 < cs < 1 are satisfied, leading to the relation that
Notably, the relation (A2) is deduced from the characteristic of the VoF method, instead of the physical fact. With (A1a,b) and (A2), we can obtain the angle relation:
From (A3), it can be deduced that the contact angle $\theta$ is fixed, given by
The theoretical contact angle for the pseudo-solid approach is 180°, which accounts for the agreement with the superhydrophobic cases.
Appendix B
Though we have explained that the theoretical contact angle in the pseudo-solid approach is naturally 180°, there still exists an important problem. There is a very thin gas film (with its thickness negligible) between the superhydrophobic cylinder and the liquid in the numerical simulations, while this gas film may not necessarily be maintained in reality. For the cases with and without the thin gas film, why do the motions of the impacting superhydrophobic cylinder (or sphere in Galeano-Rios et al. Reference Galeano-Rios, Cimpeanu, Bauman, MacEwen, Milewski and Harris2021) agree well with each other? This will be explained from the perspective of energy.
When there is no gas film between the cylinder and the liquid, the sum of surface energies of the whole system can be expressed as
where Σsg, Σlg and Σsl denote the areas of the solid–gas interface, the liquid–gas interface and the solid–liquid interface, and σsg, σlg and σsl are the three corresponding surface energy per area.
When there is a very thin gas film between the cylinder and the liquid (considering the same configurations as in (B1)), the solid–liquid interface with area Σsl mentioned above becomes the liquid–gas interface and the solid–gas interface. Thus, the sum of surface energies of the whole system becomes
According to the well-known Thomas Young's relation, the contact angle can be expressed as
When the contact angle θ of the solid (e.g. the pseudo-solid in this paper) is equal to 180°, we have the following relationship from (B3):
With (B4), it can be found that the sum of surface energies of the whole system is the same whether there is a very thin film or not (i.e. E = E’) when θ = 180°. The surface energies govern the surface tension effect on the cylinder, and the surface tension forces are related the first variations of the surface energies (Kralchevsky et al. Reference Kralchevsky, Paunov, Denkov, Ivanov and Nagayama1993). Therefore, it is demonstrated why the liquid-entry behaviour of a superhydrophobic solid is less dependent on the microscopic details of the contact with the solid surface.
Appendix C
The grid independence analysis and the domain-size independence analysis are conducted in this appendix. As shown in figure 12(a), the simulations of different minimum mesh sizes are performed for a computational domain spanning 10D in width and 20D in length (as in figure 1). It is found that the maximum grid refinement of level 12 (i.e. the minimum grid size of 10D/212) is sufficient for the numerical simulations by comparisons with the cases of levels 11 and 13. Thus, the grid refinement of level 12 is adopted in our simulations. In order to obtain the appropriate size of the computational domain, we compare the numerical results for the widths of 7.5D, 10D and 15D (with the length–width ratio fixed at 2 as in figure 1) in figure 12(b). The numerical results for the width of 10D are very close to those for the width of 15D, which shows that the width of 10D is large enough for the computational domain for the numerical simulations. Thus, the computational domain that spans 10D in width and 20D in length is adopted in our simulations.
It is worth noting that, when the physical parameters of a simulation case are very close to the limits of the parameter regimes (see the demarcation curves of the phase diagram in figure 5), the fate of the impacting cylinder can possibly change significantly with a very slight disturbance in both physical and computational parameters (e.g. from sinking to resurrecting, see the snapshots in figure 3c,d). For example, in a case with physical parameters very close to the demarcation of regimes of sinking and resurrecting, a change in the minimum mesh size (e.g. from the maximum refinement level 11 to 12) may cause a significant difference for the trajectory of the impacting cylinder, which seemingly leads to a bad result for the convergence test on grid independency. However, for the cases with physical parameters a little far from the demarcation, the grid independency is easy to verify (see e.g. figure 12). Furthermore, when exploring the critical Weber numbers for different regimes (with other physical parameters the same, see figure 5), the critical values of Wel obtained under two different refinement levels may be different but the difference is very slight (typically, much less than 1). This indicates that our numerical results are relatively reliable though the discussed impacting behaviour can be very sensitive to parameters.