Engineering:Reynolds equation

From HandWiki
Short description: Differential equation describing pressure distribution of thin viscous fluids


In fluid mechanics (specifically lubrication theory), the Reynolds equation is a partial differential equation governing the pressure distribution of thin viscous fluid films. It was first derived by Osborne Reynolds in 1886.[1] The classical Reynolds Equation can be used to describe the pressure distribution in nearly any type of fluid film bearing; a bearing type in which the bounding bodies are fully separated by a thin layer of liquid or gas.

General usage

The general Reynolds equation is: [math]\displaystyle{ \frac{\partial}{\partial x}\left(\frac{\rho h^3}{12\mu}\frac{\partial p}{\partial x}\right)+\frac{\partial}{\partial y}\left(\frac{\rho h^3}{12\mu} \frac{\partial p}{\partial y}\right) = \frac{\partial}{\partial x} \left(\frac{\rho h \left( u_a + u_b \right)}{2}\right)+\frac{\partial}{\partial y} \left(\frac{\rho h \left( v_a + v_b \right)}{2}\right)+\rho\left(w_a-w_b\right) - \rho u_a\frac{\partial h}{\partial x} - \rho v_a \frac{\partial h}{\partial y} + h\frac{\partial \rho}{\partial t} }[/math]

Where:

  • [math]\displaystyle{ p }[/math] is fluid film pressure.
  • [math]\displaystyle{ x }[/math] and [math]\displaystyle{ y }[/math] are the bearing width and length coordinates.
  • [math]\displaystyle{ z }[/math] is fluid film thickness coordinate.
  • [math]\displaystyle{ h }[/math] is fluid film thickness.
  • [math]\displaystyle{ \mu }[/math] is fluid viscosity.
  • [math]\displaystyle{ \rho }[/math] is fluid density.
  • [math]\displaystyle{ u, v, w }[/math] are the bounding body velocities in [math]\displaystyle{ x, y, z }[/math] respectively.
  • [math]\displaystyle{ a, b }[/math] are subscripts denoting the top and bottom bounding bodies respectively.

The equation can either be used with consistent units or nondimensionalized.

The Reynolds Equation assumes:

  • The fluid is Newtonian.
  • Fluid viscous forces dominate over fluid inertia forces. This is the principle of the Reynolds number.
  • Fluid body forces are negligible.
  • The variation of pressure across the fluid film is negligibly small (i.e. [math]\displaystyle{ \frac{\partial p}{\partial z} = 0 }[/math])
  • The fluid film thickness is much less than the width and length and thus curvature effects are negligible. (i.e. [math]\displaystyle{ h \ll l }[/math] and [math]\displaystyle{ h \ll w }[/math]).

For some simple bearing geometries and boundary conditions, the Reynolds equation can be solved analytically. Often however, the equation must be solved numerically. Frequently this involves discretizing the geometric domain, and then applying a finite technique - often FDM, FVM, or FEM.

Derivation from Navier-Stokes

A full derivation of the Reynolds Equation from the Navier-Stokes equation can be found in numerous lubrication text books.[2][3]

Solution of Reynolds Equation

In general, Reynolds equation has to be solved using numerical methods such as finite difference, or finite element. In certain simplified cases, however, analytical or approximate solutions can be obtained.[4]

For the case of rigid sphere on flat geometry, steady-state case and half-Sommerfeld cavitation boundary condition, the 2-D Reynolds equation can be solved analytically. This solution was proposed by a Nobel Prize winner Pyotr Kapitsa. Half-Sommerfeld boundary condition was shown to be inaccurate and this solution has to be used with care.

In case of 1-D Reynolds equation several analytical or semi-analytical solutions are available. In 1916 Martin obtained a closed form solution[5] for a minimum film thickness and pressure for a rigid cylinder and plane geometry. This solution is not accurate for the cases when the elastic deformation of the surfaces contributes considerably to the film thickness. In 1949, Grubin obtained an approximate solution[6] for so called elasto-hydrodynamic lubrication (EHL) line contact problem, where he combined both elastic deformation and lubricant hydrodynamic flow. In this solution it was assumed that the pressure profile follows Hertz solution. The model is therefore accurate at high loads, when the hydrodynamic pressure tends to be close to the Hertz contact pressure.[7]

Applications

The Reynolds equation is used to model the pressure in many applications. For example:

  • Ball bearings
  • Air bearings
  • Journal bearings
  • Squeeze film dampers in aircraft gas turbines
  • Human hip and knee joints
  • Lubricated gear contacts

Reynolds Equation adaptations - Average Flow Model

In 1978 Patir and Cheng introduced an average flow model,[8][9] which modifies the Reynolds equation to consider the effects of surface roughness on lubricated contacts. The average flow model spans the regimes of lubrication where the surfaces are close together and/or touching. The average flow model applied "flow factors" to adjust how easy it is for the lubricant to flow in the direction of sliding or perpendicular to it. They also presented terms for adjusting the contact shear calculation. In these regimes, the surface topography acts to direct the lubricant flow, which has been demonstrated to affect the lubricant pressure and thus the surface separation and contact friction.[10]

Several notable attempts have been made to taken additional details of the contact into account in the simulation of fluid films in contacts. Leighton et al.[10] presented a method for determining the flow factors needed for the average flow model from any measured surface. Harp and Salent[11] extended the average flow model by considering the inter-asperity cavitation. Chengwei and Linqing[12] used an analysis of the surface height probability distribution to remove one of the more complex terms from the average Reynolds equation, [math]\displaystyle{ d\bar{h_T}/dh }[/math] and replace it with a flow factor referred to as contact flow factor, [math]\displaystyle{ \phi_h }[/math]. Knoll et al. calculated flow factors, taking into account the elastic deformation of the surfaces. Meng et al.[13] also considered the elastic deformation of the contacting surfaces.

The work of Patir and Cheng was a precursor to the investigations of surface texturing in lubricated contacts. Demonstrating how large scale surface features generated micro-hydrodynamic lift to separate films and reduce friction, but only when the contact conditions support this.[14]

The average flow model of Patir and Cheng,[8][9] is often coupled with the rough surface interaction model of Greenwood and Tripp[15] for modelling of the interaction of rough surfaces in loaded contacts.[10][16]

References

  1. Reynolds, O. (1886). "On the Theory of Lubrication and Its Application to Mr. Beauchamp Tower's Experiments, Including an Experimental Determination of the Viscosity of Olive Oil". Philosophical Transactions of the Royal Society of London (Royal Society) 177: 157–234. doi:10.1098/rstl.1886.0005. https://archive.org/details/philtrans08214427. 
  2. Hamrock, Bernard J.; Schmid, Steven R.; Jacobson, Bo O. (2004). Fundamentals of Fluid Film Lubrication. Taylor & Francis. ISBN 978-0-8247-5371-9. https://books.google.com/books?id=aJ0engEACAAJ. 
  3. Szeri, Andras Z. (2010). Fluid Film Lubrication. Cambridge University Press. ISBN 978-0-521-89823-2. https://books.google.com/books?id=vaeMmAEACAAJ. 
  4. "Reynolds Equation: Derivation and Solution". tribonet.org. 12 November 2016. http://www.tribonet.org/wiki/reynolds-equation/. 
  5. Akchurin, Aydar (18 February 2016). "Analytical Solution of 1D Reynolds Equation". tribonet.org. http://www.tribonet.org/cmdownloads/analytical-solution-of-1d-reynolds-equation/. 
  6. Akchurin, Aydar (22 February 2016). "Semi-Analytical Solution of 1D Transient Reynolds Equation(Grubin's Approximation)". tribonet.org. http://www.tribonet.org/cmdownloads/semi-analytical-solution-of-1d-transient-reynolds-equationgrubins-approximation/. 
  7. Akchurin, Aydar (4 January 2017). "Hertz Contact Calculator". tribonet.org. http://www.tribonet.org/cmdownloads/hertz-contact-calculator/. 
  8. 8.0 8.1 Patir, Nadir; Cheng, H. S. (1978). "An Average Flow Model for Determining Effects of Three-Dimensional Roughness on Partial Hydrodynamic Lubrication". Journal of Lubrication Technology 100 (1): 12. doi:10.1115/1.3453103. ISSN 0022-2305. 
  9. 9.0 9.1 Patir, Nadir; Cheng, H. S. (1979-04-01). "Application of Average Flow Model to Lubrication Between Rough Sliding Surfaces". Journal of Lubrication Technology 101 (2): 220–229. doi:10.1115/1.3453329. ISSN 0022-2305. https://doi.org/10.1115/1.3453329. 
  10. 10.0 10.1 10.2 Leighton (2016). "Surface-specific flow factors for prediction of friction of cross-hatched surfaces" (in en). Surface Topography: Metrology and Properties 4 (2): 025002. doi:10.1088/2051-672x/4/2/025002. 
  11. Harp, Susan R.; Salant, Richard F. (2000-10-17). "An Average Flow Model of Rough Surface Lubrication With Inter-Asperity Cavitation". Journal of Tribology 123 (1): 134–143. doi:10.1115/1.1332397. ISSN 0742-4787. https://doi.org/10.1115/1.1332397. 
  12. Wu, Chengwei; Zheng, Linqing (1989-01-01). "An Average Reynolds Equation for Partial Film Lubrication With a Contact Factor". Journal of Tribology 111 (1): 188–191. doi:10.1115/1.3261872. ISSN 0742-4787. https://doi.org/10.1115/1.3261872. 
  13. Meng, F-M; Wang, W-Z; Hu, Y-Z; Wang, H (2007-07-01). "Numerical analysis of combined influences of inter-asperity cavitation and elastic deformation on flow factors". Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 221 (7): 815–827. doi:10.1243/0954406jmes525. ISSN 0954-4062. https://journals.sagepub.com/doi/10.1243/0954406JMES525. 
  14. Morris, N; Leighton, M; De la Cruz, M; Rahmani, R; Rahnejat, H; Howell-Smith, S (2014-11-17). "Combined numerical and experimental investigation of the micro-hydrodynamics of chevron-based textured patterns influencing conjunctional friction of sliding contacts". Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology 229 (4): 316–335. doi:10.1177/1350650114559996. ISSN 1350-6501. https://journals.sagepub.com/doi/full/10.1177/1350650114559996. 
  15. Greenwood, J. A.; Tripp, J. H. (June 1970). "The Contact of Two Nominally Flat Rough Surfaces". Proceedings of the Institution of Mechanical Engineers 185 (1): 625–633. doi:10.1243/pime_proc_1970_185_069_02. ISSN 0020-3483. http://dx.doi.org/10.1243/pime_proc_1970_185_069_02. 
  16. Leighton, M; Nicholls, T; De la Cruz, M; Rahmani, R; Rahnejat, H (2016-12-12). "Combined lubricant–surface system perspective: Multi-scale numerical–experimental investigation". Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology 231 (7): 910–924. doi:10.1177/1350650116683784. ISSN 1350-6501. https://journals.sagepub.com/doi/full/10.1177/1350650116683784.