非线性边界滑移挤压膜流动

用一种包含初始滑移长度和临界剪切率的非线性边界滑移模型研究了两个球体间的挤压流体膜问题．研究发现初始滑移长度对低剪切率下的滑移行为起主要作用,而临界剪切率决定了高剪切率下的边界滑移程度．球体表面挤压流体膜的边界滑移量是与半径坐标相关的高度非线性函数．在挤压膜的中心点和远离中心点处由于低剪切率滑移量等于初始滑移长度,然而在高剪切率区域滑移长度迅速增加．球体挤压膜的流体动压力随着初始滑移长度的增加和临界剪切率的减小而减小,并且临界剪切率对流体动力的影响要比初始滑移长度大的多,当临界剪切率很小的情况下,流体动压随着最小膜厚的减小几乎不再增加．所用模型给出的理论预报和实验非常吻合．     

考虑摩擦接触时简谐波透过微缝隙的传播

研究了简谐波透过微缝隙的传播问题．考虑界面摩擦接触,利用Fourier分析技术和半空间解修正法,将边界非线性波动问题简化为一组代数方程并进行了详细的数值计算．给出了分离区、滑移区和粘着区的分布以及界面应力分布和能量分配;讨论了缝隙宽度、摩擦因数、入射角等因素对波动传输的影响．结果表明:由于界面发生局部接触和滑移使得反射波和透射波中含有高频谐波．缝隙宽度与各阶谐波反射、透射系数的关系可为材料及结构的无损检测提供理论依据．     

温差影响下的局部滑移接触行为的研究

针对不同温度装配件间接触界面的局部滑移问题,建立了三维稳态热弹性局部滑移接触的半解析求解模型.基于热弹性理论与热传导方程,构建了半空间受热流载荷和力载荷作用下的频响函数并建立了相应的影响系数矩阵.借助离散卷积-快速Fourier变换等数学工具,实现了针对高温压头与热弹性半空间局部滑移接触问题的高效求解.接触界面间的热量传递满足Fourier热传导定律,并且黏/滑状态由Coulomb定律确定.基于该半解析模型分析了不同荷载及温差对表面法向压力分布、摩擦力分布、刚体位移及接触区黏/滑演化行为的影响.研究结果表明,当法向荷载和切向荷载一定时,温差的上升会导致接触区域的减小,引起接触面法向压力及摩擦力的峰值增大,并且会显著影响黏着区与滑移区的分布情况.     

压电体光滑接触界面有局部分离时的滑移脉冲波传播

应用关于压电材料的Stroh方法以及Fourier分析和奇异积分方程技术,研究了压电体光滑接触界面有局部分离时的滑移脉冲波传播问题,分析了其存在的判据．压电体由单向压应力作用而光滑接触并处于一定强度的电场中．待求问题最终转化为含Cauchy核的奇异积分方程组,并给出其解析解．数值计算结果表明:光滑接触界面上有局部分离的滑移脉冲波普遍存在,且界面法向面力和法向电位移在局部分离区两端有奇异性;对某些特殊的材料组合,外载荷不影响滑移脉冲波的存在性．     

计及Hall和离子滑移电流的旋转流体在与时间相关初值条件下的MHDCouette流动

考虑Hall和离子滑移电流的影响,在旋转系统中研究导电流体非稳定的MHD Couette流动．在小数值磁场Reynolds数假定下,推导出基本的控制方程,使用著名的Laplace变换技术,数值地求解该基本方程．分两种情况:磁场相对于流体或者移动平板固定时,得到速度和表面摩擦力统一的闭式表达式．用图形讨论了问题的不同参数,对速度和表面摩擦力的影响．所得结果显示,主流速度u和次生速度v随着Hall电流而增大．离子滑移电流的增大,也会导致主流速度u的增大,但会使次生速度v减小．还给出了旋转、Hall和离子滑移参数的综合影响,确定了次生运动对流体流动的贡献．     

基于表面阻抗张量的界面滑移波动态失稳分析

基于Stroh公式和表面阻抗张量理论,提出了研究界面滑移波动态失稳问题的一种新的方法．该方法将表面阻抗张量概念推广到复波速域,并将摩擦接触界面上的边界条件以表面阻抗张量表示．最终将边值问题化归为求解一个复多项式在单位圆内的根．以弹性半空间与刚体平面相对稳态摩擦滑移为例进行了详细的分析,导出了一个4次复特征方程并讨论了方程在单位圆内解的特性,给出了滑移界面波失稳条件的显式解析表达式．     

转动与摆动复合运动下脂润滑关节轴承的数值分析

建立了球面轴承的三维润滑模型,该模型将内圈的转动运动、轴颈倾斜引起的内圈倾斜和内圈的摆动运动等因素纳入考虑,推导出球坐标下适用于非Newton(牛顿)流体润滑的Reynolds(雷诺)方程.应用该模型,并考虑使用润滑脂的Ostwald流变模型,对向心关节轴承的润滑问题进行了数值计算,研究了在不同的幂律指数、内圈倾斜角度和摆动角速度下,脂润滑膜的压力分布、最大压力、承载力和流量.结果表明:在合适的操作条件下,脂润滑能产生明显的流体动压效应;在其它参数不变时,幂律指数对脂润滑膜的最大压力和承载能力影响显著,相对于Newton流体,剪切稠化流体可提高润滑膜的最大压力和承载能力,并增加周向流量,而剪切稀化流体的影响效果则相反;内圈倾斜角度对脂润滑膜最大压力和承载能力的影响较小,内圈摆动角速度的影响则较为明显.     

纳米摩擦学的分子动力学模拟研究

用分子动力学模拟研究了纳米级润滑薄膜的固液相变和界面滑移现象以及固体接触和黏着的微观机制.结果表明：纳米薄膜中液体的固化相变压力随膜厚减薄而下降,说明润滑剂可能处于类固态状态.薄膜中的界面滑移现象可能在较低的剪切速率下发生,并与液体的固化程度有较好的对应关系.光滑晶体表面在相互接近或分离过程中可能因部分原子突然跃迁而发生微观黏着,这对于理解界面摩擦的起源具有重要意义.     

蠕动流中圆球阻力系数的入口效应

本文研究了圆球在半无穷长圆管入口处的蠕动流。得到了速度分布,压力分布和流函数的无穷级数形式的分析解.采用配置法将无穷级数截断并确定出级数中各项系数.求出了均匀入口流绕静止圆球以及圆球以瞬时速度在管内静止流体中运动这两种情形下圆球的阻力系数以及圆球表面上的应力分布.结果表明,当圆球在入口处运动时会遭受到较无穷圆管内为大的阻力.本文还对配置法的收敛性进行了数值实验.试验证明,该法具有好的收敛性.     

一种改进的Level-Set函数重新初始化方法

改进了level-set界面捕捉方法中的level-set函数重新初始化方法,通过锁定界面位置来实现level-set函数重新初始化过程,较好地保持了在重新初始化过程中的界面位置的准确性,保障了流体体积的相对稳定．对剪切运动使流体界面变形的算例表明方法对level-set函数重新初始化过程中的体积保持比较有效,而重新初始化过程总的机器用时并没有显著增加．     

Analysis of fluid separation in microfluidic T-channels

The behaviour of a fluid, which may contain particle suspensions, flowing in micro-dimensional channels is governed by both viscous and surface tension forces as well as high shear rates and geometric effects such as bifurcations, constriction, and high surface-to-volume ratio. This paper discusses some of the key design factors affecting fluid behaviour in micro-engineered products containing a main channel, constriction and side channel bifurcations. Differences in fluid behaviour at the macro and micro-scales are discussed. The dynamic bulk fluid behaviour is characterised in terms of: (i) fluid properties, (ii) governing physics and (iii) microchannel geometric features.At this stage of the analysis the fluids are assumed to be Newtonian and single phase, where any particle suspension is represented through a bulk density and viscosity. Based on these assumptions Computational Fluid Dynamics (CFD) is used to investigate the effect of both product inlet and outlet boundary conditions on the bulk flow behaviour. Discussions are provided on how these boundary conditions can affect particle separation efficiency. In particular, the so called pull-design whereby the fluid is pulled out of the device at the outlet, is shown to offer better performance compared to the mode of operation where the fluid is pushed into the device at the inlet. It is also observed that increasing the pressure at the outlet of the main channel can achieve a balanced flow rate ratio which leads to a uniform separation performance among all bifurcations.     

Modeling of flow in a very small surface separation

When a fluid flows in a very small surface separation, the very thin boundary layer physically adhering to the solid surface will participate in the flow, while between the two boundary layers is a continuum fluid flow. An analysis is here presented for this multiscale flow. The continuum fluid is treated as Newtonian. The physical adsorbed boundary layer is treated as non-continuum across the layer thickness. The interfacial slippage can occur on the adsorbed layer-solid surface interface, while it is absent on the adsorbed layer-fluid interface. Three flow equations are derived respectively for the two adsorbed layers and the intermediate continuum fluid. They together govern the multiscale flow in such a small surface separation.     

The wear contact problem with partial slippage

A method of analysing the variation in the stress state and shape of wearable bodies, subjected to the action of an oscillating shear load with partial slippage of the surface in the contact area is proposed. The method is based on the introduction of two time scales: the time of a single cycle of variation in the shear load and the time corresponding to a specified number of cycles. To estimate the shape variation of the surface due to wear a linear relation between the wear rate and the contact pressure and velocity of relative slippage is used. Cases of complete or partial removal of the wear particles from the friction zone are considered. As an example, the kinetics of the variation in the stress state in the contact of an elastic indentor, having a flat base and rounded edges, with an elastic half-space of the same material, is investigated. Analytical expressions for calculating the asymptotic values of the stresses and shape of the worn surface are obtained.     

Asymptotic Stability of a Stationary Flowing Regime of an Ideal Incompressible Fluid

We give sufficient conditions for asymptotic stability of a stationary solution to a flowing problem of a homogeneous incompressible fluid through a given planar domain. We consider a planar problem for the Euler equation and boundary conditions for the curl and the normal component of the velocity; moreover, the latter is given on the whole boundary of the flow domain and the curl is given only on the inlet part of the boundary. We establish asymptotic stability of a stationary flow (in linear approximation), assuming it to have no rest points and to satisfy some smallness condition which means that the perturbations leave the flow domain before they become to affect the main flow. In particular, we prove asymptotic stability for an arbitrary stationary flow in a rectangular canal close to the Couette flow without rest points. Moreover, we show that stability of the main flow in the L ₂-norm under curl perturbations implies its stability in higher-order norms depending, for example, on the derivatives of the curl.     

Axisymmetric flow over a backward-facing step in an annular pipe

The incompressible flow of a Newtonian fluid over a backward-facing step is investigated numerically. The geometry is an annular pipe in which the radius of the inner cylinder decreases suddenly. Keeping the radial expansion ratio fixed axisymmetric flows are computed for outlet radius ratios from 0.1 to 1 (ratio of the inner to the outer outlet radius). The Reynolds number at which the flow separates from the outer cylinder decreases as the outlet radius ratio decreases for constant inlet geometry. The growth with Reynolds number of the recirculation zone on the inner outlet cylinder just behind the step is strongly reduced when the recirculation zone on the outer cylinder is established. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Far field boundary conditions for incompressible flows computation

Many far field boundary conditions are proposed in the literature to solve Navier-Stokes equations. It is necessary to distinguish the streamwise or outlet boundary conditions and the spanwise boundary conditions. In the first case the flow crosses the artificial frontier and it is required to avoid reflections that can change significantly the flow. In the second case the Navier-slip boundary condition is often used but if the frontier is not far enough the boundary is both inlet and outlet. Thus the Navier-slip boundary condition is not well suited as it imposes no flux through the frontier. The aim of this work is to compare some well-known boundary conditions, to quantify to which extend the artificial frontier can be close to the bodies in two- and three-dimensions and to take into account the flow rate through the spanwise directions.     

Transport equation with boundary conditions of hysteresis type

This paper is an advanced extension of the work reported in (Nonlinear Anal. 2005; 63 :1467–1473). A transport equation that describes the propagation of a substance in a moving fluid or gas is considered. The equation contains the transient, convection, and diffusion terms. The problem is formulated in a bounded domain provided with an inlet and an outlet for the fluid or gas flow. The crucial point of the problem setting is a hysteresis‐type condition posed on an active part of the boundary. This condition reflects the nondecreasing accumulation with saturation of the transported substance at each point of the active boundary part. We prove the existence and uniqueness of solutions to this problem, study the regularity properties of solutions, and perform numerical simulations that clarify the behavior of the model. Comparing with the results of (Nonlinear Anal. 2005; 63 :1467–1473), the advancement of this work consists in accounting for the motion of the fluid or gas and posing inlet and outlet boundary conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

The role of diffusion of aggressive media in the fracture of rubber under stress

Formulas expressing the time to rupture of rubbers in aggressive media were derived in terms of the rate of diffusion of these media and the rate of displacement of the leading boundary of the cracking zone. The effect of increasing stress (due to a reduction in the load-carrying capacity of the specimen cross section) on the long-time strength (according to the Bailey principle) was taken into account in these calculations. The calculated and experimental data are in satisfactory agreement when there is a destructive action of working media in the absence of overstressing due to crack formation and swelling.Mekhanika Polimerov, Vol. 3, No. 4, pp. 708–712, 1967     

A simple method for spherical indentation of an elastic thin layer with surface tension bonded to a rigid substrate

An alternative method is proposed to solve the spherical indentation problem of an elastic thin layer with surface tension bonded to a rigid substrate. Based on the Kerr model, we establish a simple modified governing equation incorporating the surface tension effects for describing the relationship between the pressure and downward deflection of the impressed surface of the layer. This modified governing equation holds both inside and outside the contact zone, making it possible to analyze the whole layer by a unified differential equation. Numerical results are presented for the contact pressure inside the contact zone, the surface deflection of the elastic layer and the load-contact zone width relation to illustrate the present method. The validity and accuracy of the present method are demonstrated by comparing our results with those available in the existing literature.