We study the lubrication process with incompressible fluid taking into account the dependence of the viscosity on the pressure. Assuming that the viscosity-pressure relation is given by the well-known Barus law, we derive an effective model using asymptotic analysis with respect to the film thickness. The key idea is to conveniently transform the governing system and then apply two-scale expansion technique.

1. Introduction

Fluid film bearings are machine elements usually studied in the broader context of tribology. Simply speaking, they consist of two surfaces in relative motion, separated by a thin fluid film, that lubricates the device and minimizes the friction and, consequently, the wear of the device. In our case, the fluid is an incompressible liquid and the two surfaces are rigid. Such elements are very important in mechanical engineering since they provide the reliability of the system and are crucial factor in limiting the dissipation of energy, that is, increasing the efficiency. If a fluid film bearing is well designed, the wear is not an issue, since two surfaces are completely separated by lubricant. It is therefore important to understand the behavior of the fluid film in the bearing. For a systematic treatment of the fundamentals of fluid film lubrication and fluid film bearings we refer the reader to [1].

The foundations of the theoretical treatment of lubrication have been laid already by Rayleigh and Stokes and in particular by the famous work of Reynolds [2]. Assuming that the viscosity is constant and using the thickness as a small parameter, a simple asymptotic approximation can be easily derived providing a well-known Reynolds equation for the pressure of the fluid. However, if the pressures generated in the lubrication flow exceed the range where the viscosity can be considered independent of the pressure, the effects of pressure-dependent viscosity should be taken into account. As reported by numerous researchers (see, e.g., [36]), for technological applications in elastohydrodynamic lubrication and in journal bearing applications, where the lubricant is forced through a very narrow gap leading to very high pressures, it is essential to include the variation of the viscosity with pressure. (When considering fluid film lubrication of highly loaded contacts, it is necessary to consider both the viscosity-pressure dependence and the elastic deformation of the surfaces. This case is essential for elastohydrodynamic lubrication. In this paper, for simplicity, we will address the lubrication regime appearing in applications that exhibit pressures high enough to effectively change the lubricant’s viscosity, but not so high to initiate significant elastic deformation in the bearing material. Such regime is known as rigid-piezoviscous regime.)

Several models have been used to describe the viscosity-pressure relation. Barus [7] proposed an exponential law for the viscosity in the form Here represents viscosity at atmospheric pressure while is the pressure-viscosity coefficient. For lubricants, the coefficient typically varies from to (see, e.g., [8]). It is important to emphasize that Barus formula (1) has been extensively used in the engineering practice, sometimes combined with temperature dependence.

The exponential law (1) rules out the possibility of deriving the analytical solutions (even for simple flows) and, thus, it has been avoided in theoretical analyses presented in the literature. The other reason lies in the fact that not much has been done in proving the well-posedness of the corresponding boundary-value problems in case of viscosity law (1). Indeed, the existence results for incompressible fluid flows with pressure-dependent viscosity have been provided only under certain technical assumptions on the viscosity which are not fulfilled by the Barus formula. We refer the reader to [911]. Recently, however, the first author of the present paper made significant progress on the subject under very general assumption on satisfied by (1) and other empiric laws. First, in [12] he proved the existence and uniqueness of the solution for the stationary Reynolds equation assuming only that the function is of class and for any value of . After that, in [13], he managed to prove the well-posedness of the stationary Stokes system (describing the real physical problem) assuming that the growth of the function and its derivative is at most exponential: It is shown that the corresponding Dirichlet problem has a solution, which is unique under some technical condition which does not rule out the Barus formula. (It is well known that in classical Stokes (or Navier-Stokes) system the pressure is determined only up to constant. To assure uniqueness we need to impose some additional condition, like prescribing the value of its integral over a specified domain. According to [13], it turns out that similar, slightly more technical condition is also needed here to fix the pressure.) For the sake of reader’s convenience, we present the key steps of the existence and uniqueness proof in the appendix of this paper.

The result from [13] enables us to investigate the asymptotic behavior of the solution of the lubrication problem with viscosity obeying (1) and that is the goal of this paper. We study the situation appearing naturally in fluid film bearings: two rigid surfaces being in relative motion are separated by a thin layer of fluid; lower surface is assumed to be perfectly smooth, while the upper one is rough with roughness described by some function . Starting from stationary Stokes system with boundary conditions applicable to lubrication flow problem, we perform a rigorous asymptotic analysis with respect to the film thickness. The key idea is to conveniently transform the original problem into the Stokes system with small nonlinear perturbation. As a result, we obtain the effective pressure in the form of the explicit formula clearly acknowledging the viscosity-pressure dependence. To our knowledge, such result cannot be found in the literature and we believe that it could be instrumental for improving the known engineering practice.

2. Description of the Problem

We consider an incompressible fluid flow in a three-dimensional thin domain defined by Here is a bounded domain and is a given smooth positive function such that while is a small parameter representing the domain thickness.

We denote by , , and the lower, the upper, and the lateral boundary of (see Figure 1); namely, In view of the application we want to model, we can assume a small Reynolds number and neglect the inertial terms in the governing equations. Thus, we assume that the flow in is governed by the following system: where . The vector field denotes the fluid velocity whereas the pressure is given by the scalar field . Throughout the text, the superscript is added in order to stress the dependence of the solution on the small parameter. We assume that the pressure-viscosity dependence is given by the Barus law: The choice of boundary conditions highly depends on the devices to be considered. Here we want to study the lubrication process where two rigid surfaces are in relative motion and are separated by a thin layer of fluid. Thus, we impose the following boundary conditions: for given outer pressure and constant velocity of relative motion of two surfaces. Note that implying . Throughout the text, will denote the standard Cartesian basis. The unit outward normal on is denoted by .

As mentioned in the Introduction section, the well-posedness of the above problem has been recently established by the first author. For reader’s convenience and in order to understand the concept of the solution and its properties, we provide the essential parts of the existence and uniqueness proof in the appendix.

3. Asymptotic Analysis

3.1. Transformation Procedure

Our main goal is to find an effective law describing the asymptotic behavior of the flow governed by (5)–(7). Following [14], first we rewrite the governing system by introducing the appropriate transformation for the pressure. Taking into account (6), the momentum equation can be written as Dividing the above equation by we obtain That motivates us to introduce a new function, denoted by , such that Obviously, there exists a continuum of such functions given by with being arbitrary. The liberty in choice of the parameter is crucial for further asymptotic analysis. Noticing that equation (9) transforms into As a result, we arrive at the following transformed system satisfied by : Note that we obtained a nonlinear Stokes-like system, with nonlinearity appearing on the right-hand side in (14) that can be made as small as we need. Indeed, by choosing small enough (such assumption is entirely consistent with the existence and uniqueness result provided in the appendix) (i.e., and large enough), we have Consequently, we will be able to control the nonlinear term in a way that it does not contribute to the macroscopic model for transformed pressure .

Remark 1. It is important to emphasize that, by choosing large and negative, the pressure becomes dominated by the large constant term. However, that is only the transformed pressure (i.e., the pressure satisfying the transformed problem), not the original one. The original pressure , satisfying the governing problem, will remain in starting, physically relevant, range, no matter what value we choose for . We simply choose the value of the parameter , construct the asymptotic approximation for (satisfying transformed Stokes-like problem (14)–(18)), and then reconstruct the approximation for the original pressure from (11). In fact, it can be easily shown by the formal computation that the effective pressure obtained in such way will not depend on at all. In view of that, we can conclude that can be seen simply as an auxiliary parameter, and, as such, it has no particular physical meaning. Choosing large and negative is, indeed, crucial for further asymptotic analysis (coming after the transformation procedure), but it produces no additional constraint to the viscosity law and does not affect the applicability of the main result.

3.2. Approximation for the Transformed Problem

There are several methods that enable us to study the processes in thin domains and to find the asymptotic behavior of the fluid flow. The most sophisticated and precise approach is based on the fine asymptotic analysis with respect to the small parameter , representing in our case the film thickness. Such method was first applied in the 50s for deriving, formally, lower-dimensional approximations mostly in theory of lubrication (see [1517] and the references therein). The method can be roughly described as follows.(1)The problem originally posed in thin domain is rewritten on domain with unit thickness by introducing new, dilated (fast) variable. To accomplish that, the differential operator has to be replaced with the new one, containing the derivatives with respect to the fast variable. Consequently, the negative powers of the domain thickness appear, singularly perturbing the operator.(2)On such, rescaled domain independent of the perturbation parameter , one seeks for the unknowns in the form of the asymptotic expansion in powers of .(3)Substituting the expansions in the governing equations and collecting equal powers of lead to the recursive sequence of linear problems yielding the asymptotic model.

We apply the above method to construct the asymptotic approximation for our transformed problem. We introduce the fast variable and in the sequel employ the following notation for the partial differential operators: for a scalar function and a vector function .

We seek for an asymptotic expansion of the unknowns and in the form The leading powers in the expansions for and are suggested by the a priori estimates which can be easily derived from the corresponding weak formulation. Now we substitute the above expansions into the transformed system (14)–(18) and collect the terms with equal powers of . Due to (19), from the momentum equation (14) we first deduce implying . In view of that, the next term from the momentum equation (14) gives Taking into account the boundary conditions for the velocity (see (16)) we deduce and Observe that On the other hand, from the divergence equation (15) we have Integrating from to with respect to yields As a consequence, we obtain the classical Reynolds equation for : The boundary condition satisfied by reads (see (18)) Thus, is given as the solution of (29) and (30), being a Dirichlet boundary value problem for linear elliptic equation of second order. The velocity is then determined straightforward from (25) and (26).

3.3. Effective Pressure

So far, we have just constructed an asymptotic approximation for the transformed problem (14)–(18). Now, we have to go back to the original problem and reconstruct the asymptotic solution. In view of (11), we deduce that the pressure distribution has the following form: where is given by (29) and (30). The asymptotic approximation for the velocity reads with  .

From the mathematical point of view, we have to ensure that the effective pressure is well defined. Namely, employing classical weak maximum principle combined with the result from [12], we are going to establish the precise upper bound for ensuring that the obtained effective pressure is well defined.

From (31) we conclude that the following condition should be fulfilled: If , we can simply apply the classical weak maximum principle on (29) and (30) (see Theorem 8.1. from [18]) providing that . Condition (33) is then obviously satisfied for any . Note that it is not realistic to assume that does not change sign. In that case we can use the special variant of the maximum principle proved recently in [12]. It yields (see Theorem  1 in [12]) where is given by In view of (34), we deduce that (33) will be satisfied if , that is, for It means that, by choosing such that (36) holds, we are positive that our effective pressure is well defined.

Remark 2. In [19] a two-dimensional plane flow with a pressure-dependent viscosity obeying Barus law is considered. An approximate system describing the effective flow is formally derived. The idea is to write the governing system in nondimensional form and to estimate the order of magnitude of the unknowns and their derivatives. Under some additional simplifying assumptions concerning this order of magnitudes, a corrected Reynolds equation is directly obtained taking into account the pressure dependence of the lubricant’s viscosity. Here we employ different approach. Starting from the original 3D problem, we first construct the asymptotic approximation of the transformed system (14)–(18). The transformation procedure is based on the introduction of new unknown (transformed pressure) satisfying (10). Motivated by the recent existence result (see the appendix), we observe that the nonlinear term appearing on the right-hand side in the transformed system does not contribute to the macroscopic model leading to a classical Reynolds equation for the transformed pressure. Original pressure, explicitly acknowledging the viscosity-pressure dependence, is then reconstructed using the inverse transformation based on (11). Both results clearly indicate that, in developing lubrication approximation associated to rigid-piezoviscous regime, the pressure dependence of viscosity should be recognized in the process, not a posteriori.

4. Concluding Remarks

In this section we give additional important remarks regarding the considered problem and its possible generalizations.

(i) Justification by Error Estimate. In the previous section, we derived an asymptotic model describing the lubrication process with viscosity obeying Barus formula (6). Though the derivation was formal, it provides a good platform for understanding the direct influence of the viscosity-pressure dependence on the effective flow of the lubricant. From the strictly mathematical point of view, formally derived model should be rigorously justified by proving some kind of convergence of the original solution (which we cannot find) towards the asymptotic one. Our aim here is to evaluate the difference between those two solutions in the appropriate norm and to obtain the satisfactory error estimates. By doing that, we validate our formally derived asymptotic model and justify its usage. It is important to notice that the domain is shrinking so that the integral norms are not appropriate for the error estimate. Indeed, the convergence of to zero in norm does not mean anything since any -bounded function converges to zero in such norm. Therefore, the appropriate convergence is that on the rescaled -independent domain. It is well known (see [2022]) that for the Reynolds pressure the error of approximation (expressed in the rescaled norm) is of order ; namely, where and constant independent of . In our case, due to the continuity of the inverse transformation, it immediately follows that for the error of approximation is of the same order. For the velocity, we have

(ii) Correctors for the Transformed Problem. We can continue with the computation and obtain the corrector for the Reynolds model (29) and (30) satisfied by the transformed pressure. More precisely, we can derive second-order model for in the form of Brinkman-type system, as recently proposed by the authors in [23]. However, being very similar to 2D Navier-Stokes system, it prevents us from applying the above maximum principles in order to deduce that the effective pressure is well defined (as we did for ). That will be the subject of our future work.

(iii) General Viscosity-Pressure Relation. Similar approach can be used to extend the analysis to a case of general viscosity-pressure dependence such that where , , . The idea is to introduce as and write the governing system in terms of the new unknown. As a consequence, we will obtain a Stokes-like system with small nonlinear perturbation which can be handled in a similar manner as presented above. For details we refer the reader to [14].


The Existence and Uniqueness Result

Let , , be a bounded smooth domain. We assume that the boundary is at least of class . We consider the following system: We assume that the function satisfies the following regularity and compatibility conditions: For the dependence of the viscosity on the pressure we assume that is given by the Barus law In the following we prove that the above problem has a solution, which is unique, under some technical conditions. In standard Stokes (or Navier-Stokes) system the pressure is obviously determined only up to a constant, so it is not unique unless we impose some additional condition, like prescribing the value of its integral over or prescribing the value of the pressure in some point of the domain . That is less obvious here, since the pressure appears in the viscosity formula. However, it turns out that similar condition is needed here to fix the pressure. Note that we prescribe velocity on the whole boundary . That is just to be able to clearly present the main ideas of the existence and uniqueness proof by avoiding unnecessary notational confusion. In the governing system (5)–(7) we prescribed pressure on the lateral part of the boundary in order to treat the situation naturally arising in the applications. Since we neglected the inertial term in the momentum equation, prescribing pressure boundary condition does not represent a serious obstacle (see, e.g., [24]), so the adaptation of the presented proof is just technical, but straightforward.

Theorem A.1 (existence result). Let satisfy (A.2). Assume in addition that (A.3) holds. Then the problem (A.1) has a solution .

Before proving the above assertion let us conveniently transform the governing system. For some we define two mappings. The first one is is of class and the function is strictly increasing (and thus injective), for any parameter . Furthermore , where and . We can now define the second one as Thus is an inverse of the function , while is treated only as a parameter. Obviously is well defined, strictly increasing, and smooth. Furthermore . Next we define the function that we need in the sequel. It is defined only on and continuous. Using the assumptions on it is an easy exercise to prove the following technical result.

Lemma A.2. Let satisfy the Barus law (A.3). Then

We now define the new unknown We can now rewrite our system in terms of that new unknown and it reads

Proof of Theorem A.1. The idea is to construct the sequence by taking ,   and for The classical result by Cattabriga [25] implies for any and any the existence of such . Furthermore, for some If, for given and , we impose then the solution is unique. Furthermore, assuming that and denoting as the imbedding constant , we obtain For any , . For any such that , we have Thus, for any and and for any , there exists such that for any It is important to notice that does not depend neither on nor on . For we choose and for (A.13) and (A.14) imply that That proves the boundness of the sequence . Then, up to a subsequence, Compact embedding and (A.17) implies the strong convergence That is enough to pass to the limit in (A.10) and to prove the existence of solution of the transformed system (A.9). Furthermore, Due to (A.12) we have in addition . To prove that , , with defined by (A.5), is the solution to the original system (A.1), we need to verify that is in the range of function , , for some . We thus need to verify that , , which is fulfilled (see Lemma A.2) for sufficiently large negative . Indeed, we know that . It is therefore sufficient to choose such that .

Theorem A.3 (uniqueness result). Assume that the viscosity function satisfies the Barus law (A.3). Let be the solution to the problem (A.1) whose existence is guaranteed by Theorem A.1. If in addition we prescribe the value of the pressure in some point , that is, for some , we impose then there are no other solutions in .

Proof. Suppose that the problem has two solutions in . We apply the transformation procedure as in the existence proof and take ,  . We denote by ,   . Then, obviously and the difference satisfies the Stokes system Thus, the standard a priori estimate (see, e.g., [25]) implies On the other hand, due to (A.6), (A.7), for given we can choose in definition of such that and thus . As was arbitrary, we conclude that and . That proves the uniqueness of the velocity since it is independent of choice of . On the other hand, and do depend on . However, their difference does not since Then, if we prescribe the value of the pressure in some point , that is, if for given we put , we have so that implying .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This research was supported by the Ministry of Science, Education and Sports, Croatia, Grant 037-0372787-2797. The authors would like to thank the referees for valuable comments and suggestions that helped to correct and improve their paper.