Rate of capillary rise in quartz nanochannels considering the dynamic contact angle by using molecular dynamics
Graphical abstract
Introduction
Capillary rise phenomena are ubiquitous in nature and play roles in important applications in the fields of powder engineering, geotechnical engineering, irrigation engineering and barrier engineering [[1], [2], [3], [4]]. Particularly in powder systems, capillary rise phenomena are of great importance in wetting processes and agglomeration processes for all types of powders, as well as in many physical behaviors involving the interaction between powders and water. For instance, the capillary rise effect can induce a dissolving liquid into agglomerate pores to improve the instant properties of food aggregations, which is widely applied in the food industry [5]. Another significant application for capillary rise is as a measurement technique to investigate the wettability [[6], [7], [8]] and surface free energy [9,10] of powders. The corresponding techniques indirectly calculate the above parameters by monitoring the rate of capillary rise. Therefore, insightful research on the capillary rise phenomenon and accurate calculation on the rising rate is essential for relevant industrial processes and measurement techniques of powder materials.
The wettability of powdered and granular solids as the description of particle-water interfacial interactions plays an important role within several operations such as granulation, rehydration, flotation and dissolution [4,11]. According to Young's .[12], the contact angle is a crucial indicator of the wetting properties of powders, as shown in Fig. 1(a), that can be measured in experiments. Several approaches, such as the sessile drop method, can be implemented to measure the contact angle. However, due to the particle size and surface roughness of powders, the measurement accuracy of these techniques is unsatisfactory [7]. Another measurement technique called the Capillary rise method (CRM) [13,14] is commonly used to assess the contact angle of powders or granular materials, in which CRM utilizes the capillary rise phenomenon along with Lucas-Washburn (L-W) equation [15,16]. A typical experimental setup for CRM is illustrated in Fig. 1(b). Although CRM has several advantages such as broad applicability and convenience, the measured contact angle measured during the wetting process is actually a non-equilibrium advancing contact angle, also called dynamic contact angle, and this parameter is dependent on the equilibrium value as well as the rate of capillary rise [17,18]. CRM neglects the dynamics process of contact angle, leading to an overestimated value of the equilibrium contact angle.
Molecular dynamics(MD) simulation is a valuable computational approach for elucidating the dynamic and thermodynamic properties of materials [19]. This method can calculate interparticle forces well, including van der Waals attractive, double-layer repulsive and capillary forces, which are the dominant forces between powder particles. Several MD simulations have been conducted on different capillary flows in nanometer pores made of different solid materials [[20], [21], [22], [23]]. However, these reports mainly focus on the transport behavior of capillary flow in nanochannels without considering the dynamic variation in the contact angle during the capillary rise. Moreover, to the best of our knowledge, few studies have involved the capillary rise phenomenon in powder materials even though it has a great impact on the physical and chemical behavior and processing technology of powder materials as mentioned above.
In this paper, the experimental process of CRM in the nanometer channels of α-quartz was simulated by using molecular dynamics. Three quartz models with different channel height were established to study the impact of the geometric restriction on the capillary rise phenomenon of water. The penetration rate and dynamic contact angle of water were measured during the capillary rise process. A modified L-W equation taking into account the concept of the effective viscosity, dynamic contact angle and slip length was proposed. The accuracy of the equilibrium contact angle obtained by the classical and modified L-W equation was discussed. Finally, the interface structure of water in the slits was further analyzed.
Section snippets
Theoretical analysis
Lucas and Washburn discovered that the capillary filling process is proportional to the square root of time, through a large number of capillary rise tests. On macroscopic scales, the law has been verified in various capillary imbibition systems, which are often used to characterize porous media [24]. The classic L-W equation demonstrates that capillary rise in the nanochannel is controlled by the balance force of capillary and viscosity forces [25], which can be expressed by:
Simulation method
To investigate the capillary interaction of α-quartz (the most common type of quartz) channels with water in the imbibition process, we established three models with different channel height of 1.5 and 2.5, 3.5 nm. Each model consisted of two parts of the α-quartz nanochannel and water reservoir. As an illustration, Fig. 2 shows the initial configuration of the water imbibition model with a 3.5 nm height. The nanochannel was built by two parallel α-quartz sheets with 9.8 nm depth, 2.2 nm
Results and discussion
In the capillary simulations, water in the reservoir can spontaneously flow into the quartz channels. Fig. 3 displays the meniscus displacement of water fluid in the 3.5 nm channel during the capillary rise process. The snapshots of the meniscus at 50, 200 and 700 ps in Fig. 3 show that the wetting frontier maintains the shape of a concave interface, and the angle of the meniscus and the capillary wall drops gradually as the height of the water column increases. The water penetration depth with
Conclusions
In this paper, molecular dynamics were used to investigate the capillary rise process of water in quartz with channels of different sizes. The change in the dynamic contact angle, the rate of penetration and the interface structure of pore water were studied by using MD method. A modified L-W method has been developed by introducing the effective viscosity, slip length and dynamic contact angle. The differences in equilibrium contact angles calculated by the classical and modified L-W equations
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgements
This research has been supported by the China National Key R&D Program during the 13th Five-year Plan Period (Grant No. 2018YFC1505104 and 2017YFC1503103) and National Natural Science Foundation of China (Grants No. 51778107). The work was carried LvLiang Cloud Computing Center of China, and the calculations were performed on TianHe-2.
References (63)
- et al.
Impact of hydrophobic surfaces on capillary wetting
Powder Technol.
(2018) - et al.
Wetting, disintegration and dissolution of agglomerated water soluble powders
Powder Technol.
(2011) - et al.
Measurements of the wettability of catalyst support materials using the Washburn capillary rise technique
Powder Technol.
(2015) - et al.
Wettability of mineral and metallic powders : applicability and limitations of sessile drop method and Washburn ’ s technique
Powder Technol.
(2012) - et al.
Colloids and surfaces A: physicochemical and engineering aspects wetting dynamics and contact angles of powders studied through capillary rise experiments
Colloids Surfaces A Physicochem. Eng. Asp.
(2013) - et al.
Effect of chemical modification on surface free energy components of Aerosil silica powders determined with capillary rise technique
Powder Technol.
(2013) - et al.
A novel method for determining surface free energy of powders using Washburn ’ s equation without calculating capillary factor and contact angle
Powder Technol.
(2015) Contact angle, adsorption and wettability-a review with respect to powders
Powder Technol.
(1990)- et al.
Comparing capillary rise contact angles of soil aggregates and homogenized soil
Geoderma.
(2008) - et al.
Effect of dynamic contact angle on capillary rise phenomena
Colloid Surf A Physicochem Eng Asp.
(2000)
Molecular dynamics simulations of nonylphenol ethoxylate on the Hatcher model of subbituminous coal surface
Powder Technol.
Na and cl immobilization by size controlled calcium silicate hydrate nanometer pores
Constr. Build. Mater.
Capillary dynamic under nanoconfinement: coupling the energy dissipation of contact line and confined water
Int. J. Heat. Mass Tran.
Capillary filling of confined water in nanopores: coupling the increased viscosity and slippage
Chem. Eng. Sci.
Dynamics of partial wetting
Adv. Colloid Interf. Sci.
Kinetics of liquidliquid displacement
J Colloid Interf. Sci.
Capillary dynamics drvien by molecular self-layering
Adv. Colloid Interf. Sci.
The kinetics of wetting in a capillary
J. Colloid Interface Sci.
Fast parallel algorithms for short-range molecular dynamics
J. Comput. Phys.
Dynamics of capillary rise
J. Colloid Interface Sci.
Multiphase capillary flows
Int. J Multiphas Flow.
A theory of capillary rise of a liquid in a vertical cylindrical tube and in a parallel-plate channel
J. Colloid Interface Sci.
Modeling the confined fluid flow in micro-nanoporous media under geological temperature and pressure
Int. J. Heat. Mass Tran.
The effect of particle size and moisture on the tensile strength of microcrystalline cellulose powder
Powder Technol.
The dynamic strength of partially saturated powder compacs: the effect of liquid porperties
Powder Technol.
Dry particle coating using magnetically assisted impaction coating: modification of surface properties and optimization of system and operating parameters
Powder Technol.
Effects of adsorptive water on the rupture of nanoscale liquid bridges
Appl. Clay Sci.
Confined fluids in clay interlayers: a simple method for density and abnormal pore pressure interpretation
Colloids Surfaces A Physicochem. Eng. Asp.
Rate of capillary rise in soil
J. Geotech. Geoenvironmental Eng.
Water movement in glass bead porous media 3. Theoretical analyses of capillary rise into initially dry media
Water Resour. Res.
A study of the applicability of the capillary rise of aqueous solutions in the measurement of contact angles in powder systems
Powder Technol.
Cited by (16)
Two-phase imbibition of water-oil displacement in silica nanochannels
2024, International Journal of Multiphase FlowEffect of particle size distribution on dynamic properties of cemented coral sand under SHPB impact loading
2022, Soil Dynamics and Earthquake EngineeringCitation Excerpt :Due to the significant difference of formation process, the mechanical properties of coral sand are quite different from that of terrestrial river sand [3–5]. Coral sand has the characteristics of high intra-particle porosities [6] and highly irregular particle shape [7], which makes it susceptible to particle breakage under load [8]. During the past decades, many scholars have studied the mechanical behavior of coral sand under loading.
Molecular dynamics study of spontaneous capillary flow and heat transfer in nanochannels
2022, International Communications in Heat and Mass TransferCitation Excerpt :The L-W equation was modified by the viscosity to more accurately predict the capillary flow process of water in copper and gold channels. Zhang et al. [13] investigated the capillary transport behavior of water flowing through α-quartz nanochannels considering the channel height, and the results showed that the capillary rise method with the classical L-W equation might obtain an overestimated equilibrium contact angle. By introducing the concepts of effective viscosity, slip length, and dynamic contact angle, an improved L-W equation was established.
Molecular dynamics modeling of hydraulic conductivity of soil considering variable viscosity and adsorbed water
2022, Applied Clay ScienceCitation Excerpt :Experimental data of the water viscosity in quartz nanochannels at room temperature measured using the capillary rise process in previous literature are compared with MD results (Haneveld et al., 2008; Jarlgaard et al., 2006; Persson et al., 2007; Thamdrup et al., 2007). These measurements can be used to determine the viscosity by using the Lucas-Washburn equation with a constant contact angle, neglecting the effect of dynamic contact angles in the capillary rise, which leads to an overestimate of the water viscosity (Feng et al., 2018; Zhang and Pei, 2020). This may be the reason for the difference between the simulation and experimental data.
Determining the bound water content of montmorillonite from molecular simulations
2021, Engineering GeologyCitation Excerpt :The dominant forces of clay particle-particle and water-particle, including van der Waals attractive, double-layer repulsive and capillary forces, can be accurately calculated by the MD method (Amarasinghe and Anandarajah, 2011; Ichikawa et al., 1999). Because the diameter of clay pore is in the nanometer to micron range, MD has been successfully applied to study behaviors of clay mineral-water systems such as EDL (Bourg and Sposito, 2011; Skipper et al., 1995), wettability (Zhang et al., 2016), swelling properties (Kawamura et al., 1999) and water transportation (Zhang and Pei, 2020; Hou et al., 2020). Therefore, MD can provide insight into the properties of water near the clay surface.
Molecular dynamics study on the zeta potential and shear plane of montmorillonite in NaCl solutions
2021, Applied Clay Science