Mathematical Problems in Engineering

Volume 2016 (2016), Article ID 4729063, 4 pages

http://dx.doi.org/10.1155/2016/4729063

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

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

Received 28 May 2016; Accepted 27 June 2016

Academic Editor: Mohamed Abd El Aziz

Copyright © 2016 Samer S. Al-Ashhab and Rubayyi T. Alqahtani. 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

We apply the Adomian decomposition method to a power-law problem for solutions that do not change the sign of curvature. In particular we consider solutions with positive curvature. The power series obtained via the Adomian decomposition method is used to estimate the shear stress parameter as well as the instant of time where the solution reaches its terminal point of a steady state. We compare our results with estimates obtained via numerical integrators. More importantly we illustrate that the error is predictable and can be reduced without further effort or using higher order terms in the approximating series.

#### 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 corresponds to a Newtonian fluid, while describes a dilatant or shear-thickening fluid and describes a pseudo-plastic or shear-thinning fluid. The corresponding nonlinear third-order problem is the following:subject towhere the primes are derivatives with respect to . To simplify notation let , , where these are referred to as the similarity variables of the problem. Equation (1) can be written, for positive curvature solutions , aswhere solution domains with are ruled out for . (Those domains are naturally ruled out since the problem is solved until the solution reaches its terminal state at, say, .)

Earlier studies of this problem date back to the work by Blasius [1] where he considered the case where . 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 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 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 ). 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 and 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 were considered and it was assumed that (or ). 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 , where the solution virtually reaches its terminal point, where , , within . 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 satisfieswhere the condition replaces the condition at infinity and where from (3) above. is an operator, where is the operator applied to on the left hand side of (3) as . Now expand in an infinite series of the form . The decomposition method is based on substituting this series for into (4) above and extracting the terms out one by one by a linearization process of the operator. In particular, ’s are obtained viawhere the ’s are referred to as the Adomian polynomials and are given by and where we have skipped writing from . This in turn yields the explicit expressions (with and using ): Proceeding into finding via (5) we obtain

The solution must go on to . We take 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 as well as the instant of time where the solution virtually reaches a constant velocity, namely, , where , . (More precisely and using the current notation is the minimum value of such that and and where the last two conditions hold for all .)

#### 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 (also referred to as the shear stress or vorticity parameter): namely, as well as . 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 whereas Tables 2 and 3 show results for and , respectively. We observe that for values of the estimated values of the shear stress parameter are quite accurate (especially for values of closer to 1), as long as stays close to , where the equation reduces to a linear equation.