Abstract

Two unsteady motions of incompressible Maxwell fluids between infinite horizontal parallel plates embedded in a porous medium are analytically studied to get exact solutions using the finite Fourier cosine transform. The motion is induced by the lower plate that applies time-dependent shear stresses to the fluid. The solutions that have been obtained satisfy all imposed initial and boundary conditions. They can be easily reduced as limiting cases to known solutions for the incompressible Newtonian fluids. For a check of their correctness, the steady-state solutions are presented in different forms whose equivalence is graphically proved. The effects of physical parameters on the fluid motion are graphically emphasized and discussed. Required time to reach the steady-state is also determined. It is found that the steady-state is rather obtained for Newtonian fluids as compared with Maxwell fluids. Furthermore, the effect of the side walls on the fluid motion is more effective in the case of Newtonian fluids.

1. Introduction

Exact solutions for initial and boundary value problems are important for many reasons. First of all, they can characterize the behavior of some fluids or solids in various circumstances. Secondly, they can be used as tests to verify numerical schemes which are implemented to study more complex deformation or motion problems. Although the numerical integration of the differential or integral equations can be realized by computers, the accuracy of results is easier established comparing with exact solutions. During the time, many exact solutions have been established for motions of the Newtonian and non-Newtonian fluids in different circumstances (see, for instance, [1, 2] for flows between two parallel plates). Generally, they correspond to initial and boundary value problems in which the velocity is given on the boundary although in some practical situations the shear stress is specified on the solid wall.

Relatively in recent studies, Renardy [3] mentioned that boundary conditions on stresses have to be imposed to formulate a well-posed boundary value problem for the flow of a Maxwell fluid between parallel plates. He also considered steady flows of viscoelastic fluids [4] and showed that the Jeffrey model is well-posed in a bounded channel if the extra-stress components are given on the boundary. Actually, in the Newtonian mechanics, the force is the cause and kinematics is the effect and prescribing the shear stress on the boundary is the same to prescribe the force applied to the boundary in order to move it [5]. The first exact solutions for motions of non-Newtonian fluids in which the shear stress is given on the boundary have been obtained by Ting [6] and Waters and King [7], respectively, for flows of second grade or Oldroyd-B fluids over an infinite plate. First exact solutions for unsteady motions of non-Newtonian fluids in which the shear stress is given on a portion of the boundary and the velocity on the other part of the boundary seem to be those of Bandelli and Rajagopal [8] for flows of the second grade fluids in cylindrical domains.

In the meantime, different exact solutions for such motions of the non-Newtonian fluids were determined both in rectangular and cylindrical domains (see, for instance, [915] and the references therein). General solutions for the motion of Newtonian fluids between parallel plates in which the shear stress is given on a part of the boundary were recently determined by Javaid et al. [16]. Other exact solutions for such motions of the incompressible Newtonian fluids with pressure-dependent viscosity were recently provided by Danish et al. [17]. Exact expressions for the steady-state solutions corresponding to oscillatory motions of the Maxwell fluids in long tubes of rectangular or triangular cross section were recently obtained by Wang et al. [18] and Sun et al. [19], respectively. Starting solutions for nonsteady flows of the Maxwell fluids induced by a moving wall through a porous plate channel were recently provided by Fetecau et al. [20]. An interesting study concerning the Maxwell fluid flow through porous plate was developed by Shahid [21]. To the best of our knowledge, closed-form expressions for starting solutions corresponding to mixed initial-boundary value problems of the Maxwell fluids moving in a porous medium lack in the literature.

The purpose of this note is to provide analytical expressions for two unsteady motions of incompressible Maxwell fluids through a horizontal porous plate channel. The fluid motion is generated by the lower plate that applies time-dependent shear stresses to the fluid. Obtained solutions fulfill all imposed initial and boundary conditions and are easily reduced to known solutions for the Newtonian fluids. For validation, the steady-state solutions are presented in different forms and their equivalence is graphically proved. Effects of pertinent parameters on the fluid movement are graphically underlined and interpreted. The required time to reach the steady-state is graphically determined. It is found that the steady-state is later obtained for Maxwell fluids as compared with Newtonian fluids. Furthermore, the effect of the side walls on the fluid motion is more effective on Newtonian fluids than the Maxwell fluids.

2. Statement of the Problem

Let us consider an incompressible upper-convected Maxwell (UCM) fluid between two infinite horizontal parallel plates embedded in a porous medium as illustrated in Figure 1.

The constitutive equations of such a fluid are as follows [22]:where T is the Cauchy stress tensor, S is the extra-stress tensor, is the constitutively indeterminate part of the stress due to the constraint of incompressibility, is the relaxation time, is the fluid viscosity, and L is the gradient of the velocity field . In the following, we establish exact solutions for two mixed initial-boundary value problems corresponding to unsteady motions of incompressible UCM fluids into abovementioned domain. All results will be reported to a convenient Cartesian coordinate system having the y-axis perpendicular to plates.

Since the plates are infinite in extent, all physical entities which characterize the fluid motion in such a domain are functions of y and t only [23]. Consequently, for such motions of the incompressible UCM fluids, we can look for a velocity field of the form as follows:where is the unit vector along the x-direction. We also suppose that the extra-stress tensor S, as well as the fluid velocity depends on y and t only.

Replacing the above velocity field in the constitutive equation (1) and bearing in mind the fact that the fluid was at rest up to the moment , i.e.,it can be proved that the components , and of S are zero and [20]where and are the non-null components of S.

Let us now assume that at the moment , the lower plate begins to apply to the fluid a shear stress having one of the forms as follows:where and are dimensional constants. In the case , corresponding to Newtonian fluids, the above relations take the simple forms as follows:

Boundary conditions of same forms (6), but for the velocity components, have been recently used by Wang and Liu (see Section 3 in [24]) for the helical flow of generalized Maxwell fluids.

Due to the shear, the fluid begins to move and, in the absence of a pressure gradient in the x-direction, the balance of linear momentum reduces to the partial differential equation as follows [25]:where is Darcy’s resistance which for such fluids has to satisfy the following relation:

In the above relations, is the fluid density while and k are the porosity and the permeability, respectively, of the porous medium. The continuity equation is identically satisfied.

Eliminating between equations (4) and (7) and using equation (8), one attains the following dimensional governing equation for the velocity field:in which is the kinematic viscosity and d is the distance between plates.

The appropriate initial and boundary conditions arefor the first motion, and the initial conditions (10) together with the boundary conditions are as follows:in the case of the second motion. It is worth to point out the fact that solving the ordinary differential equations (11) and (12) with respect to when the initial condition resulting from equation (3) is satisfied, and equality (5) is recovered.

As soon as the velocity field is known, the shear and normal stresses and can be determined solving the governing equation (4) with the initial conditions as follows:

However, in order to avoid repetition, we shall here determine the shear stress only. It is necessary, for instance, to provide the frictional force per unit area exerted by the fluid on the stationary plate.

3. Solution

In the following, the finite Fourier cosine transform, defined by equation (A.1) from Appendix, is used to determine analytical expressions for the dimensionless velocity and shear stress fields corresponding to the two unsteady motions of the incompressible UCM fluids through a porous plate channel.

3.1. Case

In order to determine solutions which are independent of the flow geometry, let us introduce the following nondimensional variables and functions:

Introducing the above nondimensional entities in equations (9)–(11) and giving up the star notation, the next mixed initial-boundary value problemis obtained. In equation (15), the Weissenberg number We and the dimensionless porosity parameter K are defined by the equalities as follows:where is a characteristic velocity. The Weissenberg number We represents the ratio of the relaxation time and the characteristic time .

The dimensionless forms of the ordinary differential equation (4) arewhere . The corresponding initial conditions are

Multiplying equation (15) with , where , integrating the result from zero to one, and bearing in mind the initial and boundary conditions (16) and (17) and equalities (A.1) and (A.2) from Appendix, it results that the finite Fourier cosine transform of has to satisfy the ordinary differential equation as follows:with the initial conditions

The solution of equation (21) with the initial conditions (22) is given by the relation as follows:where and . On the basis of equality (A.1),it results that the dimensionless velocity field corresponding to this motion is given byor equivalently (see the equality (A.3) from Appendix)

Direct computations clearly show that given by equality (25) satisfies both the governing equation (15) and all imposed initial and boundary conditions. The dimensionless volume flux across a plane normal to the x-axis and per unit width of this plane is given by the relation as follows:

For , as expected, Q becomes zero.

In order to determine the dimensionless frictional force per unit area exerted by the fluid on the stationary plate, the corresponding shear stress has to be known. Replacing the expression of from equality (25) in (19) and solving the obtained equation with the initial condition (20), it results that is given by the following expression:

The motion whose solutions have been determined in this subsection is unsteady and remains unsteady. It is clear that at large enough values of the time t, the fluid motion can be described by the long time solutions:

As expected, the dimensionless velocity and shear stress fields and corresponding to incompressible Newtonian fluids performing the same motion, namely, are immediately obtained making in equations (25) and (27). Of course, the last solutions are in agreement with the solutions obtained in equalities (29) and (31) in [16] where the boundary condition (8) has been taken with a changed sign. The corresponding volume flux can be immediately obtained making in the general expression of Q. In the absence of the porous effects, the corresponding solutions can be determined making into the above relations.

From equalities (28) and (30), it clearly results that the long time velocity fields and corresponding to the incompressible UCM and Newtonian fluids, respectively, are identical. Consequently, for large enough values of the time t, the dimensionless velocity field corresponding to incompressible UCM fluids can be approximated with that of Newtonian fluids. Of course, bearing in mind known results from the existing literature, this is not a surprise. However, this property is not valid for the corresponding shear stress fields. The corresponding frictional forces per unit area exerted by the fluid on the upper plate can be immediately determined making in any one of the relations (27), (29), or (31).

3.2. Case

Introducing the next nondimensional variables, functions, and parameters such asin equality (9) and again dropping out the star notation, one obtains for the velocity field corresponding to this motion a dimensionless governing equation of the same form as (15) with the initial condition (16). The corresponding boundary conditions are as follows:

Dimensionless governing equations for and and the corresponding initial conditions are also given by equations (19) and (20) where the new constant .

The solution of this new problem with initial and mixed boundary conditions can be determined following the same way as before. To avoid repetition, we here present the final results only. As usually, the dimensionless starting solutions and of such a motion that become steady in time can be presented under the forms as follows:in which their steady-state (permanent or long time) and transient components are given by the next relations as follows:where .

Direct computations show that and given by equation (34) satisfy both the corresponding governing equations and all imposed initial and boundary conditions. Sometime after the motion initiation, the fluid moves according to these solutions. After this time, when the transients disappear or can be neglected, the fluid motion is characterized by the steady-state solutions given by equations (35) and (36). It is worth pointing out that the steady-state components of and can be also presented in the simpler forms as follows:where , i is the imaginary unit and Im denotes the imaginary part of that which follows. The equivalence of the steady-state solutions given by the equalities (35) and (36) with those from the relations (39) is graphically proved in Figures 2 and 3.

The dimensionless steady-state and transient solutions corresponding to the motion of the incompressible Newtonian fluids induced by the lower plate that applies a oscillatory shear stress of form (6) to the fluid, namely,are immediately obtained making in equations (35)–(38). As expected, the expressions of and given by equations (40) and (42), respectively, are also in accordance with the results obtained by Javaid et al. (see equalities (40) and (41) in [16]) where the corresponding boundary condition (8) has been taken with an opposite sign and a different normalization was considered.

4. Some Numerical Results and Discussion

Exact solutions for two dimensionless mixed initial-boundary value problems corresponding to unsteady motions of incompressible UCM fluids through a porous plate channel have been established in a simple way using the finite Fourier cosine transform only. Their limiting expressions, when the relaxation time and therefore the Weissenberg number We tends to zero, correspond to motions of the incompressible Newtonian fluids induced by the lower plate that applies ramp type or oscillating shear stresses to the fluid. The steady-state solutions and corresponding to the second motion which becomes steady in time are presented in different forms whose equivalence is graphically proved in Figures 2 and 3. This proof, as well as the fact that known solutions from the existing literature for the incompressible Newtonian fluids are recovered as limiting cases of the present results, could be an evidence of the correctness of present results. The presentation of the starting solutions of a given motion as a sum of its steady-state and transient components is important for the experimentalists who want to eliminate the transients from their experiments.

To get some physical insight of present results, the influence of the Weissenberg number We and the porosity parameter K on the fluid motion has been brought to light in Figures 4 and 5. The required time to reach the steady-state for the motion which becomes steady in time is graphically determined in Figures 6 and 7 for two distinct values of We and K, respectively. This is the time after which the diagrams of the starting solution superpose over those of its steady-state component . In practice, as we previously mentioned, this time is important for those who want to eliminate the transients from their experiments. The time variations of the midplane steady-state velocity and shear stress fields are presented in Figures 8 and 9 for and 1 and different values of We.

In Figure 4, for comparison, the diagrams of the velocity fields and given by equalities (25) and (30), respectively, are together presented at two distinct values of the time t for decreasing values of We. The profiles of the starting solution are also presented in Figure 5 for three decreasing values of K. In all cases, the fluid velocity in absolute value smoothly decreases from maximum values on the lower plate to the zero value on the stationary plate. As expected, in absolute value, it is an increasing function with respect to the time t and diminishes for increasing values of K. Consequently, as expected, in the presence of a porous medium, the fluid motion becomes more difficult. Moreover, as it results from Figure 4, the fluid velocity in absolute value is a decreasing function with regard to the Weissenberg number We and converges to the Newtonian solution if . This behavior of the motion of the Maxwell fluid can be explained by the fact that, as it was proved by Poole [26], the Weissenberg number represents the ratio of elastic forces to viscous forces. Consequently, at the same elastic properties of the fluid, an increase in We means a diminution of viscous forces and therefore an increase in the fluid velocity.

Time convergence of the starting solution to its steady-state component is shown in Figures 6 and 7 for distinct values of We and K. From these figures, it also results that the required time to reach the steady-state for this motion that becomes permanent in time is an increasing function with regard to the Weissenberg number We and diminishes for increasing values of the porosity parameter K. This means that the steady-state is rather obtained for the oscillatory motion of the incompressible Newtonian fluids induced by sine oscillations of the shear stress on the boundary as compared with incompressible UCM fluids whose motion is generated by a shear stress of form (5) on the boundary. In the same time, it is clearly seen from Figure 7 that the presence of the porous medium hastens the appearance of the steady-state. The boundary condition on the upper plate is clearly satisfied.

The time variations of the midplane dimensionless steady-state velocity and shear stress fields , and , corresponding to incompressible UCM and Newtonian fluids, respectively, are together presented in Figures 8 and 9 for , , and three distinct values of the Weissenberg number We. In all cases, the convergence of the non-Newtonian solutions and to the corresponding Newtonian solutions if , as well as the oscillating specific features of the fluid motion, is clearly visualized. The amplitude of the oscillations, as expected, diminishes for increasing values of K both for velocity and shear stress. In addition, the values of the velocity and those of the corresponding shear stress have opposite signs at each time t.

5. Conclusions

Two mixed initial-boundary value problems characterizing unsteady motions of the UCM fluids between infinite parallel plates are analytically studied. Closed-form expressions are established for the dimensionless velocity and shear stress fields using the finite Fourier cosine transform. The influence of essential physical parameters on the fluid behavior was graphically investigated, and the obtained results can be such summarized:(i)Analytical solutions for two mixed initial-boundary value problems corresponding to nonsteady motions of UCM fluids through a porous plate channel are determined.(ii)The necessary time to touch the steady-state is graphically provided. It is an increasing function with respect to We and diminishes for increasing values of the parameter K.(iii)Presence of porous medium hastens the appearance of the steady-state of the motion.(iv)The amplitude of the oscillations corresponding to the steady-state solutions wanes in the presence of a porous medium both for the fluid velocity and the adequate shear stress.(v)Solutions corresponding to movements of Newtonian fluids induced by ramp type or oscillatory shear stresses on the boundary are acquired as limiting cases of the present results.

Appendix

Nomenclature

:Velocity field
λ:Relaxation time
u:Fluid velocity in x-direction
ρ:Fluid density
L:Velocity gradient
μ:Fluid viscosity
T:Cauchy stress tensor
ν = μ/ρ:Kinematic viscosity
S:Extra-stress tensor
τ = Sxy:Shear stress
k:Permeability of porous medium
σ = Sxx:Normal stress
p:Pressure
:Porosity of porous medium
R:Darcy’s resistance
We:Weissenberg number
ex:Unit vector along the x-axis
uFn(t):Finite Fourier cosine transform of ;

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.