Partial Differential Equations in Two or More Dimensions
5.3 Macroscopic Energy Balance
The general energy balance equation has the form
[pic]=[pic]([pic]+[pic]( [pic]
Let Esys be the total energy (internal + kinetic + potential) of a system, [pic]be the mass flow rate of the system input stream, and [pic]be the mass flow rates of the system output stream, then
[pic][pic] = [pic][pic] ([pic][pic] + [pic] ( [pic] (5.3-1)
where
Usys = system internal energy
Ek,sys = system kinetic energy
Ep,sys = system potential energy
[pic], [pic] = internal energies per unit mass of the system inlet and outlet streams
[pic], [pic] = average velocity of the system inlet and outlet streams
[pic] = rate of heat added to the system
[pic] = rate of work done by the system
The net rate of work done by the system can be written as
[pic]= [pic] + [pic]
where
[pic] = rate of shaft work = rate of work done by the system through a mechanical device (e.g., a pump motor)
[pic] = rate of flow work = rate of work done by the system fluid at the outlet minus rate of work done on the system fluid at the outlet
[pic]
Rate of work = Force([pic] = Force (velocity
Rate of flow work done on the system fluid = PinAinVin = Pin[pic]
Rate of flow work done by the system fluid = PoutAoutVout = Pout[pic]
Eq. (5.3-1) becomes
[pic][pic] = [pic][pic] ([pic][pic]
+ [pic] ( [pic] + Pin[pic]( Pout[pic] (5.3-2)
The internal energy can be combined with the flow work to give the enthalpy
[pic][pic] + Pin[pic] = [pic][pic] = [pic][pic]
In terms of enthalpies [pic] and [pic]
[pic][pic] = [pic][pic] ([pic][pic] + [pic] ( [pic] (5.3-3)
The internal energy and the enthalpy can be related to the heat capacities where
Cp = [pic], and Cv = [pic]
For constant values of Cp and Cv
h = Cp(T - Tref) and u = Cv(T - Tref)
For solid and liquid Cp ( Cv
If the system is at steady state with one inlet and one exit stream [pic] = [pic] = [pic], equation (5.3-3) is simplified to
hout ( hin + g(zout ( zin) + [pic] = [pic] ( [pic] (5.3-4)
Let ( = (“out”) ( (“in”), and q = [pic], w = [pic] be the heat added to the system and work done by the system, respectively, per unit mass flow rate. Equation (5.3-4) becomes
(h + g(z + [pic](V2 = q ( w (5.3-5)
This equation also applies to a system comprising the fluid between any two points along a streamline within a flow field. If these two points are only infinitesimal distance apart, the differential form of the energy equation is obtained
dh + gdz + VdV = (q ( (w (5.3-6)
The d() notation represents a total or “exact” differential and applies to the change in state properties that are determined only by the initial and final states of the properties. The (() notation represents an “inexact” differential and applies to the change in properties that depend upon the path taken from the initial to the final point of the properties. The forms of energy can be classified as either mechanical energy, associated with motion or position, or thermal energy, associated with temperature. Mechanical energy is considered to be an energy form of higher quality than thermal energy since it can be converted directly into useful work. Mechanical energy includes potential energy, kinetic energy, flow work, and shaft work. Internal energy and heat are thermal energy forms that cannot be converted directly into useful work. For systems that involve significant temperature changes, the mechanical energy terms are usually negligible compared with the thermal terms. In such cases the energy balance equation reduces to a “heat or enthalpy balance”, i.e. dh = (q.
Example 5.3-1. ----------------------------------------------------------------------------------
In a residential water heater, water at 60oF flows at a constant 5 GPM into the 100 gallons tank and leaves at 3 GPM. Initially the tank has 10 gallons of 75oF water in it. The tank gas heater heats the tank contents at a constant rate of 800 Btu/min. Assume perfect mixing, determine the temperature of the discharge water after 20 min. of operation.
Water: Cp = Cv = 1 Btu/(lb.oF), density = 62.4 lb/ft3. Unit conversion 1 ft3 = 7.481 gal.
[pic]
Solution ------------------------------------------------------------------------------------------
Step #1: Define the system.
Step #2: Find an equation that contains the temperature of the discharge water.
The energy balance for the system gives the desired equation.
Step #3: Apply the energy balance on the system with the reference temperature Tref = 0oF. Neglect the changes in kinetic and potential energies compared with the changes in thermal energies.
[pic]((VCpT) = (FoCpTo - (FCpT + [pic]
[pic](VT) = FoTo - FT + [pic]
[pic] = 5 - 3 = 2 => V = 10 + 2t
Step #4: Specify the initial condition for the differential equation.
At t = 0, T = 75oF
Step #5: Solve the resulting equation and verify the solution.
V [pic] + T[pic] = FoTo - FT + [pic]
(10 + 2t)[pic] + 2T = (5)(60) - 3T + [pic]
(2t + 10) [pic] = 395.91 - 5T
[pic] = [pic]( -[pic]ln[pic] = [pic]ln[pic]
395.91 - 5T = 20.91[pic] => T = 79.182 - 4.182 [pic]
at t = 20 min., T = 79.1oF
-----------------------------------------------------------------------------------------------------
Equation (5.3-5) and its differential form, equation (5.3-6), are not convenient for solving engineering problems.
(h + g(z + [pic](V2 = q ( w (5.3-5)
dh + gdz + VdV = (q ( (w (5.3-6)
We can use thermodynamics relations to convert the enthalpy term into a form that involves temperature, pressure, and density changes across the system.
du = Tds ( Pd(1/()
dh = du + d(P/() = Tds ( Pd(1/() + d(P/()
dh = Tds ( Pd(1/() + Pd(1/() + [pic] = Tds + [pic] (5.3-7)
For an idealized reversible process in which no energy dissipation occurs, the entropy change arises from heat transfer across the system boundaries
Tds = (q
In any real system, the process is irreversible and there is dissipation of energy, therefore
Tds = (q + (ef = du + Pd(1/() (5.3-8)
dh = (q + (ef + [pic]
In this equation (ef represents the thermal energy generated due to the irreversibility of the system. Substituting dh = (q + (ef + [pic] into the differential energy balance, Eq. (5.3-6), gives
[pic] + gdz + VdV + (ef = ( (w
This equation can be integrated along a streamline from the inlet to the outlet of the system to give
[pic] + g(zo ( zi) + [pic]((oVo2 ( (iVi2) + ef + w = 0 (5.3-9)
where ef = [pic] , w = [pic], and ( = kinetic energy correction factor, ( = 2 for laminar flow, ( = 1 for turbulent flow.
The kinetic correction factor is due to the fact that the velocity profile is not uniform over the cross-sectional area of flow. For uniform flow, the rate of kinetic energy entering a C.V. is given as
[pic] = [pic]VA
The kinetic energy per unit mass flow rate is then
[pic] = [pic]V2
For turbulent flow, the velocity profile is almost flat, therefore
[pic]( [pic]V2 ( ( = 1 for turbulent flow
[pic]
Figure 5.3-1 Laminar velocity profile in a pipe.
The velocity for laminar flow in a pipe is given as
vz = 2V[pic] = 2V (1 ( (2), where ( = [pic]
The rate of kinetic energy entering a C.V. is
[pic] = [pic]2(rdr
Therefore [pic] = (([pic]rdr = ((R2[pic][pic]d[pic]
[pic] = ((R2[pic]V3(1 ( (2)3(d( = 8((R2V3[pic](d(
[pic] = 8((R2V3[pic] = ((R2V3 = (A V3
Therefore [pic]= [pic]V2, and ( = 2 for laminar flow
................
................
In order to avoid copyright disputes, this page is only a partial summary.
To fulfill the demand for quickly locating and searching documents.
It is intelligent file search solution for home and business.
Related searches
- more than me or more than i
- more important or more importantly
- differential equations sample problems
- partial differential equations pdf
- solving differential equations in matlab
- differential equations in matlab
- system of partial differential equations linear
- more than or more than
- writing and solving equations in two variables
- graphing linear equations in two variables calculator
- equations in two variables worksheet
- more important or more importantly grammar