Recently, herringbone grooves journal bearing is used extensively in small high-speed rotating mechanisms, due to their outstanding higher stability and lower leakage, in comparison with a plain journal bearing.
Until now, numerical studies are carried out to investigate characteristics of herringbone groove journal bearing by many researchers. Vohr and Pan  proposed narrow groove theory for analysis of a herringbone grooved gas journal bearing in 1963, subsequently applied by Vohr and Chow  to get numerical solutions at small eccentricity. In 1971, Hamrock and Fleming  determined the optimal parameters for self-acting herringbone grooved journal bearings at the maximum radial load capacity.
In 1980, the film in an incompressible fluid was analyzed numerically by Murata et al.  based on a potential flow theory. In 1984, Kawabata, Ashino, and Tachibana  treated a case of large eccentricity, using a narrow groove theory. In 1994, Bonneau and Absi  applied a finite element method to a compressible Reynolds equation to get aerodynamic characteristics for 4 through 16 grooves with moderate eccentricity.
In 2010, the authors employed a spectral difference scheme to analyze an oil-lubricated herringbone grooved journal bearing with trapezoidal groove . The fluid film thickness with trapezoidal groove geometry is expressed in a Fourier series to the circumferential direction, and Reynolds equation is decomposed into each component of the Fourier series, so that it can be integrated with respect to time independently of each component to get a steady-state solution  . This process results in high resolution in space and high speed in computation. However, an inexpedience associated with Fourier series expansion of fluid film thickness has be faced at the case of rectangle groove—the persistent discrepancy, an overshoot (Gibbs phenomenon) between fluid film thickness and its approximation by a Fourier series at the number of terms in the series becomes indefinitely large. In the paper, an inexpedience of numerical analysis under a special case at rectangle groove is challenged to solve, and the spectral difference scheme will be extended as to be valid under a special case at rectangle groove.
2. Analytical Model
Consider a fluid-lubricated journal bearing equipped with herringbone grooves as shown in Figure 1. Let bearing length be 2l and groove be symmetric with respect to its center of bearing. The shaft itself rotates around its center Os with an angular velocity ω in the counter-clockwise direction, and revolves around the center Ob of the fixed-bearing with an angular velocity Ω in the counter-clockwise direction, the inner radius of the bearing is Rb0, the bearing clearance Cr is defined as , and the radius of the shaft corresponding to the plane without grooves is Rs0, the groove depth, the groove width, ridge width, and grooves angle are denoted by Cg, ag, ar, and β respectively.
The eccentricity of the shaft is given by , and the fluid film thickness, h, is defined as .
Here two coordinate systems and which are fixed at the rotation shaft, are introduced as shown in Figure 2, and relationship between the coordinates is given by
Figure 1. Herringbone-grooved journal bearing.
Figure 2. Coordinate systems and .
3. A Spectral Finite Difference Scheme for Reynolds Fluid Lubrication Equation
In coordinate system , assuming a steady state, and , the dimensionless Reynolds equation can be written as
the dimensionless velocities, Vs, at the surface of the rotating shaft are given by
In a spectral finite difference scheme, the Equation (2b) is decomposed into each component of the Fourier series to the circumferential -direction.
at , the groove shape is symmetric, the boundary conditions of pressure are
and it is assumed that the fluid is open to the atmosphere at , so that boundary conditions of pressure are
4. Fourier Series of Fluid Film Thickness Geometry
The dimensionless fluid film thickness of herringbone grooves journal bearing can be rewritten as
where , is fluid film thickness without eccentricity, E is dimensionless eccentricity as .
1) In the Case of Trapezoidal Groove
As for a trapezoidal groove which was assumed as shown in Figure 3, is defined as trapezoidal angle of groove.
In case of , the fluid film thickness with trapezoidal grooves, , is given by
The Fourier cosine and sine coefficient of fluid film thickness are obtained as
Figure 3. Cross-section perpendicular to the trapezoidal surface of groove.
where is number of grooves.
The “error” in the partial Fourier series at n = N in case of trapezoidal grooves is given as follows,
where is sinc function. Figure 4 plot the maximum error in case of Ng = 5 with the partial Fourier series term number increases, and obviously, the maximum error decays rapidly in case of trapezoidal grooves.
2) In a Special Case at Rectangle Groove
More generally, for the Fourier series expansion of fluid film thickness in case of rectangle groove, the nth partial Fourier series will overshoot this jump by
Figure 4. The maximum error with variation of the partial Fourier series term number in case of trapezoidal grooves.
approximately at a point of groove start or groove end (Gibbs phenomenon), the “error” in the partial Fourier series will be about 8.95% and 14.11% of groove depth larger than the jump in the original fluid film thickness in two cases, in the limits of increasing many terms and of increasingly high node densities, respectively, see Figure 5 by Bengt Fornberg .
How to reduce the Gibbs phenomenon, is an interesting and important topic in mathematics, several mathematician deal with this topic, and obtained some method for its field, e.g. filtering and spectral re-projection. Filtering is a classical tool for mitigating the Gibbs phenomenon in Fourier expansions, however filtering does not completely remove the Gibbs phenomenon. To completely remove the Gibbs phenomenon, one can re-expand the function in a carefully chosen different basis, it is spectral re-projection method, which is given by Gottlieb and Shu . In 1990, Gottlieb, together with Shu and other co-works, showed that it is possible to completely remove the Gibbs phenomenon by post-processing the Fourier expansion in regions in which the function is analytic, using a re-expansion of the partial sums fN(x) in a different basis (“Gibbs complementary”). But understanding that needs high mathematic knowledge, and it is not easy for industrial applications, an easy analysis method always requested from engineer, and so for compensation of its defect, the following method which is proposed.
Since Equation (8) are continuity and differentiable at a point , then Fourier coefficient of fluid film thickness in the case of rectangle groove can been replaced as taking the limit of Equation (8) at trapezoidal angle approaches 0, which are given as
Figure 5. The Gibbs Phenomenon for Fourier series by Bengt Fornberg & Natasha Flyer (2006).
and the maximum “error” at n = N of proposed method which is Equation (9) can be obtained
and since , thus
Figure 6 plots the maximum error in case of Ng = 5 with the partial Fourier series term number increases, and obviously, the maximum error decays rapidly in case of trapezoidal grooves.
Let shows an example of the film thickness at rectangle groove by Fourier’s series using the limit of trapezoidal angle approaches 0 in Figure 7, the “error” in the partial Fourier series in the case of rectangle groove was cut off, and the inexpedience which Fourier series of the film profile does not converge at jump points of groove start or groove end in the case of rectangle groove would be avoided.
5. Analyzed Examples
To confirm the applicability of the above method, Hirs model  which is Ng = 20, β = 21.8 deg., ar/ag = 1.0, L = Rb, Ω = 0, was picked, the relation between load capacity and eccentricity are numerical analyzed, and which compare to experiments of Hirs as shown in Figure 8, it shows that those numerical results are in good agreement with the experimental data.
Figure 6. The maximum error with variation of the partial Fourier series term number in case of rectangle grooves.
Figure 7. The film surface profile at rectangle groove was displayed using the limit of trapezoidal angle approaches 0.
Figure 8. Variations of load capacity with eccentricity increases, where, W is the dimensionless load capacity.
As analysis of a fluid-lubricated herringbone grooved journal bearing under a spectral difference scheme, for compensation of which Fourier series does not converge at jump points in a special case of rectangle groove. Fourier coefficient of fluid film thickness is proposed as taking the limit of which in a case trapezoidal groove at trapezoidal angle approaches 0. In addition, the difference of the film thickness and the number of terms of its Fourier series are investigated with the ratio of the Fourier’s series terms number to grooves number.
A spectral difference method is applied to get numerical solutions for a fluid-lubricated herringbone grooved journal bearing with trapezoidal grooves, and then the numerical analysis scheme will be extended to be suitable for a special case of rectangle groove.
Fourier coefficients of fluid film can be replaced as taking the limit of the trapezoidal angle approaches 0.
The author gratefully acknowledges Emeritus Professor Y. Mochimaru of the Tokyo Institute of Technology for advice.
ag, ar,= grooves width, ridge width
Cg, Cg/b = groove depth, dimensionless groove depth
Cr = bearing clearance
e, E = eccentricity, dimensionless eccentricity
H = dimensionless fluid film thickness
l, L = bearing length, dimensionless bearing length
Ng = number of grooves
p, P = pressure, dimensionless pressure
Pa = atmospheric pressure
r, θ, z = inertial coordinates
Rb = radius of bearing
Rs = radial component of coordinate at surface of shaft
Rs0 = radius of shaft without grooves
Rs/b = dimensionless radial component of coordinate at surface of shaft
t, τ= time, dimensionless time
us = circumferential velocity at surface of rotating shaft
Us = dimensionless circumferential velocity at surface of rotating shaft
vs = radial velocity at surface of rotating shaft
Vs = dimensionless radial velocity at surface of rotating shaft
= velocity components of lubricant fluid
W = dimensionless load capacity of bearing
β = groove angle
η = viscosity of fluid
Λ = bearing number
σ = dimensionless number
= attitude angle of shaft
= trapezoidal angle of groove
= angle between the fixed axis of abscissa (θ = 0) and the axis of eccentricity, dimensionless angle
ω = rotation velocity of shaft
Ω = swirl velocity of shaft
superscript*: non-inertial coordinate
 Bonneau, D. and Absi, J. (1994) Analysis of Aerodynamic Journal Bearing with Small Number of Herringbone Grooves by Finite Element Method. Journal of Tribology, 116, 698-704.
 Liu, J. and Mochimaru, Y. (2010) Analysis of Oil-Lubricated Herringbone Grooved Journal Bearing with Trapezoidal Cross-Section, Using a Spectral Finite Difference Method. Journal of Hydrodynamics, 22, 408-412.
 Liu, J. and Mochimaru, Y. (2012) The Effects of Trapezoidal Groove on a Self-Acting Fluid-Lubricated Herringbone Grooves Journal Bearing. ISRN Tribology, 2013, Article ID: 240239.
 Liu, J. and Mochimaru, Y. (2012) Pressure Distributions Generated along a Self-Acting Fluid-Lubricated Herringbone-Grooved Journal Bearing with Trapezoidal Groove. ISRN Tribology, 2013, Article ID: 560917.
 Fornberg, B. and Flyer, N. (2007) The Gibbs Phenomenon for Radial Basis Functions. In: Jerri, A., Ed., The Gibbs Phenomenon in Various Representations and Applications, Sampling Publishing, Potsdam, NY, 201-224
 Gottlieb, D. and Shu, C.W. (1994 ) On The Gibbs Phenomenon V: Recovering Exponential Accuracy in a Sub-Interval from a Gegenbauer Partial Sum of a Piecewise Analytic Function. NASA Contractor Report 194952, ICASE Report No. 94-61.