Abstract
The thermodynamic second law analysis is utilized to investigate the inherent irreversibility in an unsteady hydromagnetic generalized Couette flow with variable electrical conductivity in the presence of induced electric field. Based on some simplified assumption, the model nonlinear governing equations are obtained and solved numerically using semidiscretization finite difference techniques. Effects of various thermophysical parameters on the fluid velocity, temperature, current density, skin friction, the Nusselt number, entropy generation number, and the Bejan number are presented graphically and discussed quantitatively.
1. Introduction
Investigation of the flow of electrically conducting fluids in porous geometries is of particular importance due to the widespread prevalence in a number of industrial applications [1–4]. These include applications to geothermal reservoirs, nuclear reactor cooling, Magnetohydrodynamic (MHD) marine propulsion, electronic packaging, microelectronic device operations, thermal insulation, and petroleum reservoirs. Preliminary experimental evidence suggests promising future applications in the field of metallurgy, in particular in the MHD stirring of molten metals and in magnetic-levitation casting. Experimental and theoretical investigations of MHD flow in various geometries and under various conditions remain a topic of ongoing interest [5]. Similar investigations not including the effects of electrical conductivity are summarized in [6–12].
The cornerstone in the field of heat transfer and thermal design is the second law analysis and its design-related concept of entropy generation minimization. The foundation of knowledge of entropy production goes back to Clausius and Kelvin's studies on the irreversible aspects of the second law of thermodynamics. Since then, the theories based on these foundations have rapidly developed. However, the entropy production resulting from combined effects of velocity and temperature gradients has remained untreated by classical thermodynamics, which has motivated many researchers to conduct analyses of fundamental and applied engineering problems based on the second law analyses. Entropy generation is associated with thermodynamic irreversibility, which is common in all types of heat transfer processes, namely, conduction, convection, and radiation. In thermodynamical analysis, the fundamental principle remains the improvement of the relevant thermal systems to mitigate against energy losses and hence fully optimize the energy resources.
The purpose of this paper is to analyze the second law of thermodynamics with respect to inherent irreversibility in an unsteady flow of a viscous, incompressible, and electrically conducting fluid. The fluid flows through a channel with isothermal walls and a transversely imposed magnetic field. Expressions for the dimensionless velocity, temperature, current density, wall shear stress, all heat transfer, Bejan number, and entropy generation number are presented. These flow quantities are discussed qualitatively with respect to the embedded parameters. The mathematical formulation of the problem is established in Section 2. In Section 3, the semi-implicit finite difference technique is implemented for the solution process of the coupled nonlinear problem. Graphical results are presented and discussed qualitatively and quantitatively with respect to various parameters embedded in the system in Section 4.
2. Mathematical Model
The unsteady hydromagnetic generalized Couette flow of a viscous conducting incompressible fluid is considered in the presence of an imposed transverse magnetic field of strength taking into account the induced electric field . The induced magnetic field is assumed to be small compared with the applied magnetic field and is neglected. Initially, at time , it is assumed that the fluid is stationary with temperature . At , the fluid is subjected to a constant axial pressure gradient while the upper plate moves with a uniform velocity . The plates' surface temperatures are non-uniform with temperature at the lower plate and at the upper moving plate such that as shown in Figure 1.
The temperature-dependent variable electrical conductivity is given as [5] where is the fluid electrical conductivity at the upper plate and is the electrical conductivity variation index. Under the above assumptions, the dimensionless governing equations for the momentum and energy balance can be expressed as where is the dimensionless temperature, is the dimensionless velocity, is the dimensionless time, Re, Le, Ha, Ec, Pr, and are the moving upper plate Reynolds number, the electric field loading parameter, magnetic field parameter, the Eckert number, the Prandtl number, and the pressure-gradient parameter. The appropriate initial and boundary conditions in dimensionless form are given as follows: The following quantities have been utilized in order to obtain the dimensionless governing equations ((2) and (3)): where is the temperature, is the uniform velocity of the upper plate, is the channel width, is the thermal conductivity, is the time, is the dynamic viscosity coefficient, is the specific heat at constant pressure, and is the fluid pressure. Other quantities of interest are the skin-friction coefficient () and the Nusselt number () which are given as The current density and the total current generated within the MHD flow system are given as It is important to note that, for a short circuit configuration where there is no electric field loading, .
3. Entropy Analysis
Hydromagnetic generalized Couette flow is inherently irreversible. This may be due to the exchange of energy between the conducting fluid and the moving plate surface. According to [13], the local volumetric rate of entropy generation for a viscous incompressible conducting fluid in the presence of magnetic field and induced electric field is given as The first term in (8) describes the heat transfer irreversibility while the second and third terms represent irreversibility due to fluid friction and magnetic field, respectively. Using (6), the dimensionless form of local entropy generation rate in (8) is given as where is the temperature difference parameter and is the Brinkman number. The Bejan number, Be, is defined as where , represents the irreversibility due to heat transfer, represents fluid friction and magnetic field irreversibility, and is the irreversibility ratio.
The Bejan number () as shown in (10) has the range . If , then the irreversibility is dominated by the combined effects of fluid friction and magnetic fields, but if , then the irreversibility due to heat transfer dominates the flow system by the virtue of finite temperature differences.
4. Numerical Solution
Our numerical algorithm is based on semi-implicit finite difference schemes [14–19]. Implicit terms are taken at the intermediate time level , where . The discretization of the governing equations is based on a linear Cartesian mesh and uniform grid on which finite differences are taken. We approximate both the second and first spatial derivatives with second-order central differences. The equations corresponding to the first and last grid points are modified to incorporate the boundary conditions. The semi-implicit scheme for the velocity component reads: In (12), it is understood that . The equation for then becomes where The solution procedure for thus reduces to inversion of tridiagonal matrices, which is an advantage over a fully implicit scheme. The semi-implicit integration scheme for the temperature equation is similar to that for the velocity component. Unmixed second partial derivatives of the temperature are treated implicitly as follows: The equation for thus becomes The solution procedure again reduces to inversion of tridiagonal matrices. The schemes ((13) and (16)) were checked for consistency. For , these are first order accurate in time but second order in space. The schemes in [14] have , which improves the accuracy in time to second order. As in [14–18] we, however, use here so that we are free to choose larger time steps and still converge to the steady solutions.
5. Results and Discussion
Unless otherwise stated, we employ the following parameter values:
, , , , , , ,, , , and .
These will be the default values in this work. In the succeeding graphical, if any of these parameter values is not explicitly mentioned, it will be understood that such parameters take on the default values.
5.1. Transient and Steady-State Flow Profiles
We display the transient solutions in Figure 2. The figures show a transient increase in both fluid velocity, Figure 2(a), and temperature, Figure 2(b), until a steady state is reached.
(a)
(b)
5.1.1. Parameter Dependence of Solutions
The response of the velocity and temperature to varying values of the Prandtl number () is illustrated in Figure 3.
(a)
(b)
Larger values of the Prandtl number correspondingly increase the strength of the heat sources in the temperature equation, and hence this strength in turn increases the overall fluid temperature as clearly illustrated in Figure 3(b). Due to the implicit relationship between the Prandtl number and the velocity field via the magnetic field terms, we notice as expected that increases in the Prandtl number have limited effects on the fluid velocities as illustrated in Figure 3(a).
The response of the velocity and temperature to varying values of the Reynolds number () is illustrated in Figure 4.
(a)
(b)
The Reynolds number being coupled to the time-dependent terms shows limited effect on both the steady fluid velocity and temperature as shown in Figures 4(a) and 4(b).
The response of the velocity and temperature to varying values of the magnetic field parameter () is illustrated in Figure 5.
(a)
(b)
Larger values of correspondingly increase the strength of the heat sources in the temperature equation, and hence, this strength in turn increases the overall fluid temperature as clearly illustrated in Figure 5(b). On the other hand, higher values of decrease the strength of the source terms in the velocity equation, and hence this decreased strength in turn increases the overall fluid temperature as clearly illustrated in Figure 5(a).
The effects of on both the velocity and temperature profiles are similar to those of as illustrated in Figures 6(a) and 6(b).
(a)
(b)
Figures 7(a) and 7(b) show that the values of the initial temperature have limited effect on the steady-state flow profiles.
(a)
(b)
The response of the velocity and temperature to varying values of the Eckert number () are illustrated in Figure 8.
(a)
(b)
The effects of the Eckert number are similar to those for the Prandtl number. The effects of the electrical conductivity variation index on the velocity and temperature profiles are illustrated in Figure 9.
(a)
(b)
As expected, since , an increase in the electrical conductivity variation index () correspondingly decreases the values of . The source terms in the temperature equation will thus decrease in magnitude, whereas those in the velocity equation become less negative. This explains the slight decreases noticed in the fluid temperature, Figure 9(b), as well as the slight increases in the fluid velocity, Figure 9(a).
5.2. Wall Shear Stress and Wall Heat Transfer Rate
The wall shear stress and wall heat transfer rates are illustrated in Figures 10 and 11, respectively, with varying and .
The results in Figures 10 and 11 simply summarize those in the previous section on parameter dependence of solutions. In particular, the behavior of the velocity and temperature with respect to certain parameters can be summarized in terms of the behavior of their respective gradients. The results of the current density similarly reflect the corresponding results on parameter dependence as shown in Figure 12.
5.3. Entropy Generation
In this section, we plot the entropy generation rate across the channel under varying parameter conditions. Except for Figure 13, all graphs are otherwise drawn at the time .
Figures 13, 14, 15, 16, 17, and 18 show the expected results for . In particular, parameters that increase the velocity and temperature gradients also increase the entropy generation rate and vice versa. In Figures 13–18, the values of the entropy generation rate expectedly vary across the channel including the walls in response to the evolution of the velocity and temperature gradients with increasing parameter values.
5.4. The Bejan Number
In this section, we plot the Bejan number across the channel under varying parameter conditions. The analysis in this section is similar to that for the previous section with now replaced by .
Figures 19, 20, 21, 22, 23, and 24 show as expected that parameters which increase the entropy generation rate will correspondingly decrease the Bejan number and vice versa. In the vicinity of the walls, the strength of the fluid parameters will determine which mode of irreversibility dominates over the other.
6. Conclusion
We computationally investigate the inherent irreversibility in an unsteady hydromagnetic generalized Couette flow with variable electrical conductivity in the presence of induced electric field. We also notice that, due to the nature of the source terms, the fluid velocity and temperature will each decrease (resp., increase) with a corresponding increase in the parameters that decrease/increase the magnitudes of the source terms. We have also demonstrated computationally that parameters which increase the entropy generation rate will correspondingly decrease the Bejan number and vice versa. In particular, as the flow profiles vary in shape from linear to “parabolic” in response to varying parameter values, the velocity/temperature gradients correspondingly change in magnitude leading to noticeable effects in the entropy generation rates and Bejan numbers.