Fluid inertia in large amplitude oscillatory shear

Homogeneous shearing is required in sliding plate flow experiments with one plate fixed and the other oscillating. However, when fluid inertia becomes significant, the velocity gradient and the stress will not be uniform. MacDonald et al. (1969) and Schrag (1977) investigated this effect for a linear viscoelastic fluid. However, linear viscoelasticity does not describe the behavior of melts in large amplitude oscillatory shear (LAOS). Jeyaseelan et al. (1993) have shown that the Berkeley kinetic network model does accurately describe the LAOS behavior of polymer melts. In this work, the Berkeley model is solved for LAOS in sliding plate flow with fluid inertia, by numerical integration of spatially discretized forms of the governing equations. Nonlinear viscoelasticity is predicted to aggravate the effects of fluid inertia in LAOS and experiments confirm this. Specifically, fluid inertia amplifies the first harmonic and produces no even harmonics. Operating limits are presented graphically for minimizing inertial effects in LAOS experiments. Received: 2 January 1998 Accepted: 27 April 1998     

冲击下混凝土试样应变率效应和惯性效应探讨

结合混凝土试件的真三轴静载冲击实验结果,分别运用考虑应变率效应的Holmquist-Johnson-Cook (HJC)模型和考虑静水压效应的Drucker-Prager (DP)模型进行数值分析,以探讨研究混凝土试样应变率效应和惯性效应的方法。在探究混凝土的应变率效应和横向惯性效应的关系时,使用HJC模型的数值模拟结果来拟合DP准则的各个参数。结果表明:随着应变率的升高,混凝土的强度会提高,并且这种强度的提高,也有一部分原因是第一应力不变量I₁的增大所导致的。因此,混凝土试件的应变率效应和横向惯性约束具有较强的耦合作用。理论和数值分析了冲击下试样内部的横向应力分布特征与应变率、静水压和试样尺寸的关系,结果发现:试样内部横向应力的幅值随着应变率、静水压的升高而增大,但随着试样尺寸的增大而减小。为了探讨横向惯性带来的强度提升效果,提出了一个有关冲击方向最大应力σ_x和等效应力σ_e的参数ξ,且ξ=(σ_x?σ_e)/σ_x。此参数具有尺寸效应、应变率效应和静水压效应,但是此参数与应力三轴度的关系表现出应变率无关特性,可为应变率效应的研究提供新的思路。     

控制致裂岩体应力量测的模型试验研究

本文阐明了控制致裂岩体应力量测的基本原理、试验方法。介绍了几种控制致裂的模型试验结果及其分析，证明了控制致裂理论的正确性。     

刚性圆管中血液周期振荡流的切应力分布

本文通过求解圆管内血液振荡流的基本方程，求得圆管内血液流的压力梯度与切应力之间的关系式。在此基础上，详细讲座了圆管中轴向流速和切变率谐波的变化规律，指出流速谐波和切变率谐波的幅值都将随着谐波次数的增大而逐渐减小。为了使所得结果便于应用。文章通过管轴向中心线流速与压力梯度之间的关系式，进一步给出一种利用管轴向中心线流速计算管内切应力分布的简便方法。该方法用于检测活体血管内血液振荡流的切应力分布，具有操作简单，精度较高的优点。最后，以人体颈动脉为例，讨论血液周期振荡流的切应力的分布特性。发现在任意时刻，除了邻近管壁处切应力急剧增大到一定数值之外，沿管截面切应力分布相当均匀且接近于零，呈现出与定常流不同的切应力分布特征。     

深孔控制卸压爆破机理和防突试验研究

针对当前各种局部防突措施的优缺点和局限性,对兼有松动爆破与大直径卸压钻孔爆破两种措施优点的局部防突措施—深孔控制卸压爆破进行了研究。在理论上和利用实验室模型实验对深孔控制卸压爆破中的控制孔的作用进行了研究分析。深孔控制卸压爆破的现场防突试验表明:该法效果明显、工艺简单、经济效益好、安全且适应性强,是一种易于推广的局部防突措施,有着广阔的应用前景。     

Some exact solutions of the oscillatory motion of a generalized second grade fluid in an annular region of two cylinders

The velocity field and the associated shear stress corresponding to the longitudinal oscillatory flow of a generalized second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. Initially, the fluid and cylinders are at rest and at t = 0＋ both cylinders suddenly begin to oscillate along their common axis with simple harmonic motions having angular frequencies Ω1 and Ω2. The solutions that have been obtained are presented under integral and series forms in terms of the generalized G and R functions and satisfy the governing differential equation and all imposed initial and boundary conditions. The respective solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for the similar flow of ordinary second grade fluid and Newtonian fluid are also obtained as limiting cases of our general solutions. At the end, the effect of different parameters on the flow of ordinary second grade and generalized second grade fluid are investigated graphically by plotting velocity profiles.     

Two-fluid oscillatory flow in a channel

The validity of Navier's partial slip condition is investigated by studying the oscillatory flow in a coated channel. The two-fluid model is used to solve the unsteady viscous equations exactly. Partial slip is experienced by the core fluid. It is found that Naviers condition does not hold for an unsteady core fluid.     

Convective inertia effects in wall-bounded thin film flows

We apply the boundary layer equations to inertial flow in wall bounded films that might be characterized as ‘thin’, say ɛ ≤ 0.1 where ɛ is the ratio of the characteristic lengths, yet to which the lubrication approximation of Reynolds no longer applies. Two particular flow geometries are investigated, nominally parallel plates and nominally inclined plates, both with and without spatially periodic perturbation of the stationary plate. A Galerkin-B spline formulation of the governing equations is employed, and we rely on parametric continuation to obtain solutions at higher values of the Reynolds number. In particular, we are able to demonstrate that the boundary layer equations yield accurate results for a wide range of Reynolds numbers when the aspect ratio is less than 1/10. We also find that in both nominally parallel and nominally inclined geometries the sign of the inertial force correction is determined by the film contour in the neighborhood of the exit, this result might have implications in the design of MEMS devices.     

Acoustoelastic effects on mode waves in a fluid-filled pressurized borehole in triaxially stressed formations

Based on the nonlinear theory of acoustoelasticity, considering the triaxial terrestrial stress, the fluid static pressure in the borehole and the fluid nonlinear effect jointly, the dispersion curves of the monopole Stoneley wave and dipole flexural wave propagating along the borehole axis in a homogeneous isotropic formation are investigated by using the perturbation method. The relation of the sensitivity coefficient and the velocity-stress coefficient to frequency are also analyzed. The results show that variations of the phase velocity dispersion curve are mainly affected by three sensitivity coefficients related to third-order elastic constant. The borehole stress concentration causes a split of the flexural waves and an intersection of the dispersion curves of the flexural waves polarized in directions parallel and normal to the uniaxial horizontal stress direction. The stress-induced formation anisotropy is only dependent on the horizontal deviatoric terrestrial stress and independent of the horizontal mean terrestrial stress, the superimposed stress and the fluid static pressure. The horizontal terrestrial stress ratio ranging from 0 to 1 reduces the stress-induced formation anisotropy. This makes the intersection of flexural wave dispersion curves not distinguishable. The effect of the fluid nonlinearity on the dispersion curve of the mode wave is small and can be ignored.The project supported by the National Natural Science Foundation of China (10272004) and The Special Science Foundation of the Doctoral Discipline of the Ministry of Education of China(20050001016) The English text was polished by Keren Wang.     

狭窄血管内偶应力流体流动分析

本文考察了血管狭窄对血液流动的影响,血液以偶应力流体表示,并在求解过程中采用了在管壁上流体质点无相对涡量的边界条件,结果表明,和Young的经典工作相比流动阻抗和壁切应力大于同样程度狭窄下牛顿流体的相应值,偶应力流体对狭窄的敏感性大于牛顿流体;在狭窄发展过程中,偶应力流体的流量要小于牛顿流体的流量,和牛顿流体相比,这些结果更符合生理实际。     

Hall effects on the unsteady hydromagnetic oscillatory flow of a second-grade fluid

An exact solution of an oscillatory flow is constructed in a rotating fluid under the influence of an uniform transverse magnetic field. The fluid is considered as second-grade (non-Newtonian). The influence of Hall currents and material parameters of the second-grade fluid is investigated. The hydromagnetic flow is generated in the uniformly rotating fluid bounded between two rigid non-conducting parallel plates by small amplitude oscillations of the upper plate. The exact solutions of the steady and unsteady velocity fields are constructed. It is found that the steady solution depends on the Hall parameter but is independent of the material parameter of the fluid. The unsteady part of the solution depends upon both (Hall and material) parameters. Attention is focused upon the physical nature of the solution, and the structure of the various kinds of boundary layers is examined. Several results of physical interest have been deduced in limiting cases.     

非牛顿流体在渐变管中压力和剪切应力的二次摄动解

本文利用双摄动方法求解缓慢变化管道中Johnson-Segalman(J-S)流体流动的渐近解．将管道的扩张（或收缩）角度和粘弹性参数分别作为双摄动的参数,由流函数和涡量函数的形式,推导出压力和壁面剪切应力的渐近解．在此基础上,分析了管道角度,粘弹性参数和雷诺数等参数对压力以及剪切力影响．主要结论如下：(1) 管道扩张角度增加时,流向同一位置处径向压力以及壁面剪切应力随扩张角度减小;(2) 在同一扩张管道中,径向压力随着流向位移减小,收缩管与之相反;(3) 扩张角度与雷诺数对流场起主导作用,粘弹性系数起次要作用．     

Convective Instability of Oldroyd-B Fluid Saturated Porous Layer Heated from Below using a Thermal Non-equilibrium Model

The linear stability of a viscoelastic fluid saturated densely packed horizontal porous layer heated from below and cooled from above is investigated by considering the Oldroyd-B type fluid. A generalized Darcy model, which takes into account the viscoelastic properties, is employed as momentum equation and a two-field model is used for energy equation each representing solid and fluid phases separately. Linear stability analysis suggests that, there is a competition between the processes of viscous relaxation and thermal diffusion that causes the first convective instability to be oscillatory rather than stationary. Analytical expression for the occurrence of oscillatory onset is obtained, and it is found that the necessary condition for the existence of the same is Λ < 1. Besides, the effect of viscoelastic parameters and the thermal non-equilibrium on the stability of the system is analyzed.     

Problems of non-uniqueness when interpreting the effect of fluid inertia on the complex viscosity function

A theoretical investigation is carried out into the interpretation of the effect of fluid inertia on the complex viscosity function as measured on a controlled stress rheometer. The problem of non-unique solutions to the governing equations is considered for the parallel plate geometry. The locations of these solutions are investigated by considering the critical points of the complex mapping associated with the linear viscoelastic equations of motion. It is shown that these critical points play an important role in determining where convergence problems are likely to occur when applying numerical methods of solution to the governing equations. Analytical approximations based on a series expansion about a critical point are developed as an alternative approach to a numerical solution in the neighbourhood of a critical point. In order to verify the theoretical predictions a numerical simulation of the behaviour of a single element Maxwell fluid on a controlled stress rheometer is carried out for a parallel plate geometry. Received: 27 July 1998 Accepted: 9 April 1999     

Using instrumental inertia in controlled stress rheometry

We describe a new method for characterizing the non-linear behavior of complex fluids at both small and large deformations. For creep measurements, we use the coupling between the instrumental inertia and the material‘s elasticity to follow the rheological behavior of a solution of iota carrageenan both above and below the yield stress. It is shown that this coupling selectively excites one particular frequency of the relaxation spectrum. An analytical calculation is used to quantify the non-linear behavior near the yield stress. The “free“ oscillations observed during the first few seconds allow us to choose the most appropriate mechanical model. Comparison with experiment shows that even above the yield stress, a linear model can still give independently reliable information about the changes in each element of the mechanical model. A comparison of free and forced oscillations in controlled stress rheometry shows both experimentally and theoretically the conditions under which the use of free oscillations is advantageous. Received: 4 September1997 Accepted: 13 January 1998     

Inertia effects of the fluid and cylinder in coaxial cylindrical rheometer

The inertias of the fluid and the inner cylinder in coaxial cylinder rheometer (CACR) have great influence on the unsteady flow of non-Newtonian fluid. Even for the Newtonian fluid there exist the so called “stress overshoot” phenomenon. In the present article this phenomenon was studied in detail and a method correcting the measured results for an unsteady flow in the rheometer was proposed. It is found that the inertia effect of the fluid can be ignored when the gap between cylinders is small. The project was supported by the Youth Foundation of the Fourth Military Medical University     

A note on start-up and large amplitude oscillatory shear flow of multimode viscoelastic fluids

Start up of plane Couette flow and large amplitude oscillatory shear flow of single and multimode Maxwell fluids as well as Oldroyd-B fluids have been analyzed by analytical or semi-analytical procedures. The result of our analysis indicates that if a single or a multimode Maxwell fluid has a relaxation time comparable or smaller than the rate of change of force imparted on the fluid, then the fluid response is not singular as Elasticity Number (E ). However, if this is not the case, as E , perturbations of single and multimode Maxwell fluids give rise to highly oscillatory velocity and stress fields. Hence, their behavior is singular in this limit. Moreover, we have observed that transients in velocity and stresses that are caused by propagation of shear waves in Maxwell fluids are damped much more quickly in the presence of faster and faster relaxing modes. In addition, we have shown that the Oldroyd-B model gives rise to results quantitatively similar to multimode Maxwell fluids at times larger than the fastest relaxation time of the multimode Maxwell fluid. This suggests that the effect of fast relaxing modes is equivalent to viscous effects at times larger than the fastest relaxation time of the fluid. Moreover, the analysis of shear wave propagation in multimode Maxwell fluids clearly show that the dynamics of wave propagation are governed by an effective relaxation and viscosity spectra. Finally, no quasi-periodic or chaotic flows were observed as a result of interaction of shear waves in large amplitude oscillatory shear flows for any combination of frequency and amplitudes.     

Simulation of centrifugal inertia forces on a test bed of blow-molding machines

This paper deals with the simulation of the centrifugal inertia forces on a fixed test bed that reproduces only the opening and closing motions of the half-molds of blow-molding machines. The suggested solution is realized by the addition of springs, which allows the reproduction of the moment of the inertia forces about the vertical axis of the mold. By using the approximation method based on the mean-square-value minimization, satisfactory results are obtained. The efficiency of such a solution is illustrated by numerical example.     

The third-normal stress difference in entangled melts: Quantitative stress-optical measurements in oscillatory shear

Stress-optical measurements are used to quantitatively determine the third-normal stress difference (N ₃ = N ₁ + N ₂) in three entangled polymer melts during small amplitude (<15%) oscillatory shear over a wide dynamic range. The results are presented in terms of the three material functions that describe N ₃ in oscillatory shear: the real and imaginary parts of its complex amplitude ₃ ^* = ₃ - i ₃ , and its displacement ₃ ^d . The results confirm that these functions are related to the dynamic modulus by ² ₃ ^* ()=(1-)[G ^*())– G ^*(2)] and ² ₃ ^d ()=(1- )G() as predicted by many constitutive equations, where = –N ₂/N ₁. The value of (1-) is found to be 0.69±0.07 for poly(ethylene-propylene) and 0.76±0.07 for polyisoprene. This corresponds to –N ₂/N ₁ = 0.31 and 0.24±0.07, close to the prediction of the reptation model when the independent alignment approximation is used, i.e., –N ₂/N ₁ = 2/7 – 0.28.     

Slip effects and heat transfer analysis in a viscous fluid over an oscillatory stretching surface

In this study, we investigate the heat transfer problem in a viscous fluid over an oscillatory infinite sheet with slip condition. The sheet is moved back and forth in its own plane. The derived problem involves a dimensionless parameter indicating the relative magnitude of frequency to sheet stretching rate. A system of non‐linear partial differential equations is solved numerically using the finite‐difference scheme, in which a coordinate transformation is employed to transform the semi‐infinite physical space to a bounded computational domain. The physical features of interesting parameters on the velocity and temperature distributions are shown graphically and discussed. The values of the skin‐friction coefficient and the local Nusselt number are given in tabular form. Copyright © 2008 John Wiley & Sons, Ltd.