Überlegungen zur Scherviskosität von Suspensionen aufgrund einer Dimensionsanalyse

Dimensional analysis of the motion of solid particles suspended in a fluid phase shows that the macroscopic relative shear viscosity of suspensions generally depends not only on the volume concentration and particle shape but also on two Reynolds numbers and a dimensionless sedimentation number. These dimensionless numbers are formed using parameters characterizing the structure and motion of the suspension at the microscopic level. The analysis was based on the assumptions that the dispersed particles are rigid and sufficiently large that Brownian motion may be neglected, that the continuous fluid phase is Newtonian and that the interactions between particles and between particles and fluid phase are only hydrodynamic. The Reynolds numbers describe the influence of the inertial forces at the microscopic level, and the sedimentation number the influence of gravity. The dimensionless numbers can be neglected if their values are much smaller than one. For each of the dimensionless numbers both the shear rate and the particle size influence the shear viscosity. Thus sedimentation number is large for low shear rates, whereas the Reynolds numbers are large for high shear rates. The viscosity function for one suspension can be transformed into the viscosity function for another suspension with geometrically similar particles but of a different size. The scale-up rules are derived from the requirement that the relevant dimensionless numbers must be constant. The influence of non-hydrodynamic effects at the microscopic level on the shear viscosity can be detected by deviations from the derived scale-up rules.     

Casson Fluid Flow Due to Non-Coaxial Rotation of a Porous Disk and the Fluid at Infinity Through a Porous Medium

The flow of the Casson fluid due to non-coaxial rotation of a disk and the fluid at infinity is investigated. Partial differential equations are made dimensionless and coupled. The exact solution of the resultant nonlinear initial-boundary-value problem is solved by applying the Laplace transform. The shear stresses at the disk surface and the steady state stresses are computed. The effects of dimensionless parameters on the dimensionless primary and secondary velocities are analyzed.     

Unsteady motion of a non-linear viscoelastic fluid

In this paper, we study the unsteady flow of a generalized second grade fluid. Specifically, we solve numerically the linear momentum equations for the flow of this viscoelastic shear-thinning (shear-thickening) fluid surrounding a solid cylindrical rod that is suddenly set into longitudinal and torsional motion. The equations are made dimensionless. The results are presented for the shear stresses at the wall, related to the drag force; these are physical quantities of interest, especially in oil-drilling applications.     

Frottements et pertes de pression des fluides non newtoniens dans des conduites non circulaires

The study of fluid flow in a duct requires characteristic parameters of the flow and dimensionless numbers to correlate and compare experimental results. For Newtonian fluids in simple configurations, the definition of the Reynolds number is quite standard, but for non-Newtonian fluid flows in ducts with arbitrary shape of cross section, the dependence of the apparent viscosity with the shear rate requires a generalization of this dimensionless number. This note proposes a general method valid for a large class of non-Newtonian fluids and for all duct shapes. An application is developed for a viscoelastic flow through a rectangular duct. Results obtained in the present investigation are in a good agreement with available correlations. To cite this article: M. Mahfoud et al., C. R. Mecanique 333 (2005).     

粘弹性地基上粘弹性输流管道的稳定性分析

从Winkler假设和单轴线性粘弹性本构方程出发,推导了Kelvin-Voigt粘弹性地基上三参量固体模型输流管道的运动微分方程,采用改进的有限差分法,分析了管道和地基的粘弹性参数对输流管道无量纲复频率和无量纲流速之间的变化关系的影响。     

纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.      

Large-amplitude oscillatory shear flow from the corotational Maxwell model

Using the single relaxation time corotational Maxwell fluid, we derive explicit analytical expressions for the first, third, and fifth harmonics of the alternating shear stress response in large-amplitude oscillatory shear (LAOS). We also derive corresponding expressions for the zeroth, second, and fourth harmonics of both the first and second normal stress differences. These harmonics are found to depend upon just two dimensionless groups: the Deborah and Weissenberg numbers, each of which causes non-Newtonian behavior. The form of the solution for the corotational Maxwell model in LAOS matches the forms of the analytical solutions for two molecular models for dilute solutions and one for concentrated solutions or melts. We also derive an analytical solution for the corotational Maxwell model after startup of LAOS. For this we find that both small and large amplitude cases approach a periodic limit cycle (alternance) at the same rate for both the shear stress response and for the normal stress differences. For molten high density polyethylene that is lightly filled with carbon black, we find good quantitative agreement with measured LAOS behavior when our analytical solution is superposed for multiple relaxation times.     

Two-fluid modeling in analyzing the interfacial stability of liquid film flows

The area-averaged two-fluid model formulation of a separated two-phase flow system is used to investigate interfacial stability of liquid film flows. The analysis takes into account the effects of phase change at the interface as well as the dynamic effects of the adjacent vapor flow on the interfacial stability. Wave formation and instability criteria are established in terms of the generalized fluid and flow parameters. The criteria are applied to investigate the stability of laminar liquid film flow with interfacial shear and phase change. The influence of various dimensionless parameters characterizing film thickness, gravity, phase change and interfacial shear are studied with respect to the neutral stability, temporal growth factor and the wave propagation velocity. The results of the present study indicate that the interfacial stability analysis developed within the frame of the two-fluid model formulation proves to be quite accurate as judged by comparing its results with the available experimental data and with the results of much longer and more complex analytical investigations which are valid only for the liquid film free of interfacial shear.     

Die swell measurements of second‐order fluids: numerical experiments

An analysis of the flow of a second‐order fluid is presented. Reference values for some variables are defined, and with these a non‐dimensional formulation of the governing equations. From this formulation, three dimensionless numbers appear; one is the Reynolds number, and two numbers that are called the first‐ and second‐dimensionless normal stress (NSD) coefficients. The equations of motion are solved by a finite element method using a commercially available program (Fidap), and the steady state converged solution was used to measure the die swell. The factors that influence die swell and that are studied in this work include: the die geometry for circular cross sectional dies, including tubular, converging, diverging, half‐converging/half‐tubular shapes; fluid characteristics such as Reynolds number and first‐ and second‐DNS coefficients (both positive and negative values); and flow rates, as determined by the maximum velocity in a parabolic velocity profile at the entrance to the die. The results suggest that shear and deformation histories of the fluid directly influence not only swell characteristics, but also convergence characteristics of the numerical simulation. © 1999 John Wiley & Sons, Ltd.     

页岩气高效开采的可压裂度和射孔簇间距预测

页岩储层的可压裂性是影响页岩气产量的关键因素. 本文基于断裂力学理论, 以高围压下岩石层理弱面的剪切破坏为主要研究对象, 依据岩石抗拉强度和层理弱面抗剪强度的比值关系, 首先提出了可压裂度的概念, 给出了无量纲的定性曲线图, 涵盖了岩石脆性矿物质含量, 粘性主导和韧性主导裂缝尖端流体压强、射孔簇分布间距的综合地质与工程因素. 接着提出了一种新的表征高围压下页岩可压裂度的无量纲参数, 在保证达到充分解吸附的最小压裂间距前提下, 依据该参数可计算水平井压裂中的射孔簇间距, 可作为工程参考指标. 本文将断裂力学理论结合水力压裂高效开采页岩气工程, 具有力学理论意义和工程应用前景.      

Poiseuille Flow of a Third Grade Fluid in a Porous Medium

Effects of porous medium have been investigated on the steady flow of a third grade fluid between two stationary porous plates. The continuity and momentum equations along with modified Darcy??s law are used for the development of mathematical problem. The governing nonlinear problem is solved by a homotopy analysis method. The dimensionless velocity and shear stresses at the plates are analyzed.     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

Rotating electroosmotic flow of an Eyring fluid

A perturbation analysis is presented in this paper for the electroosmotic(EO) flow of an Eyring fluid through a wide rectangular microchannel that rotates about an axis perpendicular to its own. Mildly shear-thinning rheology is assumed such that at the leading order the problem reduces to that of Newtonian EO flow in a rotating channel, while the shear thinning effect shows up in a higher-order problem.Using the relaxation time as the small ordering parameter,analytical solutions are deduced for the leading-as well as first-order problems in terms of the dimensionless Debye and rotation parameters. The velocity profiles of the Ekman–electric double layer(EDL) layer, which is the boundary layer that arises when the Ekman layer and the EDL are comparably thin, are also deduced for an Eyring fluid. It is shown that the present perturbation model can yield results that are close to the exact solutions even when the ordering parameter is as large as order unity. By this order of the relaxation time parameter, the enhancing effect on the rotating EO flow due to shear-thinning Eyring rheology can be significant.     

Self-consistent particle simulation of model-stabilized colloidal suspensions

Dynamics of model-stabilized colloidal suspensions were investigated by the self-consistent particle simulation method (SC), a new simulation algorithm that takes into account the interaction between the particles and suspending fluid. In this method, the fluid-particle interaction is introduced self-consistently by combining the finite element method (FEM) for fluid motion with Brownian dynamics (BD) for particle dynamics. To validate the reliability of the proposed algorithm, the shear dynamics of the stable particle suspensions were investigated. Relative viscosity and microstructure as a function of dimensionless shear rate at different volume fractions were in good agreement with previous observations. The robustness of the method was also verified through numerical convergence test. The effect of the fluid-particle interaction was well represented in simulations of two model problems, pressure-driven channel flow and rotating Couette flow. Plug-shaped velocity profile was observed in pressure-driven channel flow, which arised from shear thinning behavior of the stable suspension. In rotating Couette flow, shear banded nonlinear flow profile was observed. Although full hydrodynamic interaction (HI) was not rigorously taken into account, it successfully captured the macroscopic structure-induced flow field. It also takes advantage of the geometrical adaptability of FEM and computational efficiency of BD. We expect this newly developed simulation platform to be useful and efficient for probing the complex flow dynamics of particle systems as well as for practical applications in the complex flow of complex fluids.     

Steady diffusion of an ideal fluid through a pre-stressed and reinforced hollow conduit subjected to combined finite deformations

A problem dealing with the radial steady diffusion of an ideal fluid through a pre-stressed fibrous hollow cylinder subjected to finite deformations is investigated. Such a problem has relevance to several technical problems such as: (a) improving the method for performing prosthesis conduit for use with living tissue, (b) more understanding the problem of atherogenesis, and (c) ultra filtration process. The numerical results show that the presence of a pre-stress reduces considerably the sensitivity of the dimensionless trans-mural pressure to imposed dimensionless flux. This effect is confirmed with respect to the variation of the circular and axial shear.     

Convective mass transfer across fluid interfaces in straight angular pores

Steady convective mass transfer to or from fluid interfaces in pores of angular cross-section is theoretically investigated. This situation is relevant to a variety of mass transport process in porous media, including the fate of residual non-aqueous phase liquid ganglia and gas bubbles. The model incorporates the essential physics of capillarity and solute mass transfer by convection and diffusion in corner fluid filaments. The geometry of the corner filaments, characterized by the fluid–fluid contact angle, the corner half-angle and the interface meniscus curvature, is accounted for. Boundary conditions of zero surface shear (‘perfect-slip’) and infinite surface shear (‘no-slip’) at the fluid–fluid interface are considered. The governing equations for laminar flow within the corner filament and convective diffusion to or from the fluid–fluid interface are solved using finite-element methods. Flow computations are verified by comparing the dimensionless resistance factor and hydraulic conductance of corner filaments against recent numerical solutions by Patzek and Kristensen (J. Colloid Interface Sci 236, 305–317 2001). Novel results are obtained for the average effluent concentration as a function of flow geometry and pore-scale Peclet number. These results are correlated to a characteristic corner length and local pore-scale Peclet number using empirical equations appropriate for implementation in pore network models. Finally, a previously published “2D-slit” approximation to the problem at hand is checked and found to be in considerable error.     

Heat transfer,development length,and friction factor correlations for the asymptotic region of a laminar arc constrictor

A rigorous equilibrium model is formulated for the laminar flow of plasma through the heating region of a constricted arc plasma generator. The resulting equations are solved for an argon gas using an implicit finite-difference scheme, and the wall shear stress and heat transfer are computed from momentum and energy balances for a wide range of flow rates and arc currents. Friction factors based on this method are in agreement with empirical values to within approximately 10 per cent, and heat transfer calculations agree with the orders of magnitude predicted by a simplified arc model. A dimensionless number peculiar to the constricted arc, the so-called Ohmic heating parameter, is used to correlate the mean Nusselt numbers in the asymptotic arc regime. When the fluid properties are based on the mean temperature, the product of friction factor and Reynolds number is found to be a constant for currents below 200 Amp and to vary with the Ohmic heating parameter for larger currents. Correlations are also presented for the thermal development length, and the agreement with experiment is good.     

Simulation of heat and chemical reactions on Reiner Rivlin fluid model for blood flow through a tapered artery with a stenosis

In the present article, we have analyzed the effects of heat and mass transfer on Reiner Rivlin fluid model for blood flow through a tapered artery with a stenosis. The constitutive equations for a Reiner Rivlin fluid have been modelled in cylindrical coordinates. A perturbation series in dimensionless Reiner Rivlin fluid parameter (λ ₁ ≪ 1) have been used to obtain explicit forms for the velocity, temperature, concentration, resistance impedance, wall shear stress and shearing stress at the stenosis throat. The graphical results of different type of tapered arteries i.e. converging tapering, diverging tapering, non-tapered artery have been examined for different parameters of interest.     

Numerical simulation of chemically reactive Powell-Eyring liquid flow with double diffusive Cattaneo-Christov heat and mass flux theories

numerical study is reported for two-dimensional flow of an incompressible Powell-Eyring fluid by stretching the surface with the Cattaneo-Christov model of heat diffusion. Impacts of heat generation/absorption and destructive/generative chemical reactions are considered. Use of proper variables leads to a system of non-linear dimensionless expressions. Velocity, temperature and concentration profiles are achieved through a finite difference based algorithm with a successive over-relaxation (SOR) method. Emerging dimensionless quantities are described with graphs and tables. The temperature and concentration profiles decay due to enhancement in fluid parameters and Deborah numbers.     

功能梯度石墨烯增强复合材料板的自由振动

基于高阶函数剪切变形板理论,研究了可移动简支边界条件下功能梯度石墨烯增强复合材料板的自由振动。利用修正的Halpin-Tsai模型估计了复合材料的等效杨氏模量。基于Hamilton原理构建其偏微分方程。通过Navier方法计算得到了板的量纲为一的频率,并与已有成果进行了比较,验证了本文计算模型的优越性。数值计算结果表明:功能梯度石墨烯增强复合材料板的量纲为一的基频受石墨烯小块的分布模式影响很大,随着石墨烯含量的增加而单调增加;固定石墨烯小块的厚度和长宽比,量纲为一的基频随着长厚比的增加而单调增加。