Liquid bridge force between two unequal-sized spheres or a sphere and a plane

Liquid bridge force acting between wet particles is an important property in particle characterization. This paper deals with liquid bridge force between either two unequal-sized spherical particles or a sphere and a flat plate under conditions where gravitational effect arising from bridge distortion is negligible. In order to calculate the force of the liquid bridge efficiently and accurately, expressions of liquid configuration and liquid bridge force were derived by building a mechanical model, which assumes the liquid bridge to be circular in shape between either two unequal-sized spheres or a sphere and a plane. To assess the accuracy of the numerical results of the calculated liquid bridge forces, they were compared to the published experimental data.     

Interaction between two spheres placed in tandem arrangement in steady and pulsating flow

The interaction among two spheres in tandem formation are studied for a Reynolds number of 300 using both steady and pulsating inflow conditions. The purpose is to further investigate the force characteristics as well as the shedding patterns of the two spheres as the separation distance is changed from 1.5 to 12 sphere diameters. The method used for the simulations is the volume of solid (VOS) method, an approach based on the volume of fluid (VOF) method. Comparisons with other computational methods have shown VOS to accurately resolve the flow field around solid spheres. The results show that the separation distance plays a significant role in changing the flow patterns and shedding frequencies at moderate separation distances, whereas effect on drag is observed even at a separation distance of 12 diameters.     

Exact solution for capillary interactions between two particles with fixed liquid volume

The capillary interactions, including the capillary force and capillary suction, between two unequal-sized particles with a fixed liquid volume are investigated. The capillary interaction model is used within the Young-Laplace framework. With the profile of the meridian of the liquid bridge, the capillary suction, and the liquid volume as state variables, the governing equations with two-fixed-point boundary are first derived using a variable substitution technique, in which the gravity effects are neglected. The capillary suction and geometry of the liquid bridge with a fixed volume are solved with a shooting method. In modeling the capillary force, the Gorge method is applied. The effects of various parameters including the distance between two particles, the ratio of particle radii, and the liquid-solid contact angles are discussed.     

Squeeze flow of a second-order fluid between two parallel disks or two spheres

The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reynolds‘ lubrication theory, the small parameter method was introduced to approximately analyze velocity field and stress distribution between the two disks. Then a similar procedure was carried out for analyzing the normal interaction between two nearly touching, arbitrary rigid spheres to obtain the pressure distribution and the resulting squeeze force. It has been proved that the solutions can be reduced to the case of a Newtonian fluid when the non-Newtonian terms are neelected.     

A fast and practical method to pack spheres for mesh generation

Sphere packing is an attractive way to generate high quality mesh. Several algorithms have been proposed in this topic, however these algorithms are not sufficiently fast for large scale problems. The paper presents an efficient sphere packing algorithm which is much faster and appears to be the most practical among all sphere packing methods presented so far for mesh generation. The algorithm packs spheres inside a domain using advancing front method. High efficiency has resulted from a concept of 4R measure, which localizes all the computations involved in the whole sphere packing process.     

Analytic calculation of magnetic force between two current-carrying coils

Current-carrying coils are basic elements in electromagnetic equipments, for example, in high field magnets from high temperature superconducting wires or tapes. In the assembly of these systems and their current-carrying operation, unavoidable mis-alignment and shift from the original position can be induced by disturbances such as the imbalance of magnetic force due to safety problems. For two current-carrying coils with non-coplanar axes, the analytic expression of the magnetic force between the two coils is presented according to the rule of Ampere circulation and the Biot-Savart law. Based on the expression, the dependence of the magnetic force on the size and the relative position of each other is further investigated, and the variation of the magnetic force is obtained with the above parameters.     

Squeeze flow of a power-law fluid between two rigid spheres with wall slip

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｓｑｕｅｅｚｅｆｌｏｗｏｆａｆｌｕｉｄｂｅｔｗｅｅｎｔｗｏｄｉｓｋｓｏｒｓｐｈｅｒｅｓｉｓｏｆｒｅｌｅｖａｎｃｅｔｏｍａｎｙａｐｐｌｉｃａｔｉｏｎｓ,ｉｎｃｌｕｄｉｎｇｔｈｅｆｏｒｍｉｎｇｏｆｐｏｌｙｍｅｒｍａｔｅｒｉａｌｓ ,ｓｑｕｅｅｚｅｆｌｏｗｒｈｅｏｍｅｔｅｒａｎｄｌｕｂｒｉｃａｔｉｏｎｏｆｂｅａｒｉｎｇｓ.Ｔｈｅｓｑｕｅｅｚｅｆｌｏｗｉｎｔｅｒａｃｔｉｏｎｂｅｔｗｅｅｎｓｏｌｉｄｐａｒｔｉｃｌｅｓｉｓａｌｓｏｆｕｎｄａｍｅｎｔａｌｔｏｔｈｅｃｏｍｐｌｅｘｒｈｅ…     

A geometrical model for numerical simulation of capillary imbibition in sedimentary rocks

The cylindrical model is discussed and a new tube model is proposed to describe capillary imbibition kinetics in porous sedimentary rocks. The tube consists of a periodic succession of a single hollow spherical element of which the geometry is defined by the sphere radius and the sphere access radius. These two parameters are estimated experimentally for four rock types from their specific surface areas. Introducing those parameters in the model capillary imbibition kinetics, parameters are calculated and compared with the experimental ones. A direct relation between imbibition kinetics and specific surface area has been pointed out.     

Maximising the electrorheological effect for bidisperse nanofluids from the electrostatic force between two particles

A multipole re-expansion solution for two nonidentical dielectric spheres in a parallel electric field is used to determine the critical ratio of particle radii which leads to the strongest force of attraction between the spheres at various interstices and under varying dielectric properties. These critical ratios provide genuine optimal dimensions, in the sense that the force of attraction decreases for both increasing and decreasing ratios. Numerical results are compared with experimental results from the literature and discussed from the perspective of the impact on the design of electrorheological nanofluids.     

A variational model for fracture mechanics: Numerical experiments

In the variational model for brittle fracture proposed in Francfort and Marigo [1998. Revisiting brittle fracture as an energy minimization problem. J. Mech. Phys. Solids 46, 1319-1342], the minimum problem is formulated as a free discontinuity problem for the energy functional of a linear elastic body. A family of approximating regularized problems is then defined, each of which can be solved numerically by a finite element procedure. Here we re-formulate the minimum problem within the context of finite elasticity. The main change is the introduction of the dependence of the strain energy density on the determinant of the deformation gradient. This change requires new, more general existence and Γ-convergence results. The results of some two-dimensional numerical simulations are presented, and compared with corresponding simulations made in Bourdin et al. [2000. Numerical experiments in revisited brittle fracture. J. Mech. Phys. Solids 48, 797-826] for the linear elastic model.     

A domain independent integral expression for the crack extension force of a curved crack in three dimensions

An integral expression that is domain independent in curvilinear coordinates and compatible with zero divergence of Eshelby's (Phil. Trans. Roy. Soc. (London) 244 (1951) 87.) energy momentum tensor was obtained from the principle of virtual work. By applying Eshelby's definition of the force on a material defect a general expression of the crack extension force for a curved crack in three dimensions, here called the F-integral, was derived from the domain independent integral expression. The F-integral is given explicitly for a number of curved cracks and found to be in agreement with previously known solutions, when available. The influence of crack surface and crack front curvature upon the various forms of the F-integral is discussed. The F-integral presented in this work is a generalisation of the J-integral (Rice, J. Appl. Mech. 35 (1968) 379.) to curved cracks in orthogonal curvilinear coordinates.     

A new algorithm for surface tension model in moving particle methods

A new algorithm for the surface tension model was developed for moving particle methods. The algorithm is based on the link‐list search algorithm and the continuum surface tension (CST) model. The developed algorithm with the CST model was implemented to a kind of moving particle approach, the finite volume particle (FVP) method. The FVP method with the new algorithm was tested by oscillatory behaviour of a two‐dimensional droplet. The oscillatory period agrees well with analytical one, and the transient shape of the droplet is also in good agreement with that obtained by other numerical methods. The droplet impact on a liquid surface was also studied using the new algorithm. The deposition and splashing phenomena were clearly reproduced. Simulated spread radius of the splashing phenomena was consistent with a power law. Copyright © 2007 John Wiley & Sons, Ltd.     

Simulation of non-Newtonian ink transfer between two separating plates for gravure-offset printing

The inks used in gravure-offset printing are non-Newtonian fluids with higher viscosities and lower surface tensions compared to Newtonian fluids. This paper examines the transfer of a non-Newtonian ink between two parallel plates when the top plate is moved upward with a constant velocity while the bottom plate is held fixed. Numerical simulations were carried out using the Carreau model to explore the behavior of a non-Newtonian ink in gravure-offset printing. The volume of fluid (VOF) model was adopted to demonstrate the stretching and break-up behaviors of the ink. The results indicate that the ink transfer ratio is greatly influenced by the contact angle, especially the contact angle at the upper plate (α). For lower values of α, oscillatory or unstable behavior of the position of minimum thickness of the ink between the two parallel plates during the stretching period is observed. This oscillation gradually diminishes as the contact angle at the upper plate is increased. Moreover, the number of satellite droplets increases as the velocity of the upper plate is increased. The surface tension of the conductive ink shows a positive impact on the ink transfer ratio to the upper plate. Indeed, the velocity of the upper plate has a significant influence on the ink transfer in gravure-offset printing when the Capillary number (Ca) is greater than 1 and the surface tension dominates over the ink transfer process when Ca is less than 1.     

The coefficient of normal restitution for the Hertzian contact of two rough spheres colliding in a viscous fluid

We use previous theoretical results for the added mass, history and lubrication forces between two spheres colliding in a fluid with viscosity ν to investigate the effect of viscous dissipation on the coefficient of restitution during contact. We assume that the mechanical interaction is governed by Hertzian mechanical contact of small duration τ and that the minimum approach distance between particles is approximately equal to the height σ of surface micro-asperities. A non-dimensionalization of the equation of motion indicates that the contact dynamics is governed by two parameters – the ratio ϵ of the surface roughness σ and the sphere radius a, and a contact Stokes number defined as St_c = σ²/ντ. An asymptotic solution of the equation of motion in the limit of small ϵ/St_c is used to obtain an explicit expression for the coefficient of restitution during contact and the latter is compared with estimates based on numerical solutions of the non-linear equation of motion.     

两异径湿球颗粒间静态液桥力的近似解析公式及数值评测

针对异径球形湿颗粒间液桥力,以Young-Laplace公式为基础,结合环形近似法并引入等效半径,在宽松设定的条件下推导出了简化的近似解析公式。此公式形式简洁,且当两球半径相等时,可回归到Pitois等提出的径球液桥力公式。本文对此近似解析公式进行参数敏感度分析,发现随着颗粒间距的增大,液桥力对颗粒半径的敏感度降低。此外,本文近似解析公式与其他实验和理论结果进行对比发现,该近似解析公式与上述实验和数值结果相当吻合,表明该近似解析公式可以较准确地计算接触角较小时的异径球颗粒间的液桥力,因而适于湿颗粒系统的数值模拟。     

A continuum model for initiation and evolution of deformation twinning

Within continuum dislocation theory the plastic deformation of a single crystal with one active slip system under plane-strain constrained shear is investigated. By introducing a twinning shear into the energy of the crystal, we show that in a certain range of straining the formation of deformation twins becomes energetically preferable. An energetic threshold for the onset of twinning is determined. A rough analysis qualitatively describes not only the evolving volume fractions of twins but also their number during straining. Finally, we analyze the evolution of deformation twins and of the dislocation network at non-zero dissipation. We present the corresponding stress-strain hysteresis, the evolution of the plastic distortion, the twin volume fractions and the dislocation densities.     

A hybrid force model to estimate the dynamics of curved legs in granular material

Robot locomotion on rigid terrain or in fluids has been studied to a large extent. The locomotion dynamics on or within soft substrates such as granular material (GM) has not been fully investigated. This paper proposes a hybrid force model to simulate and evaluate the locomotion performance of a legged terrestrial robot in GM. The model incorporates an improved Resistive Force Theory (RFT) model and a failure-based model. The improved RFT model integrates the force components of individual leg elements over the curved leg portion submerged in GM at any moment during a full period of leg rotation. The failure-based model is applied in a bar drag model to yield the normal and the lateral forces of the individual RFT elements as functions of the locomotion depth and speed. The hybrid model is verified by the coincidence between the theoretical predictions and the experimental results. The hybrid model is used to analyze the effects of angular velocity and leg shape with high precision and can guide the design of the legs with any profiles. Our study reveals that the interactions between locomotor and substrate are determined by the locomotor structural characteristics, the nature of the substrate, and the control strategy.     

A moving particle semi‐implicit method for free surface flow: Improvement in inter‐particle force stabilization and consistency restoring

The moving particle semi‐implicit (MPS) method has been widely applied in free surface flows. However, the implementation of MPS remains limited because of compressive instability occurred when the particles are under compressive stress states. This study proposed an inter‐particle force stabilization and consistency restoring MPS (IFS‐CR‐MPS) method to overcome this numerical instability. For inter‐particle force stabilization, a hyperbolic‐shaped quintic kernel function is developed with a non‐negative and smooth second order derivative to satisfy the stability criterion under compressive stress state. Then, a contrastive study is conducted on the contradiction between the common understanding of the conventional MPS hyperbolic‐shaped kernel function and its performance. The result shows that the conventional MPS hyperbolic‐shaped kernel function can easily cause violent repulsive inter‐particle force and then lead to the compressive instability. Therefore, the first order derivative of the modified hyperbolic‐shaped quintic kernel function is recommended as the form of the contribution of the neighbor particles to achieve a more stable inter‐particle repulsive force. For consistency restoring, the Taylor series expansion and the hyperbolic‐shaped quintic kernel are combined to improve the accuracy of the viscosity and pressure calculation. The IFS‐CR‐MPS algorithm is subsequently verified by the inviscid hydrostatic pressure, jet impacting, and viscous droplet impacting problems. These results can be used for choosing kernel function and the contribution of neighbor particles in particle methods. Copyright © 2016 John Wiley & Sons, Ltd.     

径向磁化的多环嵌套永磁轴承轴向磁力解析模型

为了解决径向磁化双环永磁轴承轴向磁力偏小的问题,设计了径向磁化多环嵌套永磁轴承新结构。基于磁路及虚功原理法,结合径向磁化多环嵌套永磁轴承结构特点及线性叠加原理,建立了径向磁化多环嵌套永磁轴承轴向磁力解析模型。模型表明:径向磁化多环嵌套永磁轴承轴向磁力与磁环剩磁平方成正比,磁力随磁环间隙的增大而减小,随磁环数的增大而增大;在正常轴向工作范围内,轴向磁力随轴承轴向偏移的增大而增大。模型计算结果与有限元计算结果基本吻合。对比计算说明:径向磁化多环嵌套永磁轴承轴向磁力远大于由原多嵌套磁环所构成的若干双环永磁轴承轴向磁力之和。     

Non-slipping adhesive contact between mismatched elastic spheres: A model of adhesion mediated deformation sensor

A generalized JKR model is established for non-slipping adhesive contact between two dissimilar elastic spheres subjected to a pair of pulling forces and a mismatch strain. We discuss the full elastic solution to the problem as well as the so-called non-oscillatory solution in which tension and shear tractions along the contact interface is decoupled from each other. The model indicates that the mismatch strain has significant effect on the contact area and the pull-off process. Under a finite pulling force, a pair of adhering spheres is predicted to break apart spontaneously at a critical mismatch strain. This study suggests an adhesion mediated deformation sensing mechanism by which cells and molecules can detect mechanical signals in the environment via adhesive interactions.