Generalized Reynolds Equation: Derivation

Generalized Reynolds Equation: Derivation

Reynolds equation is one of the fundamental equations of lubrication theory and tribology. In the <a href="http://www.tribonet.org/wiki/reynolds-equation/" target="_blank" rel="noopener">simplest case</a>, it considers incompressible and Newtonian lubricant. This model is sufficient to describe lubricant friction in <a href="http://link.springer.com/article/10.1007/s11249-015-0536-z">some cases</a> (relatively low loads and sliding speeds, no temperature rise). However, in many cases, the Newtonian  cannot describe lubricant friction. This is particularly true for non-Newtonian lubricants, highly loaded contacts. In this cases, the classical equation cannot be used and a more generalized version has to be developed. The reason for a generalization is in the fact that although pressure does not change across the lubricant film, the viscosity of the lubricant chagnes due to the variation of the speed of the lubricant layers, and correspondingly variation in temperature of the lubricant. On top of it, since the temperature varies, density is also affected. There are several forms of the generalization of Reynolds equation available, which consider various aspects (see for example <a href="http://www.sciencedirect.com/science/article/pii/S0020740362800381">Dowson, 1962</a>). one of the most general forms, incorporating non-Newtonian behavior, but also allows for density variation will be presented here following <span id="scm6MainContent_lblAuthors" class="authorNames"><a href="http://tribology.asmedigitalcollection.asme.org/article.aspx?articleid=1459762">Yang Peiran and Wen Shizhu</a>.</span> [caption id="attachment_1469" align="aligncenter" width="135"]<img class="size-full wp-image-1469" src="http://www.tribonet.org/wp-content/uploads/2017/01/Coordinate-System-1.png" alt="Coordinate System" width="135" height="119" data-wp-pid="1469" /> Coordinate System[/caption] We start with the equilibrium equations equations under thin film approximation: [math] \begin{eqnarray} \label{complete_sys1} \frac{\partial p}{\partial x} = \frac{\partial \tau_x}{\partial z}\\ \label{complete_sys2} \frac{\partial p}{\partial y} = \frac{\partial \tau_y}{\partial z} \\ \end{eqnarray} [/math] where [math] p [/math] is lubricant pressure, [math] \tau_x, \tau_y [/math] are the shear stresses in corresponding directions. In general, the shear stresses are linked to the shear strains and viscosity according to the following relations: [math] \begin{eqnarray} \label{complete_sys1} \tau_x = \eta_x^{*} \frac{\partial u}{\partial z}\\ \label{complete_sys2} \tau_y = \eta_y^{*} \frac{\partial v}{\partial z}\\ \end{eqnarray} [/math] where [math] \eta_x^{*},\eta_y^{*} [/math] are called equivalent viscosities. Following boundary conditions are imposed: [math] \begin{eqnarray} \label{complete_sys1} u = u_1, v=v_1, w = w_1, z = h_1, \\ \label{complete_sys2} u = u_2, v=v_2, w = w_2, z = h_2, \end{eqnarray} [/math] Substitution of the equations (2) into (1) and integrating with boundary conditions (3) gives following: [math] \begin{eqnarray} \label{complete_sys1} u(z) = u_1 + \frac{\partial p}{\partial x} [ \int_{h_1}^{z} \frac{z'}{\eta_x^{*}} dz' - h\frac{\eta_{ex}}{\eta_{ex}^{'}} \int_{h_1}^{z} \frac{1}{\eta_x^{*}} dz'] + \frac{u_2-u_1}{h}\eta_{ex}\int_{h_1}^{z} \frac{1}{\eta_x^{*}} dz', \\ \label{complete_sys2} v(z) = v_1 + \frac{\partial p}{\partial y} [ \int_{h_1}^{z} \frac{z'}{\eta_y^{*}} dz' - h\frac{\eta_{ey}}{\eta_{ey}^{'}} \int_{h_1}^{z} \frac{1}{\eta_y^{*}} dz'] + \frac{v_2-v_1}{h}\eta_{ey}\int_{h_1}^{z} \frac{1}{\eta_y^{*}} dz' \end{eqnarray} [/math] where [math] h=h_2-h_1, \frac{h}{\eta_{ex}}=\int_{h_1}^{h_2} \frac{1}{\eta_x^{*}} dz', \frac{h^2}{\eta_{ex}^{'}}=\int_{h_1}^{h_2} \frac{1}{\eta_x^{*}} dz',  \frac{h}{\eta_{ey}}=\int_{h_1}^{h_2} \frac{1}{\eta_y^{*}} dz', \frac{h^2}{\eta_{ey}^{'}}=\int_{h_1}^{h_2} \frac{1}{\eta_y^{*}} dz'. [/math] Reynolds equation is obtained by integrating the continuity equation acrros the film thickness. In general, the continuity equation is given in the following form: [math] \frac{\partial \rho}{\partial t} + \frac{\partial \rho u}{\partial x} + \frac{\partial \rho v}{\partial y} + \frac{\partial \rho w}{\partial z} =0 [/math]. It can be integrated from [math] h_1 [/math] to [math] h_2 [/math] and after manipulations gives the following generalized Reynolds equation: <p style="text-align: center;">[math]</p> <p style="text-align: center;">   \frac {\partial }{\partial x} (h^3 ({\frac {\rho}{\eta}})_{ex} \frac {\partial p}{\partial x}) +   \frac {\partial }{\partial y}  (h^3 ({\frac {\rho}{\eta}})_{ey}  \frac {\partial p}{\partial y}) =12\frac {\partial }{\partial x} (\rho_x^{*} \bar u h) + 12\frac {\partial }{\partial y} (\rho_y^{*} \bar v h) + 12\frac {\partial }{\partial t} (\rho_e h)</p> <p style="text-align: center;">[/math]</p> <p style="text-align: justify;">Here,</p> <p style="text-align: justify;">[math] \begin{eqnarray} \label{complete_sys1}</p> <p style="text-align: justify;">\bar u = \frac{u_1+u_2}{2},\bar v = \frac{v_1+v_2}{2}, \rho_e = \frac{1}{h} \int_{h1}^{h2} \rho dz, \\</p> <p style="text-align: justify;">\label{complete_sys2}</p> <p style="text-align: justify;">({\frac {\rho}{\eta}})_{ex}=12(\frac {\eta_{ex}\rho_{ex}^{'}}{\eta_{ex}^{'}})-\rho_{ex}^{''}, ({\frac {\rho}{\eta}})_{ey} =12(\frac {\eta_{ey}\rho_{ey}^{'}}{\eta_{ey}^{'}})-\rho_{ey}^{''} \\</p> <p style="text-align: justify;">\label{complete_sys3}</p> <p style="text-align: justify;">\rho_{ex}^{*} =  [\rho_{ex}^{'}\eta_{ex}(u_2-u_1)+\rho_e u_1]/\bar u, \rho_{ey}^{*} =  [\rho_{ey}^{'}\eta_{ey}(v_2-v_1)+\rho_e v_1]/\bar v\\</p> <p style="text-align: justify;">\label{complete_sys4}</p> <p style="text-align: justify;">\rho_{ex}^{'} = \frac{1}{h^2} \int_{h_1}^{h_2} \rho \int_{h_1}^{z} \frac{dz'}{\eta_x^*}dz, \rho_{ey}^{'} = \frac{1}{h^2} \int_{h_1}^{h_2} \rho \int_{h_1}^{z} \frac{dz'}{\eta_y^*}dz \\</p> <p style="text-align: justify;">\label{complete_sys5}</p> <p style="text-align: justify;">\rho_{ex}^{''} = \frac{1}{h^3} \int_{h_1}^{h_2} \rho \int_{h_1}^{z} \frac{z^'dz'}{\eta_x^*}dz,\rho_{ey}^{''} = \frac{1}{h^3} \int_{h_1}^{h_2} \rho \int_{h_1}^{z} \frac{z^'dz'}{\eta_y^*}dz</p> <p style="text-align: justify;">\end{eqnarray} [/math]</p> <p style="text-align: justify;">The generalized Reynolds equation obtained this way can be solved numerically, for example using finite elements (<a href="https://books.google.nl/books/about/A_Full_system_Finite_Element_Approach_to.html?id=ouzOSAAACAAJ&amp;redir_esc=y">a full system method of Habchi</a>) or finite difference method.</p> &nbsp; &nbsp; <span style="border-radius: 2px; text-indent: 20px; width: auto; padding: 0px 4px 0px 0px; text-align: center; font: bold 11px/20px 'Helvetica Neue',Helvetica,sans-serif; color: #ffffff; background: no-repeat scroll 3px 50% / 14px 14px #bd081c; position: absolute; opacity: 1; z-index: 8675309; display: none; cursor: pointer;">Save</span> <span style="border-radius: 2px; text-indent: 20px; width: auto; padding: 0px 4px 0px 0px; text-align: center; font: bold 11px/20px 'Helvetica Neue',Helvetica,sans-serif; color: #ffffff; background: no-repeat scroll 3px 50% / 14px 14px #bd081c; position: absolute; opacity: 1; z-index: 8675309; display: none; cursor: pointer;">Save</span> <span style="border-radius: 2px; text-indent: 20px; width: auto; padding: 0px 4px 0px 0px; text-align: center; font: bold 11px/20px 'Helvetica Neue',Helvetica,sans-serif; color: #ffffff; background: no-repeat scroll 3px 50% / 14px 14px #bd081c; position: absolute; opacity: 1; z-index: 8675309; display: none; cursor: pointer;">Save</span> <span style="border-radius: 2px; text-indent: 20px; width: auto; padding: 0px 4px 0px 0px; text-align: center; font: bold 11px/20px 'Helvetica Neue',Helvetica,sans-serif; color: #ffffff; background: no-repeat scroll 3px 50% / 14px 14px #bd081c; position: absolute; opacity: 1; z-index: 8675309; display: none; cursor: pointer;">Save</span> <span style="border-radius: 2px; text-indent: 20px; width: auto; padding: 0px 4px 0px 0px; text-align: center; font: bold 11px/20px 'Helvetica Neue',Helvetica,sans-serif; color: #ffffff; background: no-repeat scroll 3px 50% / 14px 14px #bd081c; position: absolute; opacity: 1; z-index: 8675309; display: none; cursor: pointer;">Save</span> <span style="border-radius: 2px; text-indent: 20px; width: auto; padding: 0px 4px 0px 0px; text-align: center; font: bold 11px/20px 'Helvetica Neue',Helvetica,sans-serif; color: #ffffff; background: no-repeat scroll 3px 50% / 14px 14px #bd081c; position: absolute; opacity: 1; z-index: 8675309; display: none; cursor: pointer;">Save</span> <span style="border-radius: 2px; text-indent: 20px; width: auto; padding: 0px 4px 0px 0px; text-align: center; font: bold 11px/20px 'Helvetica Neue',Helvetica,sans-serif; color: #ffffff; background: no-repeat scroll 3px 50% / 14px 14px #bd081c; position: absolute; opacity: 1; z-index: 8675309; display: none; cursor: pointer;">Save</span>