# Squeezing Force of the Magnetorheological Fluid Isolating Damper for Centrifugal Fan in Nuclear Power Plant

Authors: Jin Huang; Ping Wang; Guochao Wang
Journal: Science and Technology of Nuclear Installations
 (2012)

---

## Abstract

Magnetorheological (MR) disk-type isolating dampers are the semi-active control devices that use MR fluids to produce controllable squeezing force. In this paper, the analytical endeavor into the fluid dynamic modeling of an MR isolating damper is reported. The velocity and pressure distribution of an MR fluid operating in an axisymmetric squeeze model are analytically solved using a biviscosity constitutive model. Analytical solutions for the flow behavior of MR fluid flowing through the parallel channel are obtained. The equation for the squeezing force is derived to provide the theoretical foundation for the design of the isolating damper. The result shows that with the increase of the applied magnetic field strength, the squeezing force is increased.


---

## Body

## 1. Introduction

With the continuous construction of nuclear power plants in China, more and more attentions are paid to safety and reliability of the auxiliary equipments in nuclear power plant. There is a pressing need to researching earthquake resistance of auxiliary equipment in the nuclear power plant. The centrifugal fan in the nuclear power plant must have good aseismatic performance that is the important condition for the safe operation of nuclear power plants in the earthquake. However, there have been few reports on vibration insulator for centrifugal fan [1, 2]. Earthquake response attenuation of centrifugal fans can be improved by connecting them with isolator. In this paper, we propose an isolating damper based on magnetorheological (MR) for centrifugal fan.MR fluids are materials of micron-sized and magnetized particles in a carrier fluid. In the absence of an applied magnetic field, MR fluids flow freely. The fluids exhibit Newtonian-like behavior. Upon application of a magnetic field, these fluids exhibit viscoplastic behavior with yield strength [3]. Altering the strength of an applied magnetic field will precisely control the shear yield stress of the fluid. Based on the mechanical characteristics, the fluids can be used in the magnetically controlled devices such as brakes [4] and dampers [5, 6].The MR isolating damper is one such device that provides controllable squeezing force. Altering the strength of an applied magnetic field will change the squeezing force of the MR isolating damper [7]. The research of MR dampers and their applications have been done in many different ways. In the field of the flow model description for MR damper, Boelter and Janocha [8] analyzed the working mode (shear mode, flow mode, squeeze mode) in MR damper. Mcmanus et al. [9] studied the squeeze flow mode in accompany with shear in MR damper, analyzed the damper force-velocity characteristics under different values of magnetic field strength. For the research of application in building structures, Motra et al. [10] researched the response attenuation of seismically excited adjacent buildings connected by a MR damper. P. Y. Lin and T. K. Lin [11] introduced a bridge isolation system that combines the rolling pendulum system (RPS) and the MR damper. Dragasius et al. [12] researched the resistance force generated by linear hydrocylinder-type magnetorheological fluid (MRF) damper acting in different regimes. It was defined that the force increases practically linearly when increasing strength of the magnetic field and when increasing piston speed. Erkus and Johnson [13] investigated the dissipativity and performance characteristics of the semiactive control of the base-isolated benchmark structure with MR dampers. Rajamani and Larparisudthi [14] proposed a vibration control model for MR isolating damper.In this paper, biviscosity model is used to describe the constitutive characteristics of MR fluids subject to an applied magnetic field. The operational principle of the MR isolating damper is introduced. Analytical solutions for the axisymmetric squeeze flow behavior of a biviscosity fluid are obtained. The velocity equation and the location of the unyield flow region are obtained. The expression for the squeezing force is derived to provide the theoretical foundation for the design of the isolating damper. The result shows that with the increase of the applied magnetic field strength, the dynamic yield stress of the MR fluid goes up rapidly, and the squeezing force is increased.


## 2. Operational Principle

The schematic configuration of the proposed circular plate MR isolating damper is shown in Figure1. The MR isolating damper consists of an MR fluid, moving plate, working gap, electromagnet coil, flux guide, and housing. The MR fluid fills the working gap between the moving plate and housing. The electromagnet coil in the housing provides the magnetic field in the working gap. During relative vertical motion between the moving plate and housing, MR fluid is squeezed in the working gap. Thus, the pressure drop due to flow resistance of MR fluid in the working gap is induced. The MR isolating damper produces a controllable squeezing force due to the yield stress of the MR fluid, if a certain level of magnetic field is applied through the working gap.Figure 1
Operational principle of a circular plate MR isolating damper.


## 3. Modeling and Analysis

### 3.1. Analysis of Flow Velocity in Working Gap

Figure2 shows a schematic of the axisymmetric squeeze flow model. The gap between the disks is filled with the MR fluid. The bottom disk is fixed, while the top disk approaches the bottom disk with a constant velocity U. Gravity is neglected, so that the z=h/2 plane is a plane of symmetry for the flow. The force acting on bottom disk is
(1)F=∫0R2πrpdr,
where, p is the pressure, R is radius of the disks.Figure 2
Squeeze model of MR fluid between two parallel disks.To circumvent the “squeeze-flow paradox” of the Bingham model [15] and to facilitate an analytical analysis, the biviscosity constitutive model is employed in this study to describe the behavior of the MR fluid, see in Figure 3. It is known that the biviscosity model can provide very useful and convenient ways for the calculation of materials’ behavior with yield, particularly from an analytical perspective.Figure 3
Biviscosity model of MR fluid.The constitutive equation of the biviscosity model can be represented by the following two expressions [16]:(2a)τ=τ0+ηdurdz,|τ|≥τy(H),(2b)τ=η(H)durdz,|τ|≤τy(H),where τ is the shear stress, τy(H) is the dynamic yield stress developed in response to an applied magnetic field H, τ0 is the intercept stress, η is the viscosity of the fluid when stress is higher than the dynamic yield stress, η(H) is the viscosity coefficient subjected to a stress lower than the dynamic yield stress, it’s a function of the magnetic field strength. Viscosity ratio ε=η/η(H) is an important parameter for the biviscosity constitutive model, τ0=(1-ε)τy(H). Note that when ε→0, (2a) and (2b) are the Bingham constitutive model, whereas when ε=1 they are the Newtonian constitutive model.It is evident from the constitutive equation that the biviscosity constitutive model can be divided into two regions based on the dynamic yield stress. The MR fluid to be yielded when the magnitude of the fluid’s internal stress is greater than the yield stressτy(H), the fluid exhibits typical Newtonian behavior. The MR fluid is unyield when the magnitude of the fluid’s internal stress is smaller than τy(H), the fluid flows slowly with highly viscous.For the disc-type squeeze flow model in Figure2, assuming steady state condition (∂/∂t=0), symmetry (∂/∂θ=0), no tangential and axial velocity (uθ=uz=0) but only radial, no body force, and no pressure gradient in thickness direction, the Navier-Stokes equation in the r direction for the cylindrical coordinates (r,θ,z) is simplified to
(3)ηd2urdz2=-m,
where ur is the radial velocity, m=-dp/dr, dp/dr is the pressure gradient in radial direction.By integrating (3), the velocity profile can be easily obtained as
(4)ur=-m2ηz2+A1z+A2,
where A1 and A2 symbolize the integral constants.The MR fluid exhibits Newtonian behavior in the absence of an applied magnetic field, applying boundary conditions ofur=0 at z=0 and z=h, the flow velocity ur0 can be obtained as follows:
(5)ur0=m2ηz(h-z).According to the assumptions, the momentum equation in ther-direction is
(6)dτzrdz=-m.By integrating (6) along the thickness direction and applying boundary condition of τzr=0 at z=h/2, the following equation for the shear stress can be obtained:
(7)τzr=m(h2-z).Assume thatzy1 and zy2 denote the position of the unyield region boundaries measured from the bottom wall, respectively, see Figure 2. The yield surface is the locus of points where the shear stress is equal to the yield stress (|τ|=τy(H)), the location of the unyield region can be determined by satisfying the conditions such that(8a)zy1=h2-τy(H)m,(8b)zy2=h2+τy(H)m.As seen in Figure2, the flow is composed of three regions: two yield regions (0≤z≤zy1,zy2≤z≤h) contacting the bottom and top walls and an unyield region (zy1≤z≤zy2).In yield region0≤z≤zy1, the two boundary conditions are(9a)ur=0atz=0,(9b)durdz=0atz=zy1.Using the boundary conditions (9a) and (9b), determining the integral constant in (4), the velocity profile in yield region 0≤z≤zy1can be obtained as
(10)ur1=m2η(hz-z2)-τy(H)ηz,0≤z≤zy1.In unyield regionzy1≤z≤zy2, the two boundary conditions are(11a)ur=ur1atz=zy1,(11b)durdz=0atz=h2.The velocity profile in unyield regionzy1≤z≤zy2 can be obtained by using the boundary conditions (11a) and (11b) as
(12)ur2=mε2η(hz-z2)+m1-ε2η(hzy1-zy12)-τy(H)ηzy1,zy1≤z≤zy2.Similar to the previous analysis, by using the boundary conditions ofur=0 at z=h and dur/dz=0 at z=zy2, the velocity profile in yield region zy2≤z≤h can be obtained as
(13)ur3=m2η(hz-z2)-τy(H)η(h-z),zy2≤z≤h.


### 3.2. Pressure Distribution

The pressure gradientdp/dr can be determined by using the law of conservation of mass
(14)DDt∫VρdV=0,
where, ρ is the density of MR fluid, V is the volume of MR fluid. Equation (14) indicates that the derivative of liquid mass is zero
(15)∫V∂ρ∂tdV+∫Sρujnjds=0,
where, s is the integral area, uj and nj denote the velocity tensor and direction cosines.Whenzy1=h/2-τy(H)/m≤0, the MR fluid flows slowly with highly viscous, the pressure gradient can be expressed as
(16)dpdr=6η(H)rh3U.Whenzy1=h/2-τy(H)/m>0, the MR fluid flows with a floating core, submitting (10), (12), (13) into (15) is
(17)πr2U+2πr(∫0zy1ur1dz+∫zy1zy2ur2dz+∫zy2hur3dz)=0,
where
(18)∫0zy1ur1dz=mh4ηzy12-m6ηzy13-τy(H)2ηzy12,∫zy1zy2ur2dz=2∫zy1h/2ur2dz=mhε2η(h24-zy12)-mε3η(h38-zy13)+1-εηmhzy1(h2-zy1)-2mzy1(h2-zy1),-2τy(H)zy1η(h2-zy1),∫zy2hur3dz=∫0zy1ur1dz.So(19)h3m3+[6Urη-(1-ε)3h2τy(H)]m2+4(1-ε)τy3(H)=0.


### 3.3. Squeezing Force

The force acting on the bottom disc can be obtained by integrating the pressure along the radial direction:(20)F=∫dF=∫0R2πrpdr.The subsection integral is used to get the force as(21)F=[πpr2]0R-∫0Rπr2dpdrdr.The boundary conditions of pressure are(22)dpdr=0atr=0,p=0atr=R.
Substituting (21) into (20), the force acting on the bottom disc can be obtained as
(23)F=-π∫0Rr2mdr.


## 3.1. Analysis of Flow Velocity in Working Gap

Figure2 shows a schematic of the axisymmetric squeeze flow model. The gap between the disks is filled with the MR fluid. The bottom disk is fixed, while the top disk approaches the bottom disk with a constant velocity U. Gravity is neglected, so that the z=h/2 plane is a plane of symmetry for the flow. The force acting on bottom disk is
(1)F=∫0R2πrpdr,
where, p is the pressure, R is radius of the disks.Figure 2
Squeeze model of MR fluid between two parallel disks.To circumvent the “squeeze-flow paradox” of the Bingham model [15] and to facilitate an analytical analysis, the biviscosity constitutive model is employed in this study to describe the behavior of the MR fluid, see in Figure 3. It is known that the biviscosity model can provide very useful and convenient ways for the calculation of materials’ behavior with yield, particularly from an analytical perspective.Figure 3
Biviscosity model of MR fluid.The constitutive equation of the biviscosity model can be represented by the following two expressions [16]:(2a)τ=τ0+ηdurdz,|τ|≥τy(H),(2b)τ=η(H)durdz,|τ|≤τy(H),where τ is the shear stress, τy(H) is the dynamic yield stress developed in response to an applied magnetic field H, τ0 is the intercept stress, η is the viscosity of the fluid when stress is higher than the dynamic yield stress, η(H) is the viscosity coefficient subjected to a stress lower than the dynamic yield stress, it’s a function of the magnetic field strength. Viscosity ratio ε=η/η(H) is an important parameter for the biviscosity constitutive model, τ0=(1-ε)τy(H). Note that when ε→0, (2a) and (2b) are the Bingham constitutive model, whereas when ε=1 they are the Newtonian constitutive model.It is evident from the constitutive equation that the biviscosity constitutive model can be divided into two regions based on the dynamic yield stress. The MR fluid to be yielded when the magnitude of the fluid’s internal stress is greater than the yield stressτy(H), the fluid exhibits typical Newtonian behavior. The MR fluid is unyield when the magnitude of the fluid’s internal stress is smaller than τy(H), the fluid flows slowly with highly viscous.For the disc-type squeeze flow model in Figure2, assuming steady state condition (∂/∂t=0), symmetry (∂/∂θ=0), no tangential and axial velocity (uθ=uz=0) but only radial, no body force, and no pressure gradient in thickness direction, the Navier-Stokes equation in the r direction for the cylindrical coordinates (r,θ,z) is simplified to
(3)ηd2urdz2=-m,
where ur is the radial velocity, m=-dp/dr, dp/dr is the pressure gradient in radial direction.By integrating (3), the velocity profile can be easily obtained as
(4)ur=-m2ηz2+A1z+A2,
where A1 and A2 symbolize the integral constants.The MR fluid exhibits Newtonian behavior in the absence of an applied magnetic field, applying boundary conditions ofur=0 at z=0 and z=h, the flow velocity ur0 can be obtained as follows:
(5)ur0=m2ηz(h-z).According to the assumptions, the momentum equation in ther-direction is
(6)dτzrdz=-m.By integrating (6) along the thickness direction and applying boundary condition of τzr=0 at z=h/2, the following equation for the shear stress can be obtained:
(7)τzr=m(h2-z).Assume thatzy1 and zy2 denote the position of the unyield region boundaries measured from the bottom wall, respectively, see Figure 2. The yield surface is the locus of points where the shear stress is equal to the yield stress (|τ|=τy(H)), the location of the unyield region can be determined by satisfying the conditions such that(8a)zy1=h2-τy(H)m,(8b)zy2=h2+τy(H)m.As seen in Figure2, the flow is composed of three regions: two yield regions (0≤z≤zy1,zy2≤z≤h) contacting the bottom and top walls and an unyield region (zy1≤z≤zy2).In yield region0≤z≤zy1, the two boundary conditions are(9a)ur=0atz=0,(9b)durdz=0atz=zy1.Using the boundary conditions (9a) and (9b), determining the integral constant in (4), the velocity profile in yield region 0≤z≤zy1can be obtained as
(10)ur1=m2η(hz-z2)-τy(H)ηz,0≤z≤zy1.In unyield regionzy1≤z≤zy2, the two boundary conditions are(11a)ur=ur1atz=zy1,(11b)durdz=0atz=h2.The velocity profile in unyield regionzy1≤z≤zy2 can be obtained by using the boundary conditions (11a) and (11b) as
(12)ur2=mε2η(hz-z2)+m1-ε2η(hzy1-zy12)-τy(H)ηzy1,zy1≤z≤zy2.Similar to the previous analysis, by using the boundary conditions ofur=0 at z=h and dur/dz=0 at z=zy2, the velocity profile in yield region zy2≤z≤h can be obtained as
(13)ur3=m2η(hz-z2)-τy(H)η(h-z),zy2≤z≤h.


## 3.2. Pressure Distribution

The pressure gradientdp/dr can be determined by using the law of conservation of mass
(14)DDt∫VρdV=0,
where, ρ is the density of MR fluid, V is the volume of MR fluid. Equation (14) indicates that the derivative of liquid mass is zero
(15)∫V∂ρ∂tdV+∫Sρujnjds=0,
where, s is the integral area, uj and nj denote the velocity tensor and direction cosines.Whenzy1=h/2-τy(H)/m≤0, the MR fluid flows slowly with highly viscous, the pressure gradient can be expressed as
(16)dpdr=6η(H)rh3U.Whenzy1=h/2-τy(H)/m>0, the MR fluid flows with a floating core, submitting (10), (12), (13) into (15) is
(17)πr2U+2πr(∫0zy1ur1dz+∫zy1zy2ur2dz+∫zy2hur3dz)=0,
where
(18)∫0zy1ur1dz=mh4ηzy12-m6ηzy13-τy(H)2ηzy12,∫zy1zy2ur2dz=2∫zy1h/2ur2dz=mhε2η(h24-zy12)-mε3η(h38-zy13)+1-εηmhzy1(h2-zy1)-2mzy1(h2-zy1),-2τy(H)zy1η(h2-zy1),∫zy2hur3dz=∫0zy1ur1dz.So(19)h3m3+[6Urη-(1-ε)3h2τy(H)]m2+4(1-ε)τy3(H)=0.


## 3.3. Squeezing Force

The force acting on the bottom disc can be obtained by integrating the pressure along the radial direction:(20)F=∫dF=∫0R2πrpdr.The subsection integral is used to get the force as(21)F=[πpr2]0R-∫0Rπr2dpdrdr.The boundary conditions of pressure are(22)dpdr=0atr=0,p=0atr=R.
Substituting (21) into (20), the force acting on the bottom disc can be obtained as
(23)F=-π∫0Rr2mdr.


## 4. Results and Discussion

A typical MR fluid is used in this paper. Figure4 shows the yield strength of MR fluid under different magnetic field strength, measured by experiment. For the purpose of illustration, the following parameters are given: h=1 mm, R=100 mm, U=5 m/s, η=0.092Pa·s, and ε=0.001.Figure 4
Yield strength versus magnetic field strength.The relationship between the pressure gradient and the radius with apparent yield stress value of 30 kPa is shown in Figure5 under the application of magnetic field strength of 100 kAmp/m by solving (16) and (19). The absolute value of the pressure gradient increases with the radius in linear, approximately.Figure 5
The pressure gradient versus radius.When the strength of magnetic field is 100 kAmp/m, the yield surface determined by (8a) and (8b) with apparent yield stress value of 30 kPa along the radial direction is shown in Figure 6. When r≤17.7mm, the MR fluid is unyield totally, flowing slowly with highly viscous. When r>17.7mm, the MR fluid is yielded partially. The location of yield surface tends to move toward the symmetry plane of working gap as the radius increases.Figure 6
The yield surface along the radial direction.The velocity profiles can be obtained by (10), (12), and (13) with apparent yield stress value of 30 kPa, shown in Figure 7. The velocity profiles for the squeeze flow of MR fluid appear to be parabolic curves with respect to z, and their curvature tends to decrease with radius, the unyield region decreases accordingly.Figure 7
Velocity profiles for different location of radius.The squeezing force versus magnetic field strength is shown in Figure8. In the absence of an applied magnetic field, the force is 216 kN. The forces are 250 kN, 278 kN, and 302 kN at the strength of magnetic field of 50 kAmp/m, 100 kAmp/m, and 200 kAmp/m, respectively. The results indicate that with the increase of the velocity and the applied magnetic field, the squeezing force is increased.Figure 8
Squeezing force versus magnetic field strength.


## 5. Conclusions

The flow behaviors of MR fluid in circular plate MR isolating damper are investigated theoretically in this paper. The equations for the velocity and the squeezing force are derived to provide the theoretical foundation for the design of the MR isolating damper. The unyield region of MR fluid tends to move toward the symmetry plane of working gap as the radius increases. With the increase of the applied magnetic field strength, the squeezing force increased.



---

*Source: 175703-2012-10-09.xml*