Research Article  Open Access
Mathematical Model and Analysis of the WaterLubricated Hydrostatic Journal Bearings considering the Translational and Tilting Motions
Abstract
The waterlubricated bearings have been paid attention for their advantages to reduce the power loss and temperature rise and increase load capacity at high speed. To fully study the complete dynamic coefficients of two waterlubricated, hydrostatic journal bearings used to support a rigid rotor, a fourdegreeoffreedom model considering the translational and tilting motion is presented. The effects of tilting ratio, rotary speed, and eccentricity ratio on the static and dynamic performances of the bearings are investigated. The bulk turbulent Reynolds equation is adopted. The finite difference method and a linear perturbation method are used to calculate the zeroth and firstorder pressure fields to obtain the static and dynamic coefficients. The results suggest that when the tilting ratio is smaller than 0.4 or the eccentricity ratio is smaller than 0.1, the static and dynamic characteristics are relatively insensitive to the tilting and eccentricity ratios; however, for larger tilting or eccentricity ratios, the tilting and eccentric effects should be fully considered. Meanwhile, the rotary speed significantly affects the performance of the hydrostatic, waterlubricated bearings.
1. Introduction
Hydrostatic journal bearings are applied widely in spindlebearing systems owning to their favorable performance characteristics. However, with the requirement of higher machining speed, the limitations of the conventional oil film bearings are apparently due to their remarkable power loss as well as the temperature rise. Therefore, the waterlubricated bearings were developed and have been studied to fulfill the targets of lower power loss, lower temperature rise, and heavier load capacity at high speed.
Many studies related to waterlubricated bearings have been reported in the literatures in the past few years. Liu et al. compared the oillubricated and waterlubricated hybrid sliding bearings, and the results show that the latter benefits more from improved processing precision and efficiency [1]. Yuan et al. study the static and dynamic characteristics of waterlubricated hybrid journal bearings compensated by short capillaries [2]. Yoshimoto et al. investigated the static characteristics of waterlubricated hydrostatic conical bearings with spiral grooves for high speed spindles [3]. Gao et al. analyzed the effects of eccentricity ration on pressure distribution of waterlubricated plain journal bearings by computational fluid dynamics (CFD) [4]. In summary, extensive researches have been conducted in the area of water lubricated bearings in various aspects: numerical methods [1]; performance of the bearings with various geometries [2, 3, 5, 6]; effects of various kinds of restrictors upon performance of a bearing [7], et al. However, their studies were restricted to the static or dynamic characteristics of plain or grooved journal bearings considering only translational motion of the journal.
In actual practice, the bearings and the journals may not be properly aligned as a result of improper assembly or noncentral loading. As a consequence, not only should the translational motion of the rotor be studied, but also the tilting motion of the rotor should be investigated. As a result, the complete stiffness and damping coefficients of a journal bearing, which is important to the vibration of a rotor, should be taken into consideration in four degrees of freedom, including the translation in , direction and tilting about the  and axis. Numerous studies concerning the tilting motion of journal are available in the literatures [8–10]. Recently, Jang et al. [11–13] thoroughly investigated the dynamic characteristics of the journal and groove thrust bearings used to support a HDD spindle considering both of the translational motion and tilting motion of the journal. Results show that the tilting motions have an important role in the dynamic characteristics of the proposed bearings.
Unlike the conventional oil film bearings, waterlubricated bearings utilized in spindle are different in working conditions and characteristics. However, for hydrostatic waterlubricated journal bearings, we are not aware of any previous investigations to study their dynamic characteristics considering the translational and titling motions. Therefore, in this work, we aim to fully study the complete dynamic coefficients for two waterlubricated, hydrostatic journal bearings used to support a rigid rotor. The dynamic characteristics will be categorized into four groups: coefficients of force to displacement, coefficients of force to angle, coefficients of moment to angle, and coefficients of moment to displacement. In the present study, in order to fully study the variations of the complete static and dynamic characteristics of the proposed waterlubricated bearings, the influences of the tilting ratio, rotary speed, and eccentricity ratio on the bearings have been studied.
2. Mathematical Models
Figure 1 shows the schematic representations for a rigid rotor supported by a pair of identical waterlubricated journal bearings, as well as the geometry of a hydrostatic, waterlubricated journal bearing. Pressed water enters the bearing across an orifice restrictor, flows into the film lands, and then exits the bearing. The recess pressure is regarded as uniform. As shown in Figure 1, the rigid rotor moved in an inertial reference frame ; the rotor tilts about by an angle and tilts about by an angle . The rigid rotorbearing model is shown in Figure 2.
2.1. Reynolds Equation
For an isoviscous, incompressible fluid, the Reynolds equation governing the turbulent bulk flow in nondimensional form is given as where , , , and the details of the turbulent coefficients , can be obtained as follows [6, 14–16]:
With the increase of rotary speed or film depth, the flow is likely to become turbulent from laminar state. The turbulent coefficients , , and dependent on the fluid velocity field are obtained as follows [15]:
At the bearing exit plane, the pressure takes a constant value equal to the ambient pressure.
2.2. Continuity Equation
The dimensionless continuity equation at the recess is defined by the global balance between the flow through the orifice restrictor and the recess outflow into the film lands: where, , is the recess pressure, is the circumferential recess boundary, and is the axial recess boundary.
2.3. Perturbation Analysis
The journal center rotates about its steadystate position with a small whirl which is generated from the variations due to the translations of the rotor mass center and variations due to the tilting angles [11–13]. For small amplitude motions, the dimensionless film thickness and pressure fields are expressed as the sum of a zerothorder field and firstorder field, describing the steadystate condition and perturbed motion, respectively.
The dimensionless perturbed film expression considering the tilting angles is [12] where where is the distance between the grid node of each bearing and the mass center of the rotor supported on the bearings, for the front journal bearing, and for the rear journal bearing.
The dimensionless perturbed pressure expression is
Substitution of the perturbed equations (5)–(7) into the Reynolds equation yields the zeroth and firstorder expressions:
The perturbed quality into the orifice diameter can be obtained by Taylor expansion:
Substitution of the perturbed equations (5)–(7) into the continuity equation yields the zeroth and firstorder expressions:
The finite difference method (FDM) and a successive overrelaxation (SOR) scheme are implemented to solve (1)–(10) to find the pressure distribution. When the steady and perturbed pressure distributions are obtained, the static and dynamic coefficients can be solved.
2.4. Static Characteristics
The quality, frictional power loss, and pump power are calculated by integration of the pressure field on the bearing surfaces:
San Andres et al. carried out a systematic research on the waterlubricated hydrostatic journal bearings both theoretically and experimentally [17]. Numerical and experimental results show that predictions of the bearing performance characteristics like flow rate, load capacity, and rotor dynamic force coefficients are not affected by the small temperature variations () in the water hydrostatic journal bearings. As a result, an adiabatic, isothermal assumption is made in this study. All the heat produced in the bearings by friction is considered absorbed by water film. The average elevated temperature is given by where, is the density of water and is the specific heat capacity of water.
2.5. Dynamic Characteristics
The dynamic coefficients of the journal bearings can be calculated by integrating the perturbed pressure across the fluid film. There are 16 stiffness and damping coefficients, and the dynamic coefficients can be grouped into four categories: coefficients of force to displacement, coefficients of force to angle, coefficients of moment to angle, and coefficients of moment to displacement:
2.6. The Numerical Solution Procedure
The finite difference method (FDM) and a successive overrelaxation (SOR) scheme are implemented to solve (1)–(10) governing the flow in the film to find the pressure distribution. The overrelaxation factor always lies between (1~2). By trial and error, one can determine an optimum value of the relaxation factor for the fastest convergence. Normally, 1.7 is a good starting point for determining the relaxation factor [18]. In the numerical procedure implemented, any negative pressures calculated in the cavitation zone are arbitrarily set equal to zero (or ambient) pressure [18]. The pressure iterations are continued until the following convergence criterion is satisfied: where is the convergence criteria, which is set as . Once the pressure field is established for the water film, other performance parameters follow from the pressure distribution. The steady and perturbed pressure distributions obtained are subsequently integrated to yield the desired static and dynamic coefficients.
The fluid film was discretized by rectangular grid with unequal intervals, which are 0.033, 0.088, and 0.1154 for the film land, recess, and return groove, respectively, in the circumferential direction and 0.02 and 0.04 for the film land and recess area in the axial direction. The total number of the grid is 31 × 102. It takes about 35 seconds for the static performance to achieve convergence. A few validation tests were made with a coarser grid of 31 × 86 and finer grid of 37 × 158 with different intervals, and in no case did the predicted results of the static characteristics differ by more than 0.1 percent from those obtained by the initial grid.
3. Results and Discussion
3.1. Comparisons of Present Solution with Experimental Results
The present numerical solution has been correlated and validated with the experimental results available in the literature. A fiverecess, turbulentflow, waterlubricated hydrostatic bearing operating at a high rotational speed is tested by San Andres et al. [17]. Table 1 shows the bearing description and operating conditions and Figure 3 shows the schematic of the bearing. Figure 4 shows a comparison of measured load capacity and flow rate and the numerical predictions in respect to eccentricity ratio. As shown in Figure 4, a good agreement is observed between the present numerical predictions and the experimental results available in the reference. The calculated load capacity results correlate well with the experimental results, with a maximum difference of 11.5%. The average of the flow rate for the experimental results is about 1.2 kg/s, while that for the numerical results is approximately 1.45 kg/s.

(a) Comparison of load capacity
(b) Comparison of mass flow rate
3.2. Effects of Tilting Ratio on the Static Performances of the Bearings
The geometric parameters for the bearings have been presented in Figure 1, the supply pressure is 1.5 MPa, the orifice diameter is 0.6 mm, the pressurized water is supplied at a temperature of 30°C, and the viscosity is 0.00087 .
Figure 1 shows the coordinate system for the tilting motion of the rotor. As the film variation of each bearing varies with different values for the tilting angles and bearing span, the tilting ratio is proposed here to investigate the effect of misalignment. The tilting ratio is defined as follows: where is the distance between the titling center and the left edge of the front bearing or the right edge of the rear bearing, which is shown in Figure 1.
In order to obtain a better physical insight into the effect of misalignment, the static characteristics have been presented with eccentricity ratio equal to zero. Figures 5(a)–5(c) show the calculated quality, power loss, and temperature rise for different values of tilting ratios for each journal bearing operating at 10000 rpm. It is observed that the quality of each bearing undergoes a reduction by only 6.3% when the tilting ratio increases from zero to 0.75. The total power loss of each bearing keeps almost unchanged when the tilting ratio is not greater than 0.4. However, when the tilting ratio continues to increase, the power loss gets a slight increase by approximately 5%. Meanwhile, the temperature rise increases with increased tilting ratio. The reason for this is that increasing tilting angles will increase the film variation of each bearing, which promote the hydrodynamic and turbulent effect. It is noticed that the maximum temperature rise for the investigated waterlubricated bearing is 1.7°C for the rotorbearing system investigated here, which is far lower than the conventional oil film bearings. This in turn verifies the validation of the adiabatic assumption. In conclusion, when the tilting ratio is smaller than 0.4, the influence of tilting ratio on the static performance of a waterlubricated hydrostatic journal bearing can be ignored.
(a) The effect of tilting angle on the quality
(b) The effect of tilting angle on the power loss
(c) The effect of tilting angle on the temperature rise
3.3. Effects of Tilting Ratio on the Dynamic Characteristics of the Bearings
Figures 6 and 7 show the variation of the dynamic coefficients of the two identical bearings in respect to the tilting ratio with a rotary speed of 10000 rpm. The eccentricity ratio is assumed to be zero to exclude the effect of the clearance change due to the eccentricity ratio. The solid lines represent the dynamic coefficients of the front journal bearing while the dotted lines represent those of the rear bearing.
(a) The stiffness of force to displacement versus tilting ratio
(b) The stiffness of force to angle versus tilting ratio
(c) The stiffness of moment to angle versus tilting ratio
(d) The stiffness of moment to displacement versus tilting ratio
(a) The damping of force to displacement versus tilting ratio
(b) The damping of force to angle versus tilting ratio
(c) The damping of moment to angle versus tilting ratio
(d) The damping of moment to displacement versus tilting ratio
According to the results, the dynamic coefficients of force to displacement and moment to angle for the front bearing are approximately equal to those for the rear bearing due to the fact that the tilting center almost coincides with the bearing span center. However, the coupled dynamic coefficients of force to angle and moment to displacement for the two bearings have close proximity magnitudes but in the opposite direction due to the fact that the two bearings are arranged at both sides of the mass center which is the origin of the tilting motions. It is observed that, with the increase of the tilting ratio, the stiffness and damping coefficients of the bearings keep almost independent of the tilting ratio when it is not larger than 0.4; however, when the tilting ratio continues to increase, the effect on dynamic coefficients is significant, with a maximum variation rate of 16.7% for the stiffness and 49.2% for the damping coefficients. Furthermore, the larger the tilting ratio is, the greater the differences among the coefficients are. This can be ascribed to the variations of the film thickness induced by the tilting angles and bearing span, which can be calculated according to (5). The results indicate that, for a small tilting ratio, the effect of the misalignment on the performance of hydrostatic waterlubricated journal bearings can be ignored; but for a relatively greater tilting ratio, the effect of misalignment should be taken into consideration.
3.4. Effects of Rotary Speed on the Static Performances of the Bearings
Figure 8 depicts the variation of the static characteristics for each journal bearing in respect to the rotary speed. In the case of a waterlubricated, hydrostatic bearing operating with eccentricity ratio and tilting angles equal to zero, it may be observed that the static characteristics of the front bearing are the same as those of the rear bearing. As shown in Figure 8(a), the quality of each bearing apparently remains unchanged with increased rotary speed. As shown in Figures 8(b) and 8(c), the gross power loss for each bearing increases dramatically from 123 W to 1232 W when the rotary speed increases from 5000 rpm to 30000 rpm because the frictional power loss is closely related to the rotary speed according to (12). Furthermore, the temperature rise undergoes a sharp increase from 0.7°C to 8°C. It should be pointed out that most of the predicted temperature rises are higher than the actual values presumably due to the adiabatic assumption imposed on the analysis. Considering the values of the temperature rise, the adiabatic flow assumption is fully justified for the bearing studied. However, when the rotary speed continues to increase, the energy equation should be included to predict the temperature rise precisely.
(a) The effect of rotary speed on the quality
(b) The effect of rotary speed on the power loss
(c) The effect of rotary speed on the temperature rise
3.5. Effects of Rotary Speed on the Dynamic Characteristics of the Bearings
Figures 9 and 10 show the dynamic coefficients of each waterlubricated bearing in respect to the rotary speed. The eccentricity ratio and tilting angles are assumed as zero to exclude their influences on the film thickness. According to the results, a higher rotary speed generates larger coupled stiffness of force to displacement and moment to angle, but the influence of rotary speed on the relative direct coefficients is small. The direct stiffness coefficients of force to angle and moment to displacement are relatively insensitive to the variation of the rotary speed. The magnitudes of the crosscoupled stiffness of force to displacement and moment to angle are comparable to those of the direct stiffness, which demonstrates the importance of hydrodynamic effects. Unlike the stiffness, the damping coefficients are independent of the rotary speed in the aligned condition. The reason for this is that the perturbed pressure due to the perturbed velocity is not related to the rotary speed of the rotor.
(a) The stiffness of force to displacement versys rotary speed
(b) The stiffness of force to angle versus rotary speed
(c) The stiffness of moment to angle versus rotary speed
(d) The stiffness of moment to displacement versus rotary speed
(a) The damping of force to displacement versus rotary speed
(b) The damping of force to angle versus rotary speed
(c) The damping of moment to angle versus rotary speed
(d) The damping of moment to displacement versus rotary speed
3.6. Effects of Eccentricity Ratio on the Static Performances of the Bearings
The influence of eccentricity ratio in aligned condition on static performances of each bearing is as shown in Figure 11. The tilting ratio is assumed as zero to exclude the influence of tilting effect. Figure 11(a) indicates that the value of quality for each bearing decreases slowly with the eccentricity ratio. As is shown in Figure 11(b), the power loss is almost constant at first and then increases with increased eccentricity ratio. The maximum temperature rise across the bearing length is about 1.6°C at the eccentricity ratio equal to 0.5. This is expected since a smaller clearance produces a larger frictional power loss and a smaller flow rate. However, it should be pointed out that in the aligned condition the static characteristics of the hydrostatic waterlubricated journal bearings vary slightly with eccentricity ratio.
(a) The effect of rotary speed on the quality
(b) The effect of rotary speed on the power loss
(c) The effect of rotary speed on the temperature rise
3.7. Effects of Eccentricity Ratio on the Dynamic Characteristics of the Bearings
Figures 12 and 13 show a variation of stiffness and damping coefficients of each bearing versus the eccentricity ratio. It may be noticed that the dynamic coefficients are almost constant as the eccentricity ratio increases from 0 to 0.1. The direct stiffness of force to displacement and moment to angle decreases gradually with increased eccentricity ratio. Generally, the coupled stiffness varies significantly with the eccentricity ratio. The coupled coefficients of force to angle and moment to displacement for the front bearing have the same magnitude as those for the rear bearing, but in the opposite direction. The damping coefficients are also relatively insensitive to the eccentricity ratio when it is not larger than 0.1. However, in larger eccentric condition, the damping coefficients vary with eccentricity ratio, and the larger the eccentricity ratio is, the greater the coefficients change. In summary, for the small eccentric condition (≤0.1), the influence of eccentricity ratio on the full dynamic coefficients for the hydrostatic, waterlubricated journal bearing operating in aligned condition can be ignored; however, for a larger eccentric condition, the influence should be fully discussed.
(a) The stiffness of force to displacement versus eccentricity ratio
(b) The stiffness of force to angle versus eccentricity ratio
(c) The stiffness of moment to angle versus eccentricity ratio
(d) The stiffness of moment to displacement versus eccentricity ratio
(a) The damping of force to displacement versus eccentricity ratio
(b) The damping of force to angle versus eccentricity ratio
(c) The damping of moment to angle versus eccentricity ratio
(d) The damping of moment to displacement versus eccentricity ratio
4. Conclusion
This paper investigated the complete dynamic coefficients for two hydrostatic, waterlubricated journal bearings used to support a rigid rotor considering the translational and tilting motion. The bulk turbulent flow model and FDM method is used to numerically predict the performance of the bearings. The results show that the proposed waterlubricated hydrostatic bearing has the potential to fulfill the target of lower power loss, temperature rise, and larger load capacity at high speed. On the basis of the results presented, the following conclusions can be drawn.(1)For a small tilting ratio (<0.4), the influence of tilting ratio on the static and dynamic characteristics of a waterlubricated hydrostatic journal bearing is relatively small; however, when the tilting ratio continues to increase, the power loss and temperature rise increase gradually while the quality decreases, and the effect of tilting ratio on the dynamic coefficients should be taken into consideration.(2)The quality of the bearings is relatively insensitive to the rotary speed; however, the power loss and temperature rise increase sharply with the rotary speed in an aligned condition. The direct stiffness coefficients vary significantly with the rotary speed due to the hydrodynamic effect while the damping coefficients are almost constant.(3)For a relatively smaller eccentric condition (≤0.1), the static and dynamic characteristics of the hydrostatic waterlubricated journal bearings vary slightly with eccentricity ratio. However, for a larger eccentric condition, the dynamic characteristics increase or decrease significantly with the eccentricity ratio.
Nomenclature
:  Design film thickness 
:  Orifice diameter 
:  Film thickness 
:  Distance between the mass center and front journal bearing center 
:  Distance between the mass center and rear journal bearing center 
:  Distance between the mass center and the left edge of front journal bearing 
:  Distance between the mass center and the right edge of rear journal bearing 
:  Radius of the bearing 
:  Diameter of a journal bearing 
:  The bearing force 
:  The length of a journal bearing 
:  Supply pressure 
:  The external load 
:  Rotary speed 
:  Flow coefficient 
:  Density 
:  Viscosity 
:  Orifice design coefficient. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 F. Liu, B. Lin, and X. F. Zhang, “Numerical design method for waterlubricated hybrid sliding bearings,” International Journal of Precision Engineering and Manufacturing, vol. 9, no. 1, pp. 47–50, 2008. View at: Google Scholar
 X. Y. Yuan, G. Y. Zhang, B. Li, and X. Miao, “Theoretical and experimental results of waterlubricated, highspeed, shortcapillarycompensated hybrid journal bearings,” in Proceedings of IJTC, STLE/ASME International Joint Tribology Conference, pp. 391–398, San Antonio, Tex, USA, October 2006. View at: Publisher Site  Google Scholar
 S. Yoshimoto, T. Kume, and T. Shitara, “Axial load capacity of waterlubricated hydrostatic conical bearings with spiral grooves for high speed spindles,” Tribology International, vol. 31, no. 6, pp. 331–338, 1998. View at: Publisher Site  Google Scholar
 G. Y. Gao, Z. W. Yin, D. Jiang et al., “Numerical analysis of plain journal bearing under hydrodynamic lubrication by water,” Tribology International, vol. 75, pp. 31–38, 2014. View at: Google Scholar
 S. Nakano, T. Kishibe, T. Inoue, and H. Shiraiwa, “An advanced microturbine system with waterlubricated bearings,” International Journal of Rotating Machinery, vol. 2009, Article ID 718107, 12 pages, 2009. View at: Publisher Site  Google Scholar
 H. H. Feng, “Dynamic characteristics of a rigid spindle supported by waterlubricated bearings,” Applied Mechanics and Materials, vol. 401–403, pp. 121–124, 2013. View at: Google Scholar
 J. Corbett, R. J. Almond, D. J. Stephenson, and Y. B. P. Kwan, “Porous ceramic water hydrostatic bearings for improved for accuracy performance,” CIRP Annals: Manufacturing Technology, vol. 47, no. 1, pp. 467–470, 1998. View at: Publisher Site  Google Scholar
 O. Ebrat, Z. P. Mourelatos, N. Vlahopoulos, and K. Vaidyanathan, “Calculation of journal bearing dynamic characteristics including journal misalignment and bearing structural deformation,” Tribology Transactions, vol. 47, no. 1, pp. 94–102, 2004. View at: Publisher Site  Google Scholar
 S. C. Jain, S. C. Sharma, and T. Nagaraju, “Misaligned journal effects in liquid hydrostatic nonrecessed journal bearings,” Wear, vol. 210, no. 12, pp. 67–75, 1997. View at: Publisher Site  Google Scholar
 S. C. Sharma, V. M. Phalle, and S. C. Jain, “Combined influence of wear and misalignment of journal on the performance analysis of threelobe threepocket hybrid journal bearing compensated with capillary restrictor,” Journal of Tribology, vol. 134, no. 1, Article ID 011703, 2012. View at: Publisher Site  Google Scholar
 G. H. Jang and J. W. Yoon, “Dynamic characteristics of a coupled journal and thrust hydrodynamic bearing in a HDD spindle system due to its groove location,” Microsystem Technologies, vol. 8, no. 45, pp. 261–270, 2002. View at: Publisher Site  Google Scholar
 G. H. Jang and S. H. Lee, “Determination of the dynamic coefficients of the coupled journal and thrust bearings by the perturbation method,” Tribology Letters, vol. 22, no. 3, pp. 239–246, 2006. View at: Publisher Site  Google Scholar
 H. W. Kim, G. H. Jang, and S. H. Lee, “Complete determination of the dynamic coefficients of coupled journal and thrust bearings considering five degrees of freedom for a general rotorbearing system,” Microsystem Technologies, vol. 17, no. 5–7, pp. 749–759, 2011. View at: Publisher Site  Google Scholar
 G. G. Hirs, “A bulkflow theory for turbulence in lubricant films,” Journal of Lubrication Technology, vol. 95, no. 2, pp. 137–146, 1973. View at: Publisher Site  Google Scholar
 J. Frêne, M. Arghir, and V. Constantinescu, “Combined thinfilm and NavierStokes analysis in high Reynolds number lubrication,” Tribology International, vol. 39, no. 8, pp. 734–747, 2006. View at: Publisher Site  Google Scholar
 R. Bassani, E. Ciulli, B. Piccigallo, M. Pirozzi, and U. Staffilano, “Hydrostatic lubrication with cryogenic fluids,” Tribology International, vol. 39, no. 8, pp. 827–832, 2006. View at: Publisher Site  Google Scholar
 L. San Andres, D. Childs, and Z. Yang, “Turbulentflow hydrostatic bearings: analysis and experimental results,” International Journal of Mechanical Sciences, vol. 37, no. 8, pp. 815–829, 1995. View at: Publisher Site  Google Scholar
 M. M. Khonsari and E. R. Booser, Applied Tribology, WileyInterscience, London, UK, 2008.
Copyright
Copyright © 2014 HuiHui Feng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.