The Simulation of Liquid Flow in the Pore Network Model of Nanoporous Media

2020 ◽  
Vol 143 (3) ◽  
Author(s):  
Yaohao Guo ◽  
Lei Zhang ◽  
Hai Sun ◽  
Yongfei Yang ◽  
Zhi Xu ◽  
...  

Abstract The fluid–solid interaction force shows significant influence on liquid flow at nanoscale. Vast experimental observations in recent literatures have shown that Darcy's law cannot be applied to nanoporous media. In this study, the slip length and effective viscosity are adapted to characterize the nanoscale effect. First, the nanoscale effect is investigated in nanotubes through computational fluid dynamic (CFD) modeling analysis. Slip boundary condition has been studied as an important discrepancy between macroscopic flow and nanoscale liquid flow. The effect of viscosity change becomes more notable with the slip length increasing. Then, the flow equation for pore network modeling is developed to capture nanoscale effect. The results show that the apparent permeability of nanoscale systems is significantly underestimated when slip effect is neglected. The size of the pore throat determines whether the slip effect needs to be considered, and critical diameter of neglecting the slip effect for circular throat is 79.17 Ls. It is necessary to take the variation of effective viscosity into account under slip boundary condition. With the pore throat size decreasing, the nanoscale effect increases. The nanoscale effect is more sensitive to pore throat size under hydrophobic conditions than hydrophilic conditions.

Author(s):  
Derek C. Tretheway ◽  
Luoding Zhu ◽  
Linda Petzold ◽  
Carl D. Meinhart

This work examines the slip boundary condition by Lattice Boltzmann simulations, addresses the validity of the Navier’s hypothesis that the slip velocity is proportional to the shear rate and compares the Lattice Boltzmann simulations to the experimental results of Tretheway and Meinhart (Phys. of Fluids, 14, L9–L12). The numerical simulation models the boundary condition as the probability, P, of a particle to bounce-back relative to the probability of specular reflection, 1−P. For channel flow, the numerically calculated velocity profiles are consistent with the experimental profiles for both the no-slip and slip cases. No-slip is obtained for a probability of 100% bounce-back, while a probability of 0.03 is required to generate a slip length and slip velocity consistent with the experimental results of Tretheway and Meinhart for a hydrophobic surface. The simulations indicate that for microchannel flow the slip length is nearly constant along the channel walls, while the slip velocity varies with wall position as a results of variations in shear rate. Thus, the resulting velocity profile in a channel flow is more complex than a simple combination of the no-slip solution and slip velocity as is the case for flow between two infinite parallel plates.


Author(s):  
Mohammad Arif ◽  
Saurabh Kango ◽  
Dinesh Kumar Shukla

Abstract In the present study, the influence of various slip zone locations on the dynamic stability of finite hydrodynamic journal bearing lubricated with non-Newtonian and Newtonian lubricants has been investigated. Linearized equation of motion with free vibration of rigid rotor has been used to find the optimum location of the slip region with maximum stability margin limit. It has been observed that bearing with interface of slip and no-slip region near the upstream side of minimum film-thickness location is effective in improving the direct and cross stiffness coefficient, critical mass parameter, and critical whirling speed. The magnitude of dynamic performance parameters with slip effect is highly dependent on the rheology of lubricant. Shear-thinning lubricants combined with slip boundary condition shows higher dynamic stability as compared to the Newtonian lubricants under the conventional boundary condition. For all considered rheology of lubricants, the dynamic stability of bearing with slip effect is improving by increasing the eccentricity ratio.


Author(s):  
Qi Zhou ◽  
Chiu-On Ng

The hydrodynamic dispersion of a neutral non-reacting solute due to steady electro-osmotic flow in a circular channel with longitudinal step changes of zeta potential and hydrodynamic slippage is analyzed in this study. The channel wall is periodically micro-patterned along the axial position with alternating slip-stick stripes of distinct zeta potentials. Existing studies on electrically driven hydrodynamic dispersion are based on flow subject to either the no-slip boundary condition on the capillary surface or the simplification of lubrication approximation. Taking wall slippage into account, a homogenization analysis is performed in this study to derive the hydrodynamic dispersion coefficient without subject to the long-wave constraint of the lubrication approximation, but for a general case where the length of one periodic unit of wall pattern is comparable with the channel radius. The flow and the hydrodynamic dispersion coefficient are calculated numerically, using the packages MATLAB and COMSOL, as functions of controlling parameters including the period length of the wall pattern, the area fraction of the slipping region (EOF-suppressing) in a periodic unit, the ratio of the two zeta potentials, the intrinsic hydrodynamic slip length, the Debye parameter, and the Péclet number. The dispersion coefficient is found to show notable, non-monotonic in certain situations, dependence on these controlling parameters. It is noteworthy that the introduction of hydrodynamic slippage will generate much richer behaviors of the hydrodynamic dispersion than the situation with no-slip boundary condition, as slippage interacts with zeta potentials in the EOF-suppressing and EOF-supporting regions (either likewise or oppositely charged).


2021 ◽  
Vol 5 (4) ◽  
pp. 44
Author(s):  
Chen Zhang ◽  
Xuming Wang ◽  
Jiaqi Jin ◽  
Lixia Li ◽  
Jan D. Miller

Most reported slip length measurements have been made at the surfaces of synthetic materials and modified synthetic materials. In contrast, few slip length measurements at the surface of unmodified natural mineral surfaces have been reported. In this regard, flow at the silica face surfaces of the phyllosilicate minerals, talc and mica, was considered. A slip boundary condition was expected at the nonpolar hydrophobic silica surface of talc leading to enhanced flow, and a no-slip boundary condition was expected at the hydrophilic silica surface of mica. Atomic force microscopy (AFM) slip length measurements were made at the talc and mica surfaces. The slip length results for the hydrophobic silica surface of talc were contrasted to the results for the hydrophilic silica surface of mica (no-slip flow). The results are discussed based on molecular dynamics simulations (MDS), as reported in the literature, and AFM images of surface nanobubbles. For nonpolar hydrophobic surfaces (such as talc), it is doubtful that the MDS interfacial water structure and the water exclusion zone (3.2 Å) account for the AFM slip flow with slip lengths as great as 95 nm. Rather, a better explanation for the AFM slip flow condition is based on reduced interfacial viscosity due to the presence of dissolved gas and the accommodation of pancake nanobubbles at the talc surface having a height dimension of magnitude similar to the slip length.


1976 ◽  
Vol 77 (4) ◽  
pp. 665-684 ◽  
Author(s):  
E. B. Dussan V.

The singularity at the contact line which is present when the usual fluidmechanical modelling assumptions are made is removed by permitting the fluid to slip along the wall. The aim of this study is to assess the sensitivity of the overall flow field to the form of the slip boundary condition. Explicit solutions are obtained for three different slip boundary conditions. Two length scales emerge: the slip length scale and the meniscus length scale. It is found that on the slip length scale the flow fields are quite different; however, when viewed on the meniscus length scale, i.e. the length scale on which almost all fluidmechanical measurements are made, all of the flow fields appear the same. It is found that the characteristic of the slip boundary condition which affects the overall flow field is the magnitude of the slip length.


2012 ◽  
Vol 53 (4) ◽  
pp. 321-336 ◽  
Author(s):  
MICCAL T. MATTHEWS ◽  
KAREN M. HASTIE

AbstractA theoretical investigation of the unsteady flow of a Newtonian fluid through a channel is presented using an alternative boundary condition to the standard no-slip condition, namely the Navier boundary condition, independently proposed over a hundred years ago by both Navier and Maxwell. This boundary condition contains an extra parameter called the slip length, and the most general case of a constant but different slip length on each channel wall is studied. An analytical solution for the velocity distribution through the channel is obtained via a Fourier series, and is used as a benchmark for numerical simulations performed utilizing a finite element analysis modified with a penalty method to implement the slip boundary condition. Comparison between the analytical and numerical solution shows excellent agreement for all combinations of slip lengths considered.


Sign in / Sign up

Export Citation Format

Share Document