Cessation of Newtonian circular and plane Couette flows with wall slip and non-zero slip yield stress

We solve analytically the cessation flows of a Newtonian fluid in circular and plane Couette geometries assuming that wall slip occurs provided that the wall shear stress exceeds a critical threshold, the slip yield stress. In steady-state, slip occurs only beyond a critical value of the angular velocity of the rotating inner cylinder in circular Couette flow or of the speed of the moving upper plate in plane Couette flow. Hence, in cessation, the classical no-slip solution holds if the corresponding wall speed is below the critical value. Otherwise, slip occurs only initially along both walls. Beyond a first critical time, slip along the fixed wall ceases, and beyond a second critical time slip ceases also along the initially moving wall. Beyond this second critical time no slip is observed and the decay of the velocity is faster. The velocity decays exponentially in all regimes and the decay is reduced with slip. The effects of slip and the slip yield stress are discussed.     

基于板理论的层级褶皱结构失效模式分析

本文基于弹性板理论和夹层板理论对二级层级褶皱结构失效模式进行了分析。通过对基本构件进行受力分析得到了载荷与结构变形之间的关系。根据6种失效模式的定义,从极限载荷或极限应力角度出发,分析了在压缩载荷和剪切载荷工况下的各种失效模式,给出了结构单胞对应的等效正应力和等效切应力表达式。由最小失效强度得到了各失效模式之间的占优关系,并构建了失效机理图来阐释这一机制。最后通过与有限元分析结果比较,分析了本文公式的精度,与数值解吻合较好。     

Ordinary state-based peridynamics for plastic deformation according to von Mises yield criteria with isotropic hardening

This study presents the ordinary state-based peridynamic constitutive relations for plastic deformation based on von Mises yield criteria with isotropic hardening. The peridynamic force density–stretch relations concerning elastic deformation are augmented with increments of force density and stretch for plastic deformation. The expressions for the yield function and the rule of incremental plastic stretch are derived in terms of the horizon, force density, shear modulus, and hardening parameter of the material. The yield surface is constructed based on the relationship between the effective stress and equivalent plastic stretch. The validity of peridynamic predictions is established by considering benchmark solutions concerning a plate under tension, a plate with a hole and a crack also under tension.     

Wall proximity effects on the effectiveness of upstream control rod

Two dimensional flow over a circular cylinder with an upstream control rod of same diameter is simulated in unbound condition and in wall bounded conditions. The cylinders are placed at various heights from the wall and the inter-distance between cylinders is also varied. The control rod is subjected to different rotation rates. It is found that, in unbound condition, rotating the control rod decreases the critical pitch length (S/D_cr) and increases the drag and Strouhal number of the main cylinder. In presence of plane wall, the shielding provided by the separated shear layers from the control rod in cavity regime is deteriorated due to deflection of shear layers which results in higher drag and large fluctuation of lift coefficient. However, in wake impingement regime, the binary vortices from the control rod are weakened due to diffusion of vorticity and hence, the main cylinder experiences a lower drag and small lift fluctuations than that of unbound condition. The critical height of vortex suppression (H/D_cr) is higher in cavity regime than that of wake impingement regime due to the single extended-bluff body like configuration. The rotation of control rod energizes the wall boundary layer and increases the critical height of vortex suppression. Increasing the rotational rate of control rod decreases the drag force and reduces the amplitude of lift fluctuation. Analysis of the wall shear stress distribution reveals that it suffers a sudden drop at moderate height where the normal Karman vortex shedding changes to irregular shedding consisting of single row of negative vortices. Modal structures obtained from dynamic mode decomposition (DMD) reveal that the flow structures behind the main cylinder are suppressed due to wall and the flow is dominated by the wake of control rod.     

Flow instability and shear localization in a drilling mud

To have a better knowledge of problems occurring with drilling fluids in complex wells, we carried out a detailed rheological analysis of a typical drilling mud at low shear rates using both conventional rheometry and MRI velocimetry. We show the existence of a viscosity bifurcation effect: Below a critical stress value, the mud tends to completely stop flowing, whereas beyond this critical stress, it reaches an apparent shear rate larger than a finite (critical) value, and no stable flows can be obtained between this critical shear rate value and zero. These results are confirmed by MRI velocity profiles, which exhibit a slope break at the interface between the solid and the liquid phases inside the Couette geometry. Moreover, this viscosity bifurcation is a transient phenomenon, the progressive development of which can be observed by MRI. A further examination of MRI data shows that, in the transient regime, the shear rate does not vary monotonously in the rheometer gap and is particularly large along the outer (rough) cylinder, which might be at the origin of the development of a region of constant shear rate in the apparent flow curve.     

Mixed Convection in Water Near the Density Extremum Along a Vertical Plate with Sinusoidal Surface Temperature Variation Embedded in a Porous Medium

A steady laminar boundary layer flowing along a vertical plate immersed in a Darcy–Brinkman porous medium saturated with water at 4°C is studied. The plate temperature varies sinusoidally along the plate between 0 and 8°C where the density of water varies parabolically and is almost symmetrical at about 4°C. Except for the existence of the buoyancy force, it is assumed that either the plate moves upwards or the ambient water moves upwards (moving stream). The results are obtained with the direct numerical solution of the boundary layer equations taking into account the temperature dependence of water thermophysical properties (ρ, μ and c _p). Results are presented for the wall temperature gradient and the wall shear stress along the plate for free convection and mixed convection. Temperature and velocity profiles are also presented.     

Granular vortices: Identification,characterization and conditions for the localization of deformation

We relate the micromechanics of vortex evolution to that of force chain buckling and, on this basis, formulate the conditions for strain localization in a continuum model of dense granular media. Using the traditional bifurcation analysis of shear bands, we show that kinematic vortex fields are in fact solutions to the boundary value problem satisfying null boundary conditions. To establish an empirical basis for our study, we first develop a method to identify the location of the core and boundary of each vortex from a given displacement field in two dimensions. We then employ this method to characterize the residual deformation field (i.e., the deviation of particle motions from the continuum deformation) in a physical experiment and a discrete element simulation of dense granular samples submitted to biaxial compression. Vortices in the failure regime are essentially confined to the shear band. Primary vortices, the clear majority, rotate in the same direction as the shear band; secondary vortices, the so-called wakes, rotate in the opposite direction. Primary vortices align in spatial succession along the central axis of the band; wakes form next to the band boundaries, in between and beside two adjacent primary vortices. Force chain buckling, the governing mechanism for shear bands, is responsible for vortex formation in the failure regime. Vortex dynamics are consistent with stick-slip dynamics. From quiescent conditions of jamming or stick, vortical motions arise from force chain buckling and associated relative particle rotations and sliding; these in turn precipitate intermittent periods of unjamming or slip, evident in the attendant drops in stress ratio and bursts in both kinetic energy and local nonaffine deformation. A kinematic vortex field inside shear bands is proposed that is consistent with the equations of continuum mechanics and the underlying instability of force chain buckling: such a field is periodic with a repeating unit cell comprising a primary vortex at the center of the band, with two trailing wakes close next to the band boundaries.     

The yield surface of viscoelastic and plastic fluids in a vane viscometer

Vane viscometers are often used to investigate the low shear rate properties of plastic fluids. The shear stress is determined by assuming that the material is held in the space between the vane blades so that it behaves like a rigid cylinder. Experimental evidence supports this assumption and the aim of the present study is to model numerically the yield process in a vane rheometer using viscoelastic and plastic fluids. The finite element method has been used to model the behavior of Herschel-Bulkley (Bingham), Casson and viscoelastic (Maxwell type) fluids. The penalty function approach for the pressure approximation and a rotating reference frame are used together with fine meshes containing more than 1300 elements. The results show that for Herschel-Bulkley (Bingham), and Casson fluids a rotating rigid cylinder of fluid is trapped inside the periphery of the vane, the shear stress is uniformly distributed over the surface of the cylinder. Finally a modified second order fluid is used to simulate the viscoelastic behaviour, anticipated to be an intermediate between the elastic deformation and the plastic flow, to provide a more realistic simulation of the yield process about a vane. In this case, contrast with the concentration of the elastic strain rate at the blade tips, a nearly uniform distribution of the plastic shear rate is still found. This implies that the plastic shear always distributes uniformly during the entire yielding process. Evidently the assumption of uniform shear on a rotating cylinder of material occluded in the blades of a vane is a valid and useful model for many types of fluid possessing a yield stress.     

Evaluation of the stress–strain curve of metallic materials by spherical indentation

A method for deducing the stress–strain uniaxial properties of metallic materials from instrumented spherical indentation is presented along with an experimental verification.An extensive finite element parametric analysis of the spherical indentation was performed in order to generate a database of load vs. depth of penetration curves for classes of materials selected in order to represent the metals commonly employed in structural applications. The stress–strain curves of the materials were represented with three parameters: the Young modulus for the elastic regime, the stress of proportionality limit and the strain-hardening coefficient for the elastic–plastic regime.The indentation curves simulated by the finite element analyses were fitted in order to obtain a continuous function which can produce accurate load vs. depth curves for any combination of the constitutive elastic–plastic parameters. On the basis of this continuous function, an optimization algorithm was then employed to deduce the material elastic–plastic parameters and the related stress–strain curve when the measured load vs. depth curve is available by an instrumented spherical indentation test.The proposed method was verified by comparing the predicted stress–strain curves with those directly measured for several metallic alloys having different mechanical properties.This result confirms the possibility to deduce the complete stress–strain curve of a metal alloy with good accuracy by a properly conducted instrumented spherical indentation test and a suitable interpretation technique of the measured quantities.     

Mechanical characteristics and origin of wall slip in pasty biosolids

We show that wall slip in pasty sewage sludge occurs in the form of an abrupt motion after a small (critical) deformation of the bulk. The analysis or our different results under prolonged creep tests at different stresses, after different times of rest following preshear or with different tool surfaces, tends to show that this phenomenon likely results from the unsticking of molecules or particles from the tool surface. Then the motion results from the flow of a thin layer of liquid (possibly mixed with some solid elements) which has formed along the interface while the rest of the material (bulk) remains unsheared. With a rough tool surface the flow at low apparent shear rates is also localized in a thin layer which in that case resembles the bulk material, so that the flow rate may decrease until reaching stoppage if the applied stress is not high enough, probably because of the progressive fluid restructuration.     

局部冲击作用下刚塑性平板的动力响应和失效模式

基于理想刚塑性的材料模型假定,同时考虑靶板弯曲塑性变形和剪切滑动,完整地分析了平板受刚性体撞击或受局部压力脉冲冲击的动力响应和失效模式。得到了撞击物侵彻深度、靶板穿透条件、塑性变形范围等特征变形破坏参数的解析表达式。讨论了撞击物速度、压力脉冲的冲量和形状等参量对响应和破坏模式的影响。比较了刚塑性理论分析结果与相关的实验和弹塑性数值计算结果。     

Elliptic yield cap constitutive modeling for high porosity sandstone

Field and laboratory investigators have observed thin, tabular zones of localized compressional deformation without shear in high porosity sandstone. These ‘compaction bands’ display greatly reduced porosity, and may affect the withdrawal of fluids from reservoirs. Studies addressing band formation as a type of strain localization predict the onset of the bands in a range of constitutive parameters roughly consistent with experiments, but are highly dependent on the constitutive relation used. In particular, the hardening modulus in shear and the slope of the yield surface in a plot of shear stress versus mean compressive stress are critical to localization predictions. Previous yield cap constitutive models employed a single deformation criterion, linking hydrostatic and shear response. In this work, we propose an elliptic yield cap model employing separate inelastic deformation parameters along each axis of the ellipse. The two deformation parameters allow the proposed surface to change in aspect ratio as it deforms, and allow a negative hardening modulus in shear without a negative hydrostatic modulus. Some cases with simplified modeling are shown for illustrative purposes, followed by a comparison with existing models. The proposed model displays similar strain behavior to the other models, but predicts localization under less restrictive conditions.     

FLOW CHARACTERISTICS AROUND AN INCLINED ELLIPTIC CYLINDER IN A TURBULENT BOUNDARY LAYER

The flow characteristics around an inclined elliptic cylinder located near a flat plate were investigated experimentally. The axis ratio of the elliptic cylinder was AR=2. The pressure distributions along the surface of the cylinder and the flat plate were measured by varying the angle of attack of the elliptic cylinder. The velocity profiles behind the cylinder were measured using hot-wire anemometry. When the angle of attack varies, the peak pressure location on the windward cylinder surface moves towards the rear edge of the cylinder, while that on the leeward surface moves towards the front edge of the cylinder. The vortex-shedding frequency also gradually decreases, defining a critical angle of attack for each gap ratio. The location of the minimum pressure on the flat plate surface moves downstream for positive angles of attack, while it moves upstream for negative angles of attack. Negative angles of attack cause a greater disturbance in the boundary layer near the wall compared to positive angles of attack. This shows that the separated wall shear layer from the boundary layer and the lower shear layer of the cylinder wake are strongly merged compared to other cases.     

切应力对牛骨蠕变变形的影响研究

骨尤其是湿骨,在恒定载荷作用下会发生蠕变变形。为了确定切应力是否影响骨的蠕变变形,采用对骨薄板试样分别施加集中载荷和均布载荷的方式,测量试样挠度实时的变化曲线。结果显示,在载荷恒定时,骨试样的挠度随时间不断增加,体现了典型的蠕变特性。集中荷载下骨的蠕变变形远大于均布荷载下骨的蠕变变形,湿骨的蠕变位移比干骨高近7倍。分析认为,对试样粘弹性性质的影响不仅有正应力的作用,也有切应力的作用;切应力产生的蠕变变形约为正应力所产生蠕变变形的0.85倍。     

Influence of curvature on drag reduction by opposition control in turbulent flow along a thin cylinder

Opposition controlled fully developed turbulent flow along a thin cylinder is analyzed by means of direct numerical simulations. The influence of cylinder curvature on the skin-friction drag reduction effect by the classical opposition control (i.e., the radial velocity control) is investigated. The curvature of the cylinder affects the uncontrolled flow statistics; for instance, skin-friction coefficient increases while Reynolds shear stress (RSS) and turbulent intensity decrease. However, the control effect in the case of a small curvature is similar to that in channel flow. When the curvature is large, the maximum drag reduction rate decreased. However, the optimal location of the detection plane is the same as that in a flat plate. Further, the drag reduction effect is achieved even on a high detection plane where the drag increases in the flat plate. Although a difference in the drag reduction effect can be observed with a change in the curvature, its mechanism considered in this analysis based on the transport of the Reynolds stress is similar to that of the flat plate.     

New universal relations for nonlinear isotropic elastic materials

A nonlinear isotropic elastic block is subjected to a homogeneous deformation consisting of simple shear superposed on triaxial extension. Two new relations are established for this deformation which are valid for all nonlinear elastic isotropic materials, and hence are universal relations. The first is a relation between the stretch ratios in the plane of shear and the amount of shear when the deformation is supported only by shear tractions. The second relation is established for a thin-walled cylinder under combined extension, inflation and torsion. Each material element of the cylinder undergoes the same local homogeneous deformation of shear superposed on triaxial extension. The properties of this deformation are used to establish a relation between pressure, twisting moment, angle of twist and current dimensions when no axial force is applied to the cylinder. It is shown that these relations also apply for a mixture of a nonlinear isotropic solid and a fluid.     

基于不同板理论的嵌于柔性层中薄膜的皱曲分析

论文对于柔性层-薄膜-柔性层三层结构系统,基于经典板理论、一阶剪切变形理论和高阶剪切变形理论,分别推导给出薄膜皱曲的控制方程.对于两个柔性层,则把它们处理成具有有限厚度的平面应变弹性体.针对上下柔性层固支边界或自由边界条件,利用线性扰动方法得到柔性层对薄膜的横向压力差,最终获得确定薄膜具有周期性正弦型皱曲的临界载荷方程...     

颗粒物质中最大静摩擦力随填充高度的变化规律

研究了圆棒在颗粒物质中受到的最大静摩擦力与颗粒填充高度的关系,并采用连续介质模型计算了圆棒在颗粒物质中受到的最大摩擦力.结果表明,当颗粒填充高度H小于圆筒直径的2倍时,最大静摩擦力F与H2成正比;当H大于2倍圆筒直径时,F与H成正比;理论计算结果与实验结果相符.