Revision for “Reynolds Equation: Derivation and Solution” created on April 27, 2017 @ 13:28:25 [Autosave]
Title  Reynolds Equation: Derivation and Solution 

Content  Reynolds equation is a partial differential equation which describes the flow of a thin lubricant film between two surfaces. It is derived from the NavierStokes equations and is one of the fundamental equations of the classical lubrication theory. It was first derived by <a href="http://discovery.nationalarchives.gov.uk/details/rd/9c1e8ebbbaae41b28c13799b8ea03dcf">Osborne Reynolds in 1886</a>.
<h1>Derivation of Reynolds Equation</h1>
<p style="textalign: justify;">The principles of the theory are derived from the observation that the lubricant can be treated as isoviscous and laminar and the fluid film is of negligible curvature. The classical Reynold's equation can be derived from the NavierStokes equations and the equation of continuity under assumptions of:</p>
<ul>
<li>constant viscosity, Newtonian lubricant</li>
<li>thin film geometry</li>
<li>negligible body force</li>
<li>noslip boundary conditions</li>
</ul>
When these assumptions are applied, we obtain following equations (following <a href="https://books.google.nl/books/about/Fluid_Film_Lubrication.html?id=U5VXXRAgKwsC&redir_esc=y">Szeri</a>):
[math] \begin{eqnarray}
\label{complete_sys1}
\frac{\partial p}{\partial x} = \mu\frac{\partial^{2} u}{\partial y^{2}} \\
\label{complete_sys2}
\frac{\partial p}{\partial z} = \mu\frac{\partial^{2} w}{\partial y^{2}} \\
\label{complete_sys3}
\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0
\end{eqnarray} [/math]
<p style="textalign: left;">It is also assumed the lubricant is incompressible here. The third of these equations is the equation of continuity and will be used later in the derivation. The first two equations now can be integrated twice with respect to [math] y [/math], since partial derivatives of pressure [math]p [/math] dont vary across [math] y [/math] (assumption of thin film geometry):</p>
<p style="textalign: left;">[math] \begin{eqnarray}
\label{complete_sys1}
u = \frac{1}{2\mu}\frac{\partial p}{\partial x}y^2 + Ay + B \\
\label{complete_sys2}
w = \frac{1}{2\mu}\frac{\partial p}{\partial z}y^2 + Cy + D
\end{eqnarray} [/math]</p>
<p style="textalign: left;">Boundary conditions should be applied:
[math] \begin{eqnarray}
\label{complete_sys1}
u = U_1, w = 0, y = 0, \\
\label{complete_sys2}
u = U_2, w = 0, y = h,
\end{eqnarray} [/math]</p>
where [math]U_1 [/math] and [math]U_2 [/math] represent the velocity of the bearing surfaces.
Evaluation of the integration constants leads to the following velocity distribution:
[math]
\begin{eqnarray}
\label{complete_sys1}
u = \frac{1}{2\mu}\frac{\partial p}{\partial x}(y^2  yh) + (1  \frac{y}{h})U_1 + \frac{y}{h}U_2 \\
\label{complete_sys2}
w = \frac{1}{2\mu}\frac{\partial p}{\partial z}(y^2  yh)
\end{eqnarray} [/math]
<p style="textalign: left;">There are only tree equations for four unknowns. This difficulty will be alleviated by integrating, in effect averaging, the equation of continuity across the film:</p>
<p style="textalign: center;">[math] \int_{0}^{h} \frac{\partial v}{\partial y} dy = \int_{0}^{h} {\frac{\partial u}{\partial x}dy + \int_{0}^{h} \frac{\partial w}{\partial z}} dy [/math]</p>
<p style="textalign: left;">Integrating yields following Reynold's equation for lubricant pressure:
[math]
\begin{eqnarray}
\label{complete_sys1}
\frac{\partial}{\partial x}(\frac{h^3}{\mu}\frac{\partial p}{\partial x}) + \frac{\partial}{\partial z}(\frac{h^3}{\mu}\frac{\partial p}{\partial z})= 6(U_1+U_2)\frac{\partial h}{\partial x} + 12\frac{\partial h}{\partial t},
\end{eqnarray} [/math]</p>
<h1>Solution of Reynolds Equation</h1>
In general, the Reynolds equation has to be solved using numerical methods such as finite difference, or finite element. Depending on the boundary conditions and the considered geometry, however, analytical solutions can be obtained under certain assumptions.
For the case of a sphere on flat geometry (for rigid bodies) and steadystate case, the 2D Reynolds equation can be solved analytically assuming Sommerfeld (also called halfSommerfeld) cavitation boundary condition. This solution was proposed by a Nobel Prize winner <a href="https://en.wikipedia.org/wiki/Pyotr_Kapitsa">Professor Kapitza</a>. The Sommerfeld boundary condition, however is not accurate and this solution has to be used as an approximate.
In case of 1D Reynolds equation, there are several analytical, semianalytical and approximate solutions available. In 1916 Martin obtained a closed form solution for a minimum film thickness and pressure for a cylinder and plane geometry under assumptions of rigid surfaces. Derivation of the solution and the corresponding MATLAB software can be found <a href="http://www.tribonet.org/cmdownloads/analyticalsolutionof1dreynoldsequation/">here</a>. Martin employed SwiftStieber cavitation boundary conditions. This solution is not accurate in case of high loads (high pressure in the lubricant), when the elastic deformation of the surfaces contributes to the film thickness. Divergence of experimental and theoretical results by Martin for high loads leaded researchers to conclusion that elastic distortion plays a significant role in lubrication. In 1949, Grubin obtained a solution for so called elastohydrodynamic lubrication (EHL) line contact problem with certain simplifications, where he combined both elastic deformation and lubricant hydrodynamic flow. Although his solution did not satisfy both elastic and hydrodynamic equations of EHL, his analysis was recognized as particularly useful. Corresponding MATLAB code can be found <a href="http://www.tribonet.org/cmdownloads/semianalyticalsolutionof1dtransientreynoldsequationgrubinsapproximation/">here</a>.
<h1>Generalization of Reynolds Equation</h1>
Various of generalized Reynolds equations were derived to weaken the assumptions used to derive the classical form. For example, compressible, nonNewtonian lubricant behavior can be <a href="http://www.tribonet.org/wiki/generalizedreynoldsequation/" target="_blank" rel="noopener noreferrer">considered</a>. In tribology, Reynolds equation is used to predict the thickness of the lubricant film, but also to predict the friction developed by the lubricant on the surfaces. Since many tribological contacts operate in highly loaded regime and thin films, the shear rates can be very high (in the order of [math] 10^710^9 [/math]). Many of the typical lubricants start to behave nonNewtonian in the contact conditions and therefore, Reynolds equation was generalized to the case of nonNewtonian lubricant.
Another generalization includes the slip boundary conditions. This form of the Reynolds equation is used to calculate film thicknesses and friction in textured surfaces or surfaces with high slip.

Excerpt 
What are the basics of Reynolds equation? Where do these equations apply?