Printer Friendly

Estimation of the Shear Stress Parameter of a Power-Law Fluid.

1. Introduction

Due to its interesting mathematical properties and huge number of applications, the problem of non-Newtonian fluid flow has attracted attention from a large number of researchers. The most commonly used model for non-Newtonian fluid mechanics is the Ostwald-de Waele model with a power-law rheology. The value of the power-law index n = 1 corresponds to a Newtonian fluid, while n > 1 describes a dilatant or shear-thickening fluid and 0 < n < 1 describes a pseudo-plastic or shear-thinning fluid. The corresponding nonlinear third-order problem is the following:

([[absolute value of (f")].sup.n-1] f")' + 1/n + 1 ff" = 0, (1)

subject to

f (0) = 0, f (0) = [epsilon], f' (v) [right arrow] 1 as [eta] [right arrow] [infinity], (2)

where the primes are derivatives with respect to [eta]. To simplify notation let y = f, x = [eta], where these are referred to as the similarity variables of the problem. Equation (1) can be written, for positive curvature solutions y" > 0, as

y"' = [beta]y [(y").sup.2-n]; [beta] = -1/n(n + 1), (3)

where solution domains with y" = 0 are ruled out for n > 2. (Those domains are naturally ruled out since the problem is solved until the solution reaches its terminal state y" = 0 at, say, x = [t.sub.f].)

Earlier studies of this problem date back to the work by Blasius [1] where he considered the case where n = 1. Existence and uniqueness of similar power-law problems have also been studied by Guedda and Hammouch [2, 3], Nachman and Taliaferro [4], and Zheng et al. [5] to mention a few. In particular, in [5] the authors studied the problem with identical boundary conditions to the ones considered here but with coefficient equal to 1 on the second term involving ff" in (1) above. Their approach, however, enabled them to obtain bounds on an important parameter, the skin friction coefficient, in addition to proving existence and uniqueness. They compared their results to numerical estimates to illustrate their accuracy and it was observed that their bounds (upper and lower for skin friction) were very close for large n so that an accurate estimate can readily be obtained by taking, for example, the average of those two bounds. We will not be able to compare our results to theirs directly, however, due to the coefficient on the second term in (1) above which somewhat changes the "dynamics" of the governing equation and its solutions (even though the deviation may be small for some smaller values of n). Howell et al. [6] studied the problem in the context of momentum and heat transfer. In [7] the authors, Chen et al., considered boundary layer flow of an electrically conducting non-Newtonian fluid in the presence of a magnetic field. Their setting required two additional terms in the equation involving f' and [(f').sup.2] explicitly, but with slightly different boundary conditions. They proceeded to utilize Crocco variables and the Adomian decomposition method (ADM) to obtain approximate solutions to the problem (in the Crocco variables domain) as well as estimate the value of the skin friction coefficient and they illustrated the accuracy of their findings. Further numerical and approximate solutions can be found in [8-11], for example.

Existence and uniqueness were established for the exact current problem at hand in [12] where all possible values of the power-law index n were considered and it was assumed that f" [greater than or equal to] 0 (or y" [greater than or equal to] 0). In [13, 14] one finds extensive discussion of the general theory of non-Newtonian boundary layer flow fluids including the physical derivation of the power-law problem. Lastly we note that in [15] a unified equation governing positive and negative curvatures for non-Newtonian power-law fluids was derived and analytical solutions were obtained.

In Section 2 we derive the first few terms in a power series expansion (using the ADM) for the solution to the problem and use those terms to approximate the solution where our approach will help in obtaining estimates not only for the shear stress parameter but also for the instant of time [t.sub.f], where the solution virtually reaches its terminal point, where y'([t.sub.f]) = 1, y" ([t.sub.f]) = 0, within [10.sup.-4]. In Section 3 we compare our results to numerical estimates obtained via numerical computer integrators.

2. The Series Solution

To find the power series solution to our problem we utilize the Adomian decomposition method: observe that the solution y satisfies

y = [epsilon]x + [alpha][x.sup.2] + [L.sup.-1] (F), (4)

where the condition y"(0) = 2[alpha] replaces the condition at infinity and where F = [beta]y[(y").sup.2-n] from (3) above. [L.sup.-1] = [[integral].sub.x0] [[integral].sub.x0] [[integral].sub.x0] (*) dx dx dx is an operator, where L = [d.sup.3]/d[x.sup.3] is the operator applied to y on the left hand side of (3) as L[y]. Now expand y in an infinite series of the form y = [[summation].sup.[infinity].sub.i=0] [y.sub.i]. The decomposition method is based on substituting this series for y into (4) above and extracting the terms out one by one by a linearization process of the operator. In particular, [y.sub.i]'s are obtained via

[y.sub.i] = [L.sup.-1] ([A.sub.i-1]), i = 1, 2, 3, ..., (5)

where the [A.sub.i]'s are referred to as the Adomian polynomials and are given by

[mathematical expression not reproducible] (6)

and where we have skipped writing ([[partial derivative].sup.2]F/[[partial derivative].sup.2] f)([y.sub.0]) x ([([y.sub.1]).sup.2]/ 2) = 0 from [A.sub.2]. This in turn yields the explicit expressions (with [y.sub.0] = [epsilon]x + [alpha][x.sup.2] and using [gamma] = 2[alpha]):

[mathematical expression not reproducible]. (7)

Proceeding into finding [y.sub.1], [y.sub.2], [y.sub.3] via (5) we obtain

[mathematical expression not reproducible]. (8)

The solution must go on to [y.sub.4], [y.sub.5], ..., and so forth. We take y [approximately equal to] [y.sub.0] + [y.sub.1] + [y.sub.2] + [y.sub.3] as an approximate solution (which turns out to yield high accuracy relatively speaking, in approximating shear stress) and use it to estimate the values of the shear stress parameter [gamma] as well as the instant of time where the solution virtually reaches a constant velocity, namely, [t.sub.f], where f'([t.sub.f]) = 1, f"([t.sub.f]) = 0. (More precisely and using the current notation [t.sub.f] is the minimum value of x > 0 such that y'(x) = 1 and y"(x) = 0 and where the last two conditions hold for all x [greater than or equal to] [t.sub.f].)

3. Numerical and Series Estimates

Using the series in the previous section as an approximation to the solution of the problem one obtains the following estimates for the initial curvature at t = 0 (also referred to as the shear stress or vorticity parameter): namely, [gamma] = 2[alpha] as well as [t.sub.f]. We compare our results (referred to here as "series") with those obtained from numerical MATLAB integrators which are anticipated to have high accuracy (referred to here as "numerical"). We expect future works will obtain results using other numerical schemes and methods. Table 1 shows results for n = 1.5 whereas Tables 2 and 3 show results for n = 2.5 and n = 3, respectively. We observe that for values of [epsilon] < 1 the estimated values of the shear stress parameter [gamma] are quite accurate (especially for values of e closer to 1), as long as n stays close to n = 2, where the equation reduces to a linear equation.

The ranges we considered here for n and [epsilon] are very practical which signifies the current results. On the other hand, while the estimation of the values of [t.sub.f] may not be as significant as that for the shear stress parameter [gamma] we list two observations: The needed accuracy for [t.sub.f] is not as acute/crucial as that for [gamma] in applications so that our estimates are considered relatively accurate. Secondly the error in [t.sub.f] is one sided; it is an underestimate for n = 1.5 while it is an overestimate for n = 2.5 and n = 3. Therefore a correction can be simply done. The same observation holds true for the values of shear stress parameter [gamma] for n = 2.5 and n = 3. In fact, for n = 2.5, adding 0.004 to each estimated value of [gamma] will ensure an error of less than 0.001. On the other hand, for n = 3, adding 0.0083 to each estimated value of [gamma] will ensure an error of less than 0.0018. This shows that our approach proves highly accurate as the predictable error provides means to significantly improve the estimates without further/added work and effort. For n = 1.5 the estimated values for [gamma] are, relatively speaking, very accurate.

It is also noted that the error increases with decreasing [epsilon] so that the error can further be reduced by looking for a quadratic or higher order interaction between the error and [epsilon] to further reduce the error but we leave the details of this analysis.

Lastly observe that for a wide range of [epsilon] [less than or equal to] 0 the error can be singled out as staying virtually constant: an observation that can help in reducing the error, with virtually no effort, to practically almost 0 once it is determined for one or two measurements for each value of interest of n.

4. Conclusions

The Adomian decomposition method was used to obtain approximate series solution for a power-law problem that has recently been studied in literature. The approximate solution was used to obtain estimates for the shear stress parameter [gamma] as well as the time to reach terminal point/solution [t.sub.f]. It was noticed that only four terms of the infinite series were sufficient to obtain high accuracy for a certain very practical range of the power-law index n. A crucial observation in this context is that the error in the measured parameters is predictable which makes it possible to significantly reduce the error with virtually no effort. This is a peculiar phenomenon which we conjecture to apply to many other power-law problems of non-Newtonian fluids.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.


[1] H. Blasius, "Grenzschichten in Flussigkeiten mit kleiner Reibung," Zeitschrift fur Angewandte Mathematik und Physik, vol. 56, pp. 1-37, 1908.

[2] M. Guedda and Z. Hammouch, "Similarity flow solutions of a non-Newtonian power-law fluid," International Journal of Nonlinear Science, vol. 6, no. 3, pp. 255-264, 2008.

[3] M. Guedda, "Boundary-layer equations for a power-law sheardriven flow over a plane surface of non-Newtonian fluids," Acta Mechanica, vol. 202, no. 1-4, pp. 205-211, 2009.

[4] A. Nachman and S. Taliaferro, "Mass transfer into boundary layers for power law fluids," Proceedings of the Royal Society of London Series A: Mathematical and Physical Sciences, vol. 365, no. 1722, pp. 313-326, 1979.

[5] L. Zheng, X. Zhang, and J. He, "Existence and estimate of positive solutions to a nonlinear singular boundary value problem in the theory of dilatant non-Newtonian fluids," Mathematical and Computer Modelling, vol. 45, no. 3-4, pp. 387-393, 2007.

[6] T. G. Howell, D. R. Jeng, and K. J. De Witt, "Momentum and heat transfer on a continuous moving surface in a power law fluid," International Journal of Heat and Mass Transfer, vol. 40, no. 8, pp. 1853-1861, 1997.

[7] X.-H. Chen, L.-C. Zheng, and X.-X. Zhang, "MHD boundary layer flow of a non-newtonian fluid on a moving surface with a power-law velocity," Chinese Physics Letters, vol. 24, no. 7, pp. 1989-1991, 2007.

[8] X. Su, L. Zheng, and J. Feng, "Approximate analytical solutions and approximate value of skin friction coefficient for boundary layer of power law fluids," Applied Mathematics and Mechanics, vol. 29, no. 9, pp. 1215-1220, 2008.

[9] G. Bognar, "Similarity solution of a boundary layer flows for non-Newtonian fluids," International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 11-12, pp. 15551566, 2009.

[10] M. C. Ece and E. Buyiik, "Similarity solutions for free convection to power-law fluids from a heated vertical plate," Applied Mathematics Letters, vol. 15, no. 1, pp. 1-5, 2002.

[11] S.-J. Liao, "A challenging nonlinear problem for numerical techniques," Journal of Computational and Applied Mathematics, vol. 181, no. 2, pp. 467-472, 2005.

[12] D. M. Wei and S. Al-Ashhab, "Similarity solutions for non-Newtonian power-law fluid flow," Applied Mathematics and Mechanics, vol. 35, no. 9, pp. 1155-1166, 2014.

[13] H. Schlichting, Boundary Layer Theory, McGraw-Hill Press, New York, NY, USA, 1979.

[14] G. Bohme, Non-Newtonian Fluid Mechanics, North-Holland Series in Applied Mathematics and Mechanics, Elsevier Science, Amsterdam, The Netherlands, 1987.

[15] S. Al-Ashhab, "A curvature-unified equation for a non-Newtonian power-law fluid flow," International Journal of Advances in Applied Mathematics and Mechanics, vol. 2, no. 3, pp. 72-77, 2015.

Samer S. Al-Ashhab and Rubayyi T. Alqahtani

Department of Mathematics and Statistics, Al Imam Mohammad Ibn Saud Islamic University, P.O. Box 90950, Riyadh 11623, Saudi Arabia

Correspondence should be addressed to Samer S. Al-Ashhab;

Received 28 May 2016; Accepted 27 June 2016

Academic Editor: Mohamed Abd El Aziz
TABLE 1: Comparison between estimated values for [gamma] and
[t.sub.f] for n = 1.5.

                                -0.25      0       0.5      0.75

[gamma] (numerical)             0.3241   0.3647   0.2748   0.1720
[gamma] (series)                0.3243   0.3664   0.2754   0.1722
Error                           0.0002   0.0017   0.0006   0.0002
Improved [gamma] ([absolute     0.3234   0.3655   0.2745   0.1713
  value of (error)] < 0.0008)
[t.sub.f] (series)               5.12    4.095    3.065     2.57
Relative error in [t.sub.f]      2.7%     3.6%     5.0%     3.2%
Improved [t.sub.f]              5.235     4.21     3.18    2.685
  (error < 1.3%)

TABLE 2: Comparison between estimated values for [gamma] and
[t.sub.f] for n = 2.5.

                               -0.25       0        0.5      0.75

[gamma] (numerical)            0.3995   0.4326    0.3539    0.2534
[gamma] (series)               0.3945   0.4278    0.3497    0.2504
Error                          -0.005   -0.0048   -0.0042   -0.003
Improved [gamma] ([absolute    0.3985   0.4318    0.3537    0.2544
  value of (error)] < 0.001)
[t.sub.f] (series)              3.74     3.026     1.985    1.415
Relative error in [t.sub.f]     4.1%     4.8%      6.1%      6.7%
Improved [t.sub.f]             3.615     2.901     1.86      1.29
  (error < 1.8%)

TABLE 3: Comparison between estimated values for [gamma] and
[t.sub.f] for n = 3.

                                 -0.25       0        0.5      0.75

[gamma] (numerical)             0.4319    0.4624    0.3880    0.2898
[gamma] (series)                0.4218    0.4527    0.3795    0.2833
Error                           -0.0101   -0.0097   -0.0085   -0.0065
Improved [gamma] ([absolute     0.4301    0.4610    0.3896    0.2916
  value of (error)] < 0.0018)
[t.sub.f] (series)               3.45      2.805     1.785     1.22
Relative error in [t.sub.f]      7.0%      8.0%      9.9%      10.7%
Improved [t.sub.f]                3.4      2.755     1.735     1.17
  (error < 4.1%)
COPYRIGHT 2016 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2016 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Al-Ashhab, Samer S.; Alqahtani, Rubayyi T.
Publication:Mathematical Problems in Engineering
Date:Jan 1, 2016
Previous Article:Analytical Analysis and Field Observation of Break Line in the Main Roof over the Goaf Edge of Longwall Coal Mines.
Next Article:Topology Optimization of 3-DOF Peristaltic Structure Robot Based on Vector Continuous Mapping Matrix.

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters