Abstract
The extended Stokes problems, which study the flow suddenly driven by relatively moving half-planes, are reexamined for the Oldroyd-B fluid. This topic has been studied (Liu, 2011) by applying the series expansion to calculate the inverse Laplace transform. The derived solution was correct but tough to perform the calculation due to the series expansion of infinite terms. Herein another approach, the contour integration, is applied to calculate the inversion. Moreover, the Heaviside unit step function is included into the boundary condition to ensure the consistence between boundary and initial conditions. Mathematical methods used herein can be applied to other fluids for the extended Stokes’ problems.
1. Introduction
In 1995, Zeng and Weinbaum’s paper studying the viscous flow driven by relatively moving half-planes is a new as well as pioneering work for theoretical fluid mechanics [1]. This problem is later named as the extended Stokes’ problem. Different from the well-known traditional Stokes’ problems in which the flow depends on only one spatial parameter [2, 3], the extended Stokes’ problems possess two spatial dependences. The additional spatial dependence makes the problem more complex than the traditional one. To solve the velocity, one needs to perform two different integral transforms to the momentum equation with the help of boundary and initial conditions given. For the Newtonian fluid, the analytic solutions for a finite-depth and an infinite-depth cases have been provided [4, 5].
As for non-Newtonian fluids, the extended problem for the Oldroyd-B fluid was analyzed [6] by expanding the rheological parameters as well as the spatial and temporal variables in a series form. This expansion provides a way to carry out the inverse transform to pursue the exact velocity profile. As the final solution is also expressed in an infinite series form, it is sometimes time-consuming to perform the calculation. To this end, the contour integration was adopted to perform the inverse calculation for the flow of a second-grade fluid [7]. Although the contour integration is a correct and universal tool for calculating the inverse Laplace transform, it is sometimes laborious, and even impossible, due to the complex mathematical structure of the transformed variable. Different fluid models usually have different transformed types. The variety of transformed types makes a fixed contour path impossible and thus leads to the difficulty of calculating the integration.
From the mathematical viewpoint, it seems easy to reduce the Oldroyd-B model to either the Maxwell fluid or the second-grade fluid by simply setting one parameter to be zero. However, it is not always true for the inversion calculation using the contour integration. For the second-grade fluid, the contour path shown in [7] succeeds in calculating the inversion while the same path cannot be applied to the Oldroyd-B fluid. More singularity points which should be excluded from the contour make the latter case much more complicated. This is why we study the topic again.
The organization of this paper is as follows. In Section 2, the fluid system including the momentum equation, boundary and initial conditions is elucidated. The transform integrals adopted are also introduced. Next, the detailed derivation is presented in Section 3, and concluding remarks are made in Section 4.
2. Problem Description
First, the constitutive equation, boundary and initial conditions considered are shown below (see [6] for details) where represents the velocity along the direction, the fluid constants , , and are the kinematic viscosity, the relaxation time, and the retardation times, and are the plate boundary conditions for the first and second problems, respectively. The Heaviside unit step function is defined as The inclusion of the Heaviside unit step function is to avoid the inconsistence between (2.2) and (2.6) at the point (also see [8] for reference). Due to the system being the odd function of , one only needs to calculate the positive- domain, and the boundary condition (2.3) is thus replaced by Now applying the Laplace transform and then the Fourier sine transform to the present PDE system, the transformed velocity is solved to be (see [6] for details) where for two kinds of problems, respectively. Above methods and results have been elucidated in [6]. The existing paper employed a series expansion to (2.12) to calculate the inverse transformation and then derived the solution in a series form. However, this expansion may result in a possible divergence. In the next section, another approach using the contour integration will be introduced for the sake of overcoming this weakness.
3. Contour Integration
As mentioned above, a different approach for calculating the inverse Laplace transform is adopted to derive the solution. Firstly, the term in the square root of the exponential function of (2.12) is written where It is clear that the sign of strongly dominate the pole positions inside the contour path which will be demonstrated later. With the help of (3.1), the inversion of (2.12) is shown where Two cases for calculating (3.5) are discussed below.
3.1. The Case
For the case , the term in the square root of (3.3) and (3.4) can be rewritten as which possess a nonnegative value implying that both of and are negative real numbers for all values of . Thus the relation between , , and is Now the contour integration shown in Figure 1 is adopted to calculate the inner integral in (3.5). The result is where denotes the contribution from the poles. It is clear that only one pole (note ) exists for the first problem and the residue is As for the second problem, there are two poles, and , and the result is where
The definition sketch for (3.12) is displayed in Figure 2. It is remarked that , , , and depend on . Substituting (3.9) into (3.6) with the help of (3.10) to (3.12), (3.5) can be calculated to obtain the exact solution of flow velocity .
3.2. The Case
For this case, it is more complex than the previous case. Firstly the sign of has to be determined by setting the term in the square root of (3.3) and (3.4) to be zero. It reads where According to the sign of (3.13) and (3.14), the integration with respect to in (3.5) has to be divided into three parts, as shown below For the first and last integrals, the poles and are negative real numbers due to the positive sign of (3.13). Therefore, the inner integral can be calculated using the contour path shown in Figure 3: where The residue is evaluated to be equivalent to (3.10) and (3.11) for the first and second problems, respectively.
As for the second term in (3.16), where and are complex conjugates, the contour path shown in Figure 4 is used to calculate the inner integral. It is noted that the contour path goes around the origin instead of . For the first problem, the inner integral is where the second term in (3.19) is contributed by the integration around the origin. The residue is calculated to be zero.
For the second problem, we have where the residue is contributed from four poles. It reads where where the relation is shown in Figure 5. Now the exact solution of can be acquired by substituting (3.17), and (3.19)–(3.23) into (3.16) to calculate the inverse Fourier transform.
4. Concluding Remarks
The extended Stokes’ problems for the Oldroyd-B fluid are revisited in this paper by applying the contour integration to calculate the inversion of velocity rather than using the series expansion in the earlier paper. Two cases classified according to the relation between rheological parameters are analyzed. Mathematical techniques used in this paper can be generalized and applied to investigate other fluids for the extended Stokes’ problems.
Acknowledgment
The financial support from the National Science Council of Taiwan with the Grant no. NSC 100-2221-E-270-002 is acknowledged.