Numerical simulations of spontaneous capillary rises with very low capillary numbers using a front-tracking method combined with generalized Navier boundary condition

We propose a front-tracking method that considers a moving contact line with a generalized Navier boundary condition (GNBC) and with a delta function distribution approximated on the grid scale. A method of evaluating the interfacial balance at the contact line in the form of a body force, which is straightforward with the front-tracking method, cannot give a natural flow field. In contrast, the proposed method using the GNBC, which includes the unbalanced Young’s force as stress on the wall, can give a very stable and reasonable flow field. The proposed front-tracking method was applied for the capillary rise of a liquid in a tube, in which the velocity-dependent contact angle dominates the dynamic characteristics. The validity of the proposed method was confirmed by comparing simulation results with experimental measurements and simple theoretical models. The results of the present simulations with adjusted non-dimensional slip parameters agreed very well with experimental measurements. Under the present simulation conditions, the linearity of the GNBC allows the correlation between the dynamic contact angle and the contact line’s velocity to follow a simple linear expression that involves the difference of the cosine with the capillary number. The non-dimensional slip parameter, which represents the dynamic nature of the moving contact line, can therefore be easily adjusted to reproduce experimental observations under small-capillary-number conditions.     

Numerical simulation of laminar flow past a circular cylinder with slip conditions

The no‐slip condition is an assumption that cannot be derived from first principles and a growing number of literatures replace the no‐slip condition with partial‐slip condition, or Navier‐slip condition. In this study, the influence of partial‐slip boundary conditions on the laminar flow properties past a circular cylinder was examined. Shallow‐water equations are solved by using the finite element method accommodating SU/PG scheme. Four Reynolds numbers (20, 40, 80, and 100) and six slip lengths were considered in the numerical simulation to investigate the effects of slip length and Reynolds number on characteristic parameters such as wall vorticity, drag coefficient, separation angle, wake length, velocity distributions on and behind the cylinder, lift coefficient, and Strouhal number. The simulation results revealed that as the slip length increases, the drag coefficient decreases since the frictional component of drag is reduced, and the shear layer developed along the cylinder surface tends to push the separation point away toward the rear stagnation point so that it has larger separation angle than that of the no‐slip condition. The length of the wake bubble zone was shortened by the combined effects of the reduced wall vorticity and wall shear stress which caused a shift of the reattachment point closer to the cylinder. The frequency of the asymmetrical vortex formation with partial slip velocity was increased due to the intrinsic inertial effect of the Navier‐slip condition. Copyright © 2011 John Wiley & Sons, Ltd.     

微液滴在PDMS软基体表面的动态摩擦行为研究

通过固液界面摩擦力测试装置研究了微液滴在PDMS软基体表面运动时的动态摩擦学行为,并对微液滴体积、滑动速度及软基体力学性能对固液界面动态摩擦行为的影响进行了分析. 结果表明:微液滴在软基体表面运动时表现出最大静摩擦力和动态摩擦力. 最大静摩擦力与微液滴黏度和速度梯度呈正比,动态摩擦力与微液滴体积、滑动速度和基体力学性能有关. 随着微液滴体积的增加,三相接触线长度增加,动态摩擦力增加;随着相对滑动速度增加,三相接触线长度及接触角滞后增加,动态摩擦力增加;随着软基体弹性模量降低,固液界面黏附力增加,固液界面运动能量耗散增加,动态摩擦力增加. 研究结果可为PDMS软基体表面微液滴的精确驱动和运动参数优化提供理论指导,也可进一步丰富固液界面摩擦理论.      

On the regularity of some thermocapillary convection models

The classical model of confined thermocapillary convection is analyzed. Its vorticity singularity, independent of the contact angle, leads to infinite pressure values at contact lines, forbidding any numerical use of the Laplace equation to calculate free surface shapes. Four models are explored to overcome this difficulty: an explicit polynomial filtering, a Navier slip at the solid boundaries, an interface viscosity model and the combination of slip and interface viscosity. Regular solutions are obtained with the first and last approaches. Only the last one is based on physical considerations and, by the introduction of physical length scales, avoids infinite pressure values at the contact line.     

Finite element methods for a class of continuum models for immiscible flows with moving contact lines

In this paper, we present a finite element method for two‐phase incompressible flows with moving contact lines. We use a sharp interface Navier–Stokes model for the bulk phase fluid dynamics. Surface tension forces, including Marangoni forces and viscous interfacial effects, are modeled. For describing the moving contact lines, we consider a class of continuum models that contains several special cases known from the literature. For the whole model, describing bulk fluid dynamics, surface tension forces, and contact line forces, we derive a variational formulation and a corresponding energy estimate. For handling the evolving interface numerically, the level‐set technique is applied. The discontinuous pressure is accurately approximated by using a stabilized extended finite element space. We apply a Nitsche technique to weakly impose the Navier slip conditions on the solid wall. A unified approach for discretization of the (different types of) surface tension forces and contact line forces is introduced. Results of numerical experiments are presented, which illustrate the performance of the solver. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis on creeping channel flows of compressible fluids subject to wall slip

Creeping channel flows of compressible fluids subject to wall slip are widely encountered in industries. This paper analyzes such flows driven by pressure in planar as well as circular channels. The analysis elucidates unsteady flows of Newtonian fluids subject to the Navier slip condition, followed by steady flows of viscoplastic fluids, in particular, Herschel–Bulkley fluids and their simplifications including power law and Newtonian fluids, that slip at wall with a constant coefficient or a coefficient inversely proportional to pressure. Under the lubrication assumption, analytical solutions are derived, validated, and discussed over a wide range of parameters. Analysis based on the derived solutions indicates that unsteadiness alters cross-section velocity profiles. It is demonstrated that compressibility of the fluids gives rise to a concave pressure distribution in the longitudinal direction, whereas wall slip with a slip coefficient that is inversely proportional to pressure leads to a convex pressure distribution. Energy dissipation resulting from slippage can be a significant portion in the total dissipation of such a flow. A distinctive feature of the flow is that, in case of the pressure-dependent slip coefficient, the slip velocity increases rapidly in the flow direction and the flow can evolve into a pure plug flow at the exit.     

The wetting problem of fluids on solid surfaces. Part 2: the contact angle hysteresis

In part 1 (Gouin, [13]), we proposed a model of dynamics of wetting for slow movements near a contact line formed at the interface of two immiscible fluids and a solid when viscous dissipation remains bounded. The contact line is not a material line and a Young-Dupré equation for the apparent dynamic contact angle taking into account the line celerity was proposed. In this paper we consider a form of the interfacial energy of a solid surface in which many small oscillations are superposed on a slowly varying function. For a capillary tube, a scaling analysis of the microscopic law associated with the Young-Dupré dynamic equation yields a macroscopic equation for the motion of the contact line. The value of the deduced apparent dynamic contact angle yields for the average response of the line motion a phenomenon akin to the stick-slip motion of the contact line on the solid wall. The contact angle hysteresis phenomenon and the modelling of experimentally well-known results expressing the dependence of the apparent dynamic contact angle on the celerity of the line are obtained. Furthermore, a qualitative explanation of the maximum speed of wetting (and dewetting) can be given.Received: 5 June 2001, Accepted: 24 May 2003, Published online: 29 July 2003PACS: 02.90, 47.50, 66.20, 68.03, 68.08     

弹塑性微凸体侧向接触相互作用能耗

传统的结合面研究多基于光滑刚性平面与等效粗糙表面接触假设,忽略了结合面上微凸体侧向接触及相邻微凸体之间的相互作用,这导致理论模型与实际结合面存在较大出入.针对承受法向静、动态力的机械结合面,从微观上研究了微凸体侧向接触及相互作用的接触能耗.将法向静、动态力分解为法向分力和切向分力,获取弹性/弹塑性/塑性阶段考虑微凸体侧接触及相互作用的加、卸载法向分力-变形和切向分力-位移的关系.通过力的合成定理,从而获取加、卸载法向合力与总变形之间的关系,由于法向分力产生的塑性变形及切向分力产生的摩擦,导致加载、卸载法向合力-总变形曲线存在迟滞回线.通过对一个加、卸载周期内的法向合力-总变形曲线积分,获得一个周期的微凸体接触能耗,包括应变能耗及摩擦能耗.仿真分析表明:微凸体在3个阶段的能耗均随变形的增大而非线性增大.微凸体侧向接触角度越大,能耗越大,且在弹性阶段最为明显.在弹性阶段,仅存在侧向的摩擦能耗,故结合面在低载荷作用下必须采用双粗糙表面假设.在塑性阶段,由于微凸体接触能耗为应变能耗,且接触角对其能耗影响甚微,故结合面在大载荷作用下可采用单平面假设对其进行研究.相对于KE和Etsion模型,本文提出的模型与Bartier的实验结果更吻合.     

Time‐dependent flow across a step: the slip with friction boundary condition

The paper studies numerically the slip with friction boundary condition in the time‐dependent incompressible Navier–Stokes equations. Numerical tests on two‐ and three‐dimensional channel flows across a step using this boundary condition on the bottom wall are performed. The influence of the friction parameter on the flow field is studied and the results are explained according to the physics of the flow. Due to the stretching and tilting of vortices, the three‐dimensional results differ in many respects from the two‐dimensional ones. Copyright © 2005 John Wiley & Sons, Ltd.     

Transition from stick to slip in Hertzian contact with “Griffith” friction: The Cattaneo–Mindlin problem revisited

Classically, the transition from stick to slip is modelled with Amonton–Coulomb law, leading to the Cattaneo–Mindlin problem, which is amenable to quite general solutions using the idea of superposing normal contact pressure distributions – in particular superposing the full sliding component of shear with a corrective distribution in the stick region. However, faults model in geophysics and recent high-speed measurements of the real contact area and the strain fields in dry (nominally flat) rough interfaces at macroscopic but laboratory scale, all suggest that the transition from ‘static’ to ‘dynamic’ friction can be described, rather than by Coulomb law, by classical fracture mechanics singular solutions of shear cracks. Here, we introduce an ‘adhesive’ model for friction in a Hertzian spherical contact, maintaining the Hertzian solution for the normal pressures, but where the inception of slip is given by a Griffith condition. In the slip region, the standard Coulomb law continues to hold. This leads to a very simple solution for the Cattaneo–Mindlin problem, in which the “corrective” solution in the stick area is in fact similar to the mode II equivalent of a JKR singular solution for adhesive contact. The model departs from the standard Cattaneo–Mindlin solution, showing an increased size of the stick zone relative to the contact area, and a sudden transition to slip when the stick region reaches a critical size (the equivalent of the pull-off contact size of the JKR solution). The apparent static friction coefficient before sliding can be much higher than the sliding friction coefficient and, for a given friction fracture “energy”, the process results in size and normal load dependence of the apparent static friction coefficient. Some qualitative agreement with Fineberg's group experiments for friction exists, namely the stick–slip boundary quasi-static prediction may correspond to the arrest of their slip “precursors”, and the rapid collapse to global sliding when the precursors arrest front has reached about half the interface may correspond to the reach of the “critical” size for the stick zone.     

水工圆形隧洞围岩衬砌摩擦滑动接触的新解法

在实际工程中,围岩和衬砌接触时,它们之间并非完全光滑,也并非可以承受任意大的摩擦力.如果围岩与衬砌之间的剪应力大于所能承受的最大静摩擦力,接触面间将发生切向滑动,定义接触面上产生最小滑动量的状态为衬砌的真实工作状态,这种接触即为摩擦滑动接触.以库仑摩擦模型模拟围岩和衬砌之间的摩擦滑动接触,在考虑支护滞后效应的前提下,利用平面弹性复变函数方法列出了应力边界条件、应力连续条件以及位移连续条件的方程,再结合最优化理论,建立了具有一般性的摩擦滑动接触解法.在利用混合罚函数法求解最优化问题的过程中,减少了设计变量的个数,极大地简化了优化模型,提升了优化过程的迭代速度以及优化结果的精度.以此为基础,获得了围岩和衬砌相互作用下圆形水工隧洞的应力解析解.该方法可以求解光滑接触和完全接触两种极限情况,具有一般性.同时,利用一种精确的计算方法得到了不同情况下满足完全接触条件摩擦系数的阈值,还分析了衬砌和围岩边界上切向应力的变化规律.     

Implementation of physical boundary conditions into computational domain in modelling of oscillatory bottom boundary layers

This paper discusses the importance of realistic implementation of the physical boundary conditions into computational domain for the simulation of the oscillatory turbulent boundary layer flow over smooth and rough flat beds. A mathematical model composed of the Reynolds averaged Navier–Stokes equation, turbulent kinetic energy (k) and dissipation rate of the turbulent kinetic energy (ε) has been developed. Control‐volume approach is used to discretize the governing equations to facilitate the numerical solution. Non‐slip condition is imposed on the bottom surface, and irrotational main flow properties are applied to the upper boundary. The turbulent kinetic energy is zero at the bottom, whereas the dissipation rate is approaching to a constant value, which is proportional to the kinematic viscosity times the second derivative of the turbulent kinetic energy. The output of the model is compared with the available experimental studies conducted in oscillatory tunnels and wave flume. It is observed that the irrotational flow assumption at the upper boundary is not realistic in case of water tunnels. Therefore, new upper boundary conditions are proposed for oscillatory tunnels. The data of wave flume show good agreement with the proposed numerical model. Additionally, several factors such as grid aspect ratio, staggered grid arrangement, time‐marching scheme and convergence criteria that are important to obtain a robust, realistic and stable code are discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

On flows of third-grade fluids with non-linear slip boundary conditions

We prove the existence and uniqueness of steady flows of incompressible fluids of grade three subject to slip and no-slip boundary conditions in bounded domains. The slip boundary condition is a non-linear generalization of the Navier slip boundary condition and permits situations in which the solid boundary undergoes non-rigid tangential motion. The existence proof is based on a fixed point method in which the boundary-value problem is decomposed into four linear problems.     

摩擦热对UHMWPE/钢摩擦副摩擦性能的影响

通过控制销盘试验中的散热条件测定了摩擦热对UHMWPE/钢摩擦副摩擦性能的影响;利用统计学相关系数理论对摩擦系数和表层温度的相关性进行定量比较,发现不良散热条件使其相关系数显著增大;根据Taylor Hobson轮廓仪测量试样表面形貌三维初始表面高度,尝试用元胞自动机方法进行了摩擦表面温度分布模拟计算,并结合磨损表面形貌分析探讨摩擦热对UHMWPE磨损机理的影响.结果表明:摩擦系统的热传导条件严重影响UHMWPE的摩擦性能;摩擦热效应使得局部接触表面温度超过熔点,摩擦热的积累效应使UHMWPE的磨损机理由粘着磨损变为熔融剪切磨损;在聚合物材料的摩擦试验及应用研究中必须考虑摩擦热的影响.     

The wetting problem of fluids on solid surfaces. Part 1: the dynamics of contact lines

The understanding of the spreading of liquids on solid surfaces is an important challenge for contemporary physics. Today, the motion of the contact line formed at the intersection of two immiscible fluids and a solid is still subject to dispute. In this paper, a new picture of the dynamics of wetting is offered through an example of non-Newtonian slow liquid movements. The kinematics of liquids at the contact line and equations of motion are revisited. Adherence conditions are required except at the contact line. Consequently, for each fluid, the velocity field is multivalued at the contact line and generates an equivalent concept of line friction but stresses and viscous dissipation remain bounded. A Young-Dupré equation for the apparent dynamic contact angle between the interface and solid surface depending on the movements of the fluid near the contact line is proposed.Received: 5 June 2001, Accepted: 24 May 2003, Published online: 29 July 2003PACS: 47.17., 47.50, 66.20, 68.08, 83.50     

Boundary conditions for contacts lines in coextrusion flows

A bicomponent coextrusion process is modelled using a 3-D finite element formulation. The layer uniformity problem in coextrusion is addressed by examining the effects of the polymer melt/polymer melt/die wall contact line boundary condition. It has been observed that the less viscous polymer layer will tend to displace the more viscous polymer layer near the die wall. The behaviour of the contact lisle is considered to be either a stick or slip boundary condition. In the stick boundary condition, the contact line does not move from its original position after the two polymer layers meet, A slip boundary condition allows the contact line to move along the die wall. The calculated interfaces which result from different contact line assumptions are determined. Results show that if a stick boundary condition is appropriate for a given fluid/fluid/solid contact line, then a very thin entrained layer of the more viscous polymer melt will be trapped between the less viscous polymer melt and the die wall. Slip boundary conditions would allow complete displacement of the contact line along the die wall. Both slip and stick boundary conditions produce similar interface profiles far away from the die wall for small viscosity ratios. In certain eases, the displacement of the more viscous material by the less viscous material will cease and a static interface structure is produced regardless of die length. Experimental work with polycarbonate melts is compared with the numerical simulations.A. Torres on leave from Investigación y Desarrollo,, C.A. (INDESCA), P.O. Box 10319, Complejo Petroquímico El Tablazo, Maracaibo, 4001, Venezuela.     

一种模拟动边界绕流的锐利界面浸入边界法

When the structural wall moves over a fixed grid, the structure coverage will change, resulting in many dead and emerging elements. To avoid the influence of malformation and reconstruction of body-fitted grids on the calculation efficiency and accuracy of the fluid-structure interaction problems with coupled boundary movement on the fixed grid, an improved numerical method for describing the interaction between an immersed rigid body and fluid based on a sharp-interface is proposed. In this method, both the fluid and solid are regarded as pure fluid domains in the whole computational domain, and the solid boundary is divided into several Lagrangian grid points. The flow parameter or velocity is reconstructed by interpolation at the interface element, which is then directly used as the boundary condition of the flow field, thus reflecting the influence of the wall boundary conditions. The method constructs the calculation structure of “virtual point, force point and vertical foot point”, and the velocity of the virtual point is obtained by bilinear interpolation. Then, the velocity of the force point is calculated by forcing the solid boundary to meet the no-slip condition, and the equations of the coupling system based on the immersion boundary method are finally solved to realize the numerical simulation of the flow with a complex moving boundary. The numerical program for this immersed boundary method is established using C++, then the accuracy and reliability of the proposed method are validated by comparison with the literature and experimental results of the basic numerical example of flow around a cylinder. Furthermore, the effects of the structural shape and the angle of attack on the trailing vortex structure, the vortex shedding frequency, and the lift/ coefficient characteristics of the flow around the elliptical cylinder have been analyzed. The anti-symmetric S-type, “P+S” Ⅰ-type and “P+S” Ⅱ-type trailing vortex shedding modes, as well as the variation laws of the vortex structure size, vortex shedding frequency and lift-drag coefficients ratio with axis ratio and angle of attack, are captured. The critical angle of attack (25°) corresponding to the maximum lift-drag ratio is determined as 25°.     

DNS of turbulent channel flow with conjugate heat transfer: Effect of thermal boundary conditions on the second moments and budgets

Budgets of turbulent heat fluxes and temperature variance obtained from the Direct Numerical Simulation of an incompressible periodic channel flow with a Reynolds number of 150 (based on friction velocity) and a Prandtl number of 0.71 are presented and analysed for four cases: locally imposed temperature at the wall (constant Dirichlet), locally imposed heat flux (constant Neumann), heat exchange coefficient (Robin) and 3D conjugate heat transfer. The dissipation rate associated with the temperature variance is strongly impacted by the thermal boundary condition. For non-conjugate cases, a straightforward analytical analysis establishes the connection between the boundary condition, the temperature variance and the wall-normal part of the thermal dissipation rate at the wall. For the conjugate case, the two-point correlations of the thermal field in the solid domain confirms the existence of very large scale thermal structures.     

Wave induced sliding at the interface between a layered elastic medium and a half-space

The sliding interface between an unrestrained elastic half-space and a grounded layered half-space excited by an incident harmonic wave is investigated. The problem is formulated considering various possible boundary conditions and boundary inequalities at the sliding interface. The Coulomb friction model without distinction between the static and kinetic coefficients of friction is considered to govern the sliding condition. Three possible bands at the interface, namely slip, stick, and separation, are considered. The interface is assumed to be preloaded under normal and shear stresses. The solution is developed by modifying the problem of welded interface, which then is reduced to a set of algebraic equations. The effects of the incident angle, layer thickness, friction coefficient and externally applied stresses on the drift velocity of the unrestrained half-space are studied numerically for a pair of materials. It is shown that the sliding interface, and hence the drift velocity of unrestrained half-space is noticeably influenced by the layered medium. These results are expected to be useful for the development of a new kind of ultrasonic drive in future.     

Steering forces on undriven, angled wheels

The steering forces at low speed and zero camber angle were measured on undriven, angled wheels using tyres with no tread. The forces were measured in a soil bin using a moist loam soil at different levels of compaction. It was found that the coefficient of side force relative to the wheel was related to slip angle by an exponential relationship. Coefficient of rolling resistance relative to the wheel was a linear function of slip angle in the region zero to 20° but was an irregular function of slip angle at higher angles. The effects of tyre size, load, inflation pressure and soil condition were modelled well using different versions of the tyre mobility number. The most successful version of mobility number was one which incorporated both soil cohesion and internal friction angle. The coefficients of the exponential and linear relationships mentioned above were predicted with varying degrees of success using mobility number.