Abstract

Many electrorheological fluids (ERFs) as fluids with microstructure demonstrate viscoplastic behaviours. Rheometric measurements indicate that some flows of these fluids may be modelled as the flows of a Herschel-Bulkley fluid. In this paper, the flow of a Herschel-Bulkley ER fluid—with a fractional power-law exponent—in a narrow clearance between two fixed surfaces of revolution with common axis of symmetry is considered. The flow is externally pressurized, and it is considered with inertia effect. In order to solve this problem, the boundary layer equations are used. The influence of inertia forces on the pressure distribution is examined by using the method of averaged inertia terms of the momentum equation. Numerical examples of externally pressurized ERFs flows in the clearance between parallel disks and concentric spherical surfaces are presented.

1. Introduction

In recent years, the study of fluids with microstructures has gained much importance because of its numerous applications in various engineering disciplines such as chemical engineering, polymer processing, plastic forming foundry engineering, and engineering of lubrication [114].

In machines and mechanisms systems of many industrial processes, the phenomena of a flow of viscoplastic fluids are used. One of these phenomena is a slide bearing lubrication [9, 10, 13]. Advances in technology and severe operational requirements of machines necessitated the development of improved lubricants to ensure a smooth and safe operation. Generally, viscosity of lubricants decreases with temperature. For operations under high speeds and heavy loads, oils containing high molecular weight polymers as viscosity index improvers are used to increase a load carrying capacity of the modified lubricants [9, 13].

Most substances used in the lubrication technology are polymer solutions, thus, the characteristics of the bearings change when such rheological substances, known as non-Newtonian fluids, are used as lubricants. Several constitutive relations applied were used to model the non-Newtonian characteristics exhibited by some lubricants [7, 11, 13, 15, 16].

Another ones of these phenomena are processes of vibration control and torque transmission. In the last years, the electrorheological fluids (abbreviated to ERFs) have acquired a great relevance for supporting vibration control and torque transmission devices, based on the characteristic dependence of their viscosity on applied electric field strength. Since their initial discovery by Winslow [17], many particle-dispersion electrorheological fluids, consisting of dielectric particles dispersed in insulating oil, have been reported. When the external electric field is imposed to an ERF, it behaves as a viscoplastic fluids [9, 1825], displaying a field-dependent yield shear stress which is widely variable. Without the electric field, the ERF has a reversible and a constant viscosity so that it flows as a Newtonian fluid. Another salient feature of the ERF is that the time required for the viscosity variation is very short (< 0.001 sec). These attractive characteristics of the ERF provide the possibility of the appearance of new engineering technology. Recently, the application of the ERF to rotor-bearing systems has also been initiated by Basavaraja et al. [2], Jung and Choi [18], Dimarogonas and Kollias [26], and Peng and Zhu [27] and vibration control by Lee and Wereley [23], Choi et al. [24], El Wahed et al. [28], and Lim et al. [29].

To describe the rheological behaviour of viscoplastic fluids in complex geometries, the Bingham model is used [8, 13, 15, 16, 21]. Recently, the nonlinear model of Shulman et al. [21] has been successfully applied [8, 13, 15, 16, 20]. The constitutive equation for this model is given as follows: where is the shear stress, is the yield shear stress, is the coefficient of plastic viscosity, is the shear strain rate, and and are the power-law exponents. By reducing the coefficients in the Shulman equation one can obtain simpler models describing the flow of a viscoplastic fluid.

The most popular model of the ERF is the Bingham model for which in (1). Recently, the nonlinear model of Herschel-Bulkley has been successfully used to describe the viscoplastic behaviour of nonNewtonian ERFs. The constitutive equation for this model is as follows [5, 30, 31]:

The purpose of this paper is a study of pressure distribution in a flow of the Herschel-Bulkley ERF—with a fractional power-law exponent [8]—in the clearance modelled as a narrow space between two surfaces of revolution shown in Figure 1. Using the method of averaged inertia [32] (similar method is used in [33]), the influence of inertia terms of the equations of motion and viscoplastic behaviour on this distribution is analysed. To solve the problem we will use the lubrication approximation theory to the flow in a narrow clearance [8, 10, 15, 16, 25].

2. Equations of Motion of the ERF

The yield shear stress for the ERF varies with respect to the electric field. According to the experimental results reported in Shulman et al. [21], Shulman and Nosov [22], Jordan and Shaw [34], and Otsubo [35], the relation between the yield shear stress and the electric field strength is given as follows: where is the magnitude of the vector but and are the applied voltage and the film thickness, respectively. Both parameters and are the experimental constants of which the range of the exponent is 1 to 2.4 (Jordan and Shaw [34], Otsubo [35]). Other experimental data (Whittle et al. [36]) suggest that the yield shear stress for some ERFs is as follows: where and are constants; for relatively high field strengths a simpler formula may be used where is an experimental constant.

The three-dimensional constitutive equation for the Herschel-Bulkley ERF has the form [15, 16, 25]: where 1 is the unit tensor, is the pressure, is second invariant of the stretching tensor (the first Rivlin-Ericksen tensor) defined by and v is velocity vector.

The general equations of motion of the Herschel-Bulkley fluid have a form:(i)equation of continuity: (ii)equation of momentum or where here is a coefficient of viscoplasticity.

Let us consider the Herschel-Bulkley ERF in a clearance between fixed surfaces of revolution. The flow configuration is shown in Figure 1. The surfaces of revolution are defined by function which denotes the radius of the median between the surfaces, plus function which denotes the distance to each surfaces from the median, measured along a normal to the median. An intrinsic curvilinear orthogonal coordinate system is also depicted in Figure 1.

The physical parameters of the lubricant flow are the velocity components , , and pressure . With regard to the axial symmetry of the flow these parameters are not dependent on the angle .

The assumption typical for the flow in a narrow clearance [2, 8, 15, 16] can be used to make order-of-magnitude arguments for (8) (11).

A further simplification comes by noting that—in accordance to the lubrication approximation—the most important changes in an annular channel occur in the normal (to the channel median) direction. This leads to the assumption that the flow is nearly parallel to the surfaces bounding the clearance, so that in an intrinsic coordinate system.

If some asymptotic transformations are made, the same as in (Falicki [8], Walicka [15, 16]), these equations can be reduced to a simpler form: where and . The signum function takes on the value for a positive argument and for a negative argument.

The order-of-magnitude arguments show that is a function of only.

For a majority of greases, molten polymers, mush metals, and ERFs, the values of a yield shear stress are contained in the limits: , where and denotes the order of magnitude, but the exponent is equal to: or , and it is not an integer number (Falicki [8], Roussel [12]).

In the flow of a fluid with the yield shear stress, there exists a quasi-solid core flow bounded by surfaces laying at Combining the expressions (16) and (18), one obtains the boundary conditions on liquid surfaces of the core flow as The boundary conditions on solid surfaces are stated as follows: On the median line there is also Moreover, in the inlet and the outlet of the clearance conditions for the pressure can be written in the form where denotes the inlet coordinate and —outlet coordinate.

3. Solution to the Equation of Motion

Taking into account (14), one can rearrange (15) writing it in the form Then, averaging the left-hand side of (23) across the clearance thickness, solving the obtained equation and taking into account the boundary conditions (20)–(22) we get the following result for the pressure distribution [8, 15, 16, 3739]: Functions and depend upon the solution of the equation characteristic to the flow of a Herschel-Bulkley fluid in the clearance; this equation has a form [8, 16, 37, 39]: where and is the flow rate.

Its analytical solution exists only for large and for small (Falicki [8], Walicka [15, 16, 37], and Walicki and Walicka [39]):

For intermediate values of the solutions of (25) may be found only numerical methods. Figures 2 and 3 illustrate the progression of solutions versus for two values of [8]. Note that for , the above values are identical with those obtained by Covey and Stanmore [40].

The functions and are given by expressions where Note that the final form of and depend—for —on the values of and :(i) for large values of and exponent being a natural number there are where (ii) for large values of and exponent being a rational number the formula for is similar but takes the form where and denotes an entire part of ;(iii) for small values of there are where For (the coefficient but ) there are [8](iv) for large values of and exponent being a natural number: where (v) for large values of and exponent being a rational number the formula for is similar and takes the form where (vi) for small values of there are where The prim denotes the differentiation with respect to .

4. Examples of Flow and Graphic Presentation of Some Results

Taking into account the results obtained in the previous section we will present the pressure distribution in the clearance of constant thickness between two parallel disks as shown in Figure 4; to this aim introduce the following dimensionless parameters: if , then and here is the modified Reynolds number, and is the ER plasticity number (de Saint-Venant number). Note that large values of the ER plasticity number correspond to the flow with small core (small ).

Nondimensional formula for pressure distribution in ER flow of the Herschel-Bulkley fluid has the form this formula may be used to model the pressure distribution in a footstep bearing.

Note that for , all the above formulae represent the Bingham ERF flow. The nondimensional pressure distributions for the Herschel-Bulkley ERF flow in a clearance between two parallel disks are presented in Figures 5, 6, 7, 8, and 9. The values of nondimensional parameters such as , , and , are defined on the experimental data taken from the measurements presented in [8, 12].

Let us consider now the pressure distribution in the clearance of constant thickness between two concentric spheres surfaces shown in Figure 10. To this aim introduce the following dimensionless parameters: Note that if then and . The nondimensional pressure distributions are also given by formula (41); at present it models the pressure distribution in a spherical bearing.

The pressure distributions for the Herschel-Bulkley ERF flow in a clearance between concentric spherical surfaces are presented in Figures 11, 12, 13, and 14.

5. Conclusions

From the general considerations, formulae and graphs presented here for the Herschel-Bulkley ERF flows in the narrow clearance of constant thickness between parallel disks and concentric spherical surfaces shown in Figures 4 and 10 one may conclude that the pressure values(i)decrease with the increase of the modified Reynolds number ,(ii)increase with the decrease of the nonlinearity index ,(iii)increase with the decrease of the de Saint-Venant ER number ,(iv)are larger between concentric spherical surfaces than these ones between parallel disks for being natural numbers and large values of ,(v)for being rational numbers it is inversely.

For small values of the de Saint-Venant ER number the influence of the Reynolds number on the pressure values is inconsiderable. Therefore, for these values of the Herschel-Bulkley ERF flows may be considered without inertia effects. It is seen from the results obtained here, that the pressure values are larger for the flows between concentric spherical surfaces than these values for the flows between parallel disks.

Generally, it may be concluded that for the Herschel-Bulkley ERF flows in the clearance between two surfaces of revolution the pressure values(i)are larger for the flow in the clearances with curvilinear generating lines than these ones for the flow in clearances with rectilinear generating lines for the nonlinearity index being natural numbers; for being rational numbers this phenomenon progresses inversely.

Note that the results obtained here for and being natural numbers are similar to those of the work by Walicki and Walicka [39].

The formulae for the pressure distribution obtained previously may be used to model its distribution in slide thrust bearings—of arbitrary curvilinear shapes—lubricated by ER lubricants.

Nomenclature

:Auxiliary  function  defined  by  formulae
(30)1  or  (36)1
: Auxiliary  function  defined  by  formula
(30)2
: Auxiliary  function  defined  by  formula
(28)2  or   , , , and   or   and ,   respectively
: Half of the fluid film thickness
: Auxiliary  function  defined  by  formula
(30)8
: Auxiliary  function  defined  by  formula
(30)6
: Auxiliary  function  defined  by  formula
(35)3
: Plasticity number
: De Saint-Venant plasticity number
: Power-law index
: Pressure
: Flow rate
: Local radius
: Inlet radius to the bearing clearance
: Outlet radius from the bearing clearance
: Modified Reynolds number
: Signum function
: Auxiliary  function  defined  by  formulae
(28)1, ,  and   or and , respectively
: Component of the shear stress tensor
: Applied voltage
: Velocity components
: ER experimental constants
: Fluid density.
: Central angle of spherical surface
: Shear stress
: Yield shear stress.