Liquid Slip Research Articles

Acoustic Fingerprinting and Nanoslip Dynamics of Biofilms

AbstractIt is reported that bacteria can generate nanomotion, but understanding the complex dynamics of bacterial colony gliding on solid interfaces has remained unresolved. Here, this work captures the real‐time development and gliding of bacterial biofilms on vibrating solids made of piezoelectric quartz. The gliding, characterized by liquid slips, is measured in form of frequency and dissipation changes of the vibrating solid. These vibrations enable the generation of distinct acoustic fingerprints (sound/ music) of the three phases of biofilm development: viscoelastic strengthening, biofilm growth and biofilm stability. In adition, the effect of extracellular matrix secretion on the rigidity of the film and its nanoslip in each of the distinct biofilm developmental phases is quantified. This work provides a real‐time, label‐free method of quantifying bacteria biofilm dynamics and paves the way for developing libraries of acoustic signatures of bacteria and their metabolic products.

Highly transparent, self-cleaning, and UV-shielding composite coating: When eco-friendly waterborne omniphobic surface cooperates with quantum dots

Solar energy is an inexhaustible carbon-free source, and solar cells are invaluable enablers for achieving carbon neutrality; however, the surface contamination on the glass of solar cells poses a formidable obstacle to maintaining a high power conversion efficiency. Existing self-cleaning coatings for combating surface contamination generally lack the highly comprehensive properties required for applications in real-world conditions. Herein, eco-friendly, highly transparent, waterborne omniphobic Si2.5WPA/AmZnO@SiO2 coatings were developed by the ingenious design of polydimethylsiloxane-grafted waterborne polyacrylate (SixWPA) and (3-aminopropyl)triethoxysilane-modified-SiO2 encapsulated ZnO quantum dots (AmZnO@SiO2). The enriched PDMS flexible chains on the coating surface overcame the energy barriers for moving liquids, and the AmZnO@SiO2 inside absorbed and scattered UV rays, which synergistically endowed the coating with self-cleaning and UV-shielding properties. Importantly, after the self-cleaning procedure, the light transmission of dust-contaminated glass and short-circuit current density of corresponding solar cells could be restored to 99.78 % and 97.77 % of initial metrics, respectively. Furthermore, the commonly ignored self-cleaning mechanism was investigated and proposed: the coating weakened the adhesion of dust particles and utilized the drag force generated during liquid slip to achieve self-cleaning. The coating also showed good UV/humidity/thermal stability and mechanical strength, providing a design blueprint for outdoor optical device protection.

Influence of sedimentary structure and pore-size distribution on upscaling permeability and flow enhancement due to liquid boundary slip: A pore-scale computational study

Low-permeability sedimentary formations, such as tight sandstones, exhibit fluid flow and transport phenomena distinct from those in conventional porous systems due to the dominance of micro- to nanometer-sized pores and variable amounts of boundary slip. The widely used traditional no-slip boundary condition often fails to accurately describe fluid behavior in these formations. A knowledge gap exists in understanding how liquid slip influences fluid dynamics in complex, heterogeneous sedimentary structures, as previous studies have primarily focused on simplified, homogeneous pore geometries. In this study, we investigated the impact of boundary slip on low-Reynolds number fluid dynamics within synthetically designed two-dimensional graded and random pore networks with varying pore-size distributions to account for heterogeneity. Our results showed that velocity variance increased with increasing heterogeneity, following a power-law relationship. The power-law exponents decreased with boundary slip, quantifying how boundary slip mitigated the impact of heterogeneity on velocity variance. We developed a theoretical model to predict asymptotic flow enhancement and derived constitutive relations to estimate the coefficient C and maximum flow enhancement (ΔE) based on the pore-to-grain size ratio and porosity. Energy dissipation increased with both heterogeneity and boundary slip, which we identified as the primary mechanism contributing to asymptotic flow enhancement. This relationship was illustrated by a 1:1 linear correlation between maximum energy dissipation and maximum flow enhancement, regardless of heterogeneity, indicating that energy dissipation due to boundary slip entirely controls the emerging fluid dynamics. The presented theoretical model and constitutive equations offer practical applications for optimizing fluid dynamics in heterogeneous formations.

Molecular dynamics simulation study of oil slippage on polymer grafted surface: The role of grafting density

AbstractThe polymer grafted materials are frequently used in antifouling and oil–water separation, while the microscopic mechanism of liquid slippage on this surface remains unclear. This study uses molecular dynamics (MD) simulations to explore the structure and transport of oil (n‐hexadecane) on polydimethylsiloxane (PDMS) grafted surface as a function of grafting density. Simulation results indicate that PDMS chain structure undergoes a transition from liquid‐like to solid‐like structure when the grafting density increases. As a result, the entanglement between the PDMS and n‐hexadecane is also weakened. The direct consequence of this structural change is that the slip length of n‐hexadecane shows a non‐monotonic dependence on the grafting density, and a maximum slip length appears at intermediate grafting density. Further analysis indicate that the slip length also increases as the applied shear rate increases, which facilitates the transport of fluid at the grafted surface. These simulation results illustrate the subtle interplay between polymer–alkane interface structure and the transport property of the fluids, which are consistent with recent experimental studies, and will be helpful for rational design of polymer grafted materials in separation technologies.

Investigation of oil/water two-phase flow behavior in laminated shale porous media considering heterogeneous structure and fluid-solid interaction

The laminated structure of continental shale oil reservoirs introduces significant macro and micro heterogeneity, posing challenges for efficient shale oil extraction. Existing microscopic experimental and simulation methods are insufficient to accurately simulate the fluid flow behavior in mud–lime interactive laminated shale reservoirs. This study delves into the complexities of laminated shale reservoirs with a focus on their heterogeneous structure and wettability. Using scanning electron microscopy, we extracted the structural features of laminated shale reservoirs and established a Multicomponent and multiphase -multi relaxation time-Lattice Boltzmann model(MCMP-MRT-LBM) that considers nanoscale effects such as wall adsorption, liquid–solid slippage, and liquid–liquid slippage. We investigated the occurrence state, the flow behavior, and the phase seepage characteristics of crude oil in laminated shale at different water saturation stages and explored the reasons behind interlayer channeling and the conditions under which crude oil is mobile within limestone and mudstone layers. Then, we examined the impact of varying factors, including heterogeneous wettability, displacement pressure gradients, and pore size, on relative permeability, water injection capacity, sweep efficiency, and crude oil recovery in laminated reservoirs. The results indicate that wettability and pore size distribution, driven by capillary and viscous forces, dictate the occurrence state of remaining oil. Wettability primarily influences fluid distribution, whereas pore size distribution substantially impacts fluid morphology. It also highlights the pivotal role of flow velocity differences between layers in causing interlayer channeling. Increasing the capillary number promotes crude oil mobility within mudstone layers, with the initiation pressure gradient in mudstone layers being 2.5 times that of the limestone layer. The examination of the oil–water relative permeability reveals the dominance of pressure gradients in affecting oil flow capabilities, while pore size significantly influences water flow capacity. By optimizing various factors, it is possible to enhance water injection capacity and sweep efficiency, resulting in a substantial 5%–10% increase in crude oil recovery from laminated shale reservoirs.

Computational Fluid Dynamics Analysis of Wet Dust Removal in High-Gravity Countercurrent Rotating Packed Bed

High-gravity wet dust removal technology has attracted much attention because of its potential to cut liquid into smaller liquid droplets and lower energy consumption. However, the complex structure and the high-speed rotation of the rotating packed bed do not allow us to analyze the flow field using conventional methods, and thus the capture mechanism of fine particles in a high-gravity environment is poorly understood. In this study, a two-dimensional computational fluid dynamics model was established to investigate the distribution of the gas–liquid two-phase flow field inside of a rotating packed bed. The characteristics of the flow field, such as the liquid form, gas–liquid contact time, and gas flow path, were investigated, and the droplet size distribution and gas–liquid slip velocity were quantitatively analyzed. The inertial capture efficiency was calculated using the Stokes number, and the dust removal efficiency distribution in the rotating packed bed was compared. The reason for the high collection efficiency of fine particles by the high-gravity wet dust removal technology was explained by numerical simulations. Two new structures were designed to improve the total dust removal efficiency.

Molecular Mechanics of Liquid and Gas Slip Flow

By taking into account the inhomogeneity introduced by the presence of a solid boundary, slip-flow theory extends the range of applicability of the venerable Navier–Stokes description to smaller scales and into the regime where confinement starts to be important. Due to the inherently atomistic nature of solid–fluid interactions at their interface, slip flow can be described, at least in principle, predictively at this level. This review aims to summarize our current understanding of slip flow at the atomistic level in dilute gases and dense liquids. The discussion extends over the similarities and differences between slip in gases and liquids, characterization and measurement of slip by molecular simulation methods, models for predicting slip, and open questions requiring further investigation.

Studying Flotation of Gold Microdispersions with Carrier Minerals and Pulp Aeration with a Steam–Air Mixture

This work is aimed at obtaining new knowledge in the field of interactions of polydisperse hydrophobic surfaces in order to increase the extraction of mineral microdispersions via flotation. The effect of high velocity and the probability of aggregating fine particles with large ones are used to increase the extraction of finely dispersed gold in this work. Large particles act as carrier minerals, which are intentionally introduced into a pulp. The novelty of this work lies in the fact that a rougher concentrate is used as the carrier mineral. For this purpose, it is isolated from three parallel pulp streams by mixing the rougher concentrate, isolated from the first stream of raw materials, with an initial feed of the second stream; accordingly, the rougher concentrate of the second stream is mixed with the initial feed of the third stream, and the finished rougher concentrate is obtained. In this mode of extracting the rougher concentrate, the content of the extracted metal increases from stream to stream, which contributes to the growth in its content in the end product. Moreover, in order to supplement forces involved in the separation of minerals with surface forces of structural origin in the third flotation stream, the pulp is aerated for a short time (about 15%–25% of the total) with air bubbles filled with a heat carrier, i.e., hot water vapor. Within this accepted flotation method, the influence that the surface currents occurring in the wetting film have on its thinning and breakthrough kinetics is proposed to be in the form of a correction to a length of a liquid slip in the hydrophobic gap. The value of the correction is expressed as a fraction of the limiting thickness of the wetting film, determined by the condition of its thickness invariability when the streams are equal in an interphase gap: outflowing (due to an action of the downforce) and inflowing (Marangoni flows and a thermo-osmotic stream). Gold flotation experiments are performed on samples of gold-bearing ore obtained from two deposits with conditions that simulate a continuous process. Technological advantages of this developed scheme and a flotation mode of gold microdispersions are shown in comparison with the basic technology. The purpose of this work is to conduct comparative tests on the basic and developed technologies using samples of gold-bearing ore obtained from the Natalka and Olimpiada deposits. Through the use of the developed technology, an increase in gold extraction of 7.99% and in concentrate quality (from 5.09 to 100.3 g/t) is achieved when the yield of the concentrate decreases from 1.86 to 1.30%, which reduces the costs associated with its expensive metallurgical processing.

Multi-component oil–water two phase flow in quartz and kerogen nanopores: A molecular dynamics study

Two-phase flow of oil and water frequently occurs within shale nanopores. Unique phenomena, like liquid–liquid slip at the oil–water interface, have been noted within the alkane-water system. However, in the case of realistic crude oil, the presence of highly polar components like asphaltene and toluene may introduce intricate migration behaviors near phase interfaces. Consequently, it is imperative to thoroughly explore the influence of these realistic crude oil components on flow behavior. This paper employs a comparative approach to examine the impact of polar components present in crude oil on two-phase nanoflow by molecular dynamic simulations. The investigation employs a mixed alkanes model (encompassing methane, propane, octane, and tetradecane) and a black oil model (comprising a range of alkanes, toluene, and asphaltene) to analyze the behavior of oil–water two-phase flow within quartz and kerogen nanopores, respectively. Four water adsorption layers near the quartz with an individual thickness of 0.23 nm and three alkane adsorption layers near the kerogen surface with an individual thickness of 0.45 nm were observed. Near the oil–water interface region, it indicated that toluene and asphaltene tend to accumulate, leading to a reduction of alkane aggregation in that region. Such aggregation of toluene and asphaltene imparts an interfacial viscosity increase of 2.67 to 3.47 times, substantially diminishing the enhancement of liquid–liquid slip. Based Hagen-Poiseuille equation, we proposed a simplified mathematical model in which the liquid–liquid slip was neglected and the oil–water interface region was omitted. Compared with the results from molecular dynamic simulations, such a simplified model is validated and effective for the description of oil–water two-phase nanoflows containing toluene and asphaltene, and ensures that the relative error of flow rate remains below 0.5 %. This study provides valuable insights into the two-phase flow of oil–water in nanopores, particularly measuring the effect of liquid–liquid slip under the different crude oil components conditions.

Dynamics of tank trucks with barters for transportation of viscous liquids

The article deals with the issues related to the safety of transportation of viscous liquids in tanks on motor transport. The special attention in the process of research is paid to the forms of partitions and their influence on the slippage of liquid and, accordingly, on the load transfer in the roll and pitch planes. On the basis of the obtained results, an optimized design of longitudinal baffles is proposed, which allows to provide the reduction of transverse force and roll moment. Keywords tank, partition, roll, viscous fluid, pitching

A FRACTAL-BASED OIL TRANSPORT MODEL WITH UNCERTAINTY REDUCTION FOR A MULTI-SCALE SHALE PORE SYSTEM

The challenges of modeling shale oil transport are numerous and include strong solid-fluid interactions, fluid rheology, the multi-scale nature of the pore structure problem, and the different pore types involved. Until now, theoretical studies have not fully considered shale oil transport mechanisms and multi-scale pore structure properties. In this study, we propose a fractal-based oil transport model with uncertainty reduction for a multi-scale shale pore system. The fractal properties of the shale pore system are obtained using high-resolution scanning electron microscope (SEM) imaging combined with laboratory core sample gas permeability measurements to reduce the model uncertainty. This fractal-based oil transport model accounts for boundary slippage, fluid rheology, the adsorption layer, and different pore types. We further pinpoint the effects of the fractal properties (pore fractal dimension, tortuosity fractal dimension), the shale pore properties (pore type, pore size, total organic carbon in volume), and the fluid properties (yield stress, liquid slippage, adsorption layer) on the shale oil permeability and mobile oil saturation using the proposed model. The results reveal that the size of the inorganic pores has the largest influence on the shale oil transport properties, followed by the yield stress, tortuosity fractal dimension, and the fractal dimension of the inorganic pores.

A rapid performance prediction method for Two-Phase liquid metal MHD generators based on Quasi-One-Dimensional model

Liquid metal magnetohydrodynamic (LMMHD) power generation is a novel and highly efficient technology with broad application prospects. In our previous study, a novel LMMHD enhanced Closed-Brayton-Cycle power generation system for hypersonic vehicles was proposed, and its excellent power generation performance was demonstrated by thermodynamic analysis. The most important component of LMMHD system is the power generation channel, and the rapid and accurate performance prediction and optimization design methods of which are the prerequisites for further research. In this study, a quasi-one-dimensional modeling method for LMMHD power generation channels based on a two-fluid model is developed. The performance variation laws with the geometric structure, magnetic field and boundary conditions, as well as the potential loss mechanisms are investigated. Results indicate that the most optimal or appropriate design can be achieved by optimizing the channel geometry and magnetic field, such as the optimal aspect ratio of 1 and the optimal Hartmann number of 2000. But it is difficult to balance the power output and efficiency through the boundary condition optimization. The gas–liquid slip can lead to approximately 20–40% of energy losses, making the elimination and suppression of slip effects a significant research focus in the relevant field.

Mixing tank hydrodynamics and mixing efficiency for propeller impellers

In this study, the distribution of axial and radial velocities with the use of a PIV laser system was determined and the mixing power in a mixing tank with narrow and wide propeller impellers was measured. The pumping capacity and secondary circulation values in the mixing tank were calculated and on this basis, the energy efficiency of the operation of both types of propeller impellers was determined and compared with the efficiency of turbine impellers with inclined blades. It turned out that propeller impellers with wide blades are most effective as regards the amount of energy that should be supplied to the mixing tank to achieve the required mixing time. It was found that the share of pumping capacity in the total secondary circulation is almost twice as high as the share of induced circulation. Correlation equations for the dependence of the mixing power and dimensionless flow numbers on the pitch of the propeller blade were also formulated. All the results obtained were compared with the works of other authors. The paper also found that propeller impellers with a variable angle of inclination of the blade α depending on the current radius do not meet the theoretical assumption of a constant axial velocity of the liquid flowing out from of the impeller zone. For small radii of the impeller the explanation for this observation is the presence of the hub and shaft in the axis of the mixer and the effect of liquid slip for the steep agitator blades near the hub. For radii close to the radius of the stirrer, the reason is generation of induced circulation as the expense of the kinetic energy of the liquid flowing out of the impeller area.

Water-in-oil emulsion drop-wall impact: Effect of viscosity and water particles on slip length over an oleophilic surface

The study is focused on the numerical simulation of the boundary slip of a water-in-oil (W/O) emulsion along an oleophilic surface in a 2D axisymmetric formulation of the Navier–Stokes equations for an incompressible laminar flow. Drops of emulsions with a volume fraction of the dispersed phase of 0.04–0.34 collide with a smooth glass wall at initial velocities of 1.2–3.6 m/s; the motion of the interface between the liquid and gas phases is resolved by the level set method. The influence of the wall wettability by oily liquids is demonstrated by varying the static contact angle up to 0.5° and by considering the dynamic one as a function of the empirically measured contact line velocity of the liquid drop according to Tanner's law. Numerical simulation of emulsion drop spreading dynamics reveals that the boundary slip of a heterogeneous liquid along a smooth oleophilic wall is different from that of a homogeneous. The temporal changes in the average viscous stress allowed the illustration of the interrelated effects of emulsion viscosity and the presence of disperse phase particles on the boundary slip of a heterogeneous liquid and its motion in the viscous boundary layer. The analytical expression evaluating the slip length for an arbitrary W/O emulsion was derived using the effective properties of the liquid. In testing the expression, the values of the maximum spreading factor of drops of an arbitrary W/O emulsion predicted by numerical simulation correspond to those determined from experiments at a moderate error of 7.4%.

Entropy Generation in a Magnetohydrodynamic Hybrid Nanofluid Flow over a Nonlinear Permeable Surface with Velocity Slip Effect

The current study is designed to model the hydrothermal feature of a hybrid nano liquid slip flows over a permeable expanding/contracting surface with entropy generation. The model incorporates Cu-Al2O3 nanoparticles with water as the host liquid to simulate the flow. Additional impacts incorporated into the novelty of the model are viscous dissipation and Joule heating. The model is transformed appropriately to its dimensionless form using similarity quantities and the solution is numerically obtained using the spectral quasi-linearization method (SQLM). The impact of pertinent factors on the flow characteristics is communicated through graphs for the hybrid nano-suspension to discuss the hydrothermal variations. The friction factor and the rate of heat transport are also discussed with sensible judgment through tables. To ensure the code validity, a comparison with earlier studies is conducted and excellent consensus is accomplished. The result explored that diminution in the irreversibility ratio is witnessed for rising magnetic field strength along the free stream, distance away from the permeable surface as the heat dissipation to the surrounding decelerates. Also, the augmented nonlinearity parameter intensified the heat transfer rate for about 2.79% of the hybrid nano-suspension.

Jeffery’s paradox for the rotation of a single ‘stick–slip’ cylinder

We consider the problem of determining the two-dimensional fluid velocity due to the rotation of an infinitely-long ‘stick–slip’ cylinder in an otherwise quiescent Stokes flow. Stick–slip boundary conditions are introduced as a model of a rough superhydrophobic surface, via a distribution of alternating solid–liquid (stick) and gas–liquid (slip) interfaces. This leads to a mixed boundary-value problem for Stokes flow. Complex variable techniques are employed to transform the flow problem into a Hilbert problem, which involves finding a function analytic in a plane region assuming that on some portions of the boundary its real part is known, while on others its imaginary part is given. We solve the Hilbert problem to obtain semi-analytic expressions for all the pertinent fluid-dynamic quantities. We find that in the general aperiodic case there is no solution in which the velocity of the fluid vanishes at infinity. This is a form of Jeffery’s paradox, typically associated with viscous flow due to the counter-rotation of two equal rigid cylinders. Our work provides the first example of Jeffery’s paradox due to the rotation of a single cylinder.

The key role of interfacial non-bonding interactions in regulating lubricant viscosity using nanoparticles

The feasibility of using nanoparticles (NPs) as lubricant viscosity modifiers has attracted considerable attention as they can increase/decrease the lubricant viscosity. However, due to the lack of theoretical guidance for interface design, the practical application of NP viscosity modifiers is limited. In this study, molecular dynamics simulations are performed to analyze the non-bonding interactions of four typical lubricating oil/NP solid–liquid interfaces with different interface polarities. A model of the mechanism through which the NPs regulate the shear viscosity of the lubricating oil is proposed based on the solid–liquid interface slip behavior, which is verified experimentally. It is found that the shear resistance of the interface is determined by the "liquid all-atomic Coulomb potential energy" (LACPE) and is positively correlated with the “potential well” values of the LACPE distributed in the tangent plane parallel to the wall. Adjusting the characteristics of the solid surface and the "polar units", which are inside the lubricant molecules, as well as the linking modes between the polar units, is proposed as a means to regulate the interfacial shear resistance. This work could provide guidance for interface design and the R&D of nanoscale viscosity modifiers.

High-Speed Centrifugal Spinning Polymer Slip Mechanism and PEO/PVA Composite Fiber Preparation.

Composite nanofibers with excellent physical and chemical properties are widely used in new energy, biomedical, environmental, electronic, and other fields. Their preparation methods have been investigated extensively by many experts. High-speed centrifugal spinning is a novel method used to fabricate composite nanofibers. The slip mechanism of polymer solution flows is an important factor affecting the morphology and quality of composite nanofibers prepared by high-speed centrifugal spinning. As the polymer solution flows, the liquid wall slip occurs inside the nozzle, followed by liquid-liquid interface slip and gas-liquid interface slip. The factors affecting polymer slip were investigated by developing a mathematical model in the nozzle. This suggests that the magnitude of the velocity is an important factor that affects polymer slip and determines fiber quality and morphology. Under the same rotational speed, the smaller the nozzle diameter, the greater the concentration of velocity distribution and the smaller the diameter of the produced composite nanofibers. Finally, PEO/PVA composite nanofibers were prepared using high-speed centrifugal spinning equipment at 900-5000 rpm and nozzle diameters of 0.2 mm, 0.4 mm, 0.6 mm, and 0.8 mm. The morphology and quality of the collected PEO/PVA composite nanofibers were analyzed using scanning electron microscopy (SEM) and TG experiments. Then, the optimal parameters for the preparation of PEO/PVA composite nanofibers by high-speed centrifugal spinning were obtained by combining the external environmental factors in the preparation process. Theoretical evaluation and experimental data were provided for the centrifugal composite spinning slip mechanism and for the preparation of composite nanofibers.

Dynamics of co-current gas–liquid film flow through a slippery channel

We consider a thin liquid film in a wide inclined channel being driven by gravity and co-current turbulent gas flow. The bottom plate with which the liquid is in contact with is taken to be slippery, and we impose the classic Navier slip condition at this substrate. Such a setting finds application in technological processes as well as nature (e.g., distillation, absorption, and cooling devices). The gas–liquid problem can be decoupled by making certain reasonable assumptions. Under these assumptions, we solve the gas problem to obtain the tangential and normal stresses acting at the wavy gas–liquid interface for arbitrary waviness. In modeling the liquid layer dynamics, we make use of the stresses computed in the gas problem as inputs to the interface boundary conditions. We develop the long-wave model and the weighted-integral boundary layer (WIBL) model to describe the thin film dynamics. We perform a linear stability of these reduced order models to scrutinize the effect of wall slip, liquid flow rate, and the gas shear on the stability of the flat film solution. It is found that the wall slip promotes the instability of the flat interface. Furthermore, we compute solitary wave solutions of the WIBL model by implementing Keller's pseudo-arc length algorithm on a periodic domain. We observe that the wave speed as well as the wave amplitude are attenuated on incrementing the liquid slip at the substrate. We corroborate these findings with the time-dependent computations of the nonlinear WIBL model.

Axisymmetric lattice Boltzmann model for liquid flows with super-hydrophobic cylindrical surfaces

The lattice Boltzmann (LB) method has been extensively applied to the simulation of various fluid flows. Nonetheless, the standard LB model in rectangular coordinates faces some challenges when attempts are made to use it to simulate single-phase liquid flows involving super-hydrophobic cylindrical surfaces with a large slip length. Therefore, in this paper, boundary conditions for an axisymmetric LB model are proposed for the simulation of liquid slip at both convex and concave cylindrical surfaces. Based on the Navier slip boundary conditions, an equilibrium–nonequilibrium boundary scheme is developed to obtain the liquid slip in the azimuthal direction, with a combined bounce-back and specular-reflection boundary condition being developed to obtain the liquid slip in the axial direction. These boundary schemes are verified by simulating several established cases, namely, cylindrical Couette flow, harmonic oscillatory cylindrical Couette flow, Hagen–Poiseuille flow, and annular Poiseuille flow. The results are in excellent agreement with analytical results. Furthermore, some more complex flows involving super-hydrophobic cylindrical surfaces, such as stepwise oscillatory cylindrical Couette flow, are also simulated and studied using the proposed LB model. The axisymmetric LB model presented here overcomes the inability of standard LB models to accurately and efficiently simulate single-phase liquid flows involving super-hydrophobic cylindrical surfaces with a large slip length.