ChE curriculum A Connection Between TRANSPORT PHENOMENA AND THERMODYNAMICS Ross E. Swaney and R. Byron Bird University of Wisconsin-Madison • Madison, WI 53706 he core subjects of thermodynamics and transport dissipation terms are absorbed and cancel in the result. While phenomena, taught in separate courses, tend to be in- easily presented in a graduate course, a version of the main ternalized by students as seemingly unrelated subjects, ideas for an undergraduate course should likewise be possible. Tdespite the significant overlaps in their underlying physical The general procedure involves integrating the microscopic principles. In introductory thermodynamics the macroscopic balances over a macroscopic volume V(t) and then converting energy balance usually is presented by constructing a macro- from the fixed frame to the frame of the material within V(t). scopic envelope, defining the accumulation, and identifying The conversion is accomplished by combining the Leibniz various flow terms through the envelope (e.g., Reference 1, Eq. formula 3.1-4, Reference 2, Eq. 3-63). The result is then readily applied ∂ to many useful problems. Even so, complications may arise in d = f + ⋅ ∫ V(t) f dV ∫ V(t) dV ∫ A(t) n fudA (1) correctly representing the terms that appear. Examples include dt ∂t possible work by shear forces at the boundary, proper integra- with the Gauss divergence theorem tion of work and heat fluxes when not uniform over areas, and ∇ ⋅ ⋅ proper integration of velocities for kinetic energy and internal ∫ V(t) fudV = ∫ A(t) n fudA (2) energies when not uniform over the volume or flow areas. to obtain the relation Macroscopic entropy balances can be developed similarly, ∂ although the macroscopic entropy generation term appears. d f + ∇ ⋅ ∫ V(t) f dV = ∫ V(t) dV ∫ V(t) fudV (3) This macroscopic balance tends to receive less application. dt ∂t In practical cases the generation term often evades simple determination, and application is often limited to idealized R. Byron (Bob) Bird received his B.S. in chemical engineering from the University reversible processes. of Illinois (UIUC) in 1947 and his Ph.D. in chemistry from the University of Wisconsin- Many practical problems involve fluids contained within Madison in 1950. In 1953, Professor Olaf solid boundaries, and so the heat and work flows and their Hougen invited Bob to join the faculty at UW, to take over the responsibility for the courses irreversibilities can be understood by examining a fluid con- in thermodynamics. However, shortly before tinuum with diffusion of momentum (Newton’s law) and heat Bob’s arrival in Madison, Bob Marshall was (Fourier’s law). The microscopic energy and entropy balances elevated to be associate dean of the College [3] of Engineering, and Bob had to take over developed in transport phenomena capture the details of the courses in fluid dynamics and mass behavior at each point. Additionally, these balance equations transfer. In the ensuing years, he worked with W.E. Stewart and E.N. Lightfoot to develop a textbook on transport may be formally integrated to arrive at the macroscopic bal- phenomena and never taught a course in ances. In the process, the proper forms for the various inte- thermodynamics! grated terms are naturally obtained, and irreversibilities are Ross Swaney received his Ph.D. in chemi- precisely identified and explained. Moreover, one may further cal engineering from Carnegie Mellon University in 1983. After four years with show that the fundamental relations taught in thermodynamics the process design department of Atlantic do not require assumptions about reversible changes, but hold Richeld, he joined the faculty at Wisconsin. His teaching experiences include courses in quite generally, and naturally accommodate irreversibilities. thermodynamics, fluid flow, and heat transfer With the assumption only of local equilibrium, most common operations, mass transfer operations, opera- tions and processes “summer lab,” process macroscopic system behaviors can be represented. design, process control, process modeling, Presentation of this analysis can help clarify many of the and graduate courses. His research work includes fundamentals and computational systems for macroscopic complications that may make the thermodynamics balances systems modeling and flowsheet synthesis. seem vague. It can also help students unify their understanding of the two subjects. It is particularly interesting to see how the © Copyright ChE Division of ASEE 2017 Vol. 51, No. 2, Spring 2017 83 In Eqs. (1)–(3), u is the velocity of the boundary of V(t). Below we will choose u to be the local fluid velocityv , and so obtain fundamental closed system macroscopic balances. We note that alternatively choosing V(t) to be fixed in space will lead to open system macroscopic balances after integration and application of Eq. (2) to convert the divergence terms to the surface integrals defining the flow terms. ENERGY BALANCE In the recently published textbook, Introductory Transport Phenomena,[3] the equation of change for the total energy (i.e., kinetic plus internal energy) for a pure fluid is given in Eq. 11.1-7 as ∂ 1 2 1 2 The core ρv +ρUˆ + ∇ ⋅ ρv +ρUˆ v = −∇ ⋅(q+ w)+ρv ⋅ g (4) subjects of ∂t 2 2 thermodynamics Here Uˆ is the internal energy per unit mass of fluid,g is the gravitational acceleration vec- tor, q is the heat flux vector, and w = [π⋅ v] = (pδ+ τ)⋅ v = pv +[τ ⋅ v] is called the work and transport flux vector; p is the pressure, δ is the unit tensor, π is the total stress tensor, v(r,t) is the local fluid velocity vector with respect to fixed coordinates, and τ is the viscous stress phenomena, tensor. Eq. (4) is obtained in Reference 3 by applying the law of conservation of energy taught in to a differential element fixed in space, through which a pure fluid is flowing. In the textbook, it is intimated that one can get from Eq. (4) to the usual statement of separate the first law of thermodynamics courses, ∆U = Q + W (5) tend to be although a justification of this statement is not presented. This equation states that the internal energy of a closed system changes as one goes from one equilibrium state to internalized another equilibrium state because heat is added to the system and/or work is done on the system. Although some might say that it is obvious that Eq. (5) comes from Eq. (4), we by students feel that a presentation of the proof is warranted. as seemingly We select a volume V(t) moving and deforming with the fluid, and this will be the system of interest; this element does not exchange mass with the surrounding fluid. The unrelated volume element V(t) has a surface area A(t), every element of which is moving with the local velocity v(r, t) . Integration of Eq. (4) over the volume V(t) gives subjects, ∂ 1 ρ 2 +ρ ˆ + ∇ ⋅ 1 ρ 2 +ρ ˆ ∫ V(t) v UdV ∫ V(t) v UvdV despite the ∂t 2 2 = − ∇ ⋅(q+ w)dV+ ρv ⋅ gdV (6) significant ∫ V(t) ∫ V(t) overlaps in their Application of Eq. (3) then gives d 1 ρ 2 +ρ ˆ = − ⋅ + + ρ ⋅ ∫ V(t) v UdV ∫ A(t)n (q w)dA ∫ V(t) v gdV (7) underlying dt 2 physical We may integrate this equation with respect to time from t1 to t2 to get t 2 principles. t 2 t 2 1 ρ 2 +ρ ˆ = − ⋅ + + ρ ⋅ ∫ V(t) v UdV ∫ ∫ A(t)n (q w)dAdt ∫ ∫ V(t) v gdVdt (8) t1 t1 2 t1 We now define the following quantities = ρ ˆ = U(t) ∫ V(t) UdV the total internal energy within V(t) (8.1) = 1 ρ 2 = K(t) ∫ V(t) v dV the total kinetic energy within V(t) (8.2) 2 84 Chemical Engineering Education t = − 2 ⋅ = Q ∫ ∫ A(t) n qdAdt the total heat added to V(t)between t1 and t2 (8.3) t1 t t = − 2 ⋅ + 2 ρ ⋅ W ∫ ∫ A(t) n wdAdt ∫ ∫ V(t) v gdVdt t1 t1 = the total work done on V(t)between t1 and t2 (8.4) Then Eq. (8) may be rewritten as U(t2 )+ K(t2 )−U(t1 )+ K(t1 ) = Q + W (9) We now consider the two times t1 and t2 correspond to two equilibrium states. If we restrict consideration to initial and final systems at rest, there will be no kinetic energy within V(t) Many at t and t , and the terms K(t ) and K(t ) may be omitted, leaving us with 1 2 1 2 practical ∆U = U(t2 ) − U(t1 ) = Q + W (10) in which ∆U represents the increase in internal energy in the system as one goes from one problems involve equilibrium state to another, because of heat being added to the system and work being fluids contained done on the system. For some fascinating reading about the history of Eq. (10) and the personalities involved, see the text leading up to Eq. 3.37a of Reference 4. within solid For another view of the problem at hand, we can subtract the equation of change for kinetic energy (Eq. 3.3-1 of Reference 3) boundaries, ∂ 1 2 1 2 and so the heat ρv + ∇ ⋅ ρv v = −∇ ⋅ρv +ρ∇ ⋅ v − ∇ ⋅[τ ⋅ v]+ τ : ∇v +ρv ⋅ g (11) ∂t 2 2 and work flows from Eq. (4) above and obtain the equation for internal energy (Eq. 11.2-1 of Reference 3): and their ∂ ρUˆ + ∇ ⋅ρUˆ v = −∇ ⋅ q − p∇ ⋅ v − τ : ∇v (12) ∂t ( ) irreversibilities Then, integrating Eq. (12) over the volume V(t), we apply Eq. (3) to the terms on the left can be and the Gauss divergence theorem to the first term on the right, to give d understood by ∫ ρUˆ dV = − ∫ n⋅ qdA − ∫ p∇ ⋅ vdV− ∫ τ : ∇v dV (13) dt V(t) A(t) V(t) V(t) examining a Integrating over time, we get fluid continuum t 2 t 2 t 2 t 2 ρUˆ dV = − n⋅ qdAdt − p∇ ⋅ v dVdt − τ : ∇v dVdt (14) ∫V(t) ∫ t ∫ A(t) ∫ t ∫ V(t) ∫ t ∫V(t) with diffusion 1 1 1 t1 We denote the first two terms of Eq.