- Research
- Open access
- Published:
A nonlinear equation for fluids in multiconnected domain
Boundary Value Problems volume 2015, Article number: 198 (2015)
Abstract
In this paper we propose a mathematical model for Reynolds’ equation of a compressible fluid on a multiconnected field which simulates the function of a hybrid bearing. The boundary conditions on the inner boundaries are derived from the flow-rate continuity through the supplying orifices and are expressed by means of an integro-differential nonlinear equation. We propose a numerical method to solve this mathematical model.
1 Introduction
There are many studies dedicated to hybrid bearings. Some of such papers can be found in the complete book [1]. The aim of our present paper is to generalize the studies (done, for instance, in [2–4]) of hybrid bearings. In the previous papers, the domain in which the equation of Reynolds was solved is supposed multiconnected field with inside boundaries represented by continuous lines, so-called ‘supplying lines.’
We consider a multiconnected domain in which the inside frontiers have a rectangular form, which has permitted, through passing to the limit of the dimensions and of the number of frontiers, the simulation of a supplying line, so we could compare our results with those of above mentioned researchers. The physical situation considered in this paper has the feature that the boundary conditions on the inner boundaries are of more difficult forms. These boundary conditions become nonlinear integral-differential expressions with respect to the solution of the Reynolds equation. With the aim to compare our results with those reported in the open literature, we consider that the inner boundaries are disposed in a normal plane to the bearing surfaces, but it is worth to mention that this is not a mandatory condition for our method. An iterative procedure of linearization, as proposed in [5–7], is used in the present study. The resulting algebraic system is solved on using a ‘column-wise’ technique (see, for instance, Castelli and Pirvies [2]).
In [8] the authors use the Cesaro means for the components of energy and prove the asymptotic equipartition in mean of the kinetic and strain energies. Some results from the theory of semigroups of linear operators are used in [9] in order to prove the existence and uniqueness of a weak solution. The paper [10] is dedicated to heat conduction in deformable bodies which depends on two temperatures. The first one is the conductive temperature, the second is the thermodynamic temperature, and the difference between them is proportional to the heat supply.
Also, we must outline the minimum principle in [11] and remark that the solution proposed in [12] can be used for Maxwell’s fluids and, also, for Newtonian fluids.
2 Basic notions
Consider a circular bearing supplying with a compressible fluid through feed pocket disposed in planes normal on the bearing axis. The geometrical elements which appear in our mathematical model, are:
-
L = the length of the bearing;
-
R = the radius of the shaft;
-
C = the radial clearance;
-
ε = the eccentricity;
-
c = the length of a feed pocket;
-
d = the width of a feed pocket;
-
θ = the polar angle measured from the line of centers bearing surfaces;
-
H = the thickness of fluid film, \(H=C+\varepsilon\cos\theta\).
We use a Cartesian system of orthogonal axes such that the \(Ox_{1}\) and \(Ox_{2}\) are included in the bearing shaft surface, with the origin O in the intersection of the normal plane on the rotation axis, in an end of the bearing, with the plane determined by the axes of bearing shaft and bearing baking. Also, the \(Ox_{3}\)-axis is parallel with the rotation axis, the \(Ox_{1}\) has the sense of the rotation of bearing shaft. In this system the coordinates of the centers of the feed pockets are \((x_{1m},x_{3m})\), where m is a generic index. If the number of feed pockets is \(n_{0}\), then \(m=\overline{1,n_{0}}\). The domain represented by the surface of the bearing shaft, on which are projected the feed pockets can by assimilated with a domain of band type, having width L in the direction of \(Ox_{3}\) and being periodic in the direction of \(Ox_{1}\), with period \(2\pi R\). Based on the periodicity, we can consider the mathematical model in the domain
where
with c the length of a feed pocket and d the width of a feed pocket.
3 The mathematical model
We use the notations:
-
P = the pressure in the fluid film;
-
H = the gap between the rotating shaft and the bushing;
-
ϱ = the density of the fluid;
-
μ = the dynamic viscosity of the fluid;
-
V = the relative speed of the axis bearing;
-
\(p_{a}\) = the pressure on the ends of the bearing.
With these notations, the mathematical model becomes [11]:
-
the Reynolds equation for the steady state condition of a gas journal bearing:
$$\begin{aligned} \frac{\partial}{\partial x_{1}} \biggl(PH^{3}\frac{\partial P}{\partial x_{1}} \biggr)+ \frac{\partial}{\partial x_{3}} \biggl(PH^{3}\frac{\partial P}{\partial x_{3}} \biggr)=6\mu V \frac{\partial}{\partial x_{1}} (PH ); \end{aligned}$$(1) -
the equation of state, which is supplemented with the assumption that the gas undergoes an isothermal process:
$$\begin{aligned} \varrho=\frac{P}{{\mathcal{R}}T}; \end{aligned}$$(2) -
the conditions on boundary of the domain \({\mathcal{D}}\), where (1) is applied:
$$\begin{aligned}& P (x_{1},0 )=P (x_{1},L )=p_{a}; \end{aligned}$$(3)$$\begin{aligned}& \oint_{C_{m}}\varrho\biggl( \frac{\mathbf{V} H}{2}-\frac{H^{3}}{12\mu} \operatorname {grad}P \biggr)\mathbf{n}\, dl=q_{m}, \qquad (x_{1}, x_{3} )\in C_{m}, \end{aligned}$$(4)where the integral represents the flow-rate through the surface into the surface determined by \(C_{m}\), the boundary of the domain \(D_{m}\), and the generators of length H. \(q_{m}\) is the flow-rate through the incoming orifice (n is the outside normal of the curve \(C_{m}\)).
The location of the minimum film thickness is defined by the attitude angle Φ (see Figures 1 and 2). The film thickness function \(h=h(x_{3}, x_{1})\) depends on the position of the center of the shaft according to
where \(C=R_{b}-R_{j}\) is the radial clearance of the bearing and \(\varTheta =\theta+\varPhi \).
In the equations of the mathematical model, we use the following non-dimensional variables:
With these variables, (1) and the conditions (2) and (3) become
where
where θ and z are the coordinates in the circumferential and axial direction, respectively; c is the radial clearance, h is the thickness of the fluid film on the ridge region;
-
\(p_{s}\) = the supplying pressure;
-
\(p_{m}\) = the pressure on the curve \(C_{m}\);
-
\(h_{m}\) = the film thickness on the curve \(C_{m}\);
-
a = the radius of the feed orifice;
-
γ = the adiabatic constant and its value is about 1,405;
-
\(C_{D}\) = a coefficient, depending on the form of the feed orifice, experimentally was obtained: \(C_{D}\in[0.6, 1]\).
Because of the fact that the coefficients of (5) depend only on the parameter Λ, which is constant with regard to θ and z, and with regard to h, which is a periodical function on θ, one finds, as a result, the solution of the equation is periodical:
It is easy to see that the equations of the mathematical model are nonlinear. The linearization of these equations is made, in the usual way, by introducing the function \(Q=(ph)^{2}\).
With this new variable, (5), (6), and (9) become
In order to obtain the new form of (7), we consider the domain corresponding to a feed pocket.
Suppose that the pressure is constant on the boundary of the feed pocket,
and consider a mean value of h on \(C_{m}\), that is, \(\tilde{h}_{m}\), and a mean value of \(q_{m}\) on \(C_{m}\), that is, \(\tilde{q}_{m}\). Then we have as a result a constant value for Q:
The outline \(C_{m}\) is a rectangle with the sides \(AB\|CD\|Ox_{1}\) and \(BC\|DA\|Ox_{3}\) (see Figure 3). The equations of the sides are:
-
AB: \(z=z_{1m}\), \(\theta\in [\theta_{1m},\theta_{2m} ]\);
-
BC: \(\theta=\theta_{2m}\), \(z\in [z_{1m},z_{2m} ]\);
-
CD: \(z=z_{2m}\), \(\theta\in [\theta_{1m},\theta_{2m} ]\);
-
DA: \(\theta=\theta_{1m}\), \(z\in [z_{1m},z_{2m} ]\).
With these hypotheses, (7) becomes
We approximate the derivatives of the function Q on the outline by the values in the points F, G, H, I (the mid-points of the sides AB, BC, CD, DA):
such that, for the flow, we obtain the expression
To obtain the convergence of our procedure, the function \(G_{m}\) from (8) is prolonged on the domain \(p_{m}/p_{s}>1\).
4 The solution
In order to solve the system of equations of the mathematical model, we use the method of finite difference. The domain \({\mathcal{D}}\) is covered with a grid of parallel lines with \(Ox_{1}\) and \(Ox_{3}\). The feed pocket center is a node. We have
If the dimensions of the grid are \(N\times M\), then the coordinates of a node are
where
If we use the notations \(Q (\theta_{j},z_{i} )=Q_{j,i}\), then the derivatives can be written as follows:
In the neighborhood of the feed pocket, the formulas will be modified.
The derivatives of the functions Q, in the node \((k,l)\) in the neighborhood of the outline \(C_{m}\) can be calculated as follows:
-
if \((l-1)\cdot \varDelta z\in [z_{1m},z_{2m} ]\) and \(0< k\cdot \varDelta \theta-\theta_{m1}=\varDelta \theta_{1}<\varDelta \theta\), then
$$\begin{aligned} \begin{aligned} &\biggl( \frac{\partial Q}{\partial\theta} \biggr)_{k,l} = \frac{\varDelta \theta}{ \varDelta \theta_{1}(\varDelta \theta+ \varDelta \theta_{1})}Q_{m} - \frac{\varDelta \theta- \varDelta \theta_{1}}{\varDelta \theta \varDelta \theta _{1}}Q_{k,l} - \frac{\varDelta \theta_{1}}{\varDelta \theta_{1}(\varDelta \theta+ \varDelta \theta_{1})}Q_{k-1,l}, \\ &\biggl( \frac{\partial^{2} Q}{\partial\theta^{2}} \biggr)_{k,l} = \frac {2}{\varDelta \theta_{1}(\varDelta \theta+ \varDelta \theta_{1})}Q_{m} - \frac{2}{\varDelta \theta \varDelta \theta_{1}}Q_{k,l} - \frac{2}{\varDelta \theta _{1}(\varDelta \theta+ \varDelta \theta_{1})}Q_{k-1,l}; \end{aligned} \end{aligned}$$(16) -
if \((l-1)\cdot \varDelta z\in [z_{1m},z_{2m} ]\) and \(0<(k-1)\cdot \varDelta \theta-\theta_{m2}<\varDelta \theta\), then
$$\begin{aligned} \begin{aligned} &\biggl( \frac{\partial Q}{\partial\theta} \biggr)_{k,l} = \frac{\varDelta \theta_{1}}{\varDelta \theta_{1}(\varDelta \theta+ \varDelta \theta_{1})}Q_{k+1,l} + \frac{\varDelta \theta- \varDelta \theta_{1}}{\varDelta \theta \varDelta \theta _{1}}Q_{k,l} - \frac{\varDelta \theta}{ \varDelta \theta_{1}(\varDelta \theta+ \varDelta \theta_{1})}Q_{m}, \\ &\biggl( \frac{\partial^{2} Q}{\partial\theta^{2}} \biggr)_{k,l} = \frac {2}{\varDelta \theta(\varDelta \theta+ \varDelta \theta_{1})}Q_{k+1,l} - \frac{2}{\varDelta \theta \varDelta \theta_{1}}Q_{k,l} + \frac{2}{\varDelta \theta _{1}(\varDelta \theta+ \varDelta \theta_{1})}Q_{m}. \end{aligned} \end{aligned}$$(17)
With a similar procedure we can obtain the derivatives with regard to z:
-
if \((k-1)\cdot \varDelta \theta\in [\theta_{m1},\theta_{m2} ]\) and \(0< l\cdot \varDelta z-z_{m1}=\varDelta z_{1}<\varDelta z\), then
$$\begin{aligned} \biggl( \frac{\partial Q}{\partial z} \biggr)_{k,l} = \frac{\varDelta z}{\varDelta z_{1}(\varDelta z+\varDelta z_{1})}Q_{m} - \frac{\varDelta z - \varDelta z_{1}}{\varDelta z\varDelta z_{1}}Q_{k,l} - \frac {\varDelta z_{1}}{\varDelta z(\varDelta z + \varDelta z_{1})}Q_{k,l-1}; \end{aligned}$$(18) -
if \((k-1)\cdot \varDelta \theta\in [\theta_{m1},\theta_{m2} ]\) and \(0<(l-1)\cdot \varDelta z-z_{m2}<\varDelta z\), then
$$\begin{aligned} \begin{aligned} &\biggl(\frac{\partial Q}{\partial z} \biggr)_{k,l} = \frac{\varDelta z_{1}}{\varDelta z_{1}(\varDelta z + \varDelta z_{1})}Q_{k+1,l} + \frac{\varDelta z - \varDelta z_{1}}{\varDelta z\varDelta z_{1}}Q_{k,l} - \frac {\varDelta z}{\varDelta z_{1}(\varDelta z + \varDelta z_{1})}Q_{m}, \\ &\biggl(\frac{\partial^{2} Q}{\partial z^{2}} \biggr)_{k,l} = \frac {2}{\varDelta z(\varDelta z + \varDelta z_{1})}Q_{k+1,l} - \frac{2}{\varDelta z\varDelta z_{1}}Q_{k,l} + \frac{2}{\varDelta z_{1}(\varDelta z + \varDelta z_{1})}Q_{m}. \end{aligned} \end{aligned}$$(19)
In order to determine \(Q_{m}\), we made in (13) the approximations
where
Substituting \(h_{m}\) with \(\tilde{h}_{m}\) in the expression of \(\tilde{q}_{m}\) from (13), one finds, as a result,
As an initial iteration we use the function \(Q^{(0)}\) which satisfies the boundary conditions (11) and the periodicity conditions (12), such that \(Q^{(0)}=h^{2}\). In (21) we consider the unknown variable \(x=p_{m}/p_{s}\). The coefficients of this variable can be calculated as a function of the previous iteration, such that (21) can be written in the form
where
The solution of (22) is obtained by using the chord method, because \(x\in[0,1]\). In the case \(x>1\), by using the transformation \(x=1/x'\) and (22),
After (22) or (24) is solved, we can calculate \(Q_{m}\):
In the inside points of the domain \({\mathcal{D}}\), the algebraic system of equations, which results by using the finite difference in (10), is
where \(j=\overline{1,N}\), \(i=\overline{1,M-1}\).
For \(j=1\) and \(j=N\) in (26) is used the periodicity condition:
For \(i=1\) and \(i=M\) the conditions on the outside boundary become
In the neighborhood of \(C_{m}\) the equations are:
-
if \((i-1)\cdot \varDelta z\in [z_{m1},z_{m2} ]\) and \(0< j\cdot \varDelta \theta-\theta_{m2}<\varDelta \theta\), then
$$\begin{aligned} &\frac{1}{\varDelta z^{2}}Q_{j,i-1}^{(r)}-2 \biggl\{ \frac{1}{ \varDelta \theta \varDelta \theta_{1}}+ \frac{1}{\varDelta z^{2}}+\frac{1}{2h_{j,i}}\frac{\varDelta \theta- \varDelta \theta_{1}}{\varDelta \theta \varDelta \theta_{1}} \\ &\qquad {}\times\biggl[ \biggl(\frac{\partial h }{\partial\theta} \biggr)_{j,i}+ \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr]+ \frac{1}{h_{j,i}} \biggl(\frac{\partial^{2}h}{\partial\theta^{2}} \biggr)_{j,i} \biggr\} Q_{j,i}^{(r)} \\ &\qquad {}+\frac{1}{\varDelta z^{2}}Q_{j,i-1}^{(r)}+ \biggl\{ \frac{2}{ \varDelta \theta(\varDelta \theta+\varDelta \theta_{1})}+ \frac{1}{h_{j,i}} \\ &\qquad {} \times\frac{\varDelta \theta_{1}}{ \varDelta \theta(\varDelta \theta+\varDelta \theta_{1})} \biggl[ \biggl(\frac {\partial h }{\partial\theta} \biggr)_{j,i}+ \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr] \biggr\} Q_{j-1,i}^{(r)} \\ &\quad = \biggl\{ -\frac{2}{ \varDelta \theta(\varDelta \theta+ \varDelta \theta_{1})} + \frac{1}{h_{j,i}}\frac {\varDelta \theta}{ \varDelta \theta_{1}(\varDelta \theta+ \varDelta \theta_{1})} \\ &\qquad {} \times\biggl[ \biggl(\frac{\partial h }{\partial\theta} \biggr)_{j,i} + \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr] \biggr\} Q_{m}^{(r)}; \end{aligned}$$(29) -
if \((i-1)\cdot \varDelta z\in [z_{m1},z_{m2} ]\) and \(0<(j-1)\cdot \varDelta \theta-\theta_{m2}<\varDelta \theta\), then
$$\begin{aligned} &\frac{1}{\varDelta z^{2}}Q_{j,i-1}^{(r)}-2 \biggl\{ \frac{1}{ \varDelta \theta \varDelta \theta_{1}}+ \frac{1}{\varDelta z^{2}}+\frac{1}{2h_{j,i}}\frac{\varDelta \theta- \varDelta \theta_{1}}{\varDelta \theta \varDelta \theta_{1}} \\ &\qquad {} \times\biggl[ \biggl(\frac{\partial h }{\partial\theta} \biggr)_{j,i}+ \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr]+ \frac{1}{h_{j,i}} \biggl(\frac{\partial^{2}h}{\partial\theta^{2}} \biggr)_{j,i} \biggr\} Q_{j,i}^{(r)} \\ &\qquad {}+\frac{1}{\varDelta z^{2}}Q_{j,i-1}^{(r)}+ \biggl\{ \frac{2}{ \varDelta \theta(\varDelta \theta+\varDelta \theta_{1})}- \frac{1}{h_{j,i}} \\ &\qquad {} \times\frac{\varDelta \theta_{1}}{ \varDelta \theta(\varDelta \theta+\varDelta \theta_{1})} \biggl[ \biggl(\frac {\partial h }{\partial\theta} \biggr)_{j,i}+ \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr] \biggr\} Q_{j+1,i}^{(r)} \\ &\quad = \biggl\{ - \frac{2}{\varDelta \theta(\varDelta \theta+ \varDelta \theta_{1})} - \frac{1}{h_{j,i}}\frac{\varDelta \theta}{ \varDelta \theta_{1}(\varDelta \theta+ \varDelta \theta_{1})} \biggl[ \biggl(\frac{\partial h }{\partial\theta} \biggr)_{j,i} + \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr] \biggr\} Q_{m}^{(r)}; \end{aligned}$$(30) -
if \((j-1)\cdot \varDelta \theta\in [\theta_{m1},\theta_{m2} ]\) and \(0< i\cdot \varDelta z-z_{m2}<\varDelta z\), then
$$\begin{aligned} &\frac{2}{\varDelta z(\varDelta z+\varDelta z_{1})}Q_{j,i+1}^{(r)}- 2 \biggl[\frac {1}{\varDelta \theta^{2}}+ \frac{1}{\varDelta z\varDelta z_{1}}+\frac{1}{h_{j,i}} \biggl( \frac{\partial ^{2}h }{\partial\theta^{2}} \biggr)_{j,i} \biggr] Q_{j,i}^{(r)} \\ &\qquad {}+ \biggl\{ \frac{1}{\varDelta \theta^{2}}+\frac{1}{2\varDelta \theta h_{j,i}} \biggl[ \biggl( \frac{\partial h}{\partial\theta} \biggr)_{j,i}+ \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr] \biggr\} Q_{j-1,i}^{(r)} \\ &\qquad {}+ \biggl\{ \frac{1}{\varDelta \theta^{2}}-\frac{1}{2\varDelta \theta h_{j,i}} \biggl[ \biggl( \frac{\partial h}{\partial\theta} \biggr)_{j,i}+ \frac{\varLambda }{ (Q_{j,i}^{(r-1)} )^{1/2}} \biggr] \biggr\} Q_{j+1,i}^{(r)} \\ &\quad =-\frac{2}{\varDelta z_{1}(\varDelta z+\varDelta z_{1})}Q_{m}^{(r)}. \end{aligned}$$(31)
The system of equations (26)-(31) can be written in a matrix form:
where \([A_{j} ]\), \([B_{j} ]\), \([C_{j} ]\) are square matrices of dimension \(M\times M\) and \(\{Q_{j} \}\), \(\{Q_{j-1} \}\), \(\{Q_{j+1} \}\), \(\{R_{j} \}\) are the column vectors of the unknown variable on the column j, \(j-1\), \(j+1\), and, respectively, the column vector of terms in the right-sides. Here the subscripts represent the iteration order. In (32) and in the relations that follow, we consider the iteration of order r.
In order to solve the system of equations (32) we use the method called ‘column-wise’. This method is based on the relation which can be written for each j:
Substituting (33) in (32), one finds, as a result,
By solving with regard to \(\{Q_{j} \}\) and comparing with (33) we obtain the relations
where
The initial values for (35)-(38) are
After the calculation of the matrix \([E_{j} ]\), \([D_{j} ]\), \(\{F_{j} \}\), for \(j=\overline{1,N}\), we write the relation for \(j=N+1\) and use the periodicity condition \(\{Q_{N+1} \}= \{Q_{1} \}\). We find, as a result,
From (33) and (40), written for \(j=N,N-1,\ldots,2\), we obtain the following relation between \(\{Q_{j} \}\) and \(\{Q_{1} \}\):
where
By means of (33) and (41) we obtain the relation for the calculation of \([G_{j} ]\) and \(\{S_{j} \}\), for \(j=N,N-1,\ldots,2\):
After we compute the matrix \([G_{j} ]\) and \(\{S_{j} \}\), from (41), written for \(j=1\), one finds, as a result,
The other vectors \(\{Q_{j} \}\) can be computed now from (41) for \(j=\overline{2,N}\). The solution \(Q^{(r)}\) is compared with \(Q^{(r-1)}\). Our problem is solved if
where ε is a positive number of a suitable choice.
If the last condition (47) is not valid, then the calculation is repeated for a new iteration \(r+1\).
Figure 4 compares the attitude angle function of the bearing number for a compliant journal bearing at different eccentricity: 0.1, 0.4. Figure 5 represents the three dimensional plot of the hydrodynamic pressure distribution for the compliant journal bearing running at 30,000 rpm. The hydrodynamic pressure developed along the axial direction has a parabolic profile wherein the pressure is maximum at the mid-section or at the center of the bearing and the pressure is normalized at the ends of the bearing.
5 Conclusions
Even the boundary conditions become nonlinear integral-differential expressions with respect to the Reynolds equation, we can solve the problem by using the ‘column-wise’ technique.
References
Truesdell, C, Rajagopal, KR: An Introduction in the Mechanics of Fluids. Birchauser, Boston (2009)
Castelli, V, Pirvics, J: Review of numerical methods in gas bearings film analysis. Technical Report, Columbia Univ., Research Lab. (2005)
Andres, LS: Turbulent flow, flexure-pivot hybrid bearings. J. Tribol. 118(1), 190-200 (1996)
Helene, M, Arhir, M, Frene, J: Specific aspects of hybrid bearings at high rotational speeds. Méc. Ind. 5, 651-658 (2004)
Ausman, JS: Gas-lubrification bearings. J. Basic Eng. 86, 328-336 (1965)
Nazurov, SA, Videman, JH: A modified nonlinear Reynolds equilibrium for thin viscous flows. Asymptot. Anal. 52(1-2), 1-36 (2007)
Hellsten, A, Wallin, S: Explicit algebraic Reynolds stress and non-linear viscosity models. Int. J. Comput. Fluid Dyn. 23(4), 349-361 (2009)
Marin, M, Florea, O: On temporal behaviour of solutions in thermoelasticity of porous micropolar bodies. An. Ştiinţ. Univ. “Ovidius” Constanţa Ser. Mat. 22(1), 169-188 (2014)
Marin, M, Stan, G: Weak solutions in elasticity of dipolar bodies with stretch. Carpath. J. Math. 29(1), 33-40 (2013)
Marin, M, Agarwal, RP, Mahmoud, SR: Modeling a microstretch thermoelastic body with two temperatures. Abstr. Appl. Anal. 2013, Article ID 583464 (2013). doi:10.1155/2013/583464
Marin, M: On the minimum principle for dipolar materials with stretch. Nonlinear Anal., Real World Appl. 10, 1572-1578 (2009)
Fetecau, C, Prasad, SC, Rajagopal, KR: A note on the flow induced by a constantly accelerating plate in an Oldroyd-B fluid. Appl. Math. Model. 31, 647-654 (2007)
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Marin, M., Agarwal, R.P. & Florea, O.A. A nonlinear equation for fluids in multiconnected domain. Bound Value Probl 2015, 198 (2015). https://doi.org/10.1186/s13661-015-0461-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-015-0461-y