Teach Second Law of Thermodynamics via Analysis of Flow through Packed Beds and Consolidated Porous Media

: The second law of thermodynamics is indispensable in engineering applications. It allows us to determine if a given process is feasible or not, and if the given process is feasible, how e ﬃ cient or ine ﬃ cient is the process. Thus, the second law plays a key role in the design and operation of engineering processes, such as steam power plants and refrigeration processes. Nevertheless students often ﬁnd the second law and its applications most di ﬃ cult to comprehend. The second law revolves around the concepts of entropy and entropy generation. The feasibility of a process and its e ﬃ ciency are directly related to entropy generation in the process. As entropy generation occurs in all ﬂow processes due to friction in ﬂuids, ﬂuid mechanics can be used as a tool to teach the second law of thermodynamics and related concepts to students. In this article, ﬂow through packed beds and consolidated porous media is analyzed in terms of entropy generation. The link between entropy generation and mechanical energy dissipation is established in such ﬂows in terms of the directly measurable quantities such as pressure drop. Equations are developed to predict the entropy generation rates in terms of superﬁcial ﬂuid velocity, porous medium characteristics, and ﬂuid properties. The predictions of the proposed equations are presented and discussed. Factors a ﬀ ecting the rate of entropy generation in ﬂow through packed beds and consolidated porous media are identiﬁed and


Introduction
Thermodynamics is a difficult subject to learn and teach. It is considered to be one of the most abstract disciplines of the physical sciences [1]. Students all over the world face difficulties in learning thermodynamics. More specifically, it is the second law of thermodynamics dealing with entropy and entropy production that is difficult for students to fully comprehend. The second law of thermodynamics simply states that all real (irreversible) processes are accompanied by the production of entropy in the universe. However, the students of thermodynamics often find it difficult to visualize entropy generation in real processes at a mechanistic level. The quantification of entropy generation in real processes is an equally problematic issue for students. What causes entropy generation and how do we quantify entropy generation? These are some of the fundamental questions faced by students. Moreover, there are no instruments which can be used to directly measure entropy and entropy generation in real processes.
Entropy generation occurs in flow of all real (viscous) fluids [2]. When a viscous fluid is forced to flow through any geometry, such as a pipe, viscous stresses and velocity gradients are established.
For fluid to remain in motion, work has to be done against the viscous stresses which oppose its motion. Consequently, part of the mechanical energy of fluid is dissipated into frictional heat (internal energy) during motion. The dissipation of highly ordered mechanical energy into disorderly internal energy is reflected in entropy generation. Thus, the analysis of fluid mechanics problems and measurement of the appropriate flow variables such as pressure loss could be used as a tool to demonstrate and quantify entropy generation in real processes.
The main objectives of this article are: (a) to analyze the flow of viscous fluids through packed beds of discrete particles and through consolidated porous media in terms of entropy generation; (b) to quantify entropy generation in flow through packed beds and consolidated porous media in terms of directly measureable quantities such as pressure loss as a function of flow rate of fluid; and (c) to discuss various factors which affect the rate of entropy generation in flow through packed beds and consolidated porous media.

Brief Review of the Second Law of Thermodynamics
The second law of thermodynamics could be stated in several different but equivalent ways. The classical statements of the second law are the Kelvin-Planck statement and the Clausius statement.
The Kelvin-Planck statement says that "It is impossible for a system operating in a cycle and connected to a single heat reservoir to produce a positive amount of work in the surroundings [3]". According to the Kelvin-Planck statement, it is impossible to build a heat engine shown schematically in Figure 1 where heat absorbed from a heat reservoir is completely converted into work without altering the properties of the system. Mathematically, the Kelvin-Planck statement can be expressed as [4]: W cycle ≤ 0 (single heat reservoir) (1) where the system communicates thermally only with a single heat reservoir. The sign convention for work (W) used in this article is W = positive, if it is produced by the system (flows out of the system to surroundings) and W = negative, if it is absorbed by the system (flows into the system from surroundings). Thus, no cyclic process is possible where W cycle > 0 using a single heat reservoir. The equality in Equation (1) is valid for a reversible process and the inequality is valid for an irreversible process. For a reversible process, no work is produced or destroyed, that is, W cycle = 0 whereas work is destroyed when the process is irreversible, that is, W cycle < 0 (assuming that the system communicates thermally only with a single heat reservoir).
Fluids 2019, 4, x 2 of 22 motion. Consequently, part of the mechanical energy of fluid is dissipated into frictional heat (internal energy) during motion. The dissipation of highly ordered mechanical energy into disorderly internal energy is reflected in entropy generation. Thus, the analysis of fluid mechanics problems and measurement of the appropriate flow variables such as pressure loss could be used as a tool to demonstrate and quantify entropy generation in real processes.
The main objectives of this article are: (a) to analyze the flow of viscous fluids through packed beds of discrete particles and through consolidated porous media in terms of entropy generation; (b) to quantify entropy generation in flow through packed beds and consolidated porous media in terms of directly measureable quantities such as pressure loss as a function of flow rate of fluid; and (c) to discuss various factors which affect the rate of entropy generation in flow through packed beds and consolidated porous media.

Brief Review of the Second Law of Thermodynamics
The second law of thermodynamics could be stated in several different but equivalent ways. The classical statements of the second law are the Kelvin-Planck statement and the Clausius statement.
The Kelvin-Planck statement says that "It is impossible for a system operating in a cycle and connected to a single heat reservoir to produce a positive amount of work in the surroundings [3]". According to the Kelvin-Planck statement, it is impossible to build a heat engine shown schematically in Figure 1 where heat absorbed from a heat reservoir is completely converted into work without altering the properties of the system. Mathematically, the Kelvin-Planck statement can be expressed as [4]: where the system communicates thermally only with a single heat reservoir. The sign convention for work ( ) used in this article is = positive, if it is produced by the system (flows out of the system to surroundings) and = negative, if it is absorbed by the system (flows into the system from surroundings). Thus, no cyclic process is possible where > 0 using a single heat reservoir.
The equality in Equation (1) is valid for a reversible process and the inequality is valid for an irreversible process. For a reversible process, no work is produced or destroyed, that is, = 0 whereas work is destroyed when the process is irreversible, that is, < 0 (assuming that the system communicates thermally only with a single heat reservoir).   The Clausius statement of the second law [3] says that "It is impossible for a system operating in a cycle to have its sole effect the transfer of heat from a low temperature heat reservoir to a high temperature heat reservoir". According to the Clausius statement, it is impossible to construct a device based on the scheme shown schematically in Figure 2 where the sole result of the process is the transfer of heat from a cooler body at low temperature (T L ) to a hotter body at high temperature (T H ). temperature heat reservoir". According to the Clausius statement, it is impossible to construct a device based on the scheme shown schematically in Figure 2 where the sole result of the process is the transfer of heat from a cooler body at low temperature ( ) to a hotter body at high temperature ( ).
Another powerful statement of the second law of thermodynamics is that all irreversible (real) processes are accompanied by entropy generation in the universe [5], that is: where , is the total amount of entropy generated in the universe (system + surroundings), ∆ is the entropy change of the system and ∆ is the entropy change of the surroundings. The equality in Equation (2) is valid for a reversible process and the inequality is valid for an irreversible process. According to this statement of the second law, no process is possible for which , < 0. Entropy is a measure of disorderliness of a system. According to the Boltzmann entropy equation, the entropy of a system can be expressed as [3]: where Ω is the number of possible configurations of the system and is the Boltzmann constant. The larger the number of possible configurations of the system, greater is the disorderliness of the system and higher is the entropy. Therefore, we can interpret the second law of thermodynamics in yet another way, that is, "Only those processes are possible processes which lead to an increase in the disorderliness of the universe". Another powerful statement of the second law of thermodynamics is that all irreversible (real) processes are accompanied by entropy generation in the universe [5], that is: where S G,univ is the total amount of entropy generated in the universe (system + surroundings), ∆S sys is the entropy change of the system and ∆S surr is the entropy change of the surroundings. The equality in Equation (2) is valid for a reversible process and the inequality is valid for an irreversible process. According to this statement of the second law, no process is possible for which S G,univ < 0. Entropy is a measure of disorderliness of a system. According to the Boltzmann entropy equation, the entropy of a system can be expressed as [3]: where Ω is the number of possible configurations of the system and k B is the Boltzmann constant. The larger the number of possible configurations of the system, greater is the disorderliness of the system and higher is the entropy. Therefore, we can interpret the second law of thermodynamics in yet another way, that is, "Only those processes are possible processes which lead to an increase in the disorderliness of the universe". The scheme shown in Figure 1 is impossible as it decreases the disorder of the universe, that is, S G,univ < 0. Here ∆S sys = 0 but ∆S surr < 0 as the surroundings heat reservoir loses heat. For this scheme, Equation (2) gives: where T is the absolute temperature of the heat reservoir and Q is the heat transferred. As Q cycle = W cycle from the first law of thermodynamics, Equation (4) reduces to Equation (1), that is, W cycle ≤ 0 when there is only one heat reservoir involved. Similarly, the scheme shown in Figure 2 is impossible as it decreases the disorder of the universe: Note that ∆S cold−body < 0 as it loses heat whereas ∆S hot−body > 0 as it gains heat. However, due to different temperatures of the cold and hot bodies, ∆S cold−body > ∆S hot−body . Thus, the Kelvin-Planck and Clausius statements of the second law of thermodynamics are special cases of the statement of the second law expressed in the form of Equation (2).
For a flow process (see Figure 3), the second law of thermodynamics can be written as [2]: .
where . S G is the rate of entropy generation, s is the entropy per unit mass of fluid, ρ is the fluid density, n is unit outward normal to the control surface, → V is fluid velocity vector, A is the control surface area, t is time, ϑ is the volume of the control volume, . Q i is the rate of heat transfer to control volume from ith heat reservoir at an absolute temperature of T i , the subscripts CV and Surr refer to control volume and surroundings, respectively. As noted earlier, the equality in Equation (6) is valid for a reversible (frictionless) process and the inequality is valid for an irreversible process. The surface integral ρs(n· → V)dA is the net outward flow of entropy (associated with mass) across the entire control surface. The volume integral ∂(ρs) ∂t dϑ is the rate of accumulation of entropy within the entire control volume (assumed to be fixed and non-deforming). For a control volume with one inlet and one outlet, Equation (6), under steady state condition, reduces to: where . m is the mass flow rate. The quantification of entropy generation in real processes is important from a practical point of view as entropy generation is directly related to the efficiency of the process. Higher the rate of entropy generation in a process, lower is the thermodynamic efficiency of the process. According to the Gouy-Stodola theorem [6] of thermodynamics, the loss of power or work potential in a real process, due to irreversibilities in the process, is directly proportional to the total rate of entropy generation. Thus: where . W lost is the rate of work lost (wasted) as a result of irreversibilities in the process. As an example of a flow process, consider flow through a control volume shown in Figure 3. where h is the specific enthalpy of fluid, V is the fluid velocity, g is the acceleration due to gravity, z is the elevation, .
Q is the rate of heat transfer, and . W sh is the rate of shaft work. Neglecting kinetic and potential energy changes, Equation (9) simplifies to: where Q is the heat transfer per unit mass of fluid. The second law, Equation (7), can be written as: .
where T o is the absolute temperature of the heat reservoir (surroundings). For the flow process to be reversible: .
where Q rev is the heat transfer per unit mass of fluid for a reversible process. From Equations (10) and (13): The lost work, .
W lost , is defined as: .
Equation (17) is the Gouy-Stodola theorem [6] of thermodynamics. Thus, the rate of work lost due to irreversibilities is directly proportional to the rate of entropy generation. The thermodynamic efficiency of a flow process can be defined as: In a reversible process, . W lost = 0 and = 1. In any irreversible process, . W lost > 0 and < 1. If no work is produced in the process, that is, actual . W sh is zero, then all the work potential is lost due to irreversibilities in the process and consequently, = 0.
In a reversible process, = 0 and = 1. In any irreversible process, > 0 and < 1. If no work is produced in the process, that is, actual is zero, then all the work potential is lost due to irreversibilities in the process and consequently, = 0.

Flow through Unconsolidated and Consolidated Porous Media
Porous medium is a composite material in that it consists of two phases, namely pores (voids, free space pockets) and solid-phase. The pores may be occupied by a fluid (gas, oil, water, etc.). A large variety of natural and synthetic materials are porous in nature. Examples include: underground oil reservoirs, ceramics, solid foams, sand filters, wood, and packed beds of particles used widely in chemical engineering applications. The pores of a porous medium usually form a three-dimensional inter-connected network and, therefore, fluids can flow through the porous medium. If all pores of a porous medium are inter-connected, then the porosity (ε) of the porous medium is simply the fraction of the total volume of the medium that is occupied by the pores. Thus, the fraction of the total volume that is occupied by the solid phase is (1 − ε). When some pores are isolated or disconnected or have dead ends, then the effective porosity, defined as the ratio of connected void volume to total volume of the medium, is lower than the total porosity.
Porous media could be classified as consolidated or unconsolidated. In consolidated porous medium, the solid phase is basically a single piece of material or the grains of the solid phase are cemented or fused together to form a single piece of solid phase. In unconsolidated porous media, on the other hand, the grains or the particles of the solid phase are not cemented together and, therefore, the porous medium is a multi-particle system like a packed bed of individual (un-cemented, un-glued) particles. Flow of single-phase Newtonian fluids (gas, water, oil) through packed beds and consolidated porous media is important from a practical point of view.

Pressure-Loss in Flow through Packed Beds
Consider flow of a Newtonian fluid through a horizontal cylindrical packed bed of particles with bed diameter D b and bed length L b . As flow through a packed bed is quite complex, a rigorous theoretical derivation of pressure drop-flow rate relationship is not possible. Only approximate models have been developed. In one approach, used widely to model flow through packed bed of particles, the packed bed is visualized as a bundle of identical capillary tubes [7][8][9]. In its simplest form, the capillary tube bundle model assumes that the capillary tubes are straight, cylindrical of constant cross-section (uniform radius), and parallel (all oriented in the same direction), as shown in Figure 4. The average velocity in any capillary tube is the same as that in the packed bed. The capillary diameter is equal to the hydraulic diameter of the bed, d H , defined as: where ε is the void fraction or porosity, V p is the volume of a single particle, and S p is the surface area of a single particle. In writing the above expression for d H , it is assumed that all the particles are identical and that the wetted surface of the cylindrical container walls of the bed is negligible as compared with the total wetted surface area of the particles. For a bed of identical spherical particles, V p /S p = d p /6, where d p is the particle diameter. Thus, the hydraulic diameter becomes:  where is the cross-sectional area of bed through which the fluid flows. The porosity of the bed is given as: Therefore, the average or interstitial velocity of fluid in the bed is: This expression of average velocity does not consider the tortuous path taken by fluid in the bed. Due to tortuosity of the bed, the actual average velocity of fluid in the bed is larger than that given by Equation (24). The tortuosity of a bed is defined as: where is average length of the tortuous path taken by the fluid and is the straight length of the bed. To account for the tortuosity, the actual average velocity of the fluid in the bed can be expressed as follows: In laminar flow of a Newtonian fluid through a cylindrical tube, the pressure gradient in the direction of flow is given as: where is the fluid viscosity and is the tube diameter. Replacing with hydraulic diameter and the constant factor of 32 with , the pressure-gradient in a capillary tube model of the packed- The superficial velocity of fluid (V s ) is defined is: where . ϑ is the volumetric flow rate of fluid and A b is the total cross-sectional area of the bed. Thus, V s is the velocity of fluid in the bed if no particles were present in the bed. The average velocity V of fluid in the bed, also called interstitial velocity, is defined as: where A f low is the cross-sectional area of bed through which the fluid flows. The porosity of the bed ε is given as: Therefore, the average or interstitial velocity of fluid in the bed is: This expression of average velocity does not consider the tortuous path taken by fluid in the bed. Due to tortuosity of the bed, the actual average velocity of fluid in the bed is larger than that given by Equation (24). The tortuosity τ of a bed is defined as: where L e is average length of the tortuous path taken by the fluid and L b is the straight length of the bed. To account for the tortuosity, the actual average velocity of the fluid in the bed can be expressed as follows: In laminar flow of a Newtonian fluid through a cylindrical tube, the pressure gradient in the direction of flow is given as: where µ is the fluid viscosity and D is the tube diameter. Replacing D with hydraulic diameter d H and the constant factor of 32 with C, the pressure-gradient in a capillary tube model of the packed-bed model can be expressed as: Note that flow in the bed is assumed to be laminar here. The constant C is expected to be larger than 32, as the path of fluid in the bed is not straight. The fluid follows a tortuous path in the bed and consequently, the pressure drop over a certain straight length is expected to be more than that observed over the same length if the fluid path was non-tortuous. Therefore, the constant C is taken as 32 τ where τ is the tortuosity. Upon substitution of the expressions for d H and V, and taking C = 32 τ, Equation (28) gives: Equation (29) assumes that the cross-section of the representative flow passage (capillary tube) in the porous medium is circular and constant. In reality the fluid moves through converging-diverging flow passages of non-uniform cross-sections. In converging-diverging flows, fluid experiences stretching or extensional deformation. The elongation or stretching of fluid elements results in additional dissipation of energy and pressure drop. Thus, this equation needs to be modified further as follows: where K is an empirical factor that takes into account the influence of non-constant cross-section of flow passage on pressure-gradient. The tortuosity τ is often taken to be √ 2 for random packing of uniform spheres [9]. It is generally a function of the porosity of the porous medium [10]. For example, the following model, based on the Maxwell equation for electrical conductivity of composite, is often used to describe the relationship between tortuosity and porosity [10]: When τ = √ 2 and K = 25/24, the following well-known Blake-Kozeny equation for laminar flow through randomly packed bed of uniform spheres is obtained [7]: When τ = √ 2 and K = 5/4, the following Carman-Kozeny equation, another well-known equation for laminar flow through randomly packed bed of uniform spheres, is obtained [11]: Both Blake-Kozeny and Carman-Kozeny equations are popular in the literature although they use different values of the factor K. The observed difference in K values is probably related to differences in particle-shape, surface roughness of particles, and porosity of bed. Note that these equations are restricted to only laminar flow through bed of nearly uniform spheres. The packed-bed Reynolds number Re b , defined below, should be less than 10 for flow to be laminar in the bed.
For turbulent flow through packed beds (Re b > 1000), the following Burke-Plummer equation is often used [12]: In the transition region, the following equation, obtained by superposition of Blake-Kozeny and Burke-Plummer equations, proposed by Ergun [13] is widely used: As the Ergun equation, Equation (36), is obtained by the superposition of laminar and turbulent expressions, it is valid over the full range of the packed-bed Reynolds number Re b . The Ergun equation could also be re-cast in the following form [14]: where f b is the packed-bed friction factor defined as: The Ergun equation, Equation (36) or (37), is used extensively in the literature to describe pressure loss in packed beds. However, the following points should be kept in mind when using the Ergun equation: (a) it is applicable to unconsolidated beds of nearly uniform-size spherical particles with appreciable roughness. For smooth spheres, it tends to over-predict f b in the high Re b range (Re b > 700) [14][15][16]; (b) if the bed particles are non-uniform in sizes, the Sauter mean diameter of particles should be used in the application of the Ergun equation; (c) when particle shape deviates significantly from a sphere, the Ergun equation tends to under-predict f b [17]; and (d) wall effects can be important when D b /d p < 10. When wall effects are important, the experimental data show deviation from the predictions of the Ergun equation [18]. can be important when < 10 ⁄ . When wall effects are important, the experimental data show deviation from the predictions of the Ergun equation [18]. Figure 5 shows the prediction of the packed bed friction factor as a function of packed bed Reynolds number using the Ergun equation, Equation (37). In the limits of low and high , the Ergun equation predictions overlap with the predictions of the Blake-Kozeny equation (low ) and the Burke-Plummer equation (high ), as expected.

Entropy Generation in Flow through Packed Beds
Consider steady flow of an incompressible fluid through a bed of particles (see Figure 6). We can apply the following macroscopic mechanical energy balance between the inlet and outlet of the bed [2]: where is the rate of mechanical energy dissipation due to friction in fluid, ∆ is the potential energy change per unit mass of fluid, ∆( ) is the kinetic energy change per unit mass of fluid, and ∆ is the pressure change of the fluid. Neglecting potential and kinetic energy changes and taking = 0, the mechanical energy balance gives:

Entropy Generation in Flow through Packed Beds
Consider steady flow of an incompressible fluid through a bed of particles (see Figure 6). We can apply the following macroscopic mechanical energy balance between the inlet and outlet of the bed [2]: where . F l is the rate of mechanical energy dissipation due to friction in fluid, ∆ϕ is the potential energy change per unit mass of fluid, ∆(KE) is the kinetic energy change per unit mass of fluid, and ∆P is the pressure change of the fluid. Neglecting potential and kinetic energy changes and taking . W sh = 0, the mechanical energy balance gives: The rate of mechanical energy dissipation per unit length of packed bed ( . F l ) can be expressed as: .
From the second law of thermodynamics, Equation (11), we get: Note that we are assuming flow to be adiabatic with negligible heat transfer. There is no entropy generation in the surroundings. All the entropy is generated within the fluid inside the packed bed and the rate of entropy generation is the net rate of increase in entropy of the flowing stream.
The rate of mechanical energy dissipation per unit length of packed bed ( ) can be expressed as: From the second law of thermodynamics, Equation (11), we get: Note that we are assuming flow to be adiabatic with negligible heat transfer. There is no entropy generation in the surroundings. All the entropy is generated within the fluid inside the packed bed and the rate of entropy generation is the net rate of increase in entropy of the flowing stream.
We can now relate entropy change of the fluid stream to pressure change. For pure substances, the relationship between entropy and other state variables is given as [5]: where is the absolute temperature. From the first law of thermodynamics, Equation (9), the enthalpy change is zero in the absence of heat transfer and shaft work for steady flow in a horizontal bed. Consequently, Equation (43) reduces to: Assuming incompressible flow and constant temperature, Equation (44) upon integration gives: From Equations (42) and (45), it follows that: The subscript "univ" has been removed from , as , is simply the rate of entropy generation within the fluid inside a packed bed. From Equations (41) and (46), we can also express the rate of entropy generation in a packed bed on a unit volume basis as: where is the total cross-sectional area of the bed and is the rate of entropy generation per unit volume of the bed. From Equations (38) and (47), we get: The packed bed friction factor for laminar flow ( < 10) is given as: We can now relate entropy change of the fluid stream to pressure change. For pure substances, the relationship between entropy and other state variables is given as [5]: where T is the absolute temperature. From the first law of thermodynamics, Equation (9), the enthalpy change is zero in the absence of heat transfer and shaft work for steady flow in a horizontal bed. Consequently, Equation (43) reduces to: Assuming incompressible flow and constant temperature, Equation (44) upon integration gives: From Equations (42) and (45), it follows that: .
The subscript "univ" has been removed from . S G,univ as . S G,univ is simply the rate of entropy generation within the fluid inside a packed bed. From Equations (41) and (46), we can also express the rate of entropy generation in a packed bed on a unit volume basis as: where A b is the total cross-sectional area of the bed and . S G is the rate of entropy generation per unit volume of the bed. From Equations (38) and (47), we get: The packed bed friction factor f b for laminar flow (Re b < 10) is given as: Consequently, Equation (48) reduces to: .
Thus, entropy generation rate per unit volume of bed in steady laminar flow of a Newtonian fluid is directly proportional to the fluid viscosity and to the square of the superficial fluid velocity in the bed. The entropy generation rate also depends on the particle diameter, the bed porosity ε, and the temperature.
Thus, entropy generation rate per unit volume of bed in steady turbulent flow of a Newtonian fluid is independent of the fluid viscosity and is directly proportional to the cube of the superficial fluid velocity in the bed. The entropy generation rate also depends on the fluid density, the particle diameter, the bed porosity ε, and the temperature.
To cover the full range of packed bed Reynolds number Re b , we substitute the Ergun equation, Equation (37), into Equation (48) to obtain the following equation valid over the full range of packed bed Reynolds number Re b : Note that Equation (53) is simply the superposition of Equations (50) and (52).

Pressure-Loss in Flow through Consolidated Porous Media
Laminar flow in consolidated porous media (see Figure 7) is usually described by Darcy's law given as follows [19]: where k is the permeability of the porous medium. Permeability is a measure of the ability of porous medium to conduct fluid. Higher k means lower resistance to flow and consequently, larger flow rate for the same pressure gradient. k has the dimensions of (length) 2 . It is usually expressed in darcies. For example, 1 darcy = 1 (cm/s). centipoise/(atm/cm) = 9.87 × 10 −13 m 2 .
= 150 Consequently, Equation (48) reduces to: Thus, entropy generation rate per unit volume of bed in steady laminar flow of a Newtonian fluid is directly proportional to the fluid viscosity and to the square of the superficial fluid velocity in the bed. The entropy generation rate also depends on the particle diameter, the bed porosity , and the temperature.
The packed bed friction factor in turbulent flow ( > 1000) is given as [7][8][9]: Consequently, Equation (48) reduces to: Thus, entropy generation rate per unit volume of bed in steady turbulent flow of a Newtonian fluid is independent of the fluid viscosity and is directly proportional to the cube of the superficial fluid velocity in the bed. The entropy generation rate also depends on the fluid density, the particle diameter, the bed porosity , and the temperature.
To cover the full range of packed bed Reynolds number , we substitute the Ergun equation, Equation (37), into Equation (48) to obtain the following equation valid over the full range of packed bed Reynolds number : Note that Equation (53) is simply the superposition of Equations (50) and (52).

Pressure-Loss in Flow through Consolidated Porous Media
Laminar flow in consolidated porous media (see Figure 7) is usually described by Darcy's law given as follows [19]: where is the permeability of the porous medium. Permeability is a measure of the ability of porous medium to conduct fluid. Higher means lower resistance to flow and consequently, larger flow rate for the same pressure gradient. has the dimensions of (length) 2 . It is usually expressed in darcies. For example, 1 darcy = 1 (cm/s). centipoise/(atm/cm)= 9.87 × 10 −13 m 2 . Darcy's law is equally applicable to unconsolidated porous medium such as packed bed of free (un-cemented/un-fused) particles. Thus, the Blake-Kozeny equation or the Carman-Kozeny equation are special cases of the Darcy law. Upon comparison of the Darcy law with Blake-Kozeny/ Carman-Kozeny equations, the permeability of a packed bed of uniform size spherical particles can be expressed as follows: Equation (55) gives permeability based on the Blake-Kozeny equation whereas Equation (56) gives permeability based on the Carman-Kozeny equation. Darcy's law is restricted to flows where viscous forces dominate over the inertial forces. At high flow rates, the flow in the porous medium becomes turbulent and Darcy's law is no longer valid. The pressure-gradient in the turbulent regime is much higher compared with that in the laminar regime at the same superficial velocity V s . In order to extend the Darcy law to turbulent regime, Forchheimer [19] modified the Darcy law by adding a quadratic term as follows: where the first term on the right side of the equation reflects the viscous effects and the second term reflects the inertial effects. β is called the non-Darcy flow coefficient with a dimension of 1/length. This equation could be re-written in dimensionless form as: where f pm is the friction factor for flow in consolidated porous medium and Re pm is the porous-medium Reynolds number, also referred to as Forchheimer number, defined as: At low Reynolds number (Re pm < 0.1), the first term on the right side of the dimensionless Forchheimer equation, Equation (58), can be neglected and it reduces to Darcy's law. At high Reynolds number (Re pm > 10), the second term on the right side of the dimensionless Forchheimer equation, Equation (58), can be neglected and it reduces to f pm = 1. Figure 8 shows the prediction of f pm as a function of Re pm using the Forchheimer equation, Equation (58).

=
(59) At low Reynolds number ( < 0.1), the first term on the right side of the dimensionless Forchheimer equation, Equation (58), can be neglected and it reduces to Darcy's law. At high Reynolds number ( > 10), the second term on the right side of the dimensionless Forchheimer equation, Equation (58), can be neglected and it reduces to = 1. Figure 8 shows the prediction of as a function of using the Forchheimer equation, Equation (58). The relationships between f pm and f b , and Re pm and Re b are as follows:

Entropy Generation in Flow through Consolidated Porous Media
Consider one-dimensional flow of an incompressible Newtonian fluid in consolidated porous medium, as shown in Figure 9. The rate of entropy generation in consolidated porous medium per unit volume is given by Equation (47), re-written as: The porous medium friction factor f pm for laminar flow (Re pm < 0.1) is given as: Consequently, Equation (63) reduces to: .
Thus, entropy generation rate per unit volume of the porous medium in steady laminar flow of a Newtonian fluid is directly proportional to the fluid viscosity, inversely proportional to the porous medium permeability, and directly proportional to the square of the superficial fluid velocity in the medium. The entropy generation rate also depends on the temperature. volume is given by Equation (47), re-written as: The porous medium friction factor for laminar flow ( < 0.1) is given as: Consequently, Equation (63) reduces to: Thus, entropy generation rate per unit volume of the porous medium in steady laminar flow of a Newtonian fluid is directly proportional to the fluid viscosity, inversely proportional to the porous medium permeability, and directly proportional to the square of the superficial fluid velocity in the medium. The entropy generation rate also depends on the temperature.  The porous medium friction factor in turbulent flow (Re pm > 10) is given as: Consequently, Equation (63) reduces to: .
Thus, entropy generation rate per unit volume of the porous medium in steady turbulent flow of a Newtonian fluid is independent of the fluid viscosity and is directly proportional to the cube of the superficial fluid velocity in the bed. The entropy generation rate also depends on the fluid density, the non-Darcy flow coefficient β, and the temperature.
To cover the full range of porous medium Reynolds number Re pm , we substitute the Forchheimer equation, Equation (58), into Equation (63) to obtain the following equation valid over the full range of porous medium Reynolds number Re pm : Equation (68) is simply the superposition of Equations (65) and (67).

Laminar Flow
The entropy generation rate in laminar flow through packed beds, per unit volume of the bed, as a function of superficial bed velocity is shown in Figure 10

Laminar Flow
The entropy generation rate in laminar flow through packed beds, per unit volume of the bed, as a function of superficial bed velocity is shown in Figure 10 for the following conditions: = 298.15 K, = 1 mm, = 18.5 µPa·s, and = 1.184 kg/m 3 . Equation (50) is used to generate the plots. For a given value of the bed porosity , the entropy generation rate increases with the increase in the superficial velocity of the fluid in the bed. The increase in with the increase in is larger if the bed porosity is small. Additionally, the entropy generation rate falls sharply with the increase in the bed porosity. These results are as expected. With the increase in the fluid velocity, the rate of mechanical energy dissipation into frictional heat increases due to increases in viscous stresses and velocity gradients in the fluid. Consequently, there occurs an increase in the rate of entropy generation. When the bed porosity is The increase in . S G with the increase in V s is larger if the bed porosity ε is small. Additionally, the entropy generation rate falls sharply with the increase in the bed porosity. These results are as expected. With the increase in the fluid velocity, the rate of mechanical energy dissipation into frictional heat increases due to increases in viscous stresses and velocity gradients in the fluid. Consequently, there occurs an increase in the rate of entropy generation. When the bed porosity is increased, the bed structure becomes more open to fluid flow resulting in a decrease in the resistance to fluid motion and mechanical energy dissipation and, hence, a decrease in the rate of entropy generation. Figure 11 shows the effect of fluid viscosity on the rate of entropy generation in laminar flow. The plots are generated using Equation (50). With the increase in fluid viscosity, the resistance to fluid motion increases due to an increase in the viscous stresses. With the increase in viscous stresses, the rate of mechanical energy dissipation into frictional heat (internal energy) increases. The conversion of highly ordered mechanical energy into disorderly internal energy is reflected in an increase in the rate of entropy generation. increased, the bed structure becomes more open to fluid flow resulting in a decrease in the resistance to fluid motion and mechanical energy dissipation and, hence, a decrease in the rate of entropy generation. Figure 11 shows the effect of fluid viscosity on the rate of entropy generation in laminar flow. The plots are generated using Equation (50). With the increase in fluid viscosity, the resistance to fluid motion increases due to an increase in the viscous stresses. With the increase in viscous stresses, the rate of mechanical energy dissipation into frictional heat (internal energy) increases. The conversion of highly ordered mechanical energy into disorderly internal energy is reflected in an increase in the rate of entropy generation. The effect of particle diameter on the rate of entropy generation in laminar flow through packed beds is shown in Figure 12. A sharp reduction in the rate of entropy generation occurs with the increase in the particle diameter. With the increase in the particle diameter, the number density of particles is decreased for the same bed porosity . It can be readily shown that the number of particles The effect of particle diameter on the rate of entropy generation in laminar flow through packed beds is shown in Figure 12. A sharp reduction in the rate of entropy generation occurs with the increase in the particle diameter. With the increase in the particle diameter, the number density of particles is decreased for the same bed porosity ε. It can be readily shown that the number of particles per unit volume of bed n p and the fluid-solids contact area per unit volume of bed A solids are: Equation (69) assumes that the particles are uniform spheres of diameter d p . The decrease in the number density of particles with the increase in particle diameter results in a decrease in contact area between the fluid and solids resulting in a decrease in the resistance to fluid motion. Consequently the rate of mechanical energy dissipation into internal energy and hence the rate of entropy generation decrease with the increase in particle diameter.  Figure 13 shows the effect of temperature on the rate of entropy generation in laminar flow through packed beds. The rate of entropy generation decreases with the increase in temperature keeping other factors (fluid properties, bed characteristics) constant. At a given fluid velocity , the rate of mechanical energy dissipation into internal energy is the same at different temperatures as the fluid properties are kept constant. Then why does the rate of entropy generation decrease with the increase in temperature? It so happens that the increase in entropy with the increase in internal energy is sensitive to temperature. For incompressible materials, the rate of change of entropy with respect to internal energy, that is, the derivative ⁄ , is given as: Higher the temperature, smaller is the increase in entropy for the same increase in internal energy.  Figure 13 shows the effect of temperature on the rate of entropy generation in laminar flow through packed beds. The rate of entropy generation decreases with the increase in temperature keeping other factors (fluid properties, bed characteristics) constant. At a given fluid velocity V s , the rate of mechanical energy dissipation into internal energy is the same at different temperatures as the fluid properties are kept constant. Then why does the rate of entropy generation decrease with the increase in temperature? It so happens that the increase in entropy with the increase in internal energy is sensitive to temperature. For incompressible materials, the rate of change of entropy with respect to internal energy, that is, the derivative dS/dU, is given as: Higher the temperature, smaller is the increase in entropy for the same increase in internal energy.
increase in temperature? It so happens that the increase in entropy with the increase in internal energy is sensitive to temperature. For incompressible materials, the rate of change of entropy with respect to internal energy, that is, the derivative ⁄ , is given as: Higher the temperature, smaller is the increase in entropy for the same increase in internal energy.  Figures 14-16 show the entropy generation rates for turbulent flow through packed beds under different conditions. Due to broad ranges of entropy generation rates and superficial fluid velocities, the plots are drawn using log-log scale. The plots on log-log scale are linear with a slope of 3 as expected from Equation (52). With the increases in bed porosity and particle diameter, the entropy generation rates decrease due to a decrease in the resistance to fluid motion and mechanical energy dissipation. There is no dependence of entropy generation on fluid viscosity in turbulent flow. However, fluid density now plays a role. With the increase in fluid density, the rate of entropy generation in turbulent flow increases due to an increase in the resistance to fluid motion and mechanical energy dissipation caused by inertial effects.

Entropy Generation in Flow through Consolidated Porous Media
The entropy generation rate in laminar flow through consolidated porous medium is given by Equation (65). The key factor affecting the vs. behavior is the ratio of fluid viscosity to permeability, ⁄ . Figure 18 shows the plots of vs.
for different values of ⁄ . The plots on a log-log scale are linear with slopes of 2. With the increase in ⁄ ratio, the entropy generation rate increases at any given superficial velocity . When the fluid viscosity is increased (keeping the same), the entropy generation rate increases due to an increase in mechanical energy dissipation caused by viscous stresses. When the permeability of the porous medium is decreased (keeping the same), the porous medium becomes less permeable to fluid flow resulting in larger resistance to fluid motion and hence larger rates of mechanical energy dissipation and entropy generation. Note that is typically in the range of 10 −15 to 10 −12 m 2 for porous sandstones [20].

Entropy Generation in Flow through Consolidated Porous Media
The entropy generation rate in laminar flow through consolidated porous medium is given by Equation (65). The key factor affecting the . S G vs. V s behavior is the ratio of fluid viscosity to permeability, µ/k. Figure 18 shows the plots of . S G vs. V s for different values of µ/k. The plots on a log-log scale are linear with slopes of 2. With the increase in µ/k ratio, the entropy generation rate increases at any given superficial velocity V s . When the fluid viscosity is increased (keeping k the same), the entropy generation rate increases due to an increase in mechanical energy dissipation caused by viscous stresses. When the permeability k of the porous medium is decreased (keeping µ the same), the porous medium becomes less permeable to fluid flow resulting in larger resistance to fluid motion and hence larger rates of mechanical energy dissipation and entropy generation. Note that k is typically in the range of 10 −15 to 10 −12 m 2 for porous sandstones [20]. In turbulent flow through consolidated porous medium, the key factor governing the vs. behavior is where is the non-Darcy flow coefficient. With the increase in , the entropy generation rate increases at any given , as shown in Figure 19. Note that the non-Darcy flow coefficient is inversely related to the porous medium permeability [20]. For sandstones, is typically in the range of 10 8 to 10 12 m −1 . With the increase in , the porous medium becomes less permeable to fluid flow resulting in an increase in the resistance to fluid motion. Consequently, the rates of mechanical energy dissipation and entropy generation increase with the increase in . The fluid density has a similar effect on entropy generation rate in the turbulent regime. In turbulent flow through consolidated porous medium, the key factor governing the . S G vs. V s behavior is ρβ where β is the non-Darcy flow coefficient. With the increase in ρβ, the entropy generation rate increases at any given V s , as shown in Figure 19. Note that the non-Darcy flow coefficient is inversely related to the porous medium permeability [20]. For sandstones, β is typically in the range of 10 8 to 10 12 m −1 . With the increase in β, the porous medium becomes less permeable to fluid flow resulting in an increase in the resistance to fluid motion. Consequently, the rates of mechanical energy dissipation and entropy generation increase with the increase in ρβ. The fluid density has a similar effect on entropy generation rate in the turbulent regime.
generation rate increases at any given , as shown in Figure 19. Note that the non-Darcy flow coefficient is inversely related to the porous medium permeability [20]. For sandstones, is typically in the range of 10 8 to 10 12 m −1 . With the increase in , the porous medium becomes less permeable to fluid flow resulting in an increase in the resistance to fluid motion. Consequently, the rates of mechanical energy dissipation and entropy generation increase with the increase in . The fluid density has a similar effect on entropy generation rate in the turbulent regime. Figure 20 shows the entropy generation rate in a consolidated porous medium over a broad range of fluid superficial velocity (10 −4 ≤ ≤ 10 m/s) covering laminar, transition, and turbulent flow regimes. The plot is generated from Equation (68)

Conclusions
In conclusion, a novel approach is described to teach the second law of thermodynamics via the analysis of flow through packed beds and consolidated porous media. The second law of thermodynamics and the relevant background in fluid mechanics are reviewed briefly. The link between entropy generation and mechanical energy dissipation in flow through packed beds and consolidated porous media is established in terms of the directly measureable pressure loss. Equations are developed to predict the entropy generation rates in a porous medium in terms of the flow variables, fluid properties, and structural properties of the medium. Simulation results dealing with entropy generation in a porous medium are presented and discussed. The proposed approach can be implemented at an undergraduate level either as an experimental exercise dealing with pressure loss measurements in flow through a porous medium such as a packed bed or as a simulation exercise. The material presented in the article is suited for third year engineering students after they have completed introductory courses in fluid mechanics and thermodynamics.