Abstract

An exact analysis of heat and mass transfer past an oscillating vertical plate with Newtonian heating is presented. Equations are modelled and solved for velocity, temperature, and concentration using Laplace transforms. The obtained solutions satisfy governing equations and conditions. Expressions of skin friction, Nusselt number, and Sherwood number are obtained and presented in tabular forms. The results show that increasing the Newtonian heating parameter leads to increase velocity and temperature distributions whereas skin friction decreases and rate of heat transfer increases.

1. Introduction

Generally, the problems of free convection flows are usually modeled under the assumptions of constant surface temperature, ramped wall temperature, or constant surface heat flux [17]. However, in many practical situations where the heat transfer from the surface is taken to be proportional to the local surface temperature, the above assumptions fail to work. Such types of flows are termed as conjugate convective flows, and the proportionally condition of the heat transfer to the local surface temperature is termed as Newtonian heating. This work was pioneered by Merkin [8] for the free convection boundary layer flow over a vertical flat plate immersed in a viscous fluid. However, due to numerous practical applications in many important engineering devices, several other researchers are getting interested to consider the Newtonian heating condition in their problems. Few of these applications are found in heat exchanger, heat management in electrical appliances (such as computer power supplies or substation transformer), and engine cooling (such as thin fins in car radiator). Moreover, the flow over an oscillating plate with Newtonian heating also occurs in the conjugate heat transfer around fins where the conduction within the fin and the convection in the fluid surrounding it must be simultaneously analyzed in order to obtain the vital design information and in also convection flows set up when the bounding surface absorb heat by solar radiation. The literature survey shows that much attention to the problems of free convection flow with Newtonian heating is given by numerical solvers, as we can see [916] and the references therein. However, the exact solutions of these problems are very few [1721]. Exact solutions on the other hand can provide an important check for numerical methods that are used to study such flows in more complex domains.

Furthermore, the free convection flows together with heat and mass transfer are of great importance in geophysics, aeronautics, and engineering. In several process such as drying, evaporation of water at body surface, energy transfer in a wet cooling tower, and flow in a desert cooler, heat and mass transfer occurs simultaneously. In view of such applications, several authors investigated free convection flows with simultaneous heat and mass transfer phenomenon [2225].

To the best of authors’ knowledge, so far, no study has been reported in the literature which investigates the unsteady free convection flow of an incompressible viscous fluid past an oscillating vertical plate with Newtonian heating and constant mass diffusion. The present study is an attempt in this direction to fill this space. In this study, the equations of the problem are first formulated and transformed into their dimensionless forms where the Laplace transform method is applied to find the exact solutions for velocity, temperature and concentration. Moreover, expressions for skin friction, Nusselt number, and Sherwood number are obtained and presented in tabular forms. Finally, the obtained results are plotted graphically and discussed for the pertinent flow parameters.

2. Mathematical Formulation

Consider the unsteady free convection flow of a viscous incompressible fluid past a vertical plate. The is taken along the vertical plate and the is taken normal to the plate. Initially, for time , both the plate and fluid are at stationary condition with the constant temperature and concentration . At time , the plate started an oscillatory motion in its plane with the velocity against the gravitational field, where is the amplitude of the plate oscillations. At the same time, the heat transfer from the plate to the fluid is directly proportional to the local surface temperature and the concentration level near the plate is raised from to . As the plate is considered infinite in the , therefore all physical variables are independent of and are functions of and only. The physical model and coordinate system are presented in Figure 1.

Under the above assumptions, the governing equations of the free convective flow with Boussinesq’s approximation are as follows:

The initial and boundary conditions are

The radiation heat flux under Rosseland approximation [26] is expressed by

It should be noted that by using the Rosseland approximation, we limit our analysis to optically thick fluids. It is assumed that the temperature difference within the flow are sufficiently small, and then (5) can be linearized by expanding into the Taylor series about , which after neglecting higher order terms takes the form

In view of (5) and (6), (2) reduces to

To reduce the above equations into their nondimensional forms, we introduce the following nondimensional quantities:

Substituting (8) into (1), (7), and (3), we obtain the following nondimensional Partial Differential Equations: where are the Grashof number, modified Grashof number, Prandtl number, radiation parameter, and Schmidt number, respectively. The corresponding initial and boundary conditions in nondimensional form are

Here, is the Newtonian heating parameter. We note that (12) gives when , which physically means that no heating from the plate exists [12, 16].

3. Method of Solution

The Laplace transform method solves differential equations and corresponding initial and boundary value problems. The Laplace transform has the advantage that it solves problems directly, initial value problems without determining first a general solution and nonhomogeneous differential equations without solving first the corresponding homogeneous equations. In order to obtain the exact solution of the present problem, we will use the Laplace transform technique. Applying the Laplace transform with respect to time to the system of (9), we get Here, , ,  and denote the Laplace transforms of , and , respectively.

Using the initial condition (11), we get

The corresponding transformed boundary conditions are

The solutions of (15) subject to the boundary conditions (16) are where, andis the effective Prandtl number defined by Magyari and Pantokratoras [27]. By taking the inverse Laplace transform of (17) and use formulae from Appendix, we obtain

Note that the above solution is valid only for and . Moreover, the other solutions are

Case 1. When and ,

Case 2. When and ,

Case 3. When and , Here, where is the complementary error function, and are dummy variables, and , , , , and are dummy functions.
The dimensionless expression for skin friction evaluated from (20) is given by
The dimensionless expression of Nusselt number is given by
The dimensionless expression of Sherwood number is given by

4. Graphical Results and Discussion

In order to reveal some relevant physical aspects of the obtained solutions, the numerical results for velocity, temperature, and concentration are computed and shown graphically in Figures 216, to illustrate the influence of embedded flow parameters such as time , radiation parameter , Prandtl number , Grashof number , modified Grashof number , Schmidt number , Newtonian heating parameter , and phase angle . The numerical values for skin friction, Nusselt number, and Sherwood number for these parameters are also presented in Tables 13.

The velocity profiles for different values of time are shown in Figure 2. It is observed that the velocity increases with increasing values of time . The velocity profiles in case of radiation and pure convection are shown in Figure 3. It is found from this figure that the radiation parameter has an accelerating effect on velocity. Physically, it is due to the fact that an increase in the radiation parameter for fixed values of other parameters decreases the rate of radiative heat transfer to the fluid, and consequently, the fluid velocity increases. This behavior of is quite identical with that found in Figure of Mohamed et al. [28].

The effects of Prandtl number on the velocity profiles are shown in Figure 4 for (air), (electrolytic solution), (water), and (engine oil). It is seen from this figure that an increase in the values of Prandtl number results in the decrease of velocity. In heat transfer analysis, the role of Prandtl number is to control the relative thickness of the momentum and thermal boundary layers. For small value of the heat diffuses very quickly compared to the velocity. This means that for liquid metals, the thickness of the thermal boundary layer is much bigger than the velocity boundary layer. The effects of Grashof number and modified Grashof number on velocity are shown in Figures 5 and 6. It is found that the effects of Grashof number and modified Grashof number on velocity are similar. Velocity increases with increasing values of and . Physically, it is possible because an increase in the values of Grashof number and modified Grashof number has the tendency to increase the thermal and mass buoyancy effects. This gives rise to an increase in the induced flow. Further, from these figures, it is noticed that Grashof number and modified Grashof number do not have any influence as the fluid move away from the bounding surface.

The effects of on the velocity profiles are shown in Figure 7. Four different values of Schmidt number are chosen. They physically correspond to hydrogen, water vapour, ammonia, and carbon dioxide, respectively. It is clear that the velocity decreases as the Schmidt number increases. Further, it is clear from this figure that the velocity for hydrogen is the maximum and carbon dioxide carries the minimum velocity. Figure 8 displays the effect of Newtonian heating parameter on the dimensionless velocity. It is found that as the Newtonian heating parameter increases, the density of the fluid decreases, and the momentum boundary layer thickness increases and as a result, and the velocity increases within the boundary layer.

The graphical results for the phase angle are shown in Figure 9. It is interesting to note that when the phase angle is zero which physically corresponds to no oscillation, then the fluid approaches to its maximum velocity of magnitude 1 meter per second, whereas for the phase angle , the velocity gains its minimum value of magnitude 0 meter per second. The oscillations near the plate are of great significance; however, these oscillations reduce for large values of the independent variable and approach to zero as tends to infinity. The velocity profiles are plotted in Figure 10 for different values of when (impulsive motion of the plate). It is found from this figure that the behavior of on the velocity profiles quite identical with that found in Figure 6 of Narahari and Nayan [20]. Further, all these graphical results discussed above are in good agreement with the imposed boundary conditions given by (12) and (13). Hence, this ensures the accuracy of our results.

The effects of various parameters on the temperature and concentration profiles are shown in Figures 1116. In these figures, Figure 11 exhibits the influence of dimensionless time on the temperature. It is found that the temperature profiles increase with increasing time. From Figure 12, it is noted that an increase in the radiation parameter leads to an increase in the temperature due to the fact that thermal boundary layer thickness of fluid increases. The influence of Prandtl number on temperature profiles is shown in Figure 13. It is found that the temperature decreases as the Prandtl number increases. Physically, the increase of means the decrease of thermal conductivity of fluid. From Figure 14, it is observed that an increase in the Newtonian heating parameter increases the thermal boundary layer thickness and as a result the surface temperature of the plate increases. On the other hand, it is found from Figure 15 that the influence of time on concentration profiles is similar to the velocity and temperature profiles given in Figures 2 and 11. The effects of Schmidt number on the concentration profiles are shown in Figure 16. It is seen from this figure that an increase in the value of Schmidt number makes the concentration boundary layer thin, and hence, the concentration profiles decrease.

The numerical results for skin friction, Nusselt number, and Sherwood number are shown in Tables 1, 2, and 3 for various parameters of interest. It is depicted from Table 1 that skin friction decreases with, increasing and , while it increases as and are increased. Table 2 reveals that the Nusselt number increases as and are increased and decreases when and are increased. From Table 3, it is observed that the Sherwood number increases with increasing , while reverse effect is observed for .

5. Conclusions

In this paper, exact solutions of unsteady free convection flow of an incompressible viscous fluid past an oscillating vertical plate with Newtonian heating and constant mass diffusion are obtained using Laplace transform technique. The results obtained show that the velocity and temperature are increased with increasing Newtonian heating parameter. Further, the effect of Newtonian heating parameter increases the Nusselt number but reduces the skin friction. However, the Nusselt number is decreased when the radiation parameter is increased. Also, the skin friction is decreased when the radiation parameter, phase angle and Grashof number are increased. The exact solutions obtained in this study are significant not only because they are solutions of some fundamental flows, but also they serve as accuracy standards for approximate methods, whether numerical, asymptotic, or experimental.

Appendix

See Table 4.

Nomenclature

Species concentration in the fluid
Species concentration near the plate
Species concentration in the fluid far away from the plate
Dimensionless concentration
Heat capacity at a constant pressure
Mass diffusivity
Acceleration due to gravity
Heat transfer parameter for Newtonian heating
Thermal Grashof number
Modified Grashof number
Thermal conductivity of the fluid
Prandtl number
Radiative heat flux in the
Radiation parameter
Schmidt number
Temperature of the fluid
Ambient temperature
Time
Dimensionless time
Velocity of the fluid in the
Dimensionless velocity
Coordinate axis normal to the plate
Dimensionless coordinate axis normal to the plate
Thermal conductivity
Mean absorption coefficient
Volumetric coefficient of thermal expansion
Volumetric coefficient of mass expansion
Coefficient of viscosity
Kinematic viscosity
Fluid density
Stefan-Boltzmann constant
Skin friction
Dimensionless skin friction
Dimensionless temperature
Frequency of oscillation
Phase angle
:Complementary error function.

Acknowledgments

The authors would like to acknowledge MOHE and Research Management Centre, UTM for the financial support through Vote nos. 4F109 and 04H27 for this research.