`ISRN Applied MathematicsVolume 2012, Article ID 782369, 13 pageshttp://dx.doi.org/10.5402/2012/782369`
Research Article

## An Analytical Study of Weakly Nonlinear Dynamics of a Walters’ Liquid B around a Flexible Sheet Undergoing Super Linear Stretching

1Department of Mathematics, Bangalore University, Central College Campus, Bangalore 560 001, India
2Department of Civil Engineering, University of Nottingham, Malaysia Campus, Jalan Broga, 43500 Semenyih, Malaysia
3Department of Mathematics, Government First Grade College for Women, Hassan 573 201, India

Received 21 August 2011; Accepted 7 October 2011

Academic Editor: M. F. El-Sayed

Copyright © 2012 P. G. Siddheshwar et al. 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

The paper discusses the boundary layer flow of Walters’ liquid B over a stretching sheet. The stretching is assumed to be a quadratic function of the coordinate along the direction of stretching. The study encompasses within its realm both Walters’ liquid B and second order liquid. The velocity distribution is obtained by solving the nonlinear governing differential equation. Analytical expressions are obtained for stream function and velocity components as functions of the viscoelastic and stretching related parameters. It is shown that the viscoelasticity goes hand in hand with quadratic stretching in enhancing the lifting of the liquid as we go along the sheet.

#### 1. Introduction

Polymer extrusion, drawing of copper wires, continuous stretching of plastic films, and artificial fibers, hot rolling, wire drawing, glass-fiber, metal extrusion, and metal spinning are some of the examples where the problem of a stretching sheet arises. Ever since the pioneering works of Sakiadis [1, 2], several works have appeared to consider various aspects of the problem (see Siddheshwar and Mahabaleswar [3, 4], Andersson , Rollins and Vajravelu , Vleggar , Ming-I and Cha’o-Kuang , Kelly et al. , Vajravelu and Hadjinicolaou , Liao and Pop , Magyari et al. , Liao , and Dandapat and Gupta ). The core assumption in most of the reported problems is that the stretching is linearly proportional to the axial distance. This is valid provided the stretching process is delicate and slow, leading to the assumption of constant rate of stretching. It is not difficult to see that the above assumption is quite idealistic and impractical. In the strictest sense the stretching has to be nonlinearly proportional to the axial distance. In the present paper, as a first step in the general modelling exercise, we make use of a simple quadratic stretching model.

#### 2. Mathematical Formulation

We consider a steady state two-dimensional boundary layer flow of an incompressible isothermal viscoelastic liquid, of the type Walters’ liquid B, over a quadratic stretching sheet (see Figure 1). The Walters’ liquid B represents an approximation for short or rapidly fading memory liquids and is thus an approximation to first order in elasticity. The liquid is at rest and the motion is created by pulling the sheet on both ends with equal forces parallel to the sheet and with a speed , which varies quadratically with the distance from the slit as . The resulting motion of the otherwise quiescent liquid is thus caused solely by the moving sheet. On assuming, and thereby , quite small we can make use of the boundary layer theory (see Rajagopal et al. ).

Figure 1: Schematic diagram of the stretching sheet problem.

The steady two-dimensional conservation of mass and the momentum boundary layer equation for the quadratic stretching sheet problem involving Walters’ liquid B are (see Beard and Walters ): subject to the boundary conditions: Here, and are the components of the liquid velocity in the and directions, respectively, is the limiting viscosity at small rates, and is the first moment of the distribution function of relaxation times. Further we assume is quite small, that facilitates the assumption of a weakly two-dimensional flow as considered in the paper.

As pointed out by Vleggaar , in a polymer processing application involving spinning of filaments without blowing, laminar boundary layer occurs over a relatively small length of the zone 0.0–0.5 m from the die which may be taken as the origin of Figure 1. This is in fact the zone over which the major part of the stretching takes place. In such a process the initial velocity is low (about 0.3 m/s) but not very low, enough always to assume linear stretching. Thus a good approximation of the velocity of the sheet is (at any rate for the first 10–60 cm of the spinning zone), where and are the constants velocity gradients. We have adopted the quadratic stretching model in our problem. Using the dimensionless variables Equation (2.1) take the form where is the viscoelastic parameter. The parameter represents a measure of the relative importance of elastic and viscous effects and can thus be identified with the Weissenberg number.

Introducing the stream function , we get Using (2.5) in (2.6), we get

The boundary conditions to be satisfied by can be obtained from (2.2), (2.3), and (2.5) as follows: The similarity solution to (2.6), subject to (2.7), may be taken as where prime denotes differentiation with respect to . Substituting (2.8) into (2.6) and equating the coefficients of , and , we get the following three ordinary differential equations:

Equation (2.10) turns out redundant as it can be obtained by differentiating (2.9) once with respect to . In the subsequent analysis we show that (2.9) can, in fact, be obtained from (2.11), by a suitable transformation, which in turn implies consistency. The boundary conditions, for solving (2.9) for , given by (2.7) can be obtained in the formwhere . One can easily see that (2.11) is a differential equation for and we can also verify that is a solution of (2.11), and this satisfies the derivative boundary conditions in (2.12a)–(2.12c). Thus an appropriate solution of (2.9) is which satisfies the boundary condition (2.12a)–(2.12c) provided We also note that (2.13) can be a solution of the nonlinear differential equation (2.9) if and only if We may now write from (2.13)–(2.15) as

Reverting to the symbol “s” we may easily see that Using this in (2.11) we can arrive at (2.9). This proves the “consistency” of the 3 equations (2.9)–(2.11) for.

The expression for the streamline pattern of the flow in the region around the stretching sheet can be obtained from (2.8) as follows: where is a constant. The streamline can be written in the functional form as Substituting (2.8) into (2.5), we get Having obtained the analytical expression for the stream function and the velocity components and , we now move on to discuss the results obtained in the study.

#### 3. Results and Discussion

The problem of a flexible sheet undergoing quadratic stretching is investigated for the flow it generates in its immediate neighbourhood. The stretching sheet is the sole reason for the liquid flow, and liquid viscoelasticity significantly influences the flow. The flow is studied with the help of streamline patterns and also the axial and transverse velocity distributions. The results are analyzed against the background of the classical linear stretching problem involving Newtonian liquids . Before we discuss the results of the study, we make some general observations. From (2.15) it is clear that the range of applicability of the solution is . This can further be substantiated as follows. Differentiating equation (2.9) with respect to , and subject to condition (2.12a)–(2.12c), one gets From the above equation, we see that for . In conjunction with the condition in (2.12a), this would mean Obviously for , the quadratic stretching problem ceases to exist. For we note that is complex. Hence it stands reiterated that the range of applicability of must be .

We note that negative values of give us the results of a second order liquid and positive values of those of a Walters’ liquid B model. We now discuss the results of the study on Walters’ liquid B followed by those on the second order liquid.

Figure 2 is a plot of the streamline for different values of and . Increasing value of indicates the increasing rate of quadratic stretching. We find from the figure that increasing rate of stretching restricts the dynamics in the axial direction to regions close to the slit.

Figure 2: Streamline for different values of quadratic stretching parameter .

Figure 3 is a plot of various stream lines when . As is depicted in the figure, at large axial distances the streamlines converge together and are lifted up due to quadratic stretching.

Figure 3: Streamlines for different values of .

Figure 4 is a plot of the streamline for different values of and . It is evident from the aforementioned 3 figures that the viscoelastic parameter and the quadratic stretching parameter work against each other in the lifting of the liquid as we go downstream.

Figure 4: Streamline for different values of .

We now discuss the axial and transverse velocity distributions with an observation that that is, the flow outside the boundary layer becomes uniform and is directed perpendicular to the sheet. The boundary layer thickness , defined as the distance from the sheet at which the streamwise velocity has been reduced to one percent of the velocity of the quadratic stretching surface, can be expressed as

Using (2.15) in the above equation, we get Thus, we see that the boundary layer thickness is coordinate-independent for all permissible values of . Figure 5 shows the variation of with . We see from the figure that the effect of increasing is to decrease .

Figure 5: Variation of boundary layer thickness with .

Figures 615 that are three-dimensional plots of the velocity components and reveal more than the conventional two-dimensional projections on the and planes. Figure 6 is a plot of the Crane  profile of the linear stretching problem. One can easily see from the figure that the horizontal and vertical extent of the dynamics on the stretching sheet increases as we go along axial direction.

Figure 6: Axial velocity for the linear stretching sheet problem of Crane .
Figure 7: Axial velocity for the quadratic stretching sheet problem.
Figure 8: Axial velocity for the quadratic stretching sheet problem.
Figure 9: Axial velocity for the quadratic stretching sheet problem.
Figure 10: Transverse velocity for the linear stretching sheet problem of Crane .
Figure 11: Transverse velocity for the quadratic stretching sheet problem.
Figure 12: Transverse velocity for the quadratic stretching sheet problem.
Figure 13: Transverse velocity for the quadratic stretching sheet problem.
Figure 14: Axial velocity for the quadratic stretching sheet problem.
Figure 15: Transverse velocity for the quadratic stretching sheet problem.

Figure 7 brings out the effect of the quadratic stretching of the sheet as well as the viscoelasticity of the liquid. Clearly both the above effects give rise to an extended dynamic region compared to the linear stretching problem of a Newtonian liquid. Comparing Figures 7 and 8 it is obvious that the quadratic stretching increases the vertical extent of the dynamic region. Comparing Figures 8 and 9 of quadratic stretching we find that the effect of viscoelasticity is to initiate lifting of the liquid more closer to the slit compared to that of a Newtonian liquid.

The transverse velocity profile brings out the fact that quadratic stretching greatly influences the vertical velocity compared to that in the case of linear stretching. Figure 10 depicts the -independence of while Figure 11 spells out that quadratic stretching induces the -dependence of the transverse velocity component . Figures 12 and 13 explain the nature of the influence of on for the problem of quadratic stretching. The influence of on is similar to its influence on and the same is demonstrated by Figures 12 and 13. The influence of quadratic stretching on and of a second-order liquid shows that the vertical variation is comparatively less than axial variation.

Figures 14 and 15 are the axial and transverse velocity profiles for the quadratic stretching sheet problem in a second-order liquid. The corresponding graphs for a Walters’ liquid B model are Figures 7 and 11. It is clear that the lifting is initiated closer to the slit in the case of a Walters’ liquid B model compared to the second-order liquid. In the case of the latter the flow is more strongly two dimensional than in the former case. This is clearly seen on comparing Figures 11 and 15.

#### References

1. B. C. Sakiadis, “Boundary-layer behavior on continuous solid surfaces, I. The boundary layer equations for two dimensional and axi-symmetric flow,” AIChE Journal, vol. 7, pp. 26–28, 1961.
2. B. C. Sakiadis, “Boundary-layer behavior on continuous solid surfaces, II. The boundary layer on a continuous flat surface,” AIChE Journal, vol. 7, pp. 221–225, 1961.
3. P. G. Siddheshwar and U. S. Mahabaleswar, “Effects of radiation and heat source on MHD flow of a viscoelastic liquid and heat transfer over a stretching sheet,” International Journal of Non-Linear Mechanics, vol. 40, no. 6, pp. 807–820, 2005.
4. P. G. Siddheshwar and U. S. Mahabaleswar, “A new analytical procedure for solving the non-linear differential equation arising in the stretching sheet problem,” Journal of Differential Equations. resubmitted.
5. H. I. Andersson, “MHD flow of a viscoelastic fluid past a stretching surface,” Acta Mechanica, vol. 95, no. 1–4, pp. 227–230, 1992.
6. H. I. Andersson, K. H. Bech, and B. S. Dandapat, “Magnetohydrodynamic flow of a power-law fluid over a stretching sheet,” International Journal of Non-Linear Mechanics, vol. 27, no. 6, pp. 929–936, 1992.
7. H. I. Andersson, “An exact solution of the Navier-Stokes equations for magnetohydrodynamic flow,” Acta Mechanica, vol. 113, no. 1–4, pp. 241–244, 1995.
8. D. Rollins and K. Vajravelu, “Heat transfer in a second-order fluid over a continuous stretching surface,” Acta Mechanica, vol. 89, no. 1–4, pp. 167–178, 1991.
9. J. Vleggaar, “Laminar boundary-layer behaviour on continuous, accelerating surfaces,” Chemical Engineering Science, vol. 32, no. 12, pp. 1517–1525, 1977.
10. C. Ming-I and C. Cha'o-Kuang, “Temperature field in non-Newtonian flow over a stretching plate with variable heat flux,” International Journal of Heat and Mass Transfer, vol. 31, no. 5, pp. 917–921, 1988.
11. D. Kelly, K. Vajravelu, and L. Andrews, “Analysis of heat and mass transfer of a viscoelastic, electrically conducting fluid past a continuous stretching sheet,” Nonlinear Analysis: Theory, Methods & Applications, vol. 36, no. 6, pp. 767–784, 1999.
12. K. Vajravelu and A. Hadjinicolaou, “Heat transfer in a viscous fluid over a stretching sheet with viscous dissipation and internal heat generation,” International Communications in Heat and Mass Transfer, vol. 20, no. 3, pp. 417–430, 1993.
13. S. J. Liao and I. Pop, “Explicit analytic solution for similarity boundary layer equations,” International Journal of Heat and Mass Transfer, vol. 47, no. 1, pp. 75–85, 2004.
14. E. Magyari, I. Pop, and B. Keller, “New analytical solutions of a well-known boundary value problem in fluid mechanics,” Fluid Dynamics Research, vol. 33, no. 4, pp. 313–317, 2003.
15. S.-J. Liao, “On the analytic solution of magnetohydrodynamic flows of non-Newtonian fluids over a stretching sheet,” Journal of Fluid Mechanics, vol. 488, pp. 189–212, 2003.
16. B. S. Dandapat and A. S. Gupta, “Flow and heat transfer in a viscoelastic fluid over a stretching sheet,” International Journal of Non-Linear Mechanics, vol. 24, no. 3, pp. 215–219, 1989.
17. K. R. Rajagopal, T. Y. Na, and A. S. Gupta, “Flow of a viscoelastic fluid over a stretching sheet,” Rheologica Acta, vol. 23, no. 2, pp. 213–215, 1984.
18. D. W. Beard and K. Walters, “Elastico-viscous boundary-layer flows. I. Two-dimensional flow near a stagnation point,” Proceedings of the Cambridge Philosophical Society, vol. 60, pp. 667–674, 1964.
19. L. J. Crane, “Flow past a stretching plate,” Zeitschrift für angewandte Mathematik und Physik ZAMP, vol. 21, no. 4, pp. 645–647, 1970.