Transitions from static wetting to steady dewetting and deposition of liquid-film on partially wettable polymer surface

The transitions from static to steadily moving wetting perimeter and further to deposition of a liquid-film on partially wettable surface were studied with the same system under the same conditions. A polyethylene terephthalate (PET) tape was vertically withdrawn at constant velocity from glycerol–water mixture. Elevation L of the three-phase contact line above the liquid level was measured under static, steady, and dynamic wetting. The static receding Θ_R and the apparent dynamic angles Θ_app at different withdrawal velocities U were calculated from the static relationship Θ(L). It was found that the limiting static angle Θ_R,min, at which the wetting perimeter starts moving, depends on withdrawal velocity. Extrapolation of the Θ_R,min/U dependence to U = 0 yields the quasi-static value of this parameter , that coincides with the relaxation static angle Θ_R,rlx achieved after meniscus motion ceases. This conclusion holds also for the wetting mode, where the limiting static advancing angle = Θ_A,rlx. Both the limiting and relaxation angles could be used for calculation of the effective Young's contact angle on non-ideal surface following Adam's suggestion [N.K. Adam, Adv. Chem. Ser. 43 (1964) 53.].The critical velocity U_cr anfd apparent dynamic angle Θ_app,cr, at which transition between steady dewetting and dynamic wetting occurs, were determined. The value of Θ_app,cr = 0° ± 5° agrees with our previous results [R.V. Sedev, J.G. Petrov, Colloids and Surfaces, 53 (1991) 147] implying a quasi-static shape of the moving meniscus up to U_cr. At U > U_cr, the speed V of the contact line relative to the solid wall is independent of withdrawal velocity and thickness of the deposited film. The present data confirm the earlier findings [J.G. Petrov, R.V. Sedev, Colloids and Surfaces, 13 (1985) 317, T.D. Blake, K.J. Ruschak, Nature, 282 (1979) 489] that at U = U_cr, V reaches its maximum value V_max, which is most important parameter of dewetting kinetics.Weak linear decrease of Θ_app with U was found above Ca = 2.4 × 10⁻⁵ up to the critical capillary number Ca_cr = 4.1 × 10⁻⁴. Below Ca = 10⁻⁵ the apparent receding angle depends much stronger on withdrawal velocity. The hydrodynamic (HD), and the simple and more general versions of the molecular-kinetic (MK) and molecular-hydrodynamic (MHD) theories of the wetting dynamics were used for quantitative characterization of the system in the steady dewetting regime. The effective Young's angle was used in the MK and MHD treatment of the experimental data following our previous publication [J.G. Petrov, J. Ralston, M. Schneemilch, R. Hayes, J. Phys. Chem B, 107(7) (2003) 1637]. The HD theory only qualitatively satisfies our experimental data giving physically unreasonable value of the hydrodynamic cut-off (slip) length and too small static receding angle at U = 0. The MK theory gives acceptable values of the oscillation frequency K₀ of the molecules at the contact line. Its more general version, including the viscous dissipation in the contact line vicinity, yields higher oscillation frequency. Very large distance λ between adsorption centers on the solid substrate (about five times the diameter of a glycerol molecule) was obtained with both MK and MHD theories. The too small frequency K₀ obtained with the simple MHD theory is removed by the more general version, accounting for contact line and viscous friction in the inner and intermediate zone of the moving meniscus. All theories show discrepancies between theoretically expected and experimentally estimated values of some of the parameters of wetting dynamics.     

Surface Tension and Dynamic Contact Angle of Water in Thin Quartz Capillaries

The surface tension of water has been measured in quartz capillaries with radii from 200 down to 40 nm. It appears that the surface tension does not differ from the known (bulk) values in the temperature range from 8 to 70 degrees C, within 1% experimental error. The dynamic contact angle, theta(d), vanishes when the capillary surface is covered with a wetting film left behind the receding meniscus. In the case of a dry surface, theta(d) depends on the velocity of the meniscus motion. The results obtained do not agree with presently available theoretical predictions from hydrodynamic theories of dynamic contact angles. Rather the kinetics of water vapor adsorption ahead of the moving meniscus seems to be the major controlling agent of the dynamic contact angle. Copyright 2000 Academic Press.     

Novel and global approach of the complex and interconnected phenomena related to the contact line movement past a solid surface from hydrophobized silica gel

In this work a generalized hydrodynamic theory for water flow into a mesoporous matrix from hydrophobized silica gel is suggested. Although we examine a fluid dynamics problem, the motion of the water-gas-solid contact line past a hydrophobized silica gel surface, motivation for such research derives from the investigation of a novel principle of mechanical energy dissipation, called surface dissipation, and its attached machine element, named a colloidal damper (CD). Similar to a hydraulic damper, this absorber has a cylinder-piston structure, but oil is replaced by a colloid consisting of a mesoporous matrix and a lyophobic liquid. Here, the mesoporous matrix is from silica gel modified by linear chains of alkyldimethylchlorosilanes and water is the associated lyophobic liquid. Mainly, the colloidal damper energy loss can be explained by the dynamic contact angle hysteresis in advancing (liquid displaces gas) and receding motion (gas displaces liquid); such hysteresis occurs due to the geometrical and chemical heterogeneities of the solid surface. Although this new kind of dissipation could be attractive for many applications, the subject remains almost unexplored in the scientific literature. Many different, complex, and interconnected aspects are related to this subject: capillary hydrodynamics, slippage effect, contact angle hysteresis, estimation of dissipated energy, thickness optimization of the grafted layer on the surface of the mesoporous matrix, etc. For this reason, a novel and global approach to all the complex and interconnected phenomena related to the contact line movement past a solid surface from hydrophobized silica gel is proposed. Our approach has a modest experimental basis but this is compensated for with rich references to other experimental and theoretical work oriented to the study of surface phenomena in such systems. We tried to sort the existing results and to find the right place for each in building our global view of the problem. This work is structured as follows. The measurement technique of the hysteresis loop is described. From experimental data one calculates the dissipated energy versus length of the grafted molecule on the silica gel surface. These results are justified by flow analysis. Generalized hydrodynamic theory means here that the basic structure of the Navier-Stokes equations is kept, but in order to include the relation between macroscopic flow and molecular interactions, slip is allowed on the solid wall. The nanopillar architecture of the silica gel hydrophobic coating is described. Concepts of slip and contact angle hysteresis are detailed and their connection is revealed. During adsorption, water penetrates the pore space by maintaining contact with the top of the coating molecules (region of -CH(3) groups); after that, water is forced into and partially or totally fills the space between molecules (region of -CH(2) groups). In such circumstances, at the release of the external pressure, desorption occurs. An original energetic-barriers approach is proposed to understand the filling of the nanosize canals which occur in the hydrophobic grafted layer. Employing this energetic-barriers approach, one finds the optimum length of the grafted molecule which maximizes the dissipated energy of the CD reversible cycle. Such results are useful for the appropriate design of ultrahydrophobic surfaces in general, and for the optimal design of a hydrophobic coating of a mesoporous matrix destined for CD use.     

Deriving Molecular Parameters of Interfaces Between Chemically Identical Polymers from Macroscopically Observed Phenomena

Summary: By utilizing model systems, macroscopically observable phenomena like friction or dewetting allow to identify and to quantify molecular interfacial parameters like molecular interpenetration depth, interfacial tension, or slippage length between grafted and free chemically identical polymers. We present experimental results, which permit to extract these parameters from simple contact angle measurements or by following the dewetting process in real time with a simple optical microscope. We also show how these model experiments can provide valuable insight and fundamental understanding of processes like polymeric friction and adhesion.     

低能量链端对聚苯乙烯薄膜铺展和润湿行为的影响

通过低能量功能端基的表面富集作用,研究了聚苯乙烯（ＰＳ）薄膜在聚甲基丙烯酸甲酯（ＰＭＭＡ）表面上的铺展和润湿动力学．用光学显微镜跟踪了ＰＳ薄膜的润湿行为,并对高分子熔体膜中非连续部分尺寸的增大速率进行了测定．分别用ＸＰＳ和ＡＦＭ对ＰＳ薄膜的表面组成和ＰＳ液滴的平衡接触角进行了测定．发现具有低表面能的氟碳端基在薄膜表面富集使ＰＳ薄膜的表面张力下降,并使ＰＳ液滴在ＰＭＭＡ表面上的平衡接触角减小,从而使高分子熔体膜中非连续部分尺寸的增长速率下降,得到了与液液界面铺展和润湿理论一致的实验结果．     

The Neumann and Young equations for nematic contact lines

The Neumann and Young equations for three-phase nematic contact lines have been derived using the momentum balance equation and classical liquid crystal physics theories. The novel finding is the presence of bending forces, originating from the anchoring energy of nematic interfaces, and acting on the contact line. The classical Neumann triangle or tensile force balance becomes in the presence of a nematic phase the Neumann pentagon, involving the usual three tensile forces and two additional bending forces. The Young equation that describes the static contact angle of a fluid in contact with a rigid solid is again a tensile force balance along the solid, but for nematics it also involves an additional bending force. The effects of the bending forces on contact angles and wetting properties of nematic liquid crystals are thoroughly characterized. It is found that in terms of the spreading coefficient, bending forces enlarge the partial wetting window that exists between dewetting and spontaneous spreading. Bending forces also affect the behaviour of the contact angle, such that spreading occurs at contact angles greater than zero and dewetting at values greater than pi. Finally, the contact angle range in the partial wetting regime is always less than pi.     

The Neumann and Young equations for nematic contact lines

The Neumann and Young equations for three-phase nematic contact lines have been derived using the momentum balance equation and classical liquid crystal physics theories. The novel finding is the presence of bending forces, originating from the anchoring energy of nematic interfaces, and acting on the contact line. The classical Neumann triangle or tensile force balance becomes in the presence of a nematic phase the Neumann pentagon, involving the usual three tensile forces and two additional bending forces. The Young equation that describes the static contact angle of a fluid in contact with a rigid solid is again a tensile force balance along the solid, but for nematics it also involves an additional bending force. The effects of the bending forces on contact angles and wetting properties of nematic liquid crystals are thoroughly characterized. It is found that in terms of the spreading coefficient, bending forces enlarge the partial wetting window that exists between dewetting and spontaneous spreading. Bending forces also affect the behaviour of the contact angle, such that spreading occurs at contact angles greater than zero and dewetting at values greater than pi. Finally, the contact angle range in the partial wetting regime is always less than pi.     

Contact line motion in confined liquid-gas systems: Slip versus phase transition

In two-phase flows, the interface intervening between the two fluid phases intersects the solid wall at the contact line. A classical problem in continuum fluid mechanics is the incompatibility between the moving contact line and the no-slip boundary condition, as the latter leads to a nonintegrable stress singularity. Recently, various diffuse-interface models have been proposed to explain the contact line motion using mechanisms missing from the sharp-interface treatments in fluid mechanics. In one-component two-phase (liquid-gas) systems, the contact line can move through the mass transport across the interface while in two-component (binary) fluids, the contact line can move through diffusive transport across the interface. While these mechanisms alone suffice to remove the stress singularity, the role of fluid slip at solid surface needs to be taken into account as well. In this paper, we apply the diffuse-interface modeling to the study of contact line motion in one-component liquid-gas systems, with the fluid slip fully taken into account. The dynamic van der Waals theory has been presented for one-component fluids, capable of describing the two-phase hydrodynamics involving the liquid-gas transition [A. Onuki, Phys. Rev. E 75, 036304 (2007)]. This theory assumes the local equilibrium condition at the solid surface for density and also the no-slip boundary condition for velocity. We use its hydrodynamic equations to describe the continuum hydrodynamics in the bulk region and derive the more general boundary conditions by introducing additional dissipative processes at the fluid-solid interface. The positive definiteness of entropy production rate is the guiding principle of our derivation. Numerical simulations based on a finite-difference algorithm have been carried out to investigate the dynamic effects of the newly derived boundary conditions, showing that the contact line can move through both phase transition and slip, with their relative contributions determined by a competition between the two coexisting mechanisms in terms of entropy production. At temperatures very close to the critical temperature, the phase transition is the dominant mechanism, for the liquid-gas interface is wide and the density ratio is close to 1. At low temperatures, the slip effect shows up as the slip length is gradually increased. The observed competition can be interpreted by the Onsager principle of minimum entropy production.     

Motion of a drop on a solid surface due to a wettability gradient

The hydrodynamic force experienced by a spherical-cap drop moving on a solid surface is obtained from two approximate analytical solutions and used to predict the quasi-steady speed of the drop in a wettability gradient. One solution is based on approximation of the shape of the drop as a collection of wedges, and the other is based on lubrication theory. Also, asymptotic results from both approximations for small contact angles, as well as an asymptotic result from lubrication theory that is good when the length scale of the drop is large compared with the slip length, are given. The results for the hydrodynamic force also can be used to predict the quasi-steady speed of a drop sliding down an incline.     

A Study of the Authophobic Dewetting in Thin Polymer Films

We have investigated the dewetting behaviour of a thin polystyrene film on top of cross-linked network of the same polymer. By changing the cross-linking density, the dynamic of the dewetting is modified. The behaviour of the contact angle and dewetting velocity can be related to the interfacial width between the polymers. We observe moreover that while the lost of entropy of the network connected to the penetration of the free chains will favour the dewetting; the possible presence of “connectors” between the layers will tend to stabilize the wetting. By modifying the entropy of the mesh by swelling the cross-linked system with a reservoir of the same polymer, the mesh is saturated and a stable interface without interdiffusion between the layers is obtained: in this case a constant contact angle and slippage length are obtained. Networks with higher cross-linking density are more difficult in general to swell – the mesh size is smaller - and the complete saturation is not reached.     

Does macroscopic flow geometry influence wetting dynamic?

The macroscopic flow geometry has long been assumed to have little impact on dynamic wetting behavior of liquids on solid surfaces. This study experimentally studied both spontaneous spreading and forced wetting of several kinds of Newtonian and non-Newtonian fluids to study the effect of the macroscopic flow geometry on dynamic wetting. The relationship between the dynamic contact angle, θ(D), and the velocity of the moving contact line, U, indicates that the macroscopic flow geometry does not influence the advancing dynamic wetting behavior of Newtonian fluids, but does influence the advancing dynamic wetting behavior of non-Newtonian fluids, which had not been discovered before.     

Drying and hydrophobic collapse of paraffin plates

We perform molecular dynamics simulations of the hydrophobic collapse of two paraffin plates to examine how the collapse is mediated by realistic paraffin-water attractive van der Waals forces. We explore several aspects of the drying transition between the plates, including the critical separation for drying and the critical size of the vapor bubble required for the nucleation of the drying event. We also investigate the kinetics of hydrophobic collapse and find that the hydrophobic collapse occurs in about 100 ps. We compare these results with the simulations with the plate-water van der Waals attractions turned off and with recent results on the hydrophobic collapse of multidomain proteins. Last, we discuss the relationship among the dewetting transition critical distance, van der Waals potential well depth, and water contact angle on solute surface using a simple macroscopic theory.     

Rim instability by solvent-induced dewetting

We study the condition of the occurrence of the rim instability in the solvent-induced dewetting process. Our experimental results show that the film thickness not only greatly influences the occurrence of the rim instability, but also influences the wavelength lambda as characterized by the undulation of the deformed contact line. The molecular weight of polymer does not almost influence the occurrence of the rim instability and the wavelength lambda. The wavelength lambda is proportional to the width of the rim in the rim instability region. The receding contact angle theta of polymer solutions on substrates in the dewetting process is an important factor to influence the rim instability in the solvent-induced dewetting.     

Dynamic wetting of a fluoropolymer surface by ionic liquids

The spontaneous spreading of ionic liquids on a fluoropolymer surface (Teflon AF1600) in air is investigated by high-speed video microscopy. Six ionic liquids (EMIM BF(4), BMIM BF(4), OMIM BF(4), EMIM NTf(2), BMIM NTf(2) and HMIM NTf(2)) are used as probe liquids. The dependence of the dynamic contact angle on contact line velocity is interpreted with a hydrodynamic model and a molecular-kinetic model. The usefulness of the hydrodynamic model is rather limited. There is a good correspondence between the molecular dimensions of the liquids and the physical parameters of the molecular-kinetic model. The viscous and molecular-kinetic contributions to energy dissipation are calculated, revealing that energy is dissipated in the bulk as well as at the contact line during dynamic wetting. There are wide ramifications of these results in areas ranging from lubrication and biology to minerals processing and petroleum recovery.     

Droplet spreading on microstriped surfaces

When a droplet of fluid is deposited on a surface with chemical and/or topological patterns, its static shape is highly dependent on the 2D distribution of the patterns. In the case of chemical stripes, three distinct spreading regimes have been observed as a function of wettability contrast between the two kind of stripes. For low wettability contrast, the droplet spreads with the same [corrected] velocity normal and parallel to the stripes [corrected] and the macroscopic contact angle is close to Cassie's contact angle. When the wettability contrast is intermediate/high, the resulting shape of the droplets is elongated. In the intermediate wettability contrast regime, an ideal situation shows stick and slip behavior of the contact line, during which the contact line jumps from one stripe to another. For a high wettability contrast, the confinement of the fluid between two chemical stripes leads to a 2D spreading.     

On the kinetics of the capillary imbibition of a simple fluid through a designed nanochannel using the molecular dynamics simulation approach

A molecular dynamics (MD) approach was employed to simulate the imbibition of a designed nanopore by a simple fluid (i.e., a Lennard-Jones (LJ) fluid). The length of imbibition as a function of time for various interactions between the LJ fluid and the pore wall was recorded for this system (i.e., the LJ fluid and the nanopore). By and large, the kinetics of imbibition was successfully described by the Lucas-Washburn (LW) equation, although deviation from it was observed in some cases. This lack of agreement is due to the neglect of the dynamic contact angle (DCA) in the LW equation. Two commonly used models (i.e., hydrodynamic and molecular-kinetic (MK) models) were thus employed to calculate the DCA. It is demonstrated that the MK model is able to justify the simulation results in which are not in good agreement with the simple LW equation. However, the hydrodynamic model is not capable of doing that. Further investigation of the MD simulation data revealed an interesting fact that there is a direct relationship between the wall-fluid interaction and the speed of the capillary imbibition. More evidence to support this claim is presented.     

Hydrodynamic Activation of the Batch-Polymerization of Methyl methacrylate in a Taylor-Couette Reactor

Summary: In this work, the free radical batch polymerization of methyl methacrylate (MMA) premixed with xylene as a solvent, in the presence of an initiator, 2,2-azoisobutyronitrile (AIBN), in the Taylor-Couette reactor was studied. We observed an unexpected influence of hydrodynamic process parameters, i.e. angular velocity ω, on the polymer conversion, molecular weight and viscosity of the produced polymer. The polymerization process seems to be activated by hydrodynamic process parameters. Hydrodynamic activation is a promoting effect of process parameters on polymer product properties. The hydrodynamic activation is found to depend on the reaction time and the angular velocity of the inner cylinder. In addition, our results highlight both the reaction kinetics and the hydrodynamics during the polymerization. The conversion exhibits a significant difference between tests with and without the angular velocity of the inner cylinder. The conversion and the molecular weight strongly increase with the increase of the angular velocity of the inner cylinder, whereas the viscosity is less strongly dependent. There is more increase with decreasing solvent concentration. The radial Reynolds number decreases with increasing conversion. The polymerization is faster with a low solvent concentration, and the molecular weight is higher compared to the case of high solvent concentration.