作旋转流动时非线性流体的改进模型

在圆环结构中研究拟塑性流体作圆形的Couette流动.流体的粘度依赖于对守恒方程有直接影响的剪切率,守恒方程采用谱方法求解.可以证明所采用的拟塑性模型,可以被适当地表示为典型的非线性流动.在早期研究中,为了方便数值计算,粘度表达式中只考虑了剪切率的二次项,与此不同,这里考虑了二次幂项.圆形Couette流动中弯曲的流线,造成离心的不稳定性,引起环形的漩涡,称之为Taylor漩涡.进而发现,随着拟塑性影响的增加,临界Taylor数下降.与已有圆形Couette流动的实验相比较,两者有着良好的一致性.     

弹性-应变软化粘塑性材料反平面剪切动态扩展裂纹尖端的渐近解*

本文首先给出了一种用于描述材料软化,并存在有粘塑性的材料模型.用这种模型对反平面剪切型动态扩展状态下,裂纹尖端的弹粘塑性场进行了渐近分析,给出了弹性－应变软化粘塑性材料反平面剪切动态扩展裂纹尖端的渐近解方程.分析结果表明,在裂纹尖端应变具有(ln(R/r))1/(n＋1)的奇异性,应力具有(ln(R/r))-n/(n＋1)的奇异性.从而本文揭示了应变软化粘塑性材料反平面剪切动态扩展裂纹尖端的渐近行为.     

弹性-粘塑性材料Ⅰ型动态扩展裂纹尖端场的渐近解

提出了一种新的弹性-粘塑性模型用于分析Ⅰ型动态扩展裂纹尖端的应力应变场.给出了适当的位移模式,推导了渐近方程并且给出了数值解.分析和计算表明:对于低粘性情况,裂纹尖端场具有对数奇异性;对于高粘性情况,渐近方程无解.分析比较表明该结果具有高玉臣提出的单参数解的所有优点,并且消除了粘性区随裂纹扩展而移动的不足.     

血液流动模拟方法对比分析

对现有主要血液流动模拟方法进行了对比研究.以单个直血管为研究对象,分别用牛顿单相流模型、非牛顿单相流模型、液固两相流剪切稀化模型、液固两相流颗粒动力学模型和血液两相流修正模型计算了血液流场,然后对比了五种模型模拟得到的血液动力学参数:血液流速、红细胞体积分数、流体粘度和壁面剪切应力.结果表明:血液组成和粘度方程的选择对血液速度场计算结果有明显影响;流体本构方程和升力公式对壁面剪切应力计算结果均有明显影响;液固两相流颗粒动力学模型计算得到的血液粘度远偏离正常的血液粘度,模型不适用于模拟稳态血液流动;血液两相流修正模型可以较准确模拟红细胞在血管内的径向分布及其影响.研究结果可为今后相关研究中血液模拟方法的选择提供指导.     

一类Stokes方程的最小二乘混合元方法

对定常和非定常两种类型的Stokes方程建立了一类新的最小二乘混合元方法,并进行了分析,对定常的方程,采用了对u和σ的不同指标的有限元空间进行计算(LBB条件不需要),得到了最优的H1和L2模估计.对非定常的方程,采用了传统的Raviart-Thomas混合元空间,得到了最优的L2模估计.     

高速扩展平面应力裂纹尖端的各向异性塑性场

在裂纹尖端的应力分量都只是θ的函数的条件下,利用定常运动方程,Hill各向异性屈服条件及应力应变关系,我们得到高速扩展平面应力裂纹尖端的各向异性塑性场的一般解.将这个一般解用于四种各向异性特殊情形,我们就导出这四种特殊情形的一般解.最后,本文给出X=Y=Z情形的高速扩展平面应力Ⅰ型裂纹尖端的各向异性塑性场.     

简便积分方程法分析桩

本文用两种方法来分析桩受垂直载荷作用问题.一种是:将由Mindlin集中力组成的轴对称载荷沿弹性半空间z轴的[0,L]内分布,并迭加Boussinesq的解;另一种是:除上述诸虚载荷外,还将Mindlin的垂直集中力沿z轴的[0,L]内分布.前者使边界条件为: 的桩受垂直载荷问题归结为一个Fredholm第一种积分方程;后者使边界条件(其中1,3式同)(0.1)式中的2为:0≤z≤L,U(e,z)=a-e,(e→a);W(a,z)=常数(0.2)的桩受垂直载荷问题归结为两个联立的Fredholm第一种方程式.对刚性桩而言,前者适于容许桩和其侧面附着的土有相对滑动情况;后者适于无相对滑动情形.这两种方法较现有的虚载荷分布于桩表面的诸法具有下列优点:1.所得的积分方程不是二维、奇异的;而是一维、非奇异的.2.能考虑初应力的影响.第一种方法还无须预先假定沉陷函数W;在可压缩桩中容易考虑三维应力的影响的好处.本文还给出Fredholm第一种积分方程近似解误差估计的一个定理,以及两种方法用DJS—21机计算单桩沉陷的结果.     

精确计算Bingham流体圆管层流压降及速度分布规律的新方法

Bingham(宾汉)模型情况下,多采用通用公式进行圆管层流压降的解析计算,即将Bingham模型本构方程代入粘性流体圆管层流流动通用公式进行计算,仅能得到压降的解析解.新方法结合Bingham流体本构方程与运动方程,建立有关力学平衡方程,并运用代数方程的根式解理论对圆管层流流动时的非线性方程进行求解,可直接求得Bingham流体圆管层流压降及速度流核区半径的解析解,进一步可求得圆管层流速度解析解;Bingham流体圆管层流速度的直接影响因素为流量、塑性粘度和屈服值,研究发现速度流核宽度与屈服值成正比,与流量及塑性粘度成反比,且流核的宽度越大,流核区的速度越小.     

任意分布载荷下两对边简支单向非均匀变厚度矩形板的弯曲问题*

本文使用阶梯折算法[1,2],研究了任意分布载荷下两对边简支(x=0和x=a边)单向非均匀变厚度(D=D(y))矩形板的弯曲问题,得到了问题的阶梯近似解,文末对静水压力作用下四边简支线性变厚度正方形板的弯曲作了数值计算,所得结果与[3]一致,这样便验证了阶梯折算法的准确性.     

弹性—粘塑性材料I型动态扩展裂纹尖端场的渐近解

提出了一种新的弹性－粘塑性模型用于分析I型动态扩展裂纹尖端的应力应变场。给出了适当的位移模式，推导了渐近方程并且给出了数值解。分析和计算表明：对于低粘性情况，裂纹尖端场具有对数奇异性；对于高粘性情况，渐近方程无解。分析比较表明该结果具有高压臣提出的单参数解的所有优点，并且消除了粘性区随裂纹扩展而移动的不足。     

The diffusion of a discontinuity of the shear stress at the boundary of a viscoplastic half-plane

For the problem of the diffusion of a discontinuity of the shear stress at the boundary of a half-plane, which is a special case of the general problem of the diffusion of a vortex layer, the classes of media and types of assignment of boundary conditions for which self-similar solutions exist are discussed. For a viscoplastic medium in a half-plane the problem reduces to the problem in a layer of time-variable thickness, the solution of which does not possess the property of analyticity. The long-term asymptotic of this problem are investigated. In the case where, at an accessible boundary, it is possible simultaneously to measure both the shear stress and the horizontal velocity, an algorithm is proposed for finding a quantity that is difficult to measure, A namely, the thickness of the zone of viscoplastic flow.     

A fluid mechanical model for granular flow in silos

Granular materials may display both solid and fluid like behaviour. For low densities and high strain rates as in avalanches or during the discharge of silos the behaviour is mainly governed by interparticle collisions. On the other hand, frictional contacts characterise the solid state which is represented within the framework of plasticity theory. A fluid like constitutive model describes granular materials when subjected to large deformations and high strain rates. It bases upon a modified viscoplastic model that is valid for both yielded and unyielded regions. The central idea is the distinction between fluid and solid regions by means of comparing actual shear stress and Coulomb yield stress. The application to the simultion of the discharge of silos shows the feasibility of the chosen method. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The plane problem of the extrusion of a viscoplastic medium by parallel plates

A solution of the problem of the plane parallel flow of viscoplastic medium between two parallel plates when they approach (separate) at a specified velocity is given within the framework of the Bingham model in the inertialess thin-layer approximation for arbitrary values of the coefficient of viscosity and the yield stress. Analytic expressions are obtained for the velocity and pressure fields. The boundary of the flow kernel, where the shear stress on the areas of the parallel planes of the plates is less than the yield stress and the component of the velocity, parallel to the plates, does not change in a transverse direction, is determined. A single similarity parameter which defines the kinematic and dynamic flow characteristics is found. For a specified law of motion of the plates, a general expression is obtained for the force acting on plates of finite size in terms of a dimensionless function of a single dimensionless parameter. The law of approach (separation) of the plates under a constant force is found.     

Helical flow of a Bingham fluid arising from axial Poiseuille flow between coaxial cylinders

The Bingham fluid model represents viscoplastic materials that display yielding, that is, behave as a solid body at low stresses, but flow as a Newtonian fluid at high stresses. In any Bingham flow, there may be regions of solid material separated from regions of Newtonian flow by so-called yield boundaries. Such materials arise in a range of industrial applications. Here, we consider the helical flow of a Bingham fluid between infinitely long coaxial cylinders, where the flow arises from the imposition of a steady rotation of the inner cylinder (annular Coutte flow) on a steady axial pressure driven flow (Poiseuille flow), where the ratio of the rotational flow compared to the axial flow is small. We apply a perturbation procedure to obtain approximate analytic expressions for the fluid velocity field and such related quantities as the stress and viscosity profiles in the flow. In particular, we examine the location of yield boundaries in the flow and how these vary with the rotation speed of the inner cylinder and other flow parameters. These analytic results are shown to agree very well with the results of numerical computations.     

Nonsimilar boundary layer analysis of double-diffusive convection from a vertical truncated cone in a porous medium with variable viscosity

This work presents a boundary layer analysis about variable viscosity effects on the double-diffusive convection near a vertical truncated cone in a fluid-saturated porous medium with constant wall temperature and concentration. The viscosity of the fluid is assumed to be an inverse linear function of the temperature. A boundary layer analysis is employed to derive the nondimensional nonsimilar governing equations, and the transformed boundary layer governing equations are solved by the cubic spline collocation method to yield computationally efficient numerical solutions. The obtained results are found to be in good agreement with previous papers on special cases of the problem. Results for local Nusselt and Sherwood numbers are presented as functions of viscosity-variation parameter, buoyancy ratio, and Lewis number. For a porous medium saturated with a Newtonian fluid with viscosity proportional to an inverse linear function of temperature, higher value of viscosity-variation parameter leads to the decrease of the viscosity in fluid flow, thus increasing the fluid velocity as well as the local Nusselt number and the local Sherwood number.     

Natural convection boundary layer flow of fluid with temperature-dependent viscosity from a horizontal elliptical cylinder with constant surface heat flux

This study deals with the temperature-dependent viscosity effects on the natural convection boundary layer on a horizontal elliptical cylinder with constant surface heat flux. The mathematical problem is reduced to a pair of coupled partial differential equations for the temperature and the stream function, and the resulting nonlinear equations are solved numerically by cubic spline collocation method. Results for the heat transfer characteristics are presented as functions of eccentric angle for various values of viscosity variation parameters, Prandtl numbers and aspect ratios. Results show that an increase in the viscosity variation parameter tends to accelerate the fluid flow near the surface and increase the maximum velocity, thus decreasing the velocity boundary layer thickness. As the viscosity variation parameter is increased, the surface temperature tends to decrease, thus increasing the local Nusselt number. Moreover, the local Nusselt number of the elliptical cylinder increases as the Prandtl number of the fluid is increased.     

Leray–Hopf solutions to a viscoelastoplastic fluid model with nonsmooth stress–strain relation

We consider a fluid model including viscoelastic and viscoplastic effects. The state is given by the fluid velocity and an internal stress tensor that is transported along the flow with the Zaremba–Jaumann derivative. Moreover, the stress tensor obeys a nonlinear and nonsmooth dissipation law as well as stress diffusion. We prove the existence of global-in-time weak solutions satisfying an energy inequality under general Dirichlet conditions for the velocity field and Neumann conditions for the stress tensor.     

Pendant drop formation of shear-thinning and yield stress fluids

A volume-of-fluid numerical method is used to predict the dynamics of shear-thinning liquid drop formation in air from a circular orifice. The validity of the numerical calculation is confirmed for a Newtonian liquid by comparison with experimental measurements. For particular values of Weber number and Froude number, predictions show a more rapid pinch-off, and a reduced number of secondary droplets, with increasing shear-thinning. Also a minimum in the limiting drop length occurs for the smallest Weber number as the zero-shear viscosity is varied. At the highest viscosity, the drop length is reduced due to shear-thinning, whereas at lower viscosities there is little effect of shear-thinning. The evolution of predicted drop shape, drop thickness and length, and the configuration at pinch-off are discussed for shear-thinning drops. The evolution of a drop of Bingham yield stress liquid is also considered as a limiting case. In contrast to the shear-thinning cases, it exhibits a plug flow prior to necking, an almost step-change approach to pinch-off of a “torpedo” shaped drop following the onset of necking, and a much smaller neck length; no secondary drops are formed. The results demonstrate the potential of the numerical model as a design tool in tailoring the fluid rheology for controlling drop formation behaviour.     

Nonlinear mathematical analysis for blood flow in a constricted artery under periodic body acceleration

This study analyses the pulsatile flow of blood through mild stenosed narrow arteries, treating the blood in the core region as a Casson fluid and the plasma in the peripheral layer as a Newtonian fluid. Perturbation method is employed to solve the resulting coupled implicit system of non-linear partial differential equations. The expressions for shear stress, velocity, wall shear stress, plug core radius, flow rate and longitudinal impedance to flow are obtained. The effects of pulsatility, stenosis depth, peripheral layer thickness, body acceleration and non-Newtonian behavior of blood on these flow quantities are discussed. It is noted that the plug core radius, wall shear stress and longitudinal impedance to flow increase as the yield stress and stenosis depth increase and they decrease with the increase of the body acceleration, pressure gradient, width of the peripheral layer thickness. It is observed that the plug flow velocity and flow rate increase with the increase of the pulsatile Reynolds number, body acceleration, pressure gradient and the width of the peripheral layer thickness and the reverse behavior is found when the yield stress, stenosis depth and lead angle increase. It is also recorded that the wall shear stress and longitudinal impedance to flow are considerably lower for the two-fluid Casson model than that of the single-fluid Casson model. It is found that the presence of body acceleration and peripheral layer influences the mean flow rate and mean velocity by increasing their magnitude significantly in the arteries.     

Helical flow arising from the yielded annular flow of a Bingham fluid

The Bingham fluid model was developed to represent viscoplastic materials that change from rigid bodies at low stress to viscous fluids at high stress – a process termed yielding. Such a fluid model is used in the modeling of slurries, which occur frequently in food processing and other engineering applications.