Next Article in Journal
Damage Equivalent Test Methodologies as Design Elements for Journal Bearing Systems
Previous Article in Journal
Reduction of the Adhesive Friction of Elastomers through Laser Texturing of Injection Molds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Evaluation of Linear Complementarity Problem Method in Modeling the Fluid Cavitation for Squeeze Film Damper with Off-Centered Whirling Motion

Department of Mechanical and Industrial Engineering, University of Toronto, Toronto, ON M5S3G8, Canada
*
Author to whom correspondence should be addressed.
Lubricants 2017, 5(4), 46; https://doi.org/10.3390/lubricants5040046
Submission received: 4 November 2017 / Revised: 23 November 2017 / Accepted: 24 November 2017 / Published: 28 November 2017

Abstract

:
For the application of squeeze film damper (SFD) in aero-engine, a cavitation model is evaluated by means of linear complementarity problem (LCP) method. Different from the conventional SFD study that employs circular-center orbits (CCOs), a realistic condition is explored where the shaft whirling center and bearing center are misaligned. Taking into account the fluid as incompressible and compressible, the governing equations, including film cavitation, are respectively solved by developing an algorithm using the LCP method. The numerical results are compared with experimental data and the effectiveness of the model is verified. The proposed model can provide some references to investigate the competency of this cavitation method in SFDs.
Keywords:
SFD; LCP; cavitation

Graphical Abstract

1. Introduction

Squeeze film dampers (SFDs) have been widely used in modern aircraft engines to attenuate the structure vibration and reduce the cabin noise. SFDs are lubricating elements that provide mechanical support and viscous damping in a rotordynamic system. The application of SFDs on aircraft engines reduces the vibration amplitudes of the rotor at critical speeds and increases the onset of the rotor instability speed [1]. However, current analytical SFD models are often simplified and have not been calibrated for engine applications. Very little test data are available in the literature for realistic engine conditions. Therefore, it is essential to advance the modeling technology in this field and provide references for its industrial applications.
The competence of an SFD involves several nonlinear components, such as lubricant cavitation, fluid inertia, damper geometry, seal presence, and oil feeding and draining configurations. These elements play key roles in affecting the behavior of an SFD and they have undergone investigation for many years [1,2]. Despite a large quantity of work in this field, cavitation is still not fully understood to accurately evaluate its impact on the performance of SFD. Cavitation is a phenomenon of film rupture when the fluid pressure drops to the vapor pressure or even below, which is characterized by a complex process of bubble nucleation, growth, and implosion. The occurrence of film cavitation leads to a two-phase flow problem in SFD, which drastically reduces the damper capability [3].
Direct observation of cavitation in SFD is achieved through capturing high-speed photographic images and recording stroboscopic videos [3,4]. The measurement of absolute pressure in a cavitated region is attained by applying specific strain-gauge sensors [5,6,7] or pressure transducers [8,9,10]. These experiments also offer the guidelines for understanding the consequence of cavitation in a squeezed film.
The theoretical approaches to address the cavitation problem diverge on the detection of the cavitation boundary. Conventional models apply the Jakobsson-Floberg-Olsson (JFO) theory [11,12] to numerically solve the cavitation onset and the reformation boundary through extensive iterations. The Elrod cavitation algorithm [13] accelerates the computation by introducing a technique that transforms the governing equation from the elliptic form to a parabolic form. The concept of mass conserving in the cavitated film is recently interpreted as a linear complementarity problem (LCP) and a reformulated model is illustrated for an incompressible liquid [14]. For the compressible flow, two mathematical models are presented in references [15,16] for a thin film flow between two surfaces in relative motion. Although the LCP method demonstrates general agreement with analytical results in several specific bearing configurations [14,15,16], little effort has been devoted on this manner to addressing the cavitation phenomenon in SFD and comparing it with the experiment data.
A new LCP algorithm is proposed in this paper for modeling the cavitation in SFD. The developed cavitation model includes a complementarity pair that is different from previous work [14,15,16]. The paper is organized as follows: firstly, the governing equation is presented for the flow in SFD with open-ended condition; accordingly, two mathematical models are introduced to investigate the cavitation, where the flow is treated as incompressible and compressible, respectively. The formulated equations contain complementarity variables and the equations can be numerically solved. Moreover, the cavitation algorithm is applied to simulate off-centered whirling motion with no oil feeding, fully degassed lubricant, and open-end SFD application. Finally, the predictions are compared with existing experimental data [17] to verify the model accuracy.

2. Mathamatical Model

The flow equation in a hydrodynamic bearing system is usually governed by the Reynolds equation [18]. When considering an SFD that has a small length-to-diameter ratio (i.e., L / D 0 ) with open-ended configuration, the circumferential flow in the film can be regarded as effectively small [2]. Accordingly, the Reynolds equation is reduced to
z ( ρ h 3 12 μ p z ) = ρ h t
For SFD executing offset circular orbits, the film thickness at a specific time is determined by parameters including the radial clearance, the orbit ratio, the offset ratio, and the whirling frequency. Based on the coordinate relationship shown in Figure 1, the mathematical expression of film thickness is given by [19]
h = c ( 1 δ cos α λ cos ( β α ) )
β = ω t
The solution to Equation (1) provides the general dynamic response of SFD. To solve this equation with the occurrence of film cavitation, the following method is proposed.

2.1. Incompressible Flow

The LCP method for incompressible flow was introduced in [14], where the author considered cavitation only at absolute zero pressure. In this paper, the LCP method is extended to a model where the film can be ruptured at any prescribed pressure. This new incompressible flow model is more general in cavitation conditions and it follows the principle of mass conservation.

2.1.1. Governing Equation

The lubricant density is assumed to be a constant in the full-film region, namely ρ 0 , and the pressure is considered unchanged in the cavitation zone, namely p c . Instead of solving the pressure and density directly, one can solve for the relative pressure Δ p and the relative density Δ ρ first, where
Δ p = p p c Δ ρ = ρ 0 ρ
In the full-film region, the pressure is above the cavitation pressure (i.e., p > p c ) while the density is equal to the lubricant density (i.e., ρ = ρ 0 ). One can get the following relation accordingly
{ Δ p > 0 Δ ρ = 0 Δ p Δ ρ = 0 Δ p z Δ ρ = 0
In the cavitation zone, the pressure is equal to the cavitation pressure (i.e., p = p c ), while the density is below the lubricant density (i.e., ρ < ρ 0 ), rendering that
{ Δ p = 0 Δ ρ > 0 Δ p Δ ρ = 0 Δ p z Δ ρ = 0
If the Swift-Stieber boundary condition [18] is applied for the cavitation interface and we assume that the mass transfer is null between the full-film and the cavity regions, the following equation is derived
{ Δ p = 0 Δ ρ = 0 Δ p Δ ρ = 0 Δ p z Δ ρ = 0
Consequently, the relative pressure and relative density are complementary to each other everywhere within the lubricant film. In other words, the following condition is valid regardless of film rupture
{ Δ p 0 Δ ρ 0 Δ p Δ ρ = 0 Δ p z Δ ρ = 0
Since Δ p and Δ ρ are complementarity pairs, the governing equation is proposed to be interpreted into a linear complementarity equation.
Substitution of Equation (4) into Equation (1) gives
z [ ( ρ 0 Δ ρ ) h 3 12 μ ( Δ p + p c ) z ] = ( ρ 0 Δ ρ ) h t
As the cavitation pressure is assumed to be constant, the above equation can be reduced to
z [ ( ρ 0 Δ ρ ) h 3 12 μ Δ p z ] = ( ρ 0 Δ ρ ) h t
Application of the complementarity condition (i.e., Equation (8)) to Equation (10) yields
z ( ρ 0 h 3 12 μ Δ p z ) = ( ρ 0 Δ ρ ) h t
Provided that the pressure and density at the initial state, i.e., t = 0 , one can incrementally solve the above equation at any time and find the dynamic response of the system. The backward difference scheme is used to approximate the time derivations in this study.
To formulate an LCP equation, the finite difference method is applied to discretize the physical domain (which is the axial length in this case) into a uniform grid with N elements of size Δ z ( Δ z = L / N ). A central difference scheme is used to approximate the derivatives on the left side of Equation (11), and a backward difference scheme is applied on the right side. The discretized equation is given as follows
ρ 0 ( h k ) 3 12 μ Δ p i + 1 k 2 Δ p i k + Δ p i 1 k ( Δ z ) 2 = ( ρ 0 Δ ρ i k ) h k ( ρ 0 Δ ρ i k 1 ) h k 1 Δ t
where i denotes the node in the axial location, while k denotes the step in the time frame.
Equation (12) can be further organized into a standard LCP form as
{ { Δ p } = [ M ] { Δ ρ } + { q } { Δ p } 0 { Δ ρ } 0 { Δ p } { Δ ρ } T = 0
where the detailed description is provided in Appendix A.

2.1.2. Boundary Condition

The boundary condition is also constructed with complementarity variables. Given the pressure condition at specific locations where the film is unruptured, e.g., p = p i n > p c at inlet, we have
{ Δ p i n = p i n p c Δ ρ i n = 0
Furthermore, the above condition can be represented in the form of a linear equation as
Δ p i n = Δ ρ i n + p i n p c
A similar equation for the flow at the outlet boundary can also be derived. The complete LCP equation for the flow response in an SFD is formulated by these boundary equations and the governing equation for the interior points (i.e., Equation (13)).
An LCP solver would assist in providing the solution with respect to time given the initial flow state. The solution of the system in the current time step is determined by the flow status at the previous time step, and this solution will be subsequently implemented to find the fluid response for the next time step.

2.1.3. Initial Condition

The fluid at initial status can be solved through Equation (11) by neglecting the time variation for the density variables. Accordingly, a discretized form of such equation is shown as follows, and its solution is used as the input for the first time step.
ρ 0 ( h 0 ) 3 12 μ Δ p i + 1 0 2 Δ p i 0 + Δ p i 1 0 ( Δ z ) 2 = ( ρ 0 Δ ρ i 0 ) h 0 t
Similarly, one can apply the introduced LCP method to solve the above equation.

2.2. Compressible Flow

The fluid bulk modulus is a material property that is characterizing the compressibility of a fluid. It is employed by reference [13,15] when developing the cavitation algorithms. Our work continues to use this parameter to address the compressibility of the fluid in the liquid phase, and integrate it with the LCP method to find the boundaries for both the leading and trailing edges of the cavitation.

2.2.1. Governing Equation

The fluid bulk modulus B is given in Equation (17), and it shows the relationship between the density and pressure
B = ρ p ρ
Assuming that the fluid bulk modulus is a constant for the flow in the liquid phase, one obtains the following after a direct integration of Equation (17)
ρ = ρ c e ( p p c ) / B
Equation (18) can be used to describe the fluid density in the full-film domain; however, it is not compatible for the density in the cavitated region. To address this issue, a modification is proposed by adding a variable to the above equation
ρ = ρ c e ( p p c ) / B ξ
where ξ is defined as a saturation function and it satisfies the following conditions
ξ = 0 i n f u l l f i l m r e g i o n 0 < ξ < ρ c i n c a v i t a t i o n r e g i o n
However, the unknown variables p and ξ are not in a linear combination. In order to form an LCP equation, Equation (19) is further modified as
ρ = ρ c + η ξ
where η is a new variable defined as
η = ρ c [ e ( p p c ) / B 1 ]
Equation (21) provides a linear combination of the complementarity pairs η and ξ , which is in the ideal form to be implemented into the governing equation.
Since the variable η connects the unknown pressure p in the film, which is of fundamental interest in an SFD, p z in Equation (1) needs to be transformed in terms of η . Equation (22) yields
p = p c + B ln ( 1 + η ρ c )
Accordingly, the following relation is derived
p z = p η η z = B ρ c + η η z
In the full-film region, where p p c > 0 , the introduced variables follow the condition as
{ ξ = 0 η > 0 ξ η = 0 ξ η z = 0
In the cavitation region, where p p c = 0 , the following equation is derived
{ ξ > 0 η = 0 ξ η = 0 ξ η z = 0
At the cavitation interface, if the Swift-Stieber condition is considered and the nil mass transfer is assumed, one renders the following
{ ξ = 0 η = 0 ξ η = 0 ξ η z = 0
Therefore, the introduced variables η and ξ are non-negative and complementarity to each other within the squeezed film, i.e.,
{ ξ 0 η 0 ξ η = 0 ξ η z = 0
Substitution of Equations (21) and (24) into Equation (1) gives
z [ B h 3 ( ρ c + η ξ ) 12 μ ( ρ c + η ) η z ] = ( ρ c + η ξ ) h t
Applying Equation (28) to the above equation, a reduced form is obtained as
z [ B h 3 12 μ η z ] = ( η ξ ) h t + ρ c h t
Subsequently, the finite difference method is employed to numerically solve the above equation. After discretizing, the following equation is achieved
B ( h k ) 3 12 μ η i + 1 k 2 η i k + η i 1 k ( Δ z ) 2 = ( η i k ξ i k ) h k ( η i k 1 ξ i k 1 ) h k 1 Δ t + ρ c h k h k 1 Δ t
Consequently, an LCP equation is formulated with respect to the complementarity pairs. The organized LCP equation is
{ { η } = [ M ] { ξ } + { q } { η } 0 { ξ } 0 { η } { ξ } T = 0
where the detailed description is provided in Appendix B.

2.2.2. Boundary Condition

If the full-film flow is considered for the flow at the boundary, one can obtain the following equation for a compressible liquid flow
{ η i n = ρ c [ e ( p i n p c ) / B 1 ] ξ i n = 0
Hence, Equation (33) is organized into an LCP equation form as
η i n = ξ i n + ρ c [ e ( p i n p c ) / B 1 ]
The integration of the LCP equation for the boundary and the interior points formulates the complete equation set that needs to be solved using the existing solver.

2.2.3. Initial Condition

Similarly, the equation for the initial condition can be obtained through neglecting the time variation of unknown variables in Equation (30), i.e.,
B ( h 0 ) 3 12 μ η i + 1 0 2 η i 0 + η i 1 0 ( Δ z ) 2 = ( η i 0 ξ i 0 ) h 0 t + ρ c h 0 t
The solution of Equation (43) is regarded as the flow at the initial state, and it is used as the input value to solve the flow status for a dynamic response.
The simulation result and comparison are discussed in the next section.

3. Case Study

In this section, simulation results under different running conditions are presented. The simulation scenario is simplified by disregarding the thermo-hydrodynamics and air-ingestion effects. Some results from the numerical model are validated against the experiment [17]. The experiment is conducted on an SFD executing offset circular whirls, in the absence of fluid inertia and end-seals. Pressure signals at three locations that are evenly spaced around the squeeze film at the mid-plane are simultaneously recorded in the test.
Table 1 summarizes the main characteristics of three cases corresponding to the experiment at ‘II run’, ‘IV run’, and ‘V run’ [17]. These cases are selected in the presence of fluid cavitation without the effect of entrained air. The squeeze film Reynolds number at these conditions is less than 0.68 so that the effect of fluid inertia is neglected [2]. The inlet and outlet pressure are considered as approximately equal to the mean pressure of the inlet and outlet grooves in the test rig, respectively. Furthermore, in the simulation model, a linear variation for the fluid viscosity is assumed in Case I. This viscosity function is given as Equation (40), which is developed based on the test condition in order to provide an accurate lubricant property for scenarios under varying inlet pressure
μ = p i n + 20 300
The cavitation pressure is prescribed to be 0.01 bar in all three cases, and the fluid bulk modulus is tested on a value of 10MPa. Simulation results at three different locations, i.e., α1, α2, and α3, are shown in a motion of one whirling cycle.
Figure 2 and Figure 3 show the result comparison of Case I between the developed models and the experiment. Figure 2a–c illustrate the pressure for the measured position at α1, α2, and α3, respectively, under the gauge pressure of 0.2 bars. In the positive pressure ranges, the compressible model provides excellent results that are closer to the experiment, while the incompressible model over-predicts the peak pressure. In the cavitation zones, the developed models slightly underestimate the extent of cavitation, with both the incompressible model and compressible model having the same edges for the film rupture and formation boundaries. A deviation of the simulation results from the experiment is that the developed models are not capable of providing any pressure data under the cavitation pressure. This limitation is the result of the assumed cavitation boundary, where a prescribed value is applied to define the pressure in the cavitation zone and the cavitation interface is modeled using the Swift-Stieber condition. In other words, the model implies a condition that the fluid is unable to reach any subcavity pressure without considering the tensile strength of the fluid; while in the actual operation where the SFD is whirling at high speeds, the fluid pressure rapidly oscillates and it reaches negative as a result of the tensile strength. Figure 3a–c showed the comparison of pressure at three locations under the gauge pressure of 2.0 bars. It is clear that with the increase of the gauge pressure, cavitation is reduced significantly. The compressible model appears to offer better results than the incompressible model at all of the different positions.
Figure 4 and Figure 5 plot the results for Case II and Case III, respectively. It detects quite clearly some distinctive features in the behavior of the fluid film at high offset ratio. As shown in Figure 4a and Figure 5a, the films remain unruptured at measured position 1 with the smallest pressure fluctuations under both operating conditions. At position 3, where the pressure is theoretically the highest, the simulation provides an accurate prediction of the pressure profile and the extent of the cavitation region, shown in Figure 4c and Figure 5c. A deviance is recorded on position 2, where the maximum pressure is higher than predicted, especially in Case III.
Figure 6 reports a further investigation of the model sensibility by exploring the effect of fluid compressibility. To be specific, a series of the fluid bulk modulus is implemented in the compressible model to produce comparable results. The simulation is conducted under the operating condition of Case I. Figure 6a,b demonstrate the pressure profile at one cycle for position 2 and position 3, respectively. With small fluid bulk modulus, i.e., B = 105, the fluid film is predicted to be in a status with the absence of cavitation. However, as this number increases, the peak film pressure increases; meanwhile, the cavitation zone is expanded as an accompanied phenomenon. The prediction of peak pressure has a limit level for the compressible model, regardless of the magnitude of the fluid bulk modulus. If compared with the incompressible flow model, the simulation from the compressible model presents a lower maximum pressure but a larger cavitation area, and a higher bulk modulus provides closer predictions to the ones from the incompressible model.
Figure 7 examines the model competence by adjusting the operating speed and oil supply pressure, following the scenario from the Case I study. The whirling speed is carried out at low, medium, and high speeds, i.e., 1000 RPM, 2000 RPM, and 3000 RPM, respectively. Figure 7a shows the film behavior of position α1 under difference running speeds. The pressure is fairly regular, with small amplitude at low speed, in the absence of film cavitation. As the speed increases to a higher level, the maximum pressure raises and cavitation starts to occur. A large operating speed tends to cause an increased amount of cavitation and larger pressure variations. Figure 7b plots the mid-plane pressure of the film under different oil feeding conditions, with the supply pressure studied at 0, 1, and 2 bars. At small gauge pressure, the developed film pressure is relatively low, while the region of cavitation is fairly large. As the inlet pressure increases, the film pressure increases simultaneously, and the size of cavitation shrinks. It is very likely that if the inlet pressure reaches a certain level, cavitation would disappear.

4. Conclusions

This research demonstrates a new approach to address the cavitation phenomenon in SFD. The model is developed based on the LCP method for both incompressible and compressible flow assumptions. To determine the pressure distribution in the fluid system, the finite difference method is applied. The resulting equation is efficiently solved by the LCP solver. The cavitation model is validated against the experimental data that is reported in the literature under three different operating conditions. Both incompressible and compressible flow models present overall agreement with the test data. The results show that the compressible fluid model has the advantage of predicting the pressure in the positive range. Both models give excellent performance on the prediction of the cavitation extent.
The model sensitivity is studied by evaluating the value of fluid bulk modulus. Higher fluid bulk modulus provides closer predictions to that from the incompressible flow model. The computational speeds for both models are similarly fast without facing any convergence issues. Given the lubricant bulk modulus, the compressible model is suggested for application due to its advantage on the pressure prediction in the active film section.
Furthermore, simulation is performed at different operating speeds and different conditions for oil supply pressure. The predicted results shows that high speeds and low feeding pressure are contributors in the development of film cavitation, and high supply pressure would prevent the occurrence of cavitation.
Overall, the developed model is reliable to predict vibration response in SFD applications.

Acknowledgments

This research is supported by grants from Natural Science and Engineering Research Council (NSERC) and Pratt and Whitney Canada.

Author Contributions

Tieshu Fan developed the mathematical models, performed the simulation for validation and wrote the manuscript. Kamran Behdinan was the research supervisor.

Conflicts of Interest

The authors declare there is no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
SymbolQuantity
BFluid bulk modulus
cRadial clearance of the SFD
DDiameter of the SFD
dJournal orbit center offset
eJournal eccentricity
FrRadial force
FtTangential force
hFilm thickness
LSFD length
pFluid pressure
pcCavitation pressure
pinInlet film pressure
poutOutlet film pressure
rOrbit radius
RBRadius of the bearing housing
ReSqueeze film Reynolds number
RJRadius of the journal shaft
tTime
zAxial direction
αAngular coordinate of the measuring section
β=ωt, Angular coordinate of the journal shaft
δ=d/c, offset ratio
ε=e/c, eccentricity ratio
θAngular coordinate with reference to the journal shaft
θ0Attitude angle of the journal
λ=r/c, orbit ratio
μFluid dynamic viscosity
ρFluid density
ρcLiquid density at cavitation pressure for compressible flow
ρ0Liquid density for incompressible flow
ωWhirling velocity

Appendix A

The details of Equation (13) are described as follows:
{ Δ p } = { Δ p 1 k Δ p N k } T
{ Δ ρ } = { Δ ρ 1 k Δ ρ N k } T
[ M ] = h k Δ t [ ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 ρ 0 ( h k ) 3 6 μ ( Δ z ) 2 ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 0 0 ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 ρ 0 ( h k ) 3 6 μ ( Δ z ) 2 ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 ] 1
{ q } = ρ 0 h k ( ρ 0 Δ ρ i k 1 ) h k 1 Δ t [ ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 ρ 0 ( h k ) 3 6 μ ( Δ z ) 2 ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 0 0 ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 ρ 0 ( h k ) 3 6 μ ( Δ z ) 2 ρ 0 ( h k ) 3 12 μ ( Δ z ) 2 ] 1

Appendix B

The details of Equation (32) are described as follows:
{ η } = { η 1 k η N k } T
{ ξ } = { ξ 1 k ξ N k } T
[ M ] = h k Δ t [ B ( h k ) 3 12 μ ( Δ z ) 2 B ( h k ) 3 6 μ ( Δ z ) 2 h k Δ t B ( h k ) 3 12 μ ( Δ z ) 2 0 0 B ( h k ) 3 12 μ ( Δ z ) 2 B ( h k ) 3 6 μ ( Δ z ) 2 h k Δ t B ( h k ) 3 12 μ ( Δ z ) 2 ] 1
{ q } = ( ξ k 1 η k 1 ) h k 1 + ρ c ( h k h k 1 ) Δ t [ B ( h k ) 3 12 μ ( Δ z ) 2 B ( h k ) 3 6 μ ( Δ z ) 2 h k Δ t B ( h k ) 3 12 μ ( Δ z ) 2 0 0 B ( h k ) 3 12 μ ( Δ z ) 2 B ( h k ) 3 6 μ ( Δ z ) 2 h k Δ t B ( h k ) 3 12 μ ( Δ z ) 2 ] 1

References

  1. Adiletta, G.; Della Pietra, L. The squeeze film damper over four decades of investigations. Part II: Rotordynamic analyses with rigid and flexible rotors. Shock Vib. Dig. 2002, 34, 97–126. [Google Scholar]
  2. Della Pietra, L.; Adiletta, G. The squeeze film damper over four decades of investigations. Part I: Characteristics and operating features. Shock Vib. Dig. 2002, 34, 3–26. [Google Scholar]
  3. Zeidan, F.; Vance, J. Cavitation leading to a Two Phase Fluid in a Squeeze Film Damper. STLE Tribol. Trans. 1989, 32, 100–104. [Google Scholar] [CrossRef]
  4. Walton, J.F., II; Walovit, J.A.; Zorzi, E.S.; Schrand, J. Experimental Observation of Cavitating Squeeze-Film Dampers. ASME J. Tribol 1987, 109, 290–295. [Google Scholar] [CrossRef]
  5. Arauz, G.; San Andrés, L.A. Experimental Study on the Effect of a Circumferential Feeding Groove on the Dynamic Force Response of a Sealed Squeeze Film Damper. ASME J. Tribol. 1996, 118, 900–905. [Google Scholar] [CrossRef]
  6. Jung, S.Y.; San Andrés, L.A.; Vance, J.M. Measurements of Pressure Distributions and Force Coefficients in a Squeeze Film Damper Part I: Fully Open Ended Configuration. STLE Tribol. Trans. 1991, 34, 375–382. [Google Scholar] [CrossRef]
  7. Jung, S.Y.; San Andrés, L.A.; Vance, J.M. Measurements of Pressure Distributions and Force Coefficients in a Squeeze Film Damper Part II: Partially Sealed Configuration. STLE Tribol. Trans. 1991, 34, 383–388. [Google Scholar] [CrossRef]
  8. Tichy, J.A. Measurement of Squeeze Film Bearing Forces and Pressures, Including the Effect of Fluid Inertia. ASLE Trans. 1985, 28, 520–526. [Google Scholar] [CrossRef]
  9. Vance, J.M.; Kirton, A.J. Experimental Measurement of the Dynamic Force Response of a Squeeze-Film Bearing Damper. ASME J. Eng. Ind. 1975, 97, 1282–1290. [Google Scholar] [CrossRef]
  10. Ku, C.-P.; Tichy, J.A. An Experimental and Theoretical Study of Cavitation in a Finite Submerged Squeeze Film Damper. ASME J. Tribol. 1990, 112, 725–733. [Google Scholar] [CrossRef]
  11. Folberg, L.; Jakobsson, B. The Finite Journal Bearing, Considering Vaporization; Gumperts: Mississauga, ON, Canada, 1957. [Google Scholar]
  12. Otsson, K. Cavitation in Dynamically Loaded Bearings; Scandinavian University Books: Copenhagen, Denmark, 1965. [Google Scholar]
  13. Elrod, H.G. A Cavitation Algorithm. ASME J. Tribol. 1981, 103, 350–354. [Google Scholar] [CrossRef]
  14. Giacopini, M.; Fowell, M.T.; Dini, D.; Strozzi, A. A Mass Conserving Complementarity Formulation to Study Lubricant Films in the Presence of Cavitation. ASME J. Tribol. 2010, 132, 041702. [Google Scholar] [CrossRef]
  15. Almqvist, A.; Fabricius, J.; Larsson, R.; Wall, P. A New Approach for Studying Cavitation in Lubrication. ASME J. Tribol. 2013, 136, 011706-1. [Google Scholar] [CrossRef]
  16. Bertocchi, L.; Dini, D.; Giacopini, M.; Fowell, M.T.; Baldini, A. Fluid film lubrication in the presence of cavitation: A mass-conserving two-dimensional formulation for compressible, piezoviscous and non-Newtonian fluids. Tribol. Int. 2013, 67, 61–71. [Google Scholar] [CrossRef]
  17. Adiletta, G.; Della Pietra, L. Experimental Study of a Squeeze Film Damper with Eccentric Circular Orbits. ASME J. Tribol. 2005, 128, 365–377. [Google Scholar] [CrossRef]
  18. Szeri, A.Z. Fluid Film Lubricant: Theory and Design; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  19. Della Pietra, L. Analytical and Experimental Investigation of Squeeze-Film Dampers Executing Circular Orbits. Meccanina 2000, 35, 133–157. [Google Scholar] [CrossRef]
Figure 1. Coordinate relationship in SFD.
Figure 1. Coordinate relationship in SFD.
Lubricants 05 00046 g001
Figure 2. Case I with low feeding pressure.
Figure 2. Case I with low feeding pressure.
Lubricants 05 00046 g002
Figure 3. Case I with high feeding pressure.
Figure 3. Case I with high feeding pressure.
Lubricants 05 00046 g003
Figure 4. Model comparison for Case II.
Figure 4. Model comparison for Case II.
Lubricants 05 00046 g004
Figure 5. Model comparison for Case III.
Figure 5. Model comparison for Case III.
Lubricants 05 00046 g005
Figure 6. influence of fluid bulk modulus on the pressure profile.
Figure 6. influence of fluid bulk modulus on the pressure profile.
Lubricants 05 00046 g006
Figure 7. Influence of operation speed and feeding pressure on the pressure profile.
Figure 7. Influence of operation speed and feeding pressure on the pressure profile.
Lubricants 05 00046 g007
Table 1. Parameters tested in the experiment.
Table 1. Parameters tested in the experiment.
ParameterCase ICase IICase III
c [mm]0.3850.3850.385
D [mm]140140140
L [mm]303030
pin [barg]2→0.20.51
pout [barg]000
δ0.0420.230.23
λ0.50.460.46
ρ0 [kg/m3]877877877
μ [pa·s]0.06→0.0660.0980.077
ω [rad/s]314178209
α1 [rad]1.443.363.36
α2 [rad]5.631.271.27
α3 [rad]3.535.455.45

Share and Cite

MDPI and ACS Style

Fan, T.; Behdinan, K. The Evaluation of Linear Complementarity Problem Method in Modeling the Fluid Cavitation for Squeeze Film Damper with Off-Centered Whirling Motion. Lubricants 2017, 5, 46. https://doi.org/10.3390/lubricants5040046

AMA Style

Fan T, Behdinan K. The Evaluation of Linear Complementarity Problem Method in Modeling the Fluid Cavitation for Squeeze Film Damper with Off-Centered Whirling Motion. Lubricants. 2017; 5(4):46. https://doi.org/10.3390/lubricants5040046

Chicago/Turabian Style

Fan, Tieshu, and Kamran Behdinan. 2017. "The Evaluation of Linear Complementarity Problem Method in Modeling the Fluid Cavitation for Squeeze Film Damper with Off-Centered Whirling Motion" Lubricants 5, no. 4: 46. https://doi.org/10.3390/lubricants5040046

APA Style

Fan, T., & Behdinan, K. (2017). The Evaluation of Linear Complementarity Problem Method in Modeling the Fluid Cavitation for Squeeze Film Damper with Off-Centered Whirling Motion. Lubricants, 5(4), 46. https://doi.org/10.3390/lubricants5040046

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop