Research Papers: Fundamental Issues and Canonical Flows

Modeling of Hydrophobic Surfaces by the Stokes Problem With the Stick–Slip Boundary Conditions

[+] Author and Article Information
R. Kučera

VŠB-TU Ostrava,
17 listopadu 15/2172,
Ostrava-Poruba 708 33, Czech Republic
e-mail: radek.kucera@vsb.cz

V. Šátek

VŠB-TU Ostrava,
17 listopadu 15/2172,
Ostrava-Poruba 708 33, Czech Republic
e-mail: vaclav.satek@vsb.cz

J. Haslinger

Studentská 1768,
Ostrava-Poruba 708 00, Czech Republic
e-mail: hasling@karlin.mff.cuni.cz

S. Fialová

Victor Kaplan Department of Fluid Engineering,
Brno University of Technology,
Technická 2896/2,
Brno 616 69, Czech Republic
e-mail: fialova@fme.vutbr.cz

F. Pochylý

Victor Kaplan Department of Fluid Engineering,
Brno University of Technology,
Technická 2896/2,
Brno 616 69, Czech Republic
e-mail: pochyly@fme.vutbr.cz

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received February 3, 2016; final manuscript received July 12, 2016; published online October 10, 2016. Assoc. Editor: Matevz Dular.

J. Fluids Eng 139(1), 011202 (Oct 10, 2016) (9 pages) Paper No: FE-16-1073; doi: 10.1115/1.4034199 History: Received February 03, 2016; Revised July 12, 2016

Unlike the Navier boundary condition, this paper deals with the case when the slip of a fluid along the wall may occur only when the shear stress attains certain bound which is given a priori and does not depend on the solution itself. The mathematical model of the velocity–pressure formulation with this type of threshold slip boundary condition is given by the so-called variational inequality of the second kind. For its discretization, we use P1-bubble/P1 mixed finite elements. The resulting algebraic problem leads to the minimization of a nondifferentiable energy function subject to linear equality constraints representing the discrete impermeability and incompressibility condition. To release the former one and to regularize the nonsmooth term characterizing the stick–slip behavior of the algebraic formulation, two additional vectors of Lagrange multipliers are introduced. Further, the velocity vector is eliminated, and the resulting minimization problem for a quadratic function depending on the dual variables (the discrete pressure and the normal and shear stress) is solved by the interior point type method which is briefly described. To justify the threshold model and to illustrate the efficiency of the proposed approach, three physically realistic problems are solved and the results are compared with the ones solving the Stokes problem with the Navier boundary condition.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Fialová, S. , 2016, “ Identification of the Properties of Hydrophobic Layers and its Usage in Technical Practice,” Habilitation, VUTIUM, Brno University of Technology, Brno, the Czech Republic.
Fialová, S. , and Pochylý, F. , 2014, “ Identification and Experimental Verification of the Adhesive Coefficient of Hydrophobic Materials,” Wasserwirtsch., 105(1), pp. 125–129.
Haslinger, J. , Hlaváček, I. , and Nečas, J. , 1996, “ Numerical Methods for Unilateral Problems in Solid Mechanics,” Handbook of Numerical Analysis, Vol. IV, North Holland, Amsterdam, Part 2, pp. 313–485.
Anitescu, M. , and Potra, F. A. , 1997, “ Formulating Dynamic Multi-Rigid-Body Contact Problems With Friction as Solvable Linear Complementarity Problems,” Nonlinear Dyn., 14(3), pp. 231–247. [CrossRef]
Kučera, R. , Haslinger, J. , Šátek, V. , and Jarošová, M. , 2016, “ Efficient Methods for Solving the Stokes Problem With Slip Boundary Conditions,” Math. Comput. Simul. (in press).
Kučera, R. , Machalová, J. , Netuka, H. , and Ženčák, P. , 2013, “ An Interior Point Algorithm for the Minimization Arising From 3D Contact Problems With Friction,” Optim. Method Software, 28(6), pp. 1195–1217. [CrossRef]
Hammad, K. J. , 2015, “ The Flow Behavior of a Biofluid in a Separated and Reattached Flow Region,” ASME J. Fluids Eng., 137(6), p. 061104. [CrossRef]
Corredor, F. E. R. , Bizhani, M. , and Kuru, E. , 2015, “ Experimental Investigation of Drag Reducing Fluid Flow in Annular Geometry Using Particle Image Velocimetry Technique,” ASME J. Fluids Eng., 137(8), p. 081103. [CrossRef]
Ozogul, H. , Jay, P. , and Magnin, A. , 2015, “ Slipping of the Viscoplastic Fluid Flowing on a Circular Cylinder,” ASME J. Fluids Eng., 137(7), p. 071201.
ANSYS, 2015, “fluent,” ANSYS, Inc., Canonsburg, PA, accessed Nov. 15, 2015, http://www.ansys.com/Products/Fluids/ansys+Fluent
Brdička, M. , Samek, L. , and Bruno, S. , 2011, Mechanika kontinua, Academia, Prague, Czech Republic.
Elman, H. C. , Silvester, D. J. , and Wathen, A. J. , 2005, Finite Elements and Fast Iterative Solvers With Applications in Incompressible Fluid Dynamics, Oxford University Press, Oxford, UK.
Nečas, J. , 1967, Les Méthodes Directes en Théorie des Equations Elliptiques, Masson, Paris, France.
Fujita, H. , 1994, “ A Mathematical Analysis of Motions of Viscous Incompressible Fluid Under Leak and Slip Boundary Conditions,” RIMS Kokyuroku, 888, pp. 199–216.
Bulíček, M. , and Málek, J. , 2016, “ On Unsteady Internal Flows of Bingham Fluids Subject to Threshold Slip on the Impermeable Boundary,” Recent Developments of Mathematical Fluid Mechanics, H. Amann , Y. Giga , H. Kozono , H. Okamoto , and M. Yamazaki , eds., Springer, Heidelberg, Germany, pp. 135–156.
Ayadi, M. , Baffico, L. , Gdoura, M. K. , and Sassi, T. , 2014, “ Error Estimates for Stokes Problem With Tresca Friction Conditions,” ESAIM: Math. Modell. Numer. Anal., 48(5), pp. 1413–1429. [CrossRef]
Koko, J. , 2012, “ Vectorized matlab Codes for the Stokes Problem With P1-Bubble/P1 Finite Element,” http://www.isima.fr/~jkoko/Codes/StokesP1BubbleP1.pdf, accessed July 28, 2016, ISIMA Public Engineering School in Computer Science, Auvergne, France.
Jarošová, M. , Kučera, R. , and Šátek, V. , 2015. “ A New Variant of the Path-Following Algorithm for the Parallel Solving of the Stokes Problem With Friction,” Fourth International Conference on Parallel, Distributed, Grid and Cloud Computing for Engineering, P. Iványi and B. H. V. Topping, eds., Civil-Comp Press, Stirlingshire, Scotland, UK, Paper No. 11.
Nocedal, J. , Wächter, A. , and Waltz, R. A. , 2005, “ Adaptive Barrier Strategies for Nonlinear Interior Methods,” Report No. TR RC 23563, IBM T.J. Watson Research Center, Yorktown Heights, NY.
Pochylý, F. , Fialová, S. , and Kozubková, M. , 2011, “ Journal Bearings With Hydrophobic Surface,” Vibronadežnos i Germetičnos Centrobežnyh Mašin, Technical Study—Monography, pp. 314–320.
Pochylý, F. , Fialová, S. , and Malenovský, E. , 2012, “ Bearing With Magnetic Fluid and Hydrophobic Surface of the Lining,” IOP Conf. Ser.: Earth Environ. Sci., 15(2), pp. 1–9.
MATLAB, 2016, “MATLAB: The Language of Technical Computing,” http://www.mathworks.com/products/matlab/, The Mathworks, Inc., Natick, MA.
Hron, J. , Roux, C. L. , Málek, J. , and Rajagopal, K. , 2015, “ Flows of Incompressible Fluids Subject to Navier’s Slip on the Boundary,” Comput. Math. Appl., 56(8), pp. 2128–2143. [CrossRef]
Koko, J. , 2015, “ A matlab Mesh Generator for the Two-Dimensional Finite Element Method,” Appl. Math. Comput., 250, pp. 650–664. [CrossRef]


Grahic Jump Location
Fig. 1

Geometry of the channel Ω with diameter d = 0.015 (m) and length L = 1 (m)

Grahic Jump Location
Fig. 2

(a) Comparison of κnum and κs at xc and (b) dependence of |ut(xc)| on the slip bound g

Grahic Jump Location
Fig. 3

Outflow velocity profiles: (a) the same value of g on both parts of γC and (b) different values of g:=g1 on γC,1 and fixed value g = 1 on γC,2 lead to the sticking effect on γC,2. The values κs:=κs(xc) are computed by Eqs. (2) and (23).

Grahic Jump Location
Fig. 4

Dependence of the hydraulic losses on the adhesive coefficient κnum

Grahic Jump Location
Fig. 5

(a) Rotating cylinder γD and the stick–slip boundary conditions on γC and (b) comparison of κs and κnum on γC

Grahic Jump Location
Fig. 6

Dependence of ut on (a) the slip bound g and (b) the adhesive coefficient κnum

Grahic Jump Location
Fig. 7

Example: (a) eccentric rotating cylinder and (b) zoom of the mesh

Grahic Jump Location
Fig. 8

Distribution of (a) adhesive coefficient κs and (b) tangential velocity ut on γC for g = 0.18094

Grahic Jump Location
Fig. 9

Distribution of (a) shear stress |σt| and (b) tangential velocity ut on γC for g = 1.60313

Grahic Jump Location
Fig. 10

The distribution of |σt| on γC for g = 50

Grahic Jump Location
Fig. 11

Zoomed velocity fields u

Grahic Jump Location
Fig. 12

Distribution of the pressure p on γC




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In