Reynolds equation

Reynolds equation. It is a fundamental partial differential equation in the field of tribology and fluid film lubrication, describing the pressure distribution within a thin fluid film between two surfaces in relative motion. Formulated by Osborne Reynolds in 1886, it is derived from the Navier–Stokes equations and the continuity equation under specific simplifying assumptions. The equation is central to the design and analysis of journal bearings, thrust bearings, and other lubrication systems, providing the theoretical basis for predicting load capacity, friction, and film thickness.

● Derivation

The derivation begins with the Navier–Stokes equations, which govern fluid motion, and the continuity equation for mass conservation. Key assumptions are applied to model a thin film: the fluid is Newtonian, incompressible, and laminar; inertia and body forces are negligible compared to viscous forces; and the pressure is constant across the film thickness. The no-slip condition is enforced at the boundaries with the moving surfaces, such as a rotating shaft in a bearing. Integrating the simplified momentum equations across the film and applying the boundary conditions yields an expression relating the pressure gradient to the flow. Substituting this into the integrated continuity equation results in the final form, which balances Poiseuille flow from pressure gradients with Couette flow from surface motion and squeeze film effects from changing gap height.

● Forms and simplifications

The general form is a complex partial differential equation for pressure as a function of spatial coordinates. For practical engineering, it is often simplified. The most common is the classical form for an incompressible, isoviscous fluid, which neglects fluid compressibility and thermal effects. For one-dimensional flow, such as in an infinitely long journal bearing, it reduces to an ordinary differential equation known as the infinitely long bearing approximation. The short bearing approximation, attributed to Sommerfeld and later Kingsbury, assumes pressure gradients are dominant in one direction. For gas bearings, the compressible form must be used, where density is coupled to pressure via an equation of state, often the ideal gas law. Other specialized forms include those for elastohydrodynamic lubrication, which couples the equation with the elasticity equations of the contacting solids.

● Applications

Its primary application is in the analysis and design of fluid film bearings, including plain bearings, tilting-pad bearings, and foil bearings used in turbomachinery like gas turbines and centrifugal compressors. It is essential for predicting the performance characteristics of internal combustion engine components, such as piston rings and camshaft followers. In magnetic storage devices, it models the air bearing between the read/write head and the disk surface. The equation is also applied in biomechanics to study synovial fluid lubrication in human joints like the hip joint and knee joint. Furthermore, it underpins the design of mechanical seals and gear lubrication systems in automotive engineering and aerospace engineering.

● Limitations and extensions

The classical equation has significant limitations, as it assumes a rigid, smooth surface and constant viscosity. Real contacts often involve surface roughness, which disrupts the thin film; this led to the development of average flow models by researchers like Tichy and Patir. In elastohydrodynamic lubrication, where pressures are high enough to deform the surfaces and alter fluid properties, the equation is coupled with the elasticity equations of Hertzian contact theory and a pressure-viscosity relationship. For thermohydrodynamic lubrication, the energy equation and temperature-dependent viscosity must be incorporated, as studied by institutions like the MIT. Extensions also address non-Newtonian fluid behavior, turbulence in high-speed bearings, and multiphase flow with cavitation, the latter often modeled using the Jakobsson–Floberg–Olsson boundary condition.

● Numerical solution methods

Analytical solutions exist only for highly simplified geometries, making numerical analysis essential for most practical problems. The primary method is the finite difference method, which discretizes the domain into a grid, as implemented in early codes from NASA. The finite element method, supported by software like ANSYS and COMSOL Multiphysics, offers greater flexibility for complex geometries and coupled physics, such as with structural analysis. For gas bearings with strong compressibility, techniques like the successive over-relaxation method are often employed. The emergence of computational fluid dynamics packages has integrated its solution into broader multiphysics simulation environments. Key challenges in numerical solution include handling the Reynolds boundary condition for cavitation and achieving stable convergence in highly nonlinear regimes like elastohydrodynamic lubrication.

Category:Equations of fluid dynamics Category:Tribology Category:Partial differential equations