Mathematical Problems in Engineering

Volume 2011 (2011), Article ID 709836, 16 pages

http://dx.doi.org/10.1155/2011/709836

## Extended Stokes' First Problem of an Oldroyd-B Fluid for Relatively Moving Half-Planes

Division of Mathematics, General Education Center, Chienkuo Technology University, 500 Changhua City, Taiwan

Received 20 October 2010; Accepted 31 January 2011

Academic Editor: Horst Ecker

Copyright © 2011 Chi-Min Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

An Oldroyd-B fluid suddenly disturbed by relatively moving half-planes is theoretically studied in this paper. This new problem, extended from the traditional Stokes' problem, recently attracts a great deal of attention due to its potential applications in engineering. Using integral transformations and dividing the original system into two subsystems, the exact solution is derived and shown in a series form. In addition to the general understanding of velocity developments, the effects of the relaxation time and the retardation time on the velocity profiles are examined in a series of figures. It is found that above two rheological parameters influence the induced flow in an opposite way. Other interesting characteristics are also elucidated from present results.

#### 1. Introduction

The most classical problem in fluid mechanics might be a viscous flow driven by an impulsive or oscillatory plate below. This problem usually bears Stokes’ name [1] due to his great contributions and profound influences on related studies. Today solutions of the Stokes’ problems have been applied to many fields, which include the industry manufacturing, chemical engineering, geophysical flows, and even the heat conduction problems. For a Newtonian fluid, exact solutions to velocity profiles can be obtained using integral transforms (e.g., Panton [2] and Erdogan [3]) and the solution to wall stress in terms of Fresnel integrals was recently given by C. M. Liu and I. C. Liu [4].

For the rapidly growing applications either in academic studies or engineering fields, the studies on the Stokes’ problems traditionally for a Newtonian fluid were extended to cases of non-Newtonian fluids. Since the rheological properties vary within a quite wide range, how to choose a suitable model to accurately describe fluid properties always becomes an important issue in related studies. To this end, the Maxwell model which is one of the simplest models for viscoelastic fluids was studied to understand basic fluid characteristics. Fetecau and Fetecau [5] solved the first Stokes’ problem of a Maxwell fluid by applying the Fourier transform. Later, Jordan et al. [6] revisited the same problem by applying both Laplace and Fourier transform methods. In addition to Maxwell fluids, another common fluid, the second-grade fluid, was frequently investigated as well. Erdogan [7] and Tan and Masuoka [8] solved the first problem of a second-grade fluid. The second problem was later studied by Fetecau and Fetecau [9]. More recently, a series of papers concerning the first problem of Burger’ fluids were also presented [10, 11].

The Oldroyd-B fluid model [12], which includes the Newtonian fluid, the Maxwell fluid and the second-grade fluid as special cases, was also widely studied in the past decades. Solution to the Stokes’ first problem of an Oldroyd-B fluid was given by Fetecau and Fetecau [13], and also by Tan and Masuoka [14] with the consideration of a porous boundary. Aksel et al. [15] later provided an analysis of the second problem and other unidirectional flows. Qi and Xu [16] considered a generalized Oldroyd-B model by applying the fractional calculus approach in the constitutive relationship. Their solution to the first problem was shown in a series form. Other specific flows were also well investigated by Rajagopal and Bhatnagar [17] and Hayat et al. [18].

In addition to studies contributing to various fluids based on various models, more boundary conditions were involved to meet the practical requirements. The most common modified boundary condition might be a porous plate which allows the mass injection and suction through the plate. The corresponding solutions have been obtained for different fluids [8, 14, 19, 20]. Pressure gradient [9] and heated boundary were also studied [8]. For practical applications in predicting flows induced by either earthquakes or fracture of ice sheets, Zeng and Weinbaum [21] analyzed the extended Stokes’ problems for a Newtonian fluid by imposing relatively moving half-planes as a boundary condition. They provided the exact solution for the first problem and the steady-state solution for the second problem which is accurate enough for long-term phenomena. For capturing the characteristics of flows at earlier stage, the transient solution to the second problem, however, is still lacking. To this end, C. M. Liu [22] derived the exact solution by using integral transforms and the concept of antisymmetry. His exact solution can calculate the flow not only at larger times, but also at smaller times. The finite-depth cases as well as the infinite-depth cases were investigated to calculate the velocity distributions.

In this paper, the flow of an Oldroyd-B fluid driven by relatively moving half-planes with constant speed is analyzed. This problem is named as the extended Stokes’ first problem. The organization of present paper is as follows. The constitutive equations for the present problem are firstly derived in Section 2. The detailed derivation and mathematical methods used are shown in Section 3. In addition to integral transforms, the main mathematical techniques include the division of the original fluid system into two subsystems and the expansion of the specific functions into series forms. Results are discussed in Section 4, and concluding remarks are made in Section 5.

#### 2. Constitutive Equations

The flow field considered is depicted as Figure 1. An Oldroyd-B fluid rather than a Newtonian fluid [22] is considered to occupy all the space above the plate. The fluid and the plate are initially at rest everywhere. For , the positive half-plate () located at is suddenly moving at a constant speed along the coordinate and the negative half-plate () remains still. This is the so-called extended Stokes’ first problem. Firstly, the constitutive equations for an incompressible Oldroyd-B fluid [23] are introduced where denotes the Cauchy stress tensor, the extra-stress tensor, the pressure, the isotropic unit tensor, where is the velocity and is the gradient operator, and is defined to be . The fluid constants , , and are the dynamic viscosity, the relaxation time, and the retardation times, respectively. The superposed dot represents the material time derivative. This model includes several fluid models as special cases. They are the Maxwell fluid () and the Newtonian fluid (). Since the fluid is assumed to be incompressible, it leads to As the flow only moves along the direction, we assume where represents a unit vector along the direction. Substituting (2.4) into the linear momentum equation and neglecting the body forces, it reads where denotes the density of the fluid. For the present flow, it is noted that all physical variables are independent of implying the first two terms of (2.5) are zero. From (2.1), relations between , , , , , and are By considering the initial condition , solutions to (2.8) to (2.10) can be solved By using (2.12), (2.7) and (2.8) can be simplified to be Now we combine (2.5), (2.6), and (2.13) together to obtain the final constitutive equation where denotes the kinematic viscosity. In comparison with the traditional Stokes’ problems of an Oldroyd-B fluid, the additional terms of appear. The detailed transform methods as well as other mathematical techniques for solving (2.14) and the associated conditions will be demonstrated in the following section.

#### 3. Transform Methods and Solutions

For the flow system shown in Figure 1, the governing equation, the associated boundary conditions, and initial conditions are shown below
where is the constant plate speed. To solve this flow system, the original problem is decomposed into two subproblems
where and , respectively, satisfy
and are respectively constrained by
and
The solution to , which is for the traditional Stokes’ first problem, is first introduced. Applying the Laplace transform
to (3.3), (3.5) and (3.6) with the help of (3.7) results in
Thus, the solution to is
In order to obtain an analytical solution and avoid the tedious calculations of residues and contour integrals, (3.15) is rewritten in a series form
where denotes the factorial and the Gamma function. By applying the inverse Laplace transform to (3.16), it yields
To rewrite (3.17) in a dimensionless form, substituting the following dimensionless parameters
into (3.17) results in
where all asterisks are dropped for convenience. As for the second subsystem of , since the positive- and negative-*z* plates move with the same speed but in the opposite directions, the induced flow which is antisymmetrical with respect to satisfies
which requires
Above results imply that the velocity in the whole domain can be calculated by only solving the flow in the positive- domain. For flows in the positive- domain, the boundary conditions of (3.9) and (3.11) can be replaced by
Since the velocity varies with the time and the space coordinates and , in addition to the Laplace transform, it is noted that an additional integral transform in space coordinate is needed to obtain the solvable boundary-value problem. Firstly, applying the Laplace transform to (3.4), (3.8), (3.10), and (3.22) with the help of (3.12) yields
Now one further applies the Fourier sine transform
to (3.23); it yields
where
By solving (3.25), is readily obtained
Above equation can be further rewritten in a series form
After applying the inverse Fourier sine transforms to (3.28), the result is
where is the modified Bessel function of second kind. After applying the calculation shown in the appendix, the dimensionless form for can be obtained by using the dimensionless parameters defined in (3.18) as well as a new parameter
where all asterisks for dimensionless variables are dropped. The final step of the derivation is to apply the Laplace inverse transform to (3.30). To this end, terms containing the transform variable in (3.30) are classified into two kinds. For terms only involving powers of , it requires the following rules of inverse Laplace transform to calculate the inversion
where represents the Dirac delta function, indicates the Bessel function of the first kind and the absolute value. As for terms containing the product of powers and the logarithm of , the inversion will be calculated in a convolution form with the help of
Finally, by combing the solutions of and , the velocity distribution either in the positive- or the negative- domain can be analytically calculated.

#### 4. Results and Discussion

In this section, the total solution consisting of and will be calculated and plotted in dimensionless forms. The focus of the analysis is how the rheological parameters influence the velocity profiles at various sections. Firstly, the velocity profiles at various sections for the conditions , , and are plotted in Figure 2. This figure gives a general understanding of flows driven by relatively moving half-planes. It indicates that for the positive -plane, the influence of the still half-plate will decay as grows. Similarly, at the far end of the negative -plane, the effect of the moving plane becomes much weaker. It is remarked that velocities at are calculated by the solution of the traditional Stokes’ problem, and the singularity exists along the line . To elucidate the velocity developments more clearly, the velocity profiles at , and −0.01 at different times ( and 1.0) for the conditions and are plotted in Figure 3. Solid, broad-solid and dashed curves are drawn for the velocity distributions at , and −0.01, respectively. At the early stage, flows above the moving plane () are quite different from those above the still plane () as it has not enough time to fully deliver the effects of the moving plane and the still plane to flows above each other. It is also found that velocity profiles either at or at will move with that at as time goes by.

In addition to a general understanding given by Figures 2 and 3, the effects of two important rheological parameters for an Oldroyd-B fluid, the relaxation time and the retardation time, on the velocity developments are investigated by plotting the velocity diagrams at level. Figure 4 shows the velocity developments for different relaxation times (, and 4) with a fixed retardation time . Solid and dashed curves denote the results at and , respectively. It is clearly seen that a larger value of will result in a faster development of the velocity profile. Moreover, flows above the negative- plane are much weaker than those above the positive- plane which can be also seen in Figure 2. On the other hand, the effects of changing the retardation time ( and 4) on the velocity field are depicted in Figure 5. It is of interest that, in comparison with the relaxation time, the retardation time plays an opposite role in the flow development. Namely, a larger value of will yield a slower growth of velocity profiles at both sides of half-planes.

#### 5. Concluding Remarks

The extended Stokes’ first problem of an Oldroyd-B fluid driven by relatively moving half-planes is theoretically studied in this paper. Using integral transformations and dividing the original problem into two subsystems, the total solution is derived and shown in a series form. Velocity profiles and their developments influenced by the relaxation time and the retardation time are examined. From the present results, the following conclusions are summarized: (1) the velocity profiles above both half-planes will move with that at the central plane (), (2) a larger value of the relaxation time will result in a faster development of the velocity profile, and (3) the effects of the magnitude of the retardation time on the velocity field are opposite to those of the relaxation time.

#### Appendix

In this appendix, the derivation of displaying (3.29) in a series form is presented. The key technique to acquire the inversion of (3.29) is to expand the Bessel function in a series form. First the corresponding identities of Bessel functions are introduced below [24] where is the digamma function defined as where is integer and is the Euler constant (). Since there exist two kinds of expansion depending on the index , Bessel functions in (3.29) will be classified and expanded in two different types. With the help of (A.1), terms of odd in (3.29) are collected and rewritten as and can be also shown as To perform the inverse transform from -domain to -domain, terms of which contain the variable have to be expanded in a series form as well where denotes the Gamma function. Applying (A.6) into (A.5) yields where where and only contain the powers of and has the products of logarithm and powers of . Now the terms of even in (3.29) are expressed as and can be further rewritten as Applying (A.6) into (A.10) leads to where only the powers of exist for the inverse calculation. Finally, (3.29) can be rewritten by using the summation of (A.6) and (A.11) and its dimensionless form is shown as (3.30).

#### Acknowledgment

The author is indebted to the financial support from National Science Council of Taiwan with the Grant no. NSC 97-2221-E-270-013-MY3.

#### References

- G. G. Stokes, “On the effect of the internal friction of fluids on the motion of pendulums,”
*Transactions of the Cambridge Philosophical Society*, vol. 9, pp. 8–106, 1851. View at Google Scholar - R. Panton, “The transient for Stokes’ oscillating plate: a solution in terms of tabulated functions,”
*Journal of Fluid Mechanics*, vol. 31, pp. 819–825, 1968. View at Publisher · View at Google Scholar - M. E. Erdogan, “Note on an unsteady flow of a viscous fluid due to an oscillating plane wall,”
*International Journal of Non-Linear Mechanics*, vol. 35, no. 1, pp. 1–6, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - C. M. Liu and I. C. Liu, “A note on the transient solution of Stokes' second problem with arbitrary initial phase,”
*Journal of Mechanics*, vol. 22, no. 4, pp. 349–354, 2006. View at Google Scholar · View at Scopus - C. Fetecau and C. Fetecau, “A new exact solution for the flow of a Maxwell fluid past an infinite plate,”
*International Journal of Non-Linear Mechanics*, vol. 38, no. 3, pp. 423–427, 2003. View at Publisher · View at Google Scholar - P. M. Jordan, A. Puri, and G. Boros, “On a new exact solution to Stokes' first problem for Maxwell fluids,”
*International Journal of Non-Linear Mechanics*, vol. 39, no. 8, pp. 1371–1377, 2004. View at Publisher · View at Google Scholar - M. E. Erdogan, “On unsteady motions of a second-order fluid over a plane wall,”
*International Journal of Non-Linear Mechanics*, vol. 38, no. 7, pp. 1045–1051, 2003. View at Publisher · View at Google Scholar · View at Scopus - W. Tan and T. Masuoka, “Starting solutions for some unsteady unidirectional flows of a second grade fluid,”
*International Journal of Engineering Science*, vol. 43, no. 10, pp. 781–789, 2005. View at Publisher · View at Google Scholar - C. Fetecau and C. Fetecau, “Stokes' first problem for a second grade fluid in a porous half-space with heated boundary,”
*International Journal of Non-Linear Mechanics*, vol. 40, no. 4, pp. 515–522, 2005. View at Publisher · View at Google Scholar · View at Scopus - D. Vieru, T. Hayat, C. Fetecau, and C. Fetecau, “On the first problem of Stokes for Burgers' fluids. II: the cases $\nu ={\lambda}^{2}/4$ and $\nu >{\lambda}^{2}/4$,”
*Applied Mathematics and Computation*, vol. 197, no. 1, pp. 76–86, 2008. View at Publisher · View at Google Scholar - C. Fetecau, D. Vieru, T. Hayat, and C. Fetecau, “On the first problem of Stokes for Burgers' fluids, I: case $\nu <{\lambda}^{2}/4$,”
*Nonlinear Analysis: Real World Applications*, vol. 10, no. 4, pp. 2183–2194, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. G. Oldroyd, “On the formulation of rheological equations of state,”
*Proceedings of the Royal Society. London A*, vol. 200, pp. 523–541, 1950. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - C. Fetecau and C. Fetecau, “The first problem of Stokes for an Oldroyd-B fluid,”
*International Journal of Non-Linear Mechanics*, vol. 38, no. 10, pp. 1539–1544, 2003. View at Publisher · View at Google Scholar · View at Scopus - W. Tan and T. Masuoka, “Stokes' first problem for an Oldroyd-B fluid in a porous half space,”
*Physics of Fluids*, vol. 17, no. 2, Article ID 023101, 7 pages, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - N. Aksel, M. Scholle, and C. Fetecau, “Starting solutions for some unsteady unidirectional flows of Oldroyd-B fluids,”
*Zeitschrift fur Angewandte Mathematik und Physik*, vol. 57, no. 5, pp. 815–831, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. Qi and M. Xu, “Stokes' first problem for a viscoelastic fluid with the generalized Oldroyd-B model,”
*Acta Mechanica Sinica*, vol. 23, no. 5, pp. 463–469, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - K. R. Rajagopal and R. K. Bhatnagar, “Exact solutions for some simple flows of an Oldroyd-B fluid,”
*Acta Mechanica*, vol. 113, no. 1–4, pp. 233–239, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - T. Hayat, A. M. Siddiqui, and S. Asghar, “Some simple flows of an Oldroyd-B fluid,”
*International Journal of Engineering Science*, vol. 39, no. 2, pp. 135–147, 2001. View at Publisher · View at Google Scholar · View at Scopus - C. M. Liu, H. H. Hwung, and C. H. Kong, “The unsteady viscous flow generated by an oscillating porous plate,”
*Journal of Mechanics*, vol. 24, no. 2, pp. 145–152, 2008. View at Google Scholar · View at Scopus - P. M. Jordan and P. Puri, “Stokes' first problem for a Rivlin-Ericksen fluid of second grade in a porous half-space,”
*International Journal of Non-Linear Mechanics*, vol. 38, no. 7, pp. 1019–1025, 2003. View at Publisher · View at Google Scholar · View at Scopus - Y. Zeng and S. Weinbaum, “Stokes problems for moving half-planes,”
*Journal of Fluid Mechanics*, vol. 287, pp. 59–74, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - C. M. Liu, “Complete solutions to extended Stokes' problems,”
*Mathematical Problems in Engineering*, vol. 2008, Article ID 754262, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. A. Barnes, J. F. Hutton, and K. Walters,
*An Introduction to Rheology*, Elsevier Science, New York, NY, USA, 1989. - I. S. Gradshteyn and I. M. Ryzhik,
*Table of Integrals, Series, and Products*, chapter, Academic Press, New York, NY, USA, 1980.