#### Abstract

Squeeze film dampers are novel rotor dynamic devices used to alleviate small amplitude, large force vibrations and are used in conjunction with antifriction bearings in aircraft jet engine bearings to provide external damping as these possess very little inherent damping. Electrorheological (ER) fluids are controllable fluids in which the rheological properties of the fluid, particularly viscosity, can be controlled in accordance with the requirements of the rotor dynamic system by controlling the intensity of the applied electric field and this property can be utilized in squeeze film dampers, to provide variable stiffness and damping at a particular excitation frequency. The paper investigates the effect of temperature and electric field on the apparent viscosity and dynamic (stiffness and damping characteristics) of ER fluid (suspension of diatomite in transformer oil) using the available literature. These characteristics increase with the field as the viscosity increases with the field. However, these characteristics decrease with increase in temperature and shear strain rate as the viscosity of the fluid decreases with temperature and shear strain rate. The temperature is an important parameter as the aircraft jet engine rotors are located in a zone of high temperature gradients and the damper fluid is susceptible to large variations in temperature.

#### 1. Introduction

The application of controllable fluids in engineering is an area which began to be explored in the early 1980s and fluids of interest in this study are materials that respond to an applied electric field with a remarkable change in its rheological behavior. The essential characteristic of these controllable fluids is their capability to change from a liquid state into a solid-like gel under the action of an electric field. These materials are commonly called electrorheological (ER) fluids and are normally low-viscosity colloidal suspensions. When an electric field is applied, the fluid undergoes a seemingly reversible transition to a solid in milliseconds, being able to support considerable stress under static load without yield. Willis M. Winslow (1949) was the first to discover an ER fluid and introduced the concept of controlling the apparent viscosity of an electroviscous fluid by using an electric field.

ER fluids usually consist of a base fluid, usually some type of low viscosity insulating liquid, mixed with non-conductive particles (in the range of 1–10 *μ*m diameter), such as cornstarch, gypsum and lithium salt.

Since its mechanical properties can be easily controlled within a wide range by applied voltage (almost from a pure liquid to a solid), the ER fluid could be used as an electrical and mechanical interface in various industrial areas. ER fluids can be utilized in automobile clutch, brake, and damping systems. In a vibratory damping system, the fluid is to be of variable viscosity which is best satisfied by the application of smart fluids (ER and MR fluid). The change in the field induced shear stress produces a change in their apparent viscosity and this property is successfully utilized in applications involving valves, dampers, clutches, and adaptive structures but applications in rotor dynamics are far from satisfactory and require considerable attention [1–3]. Limited information is available regarding the stiffness and damping characteristics of ER fluid and their applications in active journal bearings and squeeze film dampers. While general aspects and applications of MR and ER fluids have been dealt with to some extent [4, 5], sharing their promising characteristics, it was only in the last decade that some pieces of literature have appeared on rotor dampers [6, 7]. Since the available literature regarding the effect of temperature on the stiffness and damping characteristics of ER Fluid squeeze film dampers is very scanty, this paper investigates these aspects on a short squeeze film damper using the relation obtained between viscosity and temperature for a particular ER fluid which is a suspension of diatomite in transformer oil.

#### 2. The Herschel-Bulkley Model

The one-dimensional Herschel-Bulkley model is used to model the shear thinning of ER fluid in the post-yield region. The constitutive equation for this model [9] iswhere is the shear stress, is the shear strain rate, and is the dynamic yield stress. and are two parameters that describe the magnitude and slope of the curve. The parameter is sometimes referred to as the “consistency index.” The exponent is dimensionless and is called the “behavior index.” Both and depend on material properties and working conditions such as friction and temperature.

For steady shear flow , assume that , (1b) can be written as

The apparent viscosity thus is For shear thinning of ER suspensions with , the apparent viscosity decreases with the increment of shear rate. The temperature affects both the values and and is responsible for the variation in the apparent viscosity of the ER fluid.

#### 3. Temperature Effect on Viscosity of ER Fluid

Changes in the temperature of ER fluid due to working conditions like friction, or electrical power dissipation, may influence its performance in ER devices. The temperature decreases the viscosity of the ER fluid and consequently affects both the static and dynamic characteristics of the squeeze film damper. The viscosity-temperature relation of the ER Fluid at a particular electric field intensity can be modeled by the exponential decay type Reynolds equation or the Arrhenius equation outlined as (4) and (5), respectively, as where is the temperature and and are coefficients. This is an empirical model that usually works for a limited range of temperatures where is the activation energy and is the universal gas constant. A first-order fluid is another name for a power law fluid with exponential dependence of viscosity on temperature.

This can be also be presented as where .

#### 4. Stiffness and Damping Characteristics of Squeeze Film Dampers with Orbital Motion

These squeeze film dampers (Figure 2) find application in the air craft engine rotors and are usually mounted on rolling element bearings, with little or no inherent damping to fight the unbalance response, especially at critical speeds, and to eliminate rotor dynamic instabilities. The coordinates of the orbit are presented as [10] Ignoring the stiffness terms, the radial and tangential forces acting on the system are The direct damping coefficient is and the cross-coupled damping coefficient is The radial component of the force is considered as a stiffness force defined as the direct stiffness coefficient given by The squeeze film velocity is The Reynolds equation for the present case with zero spin velocity for the short bearing case reduces to The governing Reynolds equation for the short bearing case with zero spin velocity is

The short bearing approximation is applicable in cases in which the end seals are absent. For the short bearing approximation, the journal axis remains parallel to the bearing under the steady state operating condition. In such a case, the change in pressure in the circumferential direction is negligible in relation to the pressure change in axial direction, The Reynolds equation modifies to The pressure distribution in the film is obtained by integration as The constants in the above equation are evaluated from the boundary conditions

where the coordinate system of Figure 1 is reckoned midway along the length of the bearing (perpendicular to the plane of the paper). The two constants are evaluated to give the pressure, , as Using the film thickness Equation and ignoring the pressure distribution in the diverging wedge, that is, using the half Sommerfeld condition, one getsConsidering the equilibrium of the journal, the load capacity can be obtained from the integration of the following equations: where The load capacity, , is The radial and the tangential force generated in the damper are components of the load capacity and are defined as The direct stiffness, direct damping, and cross-coupled damping coefficients are, respectively,

#### 5. Methodology Proposed

The effect of electric field intensity and temperature on the theoretical squeeze film dynamic (damping and stiffness) characteristics can be obtained using the viscosity-temperature relations [8] using the damper specifications and relevant ER fluid data outlined in Tables 1, 2(a), and 2(b), respectively.

#### 6. Results and Discussions

The theoretical results are discussed under two headings:(a)effect of temperature on the viscosity and dynamic characteristics (stiffness and damping characteristics),(b)effect of electric field intensity and shear strain rate on the viscosity and dynamic characteristics (stiffness and damping characteristics).

##### 6.1. Effect of Temperature on the Viscosity, Stiffness, and Damping Characteristics

This section presents the effect of temperature on the viscosity and dynamic characteristics of the squeeze film damper for the specifications presented in Tables 1, 2(a), and 2(b).

The effect of temperature on the viscosity for different strain rates without the electric field is presented in Figure 2. The viscosity decreases with the temperature owing to loosening of the fibrillar structures with temperature.

Figure 3 indicates the variation of viscosity with temperature with the applied electric field at 2.5 kV/mm. The viscosity is also found to decrease with the strain rate as indicated by Figures 2 and 3 since the ER Fluid is a shear thinning fluid. However, it should be noticed that the viscosity increases by many orders under the action of the field.

The constants and of the temperature- viscosity model of the ER fluid (using either the Reynolds equation or the Arrhenius equation as outlined in (4) and (5), resp.) are presented in Table 3.

The effect of temperature on the direct stiffness, direct damping, and cross-coupled damping coefficients are presented through Figures 4(a), 4(b), 5(a), 5(b), 6(a), and 6(b) for electric field intensities of 0 kV/mm and 2.5 kV/mm. The dynamic coefficients decrease with temperature in accordance with (4) and (5), respectively.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

##### 6.2. Effect of Electric Field Intensity and Shear Strain Rate on the Viscosity, Stiffness, and Damping Characteristics of ER Squeeze Film Dampers

The apparent viscosity of the ER fluid is depicted in Figures 7(a) and 7(b) as a function of the shear strain rate for different electric field intensities, kV/mm and 2.5 kV/mm, respectively. The viscosity increases by many orders on the application of the field and decreases with increase in the shear strain rate as envisaged in these diagrams.

**(a)**

**(b)**

The effect of electric field on the dynamic characteristics of the ER damper are presented for temperatures 20°C and 60°C with and without the electric field through Figures 8, 9, and 10 as a function of the shear strain rate. All the dynamic characteristics show identical behavior and increase with the electric field through Figures 8, 9, and 10, respectively as is apparent from these plots. However, these characteristics decrease with temperature as can be seen through Figures 8, 9, and 10, respectively, for 20°C and 60°C temperatures. However, the dynamic characteristics decrease with increase in the strain rate in accordance with the viscosity behavior (shear thinning fluid).

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

#### 7. Conclusions

The effect of temperature on the dynamic (stiffness and damping) characteristics of an ER Fluid short squeeze film damper is presented in this paper. The dynamic coefficients increase with the electric field due to increase in viscosity with the electric field since it produces fibrillar chains or layers which retard the flow of the fluid and consequently increases the viscosity. However, the opposite effect happens at high temperatures and the fibrillar structures are broken at this temperature resulting in decreased viscosity. These coefficients, thus, decrease with increase in temperature of the ER fluid. They also decrease with the shear strain rate as the viscosity decreases with shear strain rate. For a particular damper configuration, these coefficients increase with decrease in clearance, increase in ratio and eccentricity ratio.

#### Nomenclature

: | Total shear stress, N/m^{2} (Pa) |

: | Dynamic yield stress, N/m^{2} (Pa) |

: | Field independent plastic viscosity, Ns/m^{2} |

: | Excitation frequency, rad/s |

: | Density, kg/m^{3} |

: | Attitude angle, degree |

: | Angular coordinate of the damper, degree |

: | Shear strain rate, s^{−1} |

: | Apparent viscosity due to the field, Pas |

: | Clearance in the damper, m |

: | Short bearing radial damping coefficient, MNs/m |

: | Short bearing cross-coupled damping coefficient, MNs/m |

: | Short bearing tangential or direct damping coefficient, MNs/m |

: | Damper diameter, m |

: | Eccentricity or whirl radius, m |

: | Electric field intensity, kV/mm |

: | Radial force in the damper, N |

: | Tangential force, N |

: | Consistency index and |

: | Short bearing direct stiffness coefficient, GN/m |

: | Length of the damper, m |

: | Fluid parameters and |

: | Eccentricity ratio = |

: | Viscosity exponent or behavior index |

: | Damper radius, m |

: | Signum function |

: | Time, s |

: | Surface velocity, m/s |

: | Radial whirl velocity, m/s |

: | Tangential whirl velocity, m/s |

, , and : | Coordinate axes. |

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.