Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Polymath Demo Solutions 2002

Download as pdf or txt
Download as pdf or txt
You are on page 1of 40

Sessions 16 and 116

POLYMATH SOLUTIONS TO THE CHEMICAL ENGINEERING DEMONSTRATION PROBLEM SET


Mathematical Software - Sessions 16 and 116*
Michael B. Cutlip, Department of Chemical Engineering, Box U-3222, University of Connecticut, Storrs, CT 06269-3222 (Michael.Cutlip@Uconn.Edu) Mordechai Shacham, Department of Chemical Engineering, Ben-Gurion University of the Negev, Beer Sheva, Israel 84105 (shacham@bgumail.bgu.ac.il) INTRODUCTION These solutions are for a set of representative numerical problems in chemical engineering developed for Session 16 and 116 at the ASEE Chemical Engineering Summer School held at the University of Colorado in Boulder, CO from July 27 - August 1, 2002. The problems in this set are intended to utilize the basic numerical methods in problems which are appropriate to a variety of chemical engineering subject areas. The package used to solve each problem is the POLYMATH Numerical Computation Package Version 5.1 which is widely used in Chemical Engineering. Inexpensive educational site licenses and academic evaluation copies of the software are available from the CACHE Corporation** with information at http://www.che.utexas.edu/cache/ polymath.html. More details on POLYMATH and on-line purchasing options for individual copies are found at http:/ /www.polymath-software.com. The POLYMATH Numerical Computation Package has four companion programs. - SIMULTANEOUS DIFFERENTIAL EQUATIONS - SIMULTANEOUS ALGEBRAIC EQUATIONS - SIMULTANEOUS LINEAR EQUATIONS - CURVE FITTING AND REGRESSION POLYMATH is a proven computational system which has been specifically created for educational use by Mordechai Shacham and Michael B. Cutlip. The latest version runs on all WindowsTM operating systems. The various POLYMATH programs allow the user to apply effective numerical analysis techniques during interactive problem solving on personal computers. Results are presented graphically for easy understanding and for incorporation into papers and reports. Students with a need to solve numerical problems will appreciate the efficiency and speed of problem solution. With POLYMATH, the user is able to focus complete attention to the problem rather that spending valuable time in learning how to use or reuse the software. NOTE - The box around the references to the solution files at the end of each problem can be clicked to launch POLYMATH with the particular problem solution ready to be solved. If this is not available, then the problem solution files can be executed from the Polymath Solution Files directory.
*The Ch. E. Summer School was sponsored by the Chemical Engineering Division of the American Society for Engineering Education. This material is copyrighted by the authors, and permission must be obtained for duplication unless for educational use within departments of chemical engineering. **A non-profit educational corporation supported by most North American chemical engineering departments and many chemical corporation. CACHE stands for computer aides for chemical engineering.

Page PD-1

Page PD-2

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 1D Solution - Steady State Material Balances on a Separation Train (a) The coefficients and the constants in the problem given as Equation Set (D-1) can be directly introduced into the POLYMATH Linear Equation Solver in matrix form as shown

The solution is
[1] [2] [3] [4] [1] [2] [3] [4] D1 B1 D2 B2 = = = = 26.25 17.5 8.75 17.5 + + + + 0.18B1 0.24B1 0.42B1 0.16B1 + + + + 0.15D2 + 0.24B2 = 10.5 0.1D2 + 0.65B2 = 17.5 0.54D2 + 0.1B2 = 28 0.21D2 + 0.01B2 = 14

and the equations are


0.07D1 0.04D1 0.54D1 0.35D1

in which these flow rates have units of mol/min. (b) The overall balances and individual component balances on distillation column #2 given as Equation Set (D-2) can be solved algebraically to give XDx = 0.114, XDs = 0.120, XDt = 0.492 and XDb = 0.274. Similarly, the overall balance and individual component balances on distillation column #3 presented as Equation Set (D-3) yield XBx = 0.210, XBs = 0.4667, XBt = 0.2467 and XBb = 0.0767. (a&b) A combined solution is possible for all equation sets using the POLYMATH Simultaneous Algebraic Equation Solver. This program will solve nonlinear equations and explicit algebraic equations. The linear equations of Equation Set (D-1) of the problem must be entered as nonlinear equations where the function is equal to zero at the solution. Thus, the POLYMATH equations for this combined solution can be expressed as: Nonlinear equations
[1] [2] [3] [4] f(D1) = 0.07*D1+0.18*B1+0.15*D2+0.24*B2-0.15*70 = 0 f(B1) = 0.04*D1+0.24*B1+0.10*D2+0.65*B2-0.25*70 = 0 f(D2) = 0.54*D1+0.42*B2+0.54*D2+0.10*B2-0.40*70 = 0 f(B2) = 0.35*D1+0.16*B1+0.21*D2+0.01*B2-0.20*70 = 0 DD = D1+B1 BB = D2+B2 XDx = (0.07*D1+0.18*B1)/DD XDs = (0.04*D1+0.24*B1)/DD XDt = (0.54*D1+0.42*B1)/DD XDb = (0.35*D1+0.16*B1)/DD

Explicit equations
[1] [2] [3] [4] [5] [6]

Problem 1D Solution

Page PD-3

[7] XBx = (0.15*D2+0.24*B2)/BB [8] XBs = (0.10*D2+0.65*B2)/BB [9] XBt = (0.54*D2+0.10*B2)/BB [10] XBb = (0.21*D2+0.01*B2)/BB

The nonlinear equations (really only linear equations in this example) need to have initial guesses for the solutions entered into POLYMATH. A screen display is given below.

The Polymath solution output file yields the following results:


Variable D1 B1 D2 B2 DD BB XDx XDs XDt XDb XBx XBs XBt XBb Value f(x) Ini Guess 26.25 1.121E-09 20 17.5 -1.965E-09 20 8.75 4.112E-09 20 17.5 -9.819E-10 20 43.75 26.25 0.114 0.12 0.492 0.274 0.21 0.4666667 0.2466667 0.0766667

The POLYMATH problem solution file for this problem is D01.pol. An alternate solution where all three equation sets are solved simultaneously is given in D01alt.pol.

Page PD-4

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 2D Solution - Molar Volume and Compressibility Factor from Van Der Waals Equation Equation (D4) can not be rearranged into a form where V can be explicitly expressed as a function of T and P. However, it can easily be solved numerically using techniques for nonlinear equations. In order to solve Equation (D4) using the POLYMATH Simultaneous Algebraic Equation Solver, it must be rewritten in the form a + ----- ( V b ) RT f(V )= P 2 V PD-(1)

where the solution is obtained when the function is close to zero, f ( V ) 0 . Additional explicit equations and data can be entered into the POLYMATH program in direct algebraic form. The POLYMATH program will reorder these equations as necessary in order to allow sequential calculation. The POLYMATH equation set for this problem is given by Nonlinear equations Explicit equations
[1] [2] [3] [4] [5] [6] [7] [8] [9] [1] f(V) = (P+a/(V^2))*(V-b)-R*T = 0 P = 56 R = 0.08206 T = 450 Tc = 405.5 Pc = 111.3 Pr = P/Pc a = 27*(R^2*Tc^2/Pc)/64 b = R*Tc/(8*Pc) Z = P*V/(R*T)

The POLYMATH input display for this problem is given below.

In order to solve a single nonlinear equation with POLYMATH, an interval for the expected solution variable,V

Problem 2D Solution

Page PD-5

in this case, must be entered into the program. This interval can usually be found by consideration of the physical nature of the problem. (a) For part (a) of this problem, the volume calculated from the ideal gas law as V = 0.66 liter/g-mol can be a basis for specifying the required solution interval. An interval for the expected solution for V can be entered as between 0.4 as the lower limit and 1.0 as the higher limit. The POLYMATH solution, which is given in Figure PD-(1) for T = 450 K and P = 56 atm, yields V = 0.5749 liter/gmol where the compressibility factor is Z = 0.8718.

NLE Solution
Variable V P R T Tc Pc Pr a b Z Value 0.5748919 56 0.08206 450 405.5 111.3 0.5031447 4.1969459 0.0373712 0.8718268 f(x) 6.395E-13 Ini Guess 0.7

Figure PD-1 Plot of f(V) versus V for van der Waals Equation and Solution Summary Table for Problem 2(a).

(b) Solution for the additional pressure values can be accomplished by changing the equations in the POLYMATH program for P and Pr to
Pr=1 P=Pr*Pc

Additionally, the bounds on the molar volume V may need to be altered to obtain an interval where there is a solution. Subsequent program execution for the various Pr s is required.

Page PD-6

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

(c) The calculated molar volumes and compressibility factors are summarized in Table PD-(1). These calculated results indicate that there is a minimum in the compressibility factor Z at approximately Pr = 2. The compressibility factor then starts to increase and reaches Z = 2.783 for Pr = 20.
Table PD-1 Compressibility Factor for Gaseous Ammonia at 450 K

P(atm) 56 111.3 222.6 445.2 1113.0 2226.0

Pr 0.503 1.0 2.0 4.0 10.0 20.0

V .574892 .233509 .0772676 .0606543 .0508753 .046175

Z 0.871827 0.703808 0.465777 0.731261 1.53341 2.78348

A graph can be prepared using the POLYMATH Data Table Program to yield Figure PD-(2).

Figure PD-2 Compressibility Factor versus Reduced Pressure

The POLYMATH problem solution files for this problem are D02a.pol and D02b.pol.

Problem 3D Solution

Page PD-7

Problem 3D Solution - Three Phase Equilibrium - Bubble Point This problem contains simultaneous nonlinear algebraic equations and also explicit algebraic equations. Equations (9) and (10) can be written for each of the two components which result in four nonlinear equations. Equations (11) and (12) can also be written as nonlinear equations. Note that the entry of a nonlinear equation into the POLYMATH Simultaneous Algebraic Equation Solver, it must be rewritten in the form where the right hand side expression should be zero at the solution. f(x) = an expression The argument x in the above equation can be any variable in the problem, and this variable does not need to be in the particular expression. This argument just served to identify the name of a problem variable to the POLYMATH software. There are many ways to arrange the nonlinear and explicit algebraic equations in this problem. The following POLYMATH Report presents the results from one such Problem 3 solution that uses the default solution algorithm and reasonable initial estimates of the solution. NLES Solution
Variable x11 x12 x21 x22 t beta p1 p2 A B gamma11 gamma21 gamma12 gamma22 k11 k21 k12 k22 z1 z2 y1 y2 Value 0.0226982 0.6867476 0.9773018 0.3132524 88.53783 0.7329991 357.05029 498.65881 1.7 0.7 33.36649 1.0046055 1.1028207 3.1342223 15.675678 0.6591518 0.5181085 2.0564573 0.2 0.8 0.3558097 0.6441903 f(x) 8.538E-10 7.962E-10 -5.656E-11 8.117E-10 -1.558E-08 0 Ini Guess 0 1 1 0 100 0.8

NLES Report (fastnewt) Nonlinear equations


[1] [2] [3] [4] [5] [6] f(x11) = x11*(beta+(1-beta)*k11/k12)-z1 = 0 f(x12) = x12*k12-x11*k11 = 0 f(x21) = x21*(beta+(1-beta)*k21/k22)-z2 = 0 f(x22) = x22*k22-x21*k21 = 0 f(t) = x11-y1+x21-y2 = 0 f(beta) = (x11-x12)+(x21-x22) = 0 p1 = 10^(7.62231-1417.9/(191.15+t)) p2 = 10^(8.10765-1750.29/(235+t)) A = 1.7 B = 0.7 gamma11 = 10^(A*x21*x21/((A*x11/B+x21)^2))

Explicit equations
[1] [2] [3] [4] [5]

Page PD-8

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

[6] gamma21 = 10^(B*x11*x11/((x11+B*x21/A)^2)) [7] gamma12 = 10^(A*x22*x22/((A*x12/B+x22)^2)) [8] gamma22 = 10^(B*x12*x12/((x12+B*x22/A)^2)) [9] k11 = gamma11*p1/760 [10] k21 = gamma21*p2/760 [11] k12 = gamma12*p1/760 [12] k22 = gamma22*p2/760 [13] z1 = 0.2 [14] z2 = 0.8 [15] y1 = k11*x11 [16] y2 = k21*x21 Comments [1] f(x11) = x11*(beta+(1-beta)*k11/k12)-z1 Rearrangement of Equation (9) for component i = 1 [2] f(x12) = x12*k12-x11*k11 Rearrangement of Equation (10) for component i = 1 [3] f(x21) = x21*(beta+(1-beta)*k21/k22)-z2 Rearrangement of Equation (9) for component i = 2 [4] f(x22) = x22*k22-x21*k21 Rearrangement of Equation (10) for component i = 2 [5] f(t) = x11-y1+x21-y2 Rearrangement of Equation (11) [6] f(beta) = (x11-x12)+(x21-x22) Rearrangement of Equation (12) [7] p1 = 10^(7.62231-1417.9/(191.15+t)) Antoine equation for component i = 1 [8] p2 = 10^(8.10765-1750.29/(235+t)) Antoine equation for component i = 2 [9] A = 1.7 Numerator constant in Equation (15) [10] B = 0.7 Numerator constant in Equation (16) [11] gamma11 = 10^(A*x21*x21/((A*x11/B+x21)^2)) Equation (15) for component i =1 in liquid phase j = 1 [12] gamma21 = 10^(B*x11*x11/((x11+B*x21/A)^2)) Equation (16) for component i =2 in liquid phase j = 1 [13] gamma12 = 10^(A*x22*x22/((A*x12/B+x22)^2)) Equation (15) for component i =1 in liquid phase j = 2 [14] gamma22 = 10^(B*x12*x12/((x12+B*x22/A)^2)) Equation (16) for component i =1 in liquid phase j = 2 [15] k11 = gamma11*p1/760 Equation (13) for component i = 1 in liquid phase j = 1 [16] k21 = gamma21*p2/760 Equation (13) for component i = 2 in liquid phase j = 1 [17] k12 = gamma12*p1/760 Equation (13) for component i = 1 in liquid phase j = 2 [18] k22 = gamma22*p2/760 Equation (13) for component i = 2 in liquid phase j = 2 [19] z1 = 0.2 Mole fraction of component i = 1 in feed [20] z2 = 0.8 Mole fraction of component i = 2 in feed [21] y1 = k11*x11 Equation (10) for mole fraction of i = 1 in vapor phase [22] y2 = k21*x21 Equation (10) for mole fraction of i = 2 in vapor phase

The POLYMATH problem solution file for this problem is D03.pol.

Problem 4D Solution

Page PD-9

Problem 4D Solution - Terminal Velocity of Falling Particles (a) For conditions similar to those of this problem, the Reynolds number will not exceed 1000 so that only Equations (D-18) and (D-19) need to be applied. The logic which selects the proper equation based on the value of Re can be employed using the if... then... else... statement within the POLYMATH Simultaneous Algebraic Equation Solver. C D = if ( Re < 0.1 ) then ( 24 Re ) else ( 24 ( 1 + 0.14 Re
0.7

))

PD-(2)

Equation (D17) should be rearranged in order to avoid possible division by zero and negative square roots as it is entered into the form of a nonlinear equation for POLYMATH. f ( vt ) = vt ( 3 CD ) 4 g ( p ) Dp The following equation set can be solved by POLYMATH.
Nonlinear equations [1] f(vt) = vt^2*(3*CD*rho)-4*g*(rhop-rho)*Dp = 0 Explicit equations [1] rho = 994.6 [2] g = 9.80665 [3] rhop = 1800 [4] Dp = 0.208e-3 [5] vis = 8.931e-4 [6] Re = Dp*vt*rho/vis [7] CD = if (Re<0.1) then (24/Re) else (24*(1+0.14*Re^0.7)/Re)
2

PD-(3)

Specifying v t, min = 0.0001 and v t, max = 0.05 leads to the results summarized below from the POLYMATH Report.
Variable vt rho g rhop Dp vis Re CD Value 0.0157816 994.6 9.80665 1800 2.08E-04 8.931E-04 3.6556385 8.8426582 f(x) 2.665E-15 Ini Guess 0.02505

(b) The terminal velocity in the centrifugal separator can be calculated by replacing the g in Equation PD-(3) by 30g. Introduction of this change to the equation set gives the following results:
Variable vt rho g rhop Dp vis Re CD Value 0.2060215 994.6 294.1995 1800 2.08E-04 8.931E-04 47.722612 1.5566185 f(x) 2.842E-13 Ini Guess 0.02505

The POLYMATH problem solution files for this problem are D04a.pol and D04b.pol.

Page PD-10

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 5D Solution - Reaction Equilibrium for Multiple Gas Phase Reactions The Equation Set (D-22) can be entered into the POLYMATH Simultaneous Algebraic Equation Solver, but the nonlinear equilibrium expressions must be written as functions which are equal to zero at the solution. A simple transformation of the equilibrium expressions of Equation Set (D-22) to the required functional form yields CCCD - KC 1 f ( C D ) = -------------CA CB CX CY f ( C X ) = -------------- K C 2 CB CC CZ - KC 3 f ( C Z ) = ------------CA CX The above equation set may be difficult to solve because the division by unknowns may make most solution algorithms diverge. Expediting the Solution of Nonlinear Equations An additional simple transformation of the nonlinear function can make many functions much less nonlinear and easier to solve by simply eliminating division by the unknowns. In this case, the Equation Set PD-(4) can be modified to f ( CD ) = CC CD KC 1 CA CB f ( CX ) = CX CY KC 2 CB CC f ( CZ ) = CZ KC 3 CA CX The POLYMATH Report file utilizing Equation Set PD-(5) with the initial conditions for part (a) of C D = C X = C Z = 0 is given below.
NLES Solution Variable CD CX CZ KC1 CY KC2 KC3 CA0 CB0 CC CA CB Value 0.7053344 0.1777924 0.3739766 1.06 0.551769 2.63 5 1.5 1.5 0.1535654 0.420689 0.2428966 f(x) 3.577E-13 3.588E-13 -2.287E-13 Ini Guess 0 0 0

PD-(4)

PD-(5)

NLES Report (safenewt)

Nonlinear equations

[1] f(CD) = CC*CD-KC1*CA*CB = 0 [2] f(CX) = CX*CY-KC2*CB*CC = 0 [3] f(CZ) = CZ-KC3*CA*CX = 0 [1] KC1 = 1.06

Explicit equations

Problem 5D Solution

Page PD-11

[2] [3] [4] [5] [6] [7] [8] [9]

CY = CX+CZ KC2 = 2.63 KC3 = 5 CA0 = 1.5 CB0 = 1.5 CC = CD-CY CA = CA0-CD-CZ CB = CB0-CD-CY

(a), (b) and (c) The POLYMATH solutions are summarized in Table PD-(2) for the three sets of initial conditions. Note that the initial conditions for problem part (a) converged to all positive concentrations. However the initial conditions for parts (b) and (c) converged to some negative values for some of the concentrations. Thus a reality check on Table PD-(2) for physical feasibility reveals that the negative concentrations in parts (b) and (c) are the basis for rejecting these solutions as not representing a physically valid situation.
Table PD-2 POLYMATH Solutions of the Chemical Equilibrium Problem

Variable CD CX CZ CA CB CC CY

Part (a) 0.7053 0.1778 0.3740 0.4207 0.2429 0.1536 0.5518

Part (b) 0.05556 0.5972 1.082 0.3624 -0.2348 -1.624 1.679

Part (c) 1.070 -0.3227 1.131 -0.7006 -0.3779 0.2623 0.8078

Alternate Constrained Solution The POLYMATH Simultaneous Algebraic Equation Solver also offers a selection of algorithms for solving the nonlinear equations. For this example, the constrained algorithm selection allows selected variables to be either (1) be positive or negative at solution, (2) positive at the solution, or (3) positive during iterations and at the solution. The specification of a positive solution is shown below in the input box for variable CD.

In this problem, it is also helpful to express the other gas concentrations as nonlinear equations so that the constraints will allow all gas concentrations to be positive at the problem solution. The POLYMATH input display that

Page PD-12

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

allows this alternate solution is shown below for the initial conditions of part (c).

The results of all three sets of initial conditions with this alternate treatment are equivalent to the original problem solution of part (a) indicating the value of constrained solutions to sets of nonlinear equations.
Variable CD CX CZ CC CA CB KC1 CY KC2 KC3 CA0 CB0 Value 0.7053344 0.1777924 0.3739766 0.1535654 0.420689 0.2428966 1.06 0.551769 2.63 5 1.5 1.5 f(x) 0 0 0 -5.551E-17 0 -2.776E-17 Ini Guess 10 10 10 10 10 10

The POLYMATH problem solution files for this problem are D05a.pol, D05b.pol, D05c.pol, D05b(alt).pol, and D05c(alt).pol.

Problem 6D Solution

Page PD-13

Problem 6D Solution - Vapor Pressure Data Representation by Polynomials and Equations (a) Data Regression with a Polynomial The POLYMATH Polynomial, Multiple Linear and Nonlinear Regression Program can be used to solve this problem by first entering the data in a similar manner to using a spreadsheet. Let us denote the column of temperature data in C as TC and the column of pressure data as P. This POLYMATH Data Table worksheet is reproduced below where the first two columns of data have been entered with appropriate titles.

The Regression tab at the bottom of the POLYMATH Data Table allows a polynomial regression option with the dependent variable column P and the independent variable column TC. This corresponds directly to Equation (D23).

Page PD-14

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

The resulting POLYMATH Report give the details of the polynomial regression including the variance.

POLYMATH Results
Problem 6 - Vapor Pressure Data Representation by Polynomials and Equations

Linear Regression Report


Model: P = a0 + a1*TC + a2*TC^2 + a3*TC^3 + a4*TC^4 Variable a0 a1 a2 a3 a4 Value 24.678757 1.6061958 0.0360443 4.131E-04 3.963E-06 95% confidence 0.7872334 0.0544632 0.0010089 4.005E-05 4.514E-07

General Order of polynomial = 4 Regression including free parameter Number of observations = 10 Statistics R^2 = 0.9999963 R^2adj = 0.9999934 Rmsd = 0.1410532 Variance = 0.3979203

A POLYMATH plot of this resulting polynomial is given below.

Successive polynomial regressions indicate that the polynomial with the minimum variance is the 4th degree.

Problem 6D Solution

Page PD-15

(b) Regression with Clausius-Clapeyron Equation Data regression with the Clausius-Clapeyron expression, Equation (D24), can be accomplished by three additional transformed variables (columns) in the POLYMATH Data Table used for part (a). Additional columns can be defined by the relationships: logP = log(P), TK= T + 273.15, and Trec = 1/TK as indicated below.

A request for linear regression (polynomial with 1st degree) when the dependent variable column is logP and the independent variable column is Trec yields the following plot and numerical results from POLYMATH.

Linear Regression Report


Model: logP = a0 + a1*Trec Variable a0 a1 Value 8.7520167 -2035.3331 95% confidence 0.5423357 153.62853

(c) Riedel Equation Data Correlation

The Riedel equation correlation requires two additional columns for

Page PD-16

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

transformed variables, log T = log ( TK ) and T 2 = TK TK .The Multiple linear option from the Regression tab with Trec, logT, and T2 as the independent variables and logP as the dependent variable can be selected as shown below.

The resulting Polymath Report yields the following results.

POLYMATH Results
Problem D6(c) - Vapor Pressure Data Representation by Polynomials and Equations

Multiple linear regression


Model: logP = a0 + a1*Trec + a2*logT + a3*T2 Variable a0 a1 a2 a3 Value 95% confidence 216.72144 156.41354 -9318.66 4857.1966 -75.748179 58.427175 4.445E-05 5.001E-05

The Graph plot shows that there is a fairly good agreement between the experimental and calculated values of

Problem 6D Solution

Page PD-17

log ( P ) .

The Residual plot showing the error distribution given in is more random than for either the polynomial or the Clapeyron equation.

(d) Nonlinear Regression with the Antoine Equation This expression, Equation (D26), cannot be linearized and so it must be regressed with nonlinear regression option of the POLYMATH Polynomial, Multiple Linear and Nonlinear Regression Program. With this option, the user must supply initial estimates. In this case, it is helpful to use the initial estimates for A and B which were determined in part (b) and use the estimate for C as 273.15. Direct entry of Equation (D26) with the initial guesses for the parameters is accomplished using the Nonlinear tab of the

Page PD-18

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Regression options in the POLYMATH Data Table.

Problem 6D Solution

Page PD-19

The POLYMATH Report of the regression gives more statistical information as shown below. Nonlinear regression (L-M)
Model: logP = A-B/(C+TC) Variable A B C Ini guess 8.752 2035 273.15 Value 5.7673466 677.09401 153.88537 95% confidence 0.1520845 48.159076 5.6870913

Nonlinear regression settings Max # iterations = 64 Precision R^2 R^2adj Rmsd Variance = = = = 0.9996879 0.9995987 0.0047228 3.186E-04 = 10 =3 =1 = 24

General Sample size # Model vars # Indep vars # Iterations

A graph of the residuals from the regressed equation can be used to verify that the errors are approximately randomly distributed.

The POLYMATH problem solution files for this problem are D06a.pol and D06bcd.pol.

Page PD-20

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 7D Solution - Unsteady State Heat Exchange in a Series of Agitated Tanks Equations (D-29) to (D-31), together with the numerical data and initial values given in the problem statement, can be entered into the POLYMATH Simultaneous Differential Equation Solver. The initial startup is from a temperature of 20C in all three tanks, thus this is the appropriate initial condition for each tank temperature. The final value or steady state value can be determined by solving the differential equations to steady state by giving a large time interval for the numerical solution. Alternately one could set the time derivatives to zero, and solve the resulting algebraic equations. In this case, it is easiest just to numerically solve the differential equations to large value of t where steady state is achieved. The POLYMATH differential equation input display with the input of the appropriate equations is shown below where the final value of t is large enough to reach steady state.

The POLYMATH Report provides an overview of the problem solution as given below. Problem D7 - Unsteady State Heat Exchange in a Series of Agitated Tanks Calculated values of the DEQ variables
Variable t T1 T2 T3 W Cp T0 UA Tsteam M initial value 0 20 20 20 100 2 20 10 250 1000 minimal value 0 20 20 20 100 2 20 10 250 1000 maximal value 200 30.952381 41.38322 51.31735 100 2 20 10 250 1000 final value 200 30.952381 41.38322 51.31735 100 2 20 10 250 1000

ODE Report (RKF45)

Problem 7D Solution

Page PD-21

Differential equations as entered by the user [1] d(T1)/d(t) = (W*Cp*(T0-T1)+UA*(Tsteam-T1))/(M*Cp) [2] d(T2)/d(t) = (W*Cp*(T1-T2)+UA*(Tsteam-T2))/(M*Cp) [3] d(T3)/d(t) = (W*Cp*(T2-T3)+UA*(Tsteam-T3))/(M*Cp) Explicit equations as entered by the user [1] W = 100 [2] Cp = 2.0 [3] T0 = 20 [4] UA = 10. [5] Tsteam = 250 [6] M = 1000

The time to reach steady state is usually considered to be the time to reach 99% of the final steady state value for the variable which is increasing and responds the most slowly. For this problem, T3 increases the most slowly, and the steady state value is found to be 51.317C. In POLYMATH, this can be easily done by displaying the output in tabular form for T1, T2, and T3 so that the approach to steady state can accurately be observed. Thus the time must be determined when T3 reaches 0.99(51.317) or 50.804 C. The POLYMATH Data Table output from the Differential Equation Solver is useful in determining this time as illustrated in the table which is partly reproduced below.

Thus the time to reach steady state for T3 is approximately 63.3 minutes as estimated from the above table. The temperatures in the three tanks can be easily plotted using the POLYMATH Graph option from the main

Page PD-22

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Differential Equation display.

Alternate Solution The POLYMATH if... then... else... statement can be used in this solution to calculate the time to reach 99% of the steady state temperature for T3. This involves creating a new differential equation for variable, named ts for example, which follows the time variable until T3 reaches 99% of the steady state temperature, and then the differential change of this variable is set to zero for all larger times. The following statement can be entered into the POLYMATH to provide this differential equation.
d(ts)/d(t) = if (T3<50.804) then (1) else (0)

Part of the POLYMATH Report for this alternate solution is shown below where the ts (steady state time) is calculated to be 83.00
Variable t T1 T2 T3 ts initial value 0 20 20 20 20 minimal value 0 20 20 20 20 maximal value 200 30.952381 41.38322 51.31735 83.004524 final value 200 30.952381 41.38322 51.31735 83.004524

The POLYMATH problem solution files for this problem are D07.pol and D07(alt).pol.

Problem 8D Solution

Page PD-23

Problem 8D Solution - Diffusion with Chemical Reaction in a One Dimensional Slab Solving Higher Order Ordinary Differential Equations POLYMATH, like most mathematical software packages, can solve only systems of first order ordinary differential equations (ODEs). Fortunately, the solution of an n-th order ODE can be accomplished by expressing the equation with a series of simultaneous first order differential equations. This is the approach that is typically used for the integration of higher order ODEs. (a) Equation (D32) is a second order ODE, but it can be converted into a system of first order equations by substituting new variables for the higher order derivatives. In this particular case, a new variable y can be defined which represent the first derivation of CA with respect to z. Thus Equation (D32) can be written as the equation set dC A --------- = y dz k dy ----- = ---------- C A D AB dz

PD-(6)

This set of first order ODEs can be entered into the POLYMATH Simultaneous Differential Equation Solver for solution, but initial conditions for both CA and y are needed. Since the initial condition of y is not known, an iterative method (also referred to as a shooting method) can be used to find the correct initial value for y which will yield the boundary condition given by Equation (D33). Shooting Method-Trial and Error The shooting method is used to achieve the solution of a boundary value problem to one of an iterative solution of an initial value problem. Known initial values are utilized while unknown initial values are optimized to achieve the corresponding boundary conditions. Either trial and error or variable optimization techniques are used to achieve convergence on the boundary conditions. For this problem, a first trial and error value for the initial condition of y, for example y0 = -150, is used to carry out the integration and calculate the error for the boundary condition designated by . Thus the difference between the calculated and desired final value of y at z = L is given by ( y 0 ) = y f, calc y f, desired PD-(7)

Note that for this example, yf ,desired = 0 and thus (y0) = yf,calc only because this desired boundary condition is zero. The equations as entered in the POLYMATH Simultaneous Differential Equation Solver for an initial trial and error solution are shown in Figure PD-(3). The calculation of err in the POLYMATH equation set which corresponds to Equation PD-(7) is only valid at the end of the ODE solution. Repeated reruns of this POLYMATH equation set with different initial conditions for y can be used in a trial and error mode to converge upon the desired boundary condition for y0 where (y0) or err 0. Some results are summarized in Table PD-(3) for various values of y0. The
Table PD-3 Trial Boundary Conditions for Equation Set PD-(6)

y0 (z = 0) yf,calc (z = L) (y0)

-120. 17.23 17.23

-130. 2.764 2.764

-140. -11.70 -11.70

-150. -26.16 -26.16

desired initial value for y0 lies between -130 and -140. This trial and error approach can be continued to obtain a

Page PD-24

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Figure PD-3 POLYMATH Equation Entry

more accurate value for y0, or an optimization technique can be applied. Newtons Method for Boundary Condition Convergence A very useful method for optimizing the proper initial condition is to consider this determination to be a problem in finding the zero of a function. In the notation of this problem, the variable to be optimized is y0 and the objective function is (y0) which is defined by Equation PD-(7). Newtons method, an effective method for optimizing a single variable, can be applied here to minimize the above objective function. According to this method, an improved estimate for y0 can be calculated using the equation y 0, new = y 0 ( y 0 ) ' ( y 0 ) PD-(8)

where ' ( y 0 ) is the derivative of at y = y 0 . The derivative, ' ( y 0 ) , can be estimated using a finite difference approximation ( y0 + y0 ) ( y0 ) ' ( y 0 ) --------------------------------------------- y0 PD-(9)

where y 0 is a small increment in the value of y 0 . It is very convenient that ( y 0 + y 0 ) can be calculated simultaneously with the numerical ODE solution for ( y 0 ) thereby allowing calculation of ' ( y 0 ) from Equation PD-(9) and a new estimate for y 0 from Equation PD-(8). Using = 0.0001 for this example, the POLYMATH equation set for carrying out the first step in Newtons method procedure is given by

Problem 8D Solution

Page PD-25

This set of equations yields the results POLYMATH Report give below where the new estimate for y0 is the final value of the POLYMATH variable ynew or -131.911. Another iteration of Newtons method can be obtained by starting with the new estimate and modifying the initial conditions for y and y1 and the value of y0 in the POLYMATH equation set. The second iteration indicates that the err is approximately 3.e-4 and that ynew is unchanged indicating that convergence has been obtained. For the value of y0 = -131.911, the numerical and analytical solutions are equal to at least six significant digits. Calculated values of the DEQ variables (FIRST ITERATION)
Variable z CA y CA1 y1 k DAB err err1 y0 L delta CAanal derr ynew initial value 0 0.2 -130 0.2 -130.013 0.001 1.2E-09 -130 -130.013 -130 0.001 1.0E-04 0.2 1 -5.227E-11 minimal value 0 0.1404279 -130 0.1404135 -130.013 0.001 1.2E-09 -130 -130.013 -130 0.001 1.0E-04 0.1382726 1 -131.91119 maximal value 0.001 0.2 2.7643829 0.2 2.7455795 0.001 1.2E-09 2.7643829 2.7455795 -130 0.001 1.0E-04 0.2 1.4464177 -5.227E-11 final value 0.001 0.1404606 2.7643829 0.1404457 2.7455795 0.001 1.2E-09 2.7643829 2.7455795 -130 0.001 1.0E-04 0.1382726 1.4464177 -131.91119

The POLYMATH problem solution files for this problem are D08a1.pol, D08a2.pol, and D08a3.pol.

Page PD-26

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 9D Solution - Reversible, Exothermic, Gas Phase Reaction in a Catalytic Reactor Introduction of the given equations and the numerical values of the parameter provided in the problem statement into the POLYMATH Simultaneous Differential Equation Solver is shown below.

(a) The requested plot for the conversion (X), reduced pressure (y) and temperature (T 10-3) can be accomplished by requesting the POLYMATH Table during the problem solution that is partially shown below The scaled

Page PD-27

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

temperature variable is calculated by adding a column to the solution Data Table.as defined below.

The resulting plot indicates that there is a rapid increase in conversion and temperature within the reactor at approximately the midpoint of the catalyst bed. The bed pressure drop is enhanced by the increased temperature and reduced pressure even though the number of moles is decreasing.

(b) The dramatic increase in conversion and temperature is due to the exothermic reaction rapidly accelerating due to the increasing temperature even though the reactant concentration falling. Equilibrium is rapidly achieved after this hot spot is achieved with the temperature and conversion only reducing slightly due to the external heat transfer which tends to slightly cool the reactor as the reacting mixture continues toward the reactor exit.

Page PD-28

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

(c) The concentration profiles reflect the net effects of reaction rate and changes in temperature and pressure within the reactor.

The POLYMATH problem solution files for this problem are D09.pol and D09a.pol.

Page PD-29

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 10D Solution - Dynamics of a Heated Tank with PI Temperature Control This problem requires the solution of Equations (D-46) and (D-48) through (D-53) which can be accomplished with the POLYMATH Simultaneous Differential Equation Solver. The step change in the inlet temperature can be introduced at t = 10 by using the POLYMATH if... then... else... statement to provide the logic for a variable to change at a particular value of t. The generation of a step change at t = 10, for example, is accomplished by the following POLYMATH program statement step=if (t<10) then (0) else (1) (a) Open Loop Performance The step down of 20C in the inlet temperature at t = 10 is implemented below in the equation set for the case where Kc = 0 which gives the open loop response.

Page PD-30

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

A plot of the temperatures T, T0 and Tm as generated by POLYMATH is given in Figure PD-(4) which also verifies the steady state operation for t < 10 min as there is no change in any of the temperature values.

Figure PD-4 Open Loop Response to Step Down in Inlet Feed Temperature at t = 10 min.

(b) Closed Loop Performance The closed loop performance of the PI controller requires the change of Kc from zero in part (a) to the baseline proportional gain of 50. This simple change results in the temperature transients shown in Figure PD-(5).

Figure PD-5 Closed Loop Response to Step Down in Inlet Feed Temperature at t = 10 min

Page PD-31

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

(c) Closed Loop Performance for Kc = 500 The increase of a factor of 10 in the proportional gain from the baseline case gives the unstable result plotted in Figure PD-(6). This is clearly an undesirable result.

Figure PD-6 Closed Loop Response to Step Down in Inlet Feed Temperature at t = 10 min for Kc = 500.

(d) Closed Loop Performance for Only Proportional Control The removal of the integral control action gives the stable result plotted in Figure PD-(7). Note that there is offset from the set point when the system returns to steady state operation. This is always the case for only proportional control, and the use of integral control allows the offset to be eliminated.

Figure PD-7 Closed Loop Response for only Proportional Control.

Page PD-32

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

(e) Closed Loop Performance with Limits on q There are many times in control when limits must be established. In this example, the limits on q can be achieved by a POLYMATH if... then... else... statement which is indicated on the complete POLYMATH Report given below. Problem 10(e) - Closed Loop Response for only Proportional Control with Limits on q Calculated values of the DEQ variables
Variable t T T0 Tm errsum WC Ti rhoVCp taud taum Kc tauI step Tr q qlim dTdt initial value 0 80 80 80 0 500 60 4000 1 5 5000 2 0 80 10000 10000 0 minimal value 0 80 79.42612 79.944186 0 500 60 4000 1 5 5000 2 0 80 -1587.6249 0 -3.9572635 maximal value 200 98.38435 98.382229 92.31428 211.07106 500 60 4000 1 5 5000 2 1 90 6.028E+04 2.6E+04 4 final value 200 89.088128 89.090896 89.093176 211.07106 500 60 4000 1 5 5000 2 1 90 1.453E+04 1.453E+04 -0.0024859

ODE Report (RKF45)


Differential equations as entered by the user [1] d(T)/d(t) = (WC*(Ti-T)+qlim)/rhoVCp [2] d(T0)/d(t) = (T-T0-(taud/2)*dTdt)*2/taud [3] d(Tm)/d(t) = (T0-Tm)/taum [4] d(errsum)/d(t) = Tr-Tm Explicit equations as entered by the user [1] WC = 500 [2] Ti = 60 [3] rhoVCp = 4000 [4] taud = 1 [5] taum = 5 [6] Kc = 5000 [7] tauI = 2 [8] step = if (t<10) then (0) else (1) [9] Tr = 80+step*(10) [10] q = 10000+Kc*(Tr-Tm) [11] qlim = if(q<0)then(0)else(if(q>=2.6*10000)then(2.6*10000)else (q)) [12] dTdt = (WC*(Ti-T)+qlim)/rhoVCp Independent variable variable name : t initial value : 0 final value : 200

Page PD-33

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

The values of q and qlim plotted in Figure PD-(8) indicate that this proportional controller has wide oscillations before settling to a steady state, and the limits imposed on qlim are evident. The corresponding plots of the system temperatures are presented in Figure PD-(9).

Figure PD-8 Closed Loop Response for only Proportional Control with Limits on q.

Figure PD-9 Closed Loop Response for only Proportional Control with Limits on q.

The POLYMATH problem solution files for this problem are D10a.pol, D10b.pol, D10c.pol, D10d.pol, and D10e.pol.

Page PD-34

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 11D Solution - Binary Batch Distillation This problem requires the simultaneous solution of Equation (D54) while the temperature is calculated from the bubble point considerations implicit in Equation (D56). A system of equations comprising of differential and implicit algebraic equations is called differential algebraic or a DAE system. There are several numerical methods for solving DAE systems. Most problem solving software packages including POLYMATH do not have the specific capability for DAE systems. Approach 1 The first approach will be to use the controlled integration technique proposed by Shacham, et al. Using this method, the nonlinear Equation (D56) is rewritten with an error term given by = 1 k1 x1 k2 x2 PD-(10)

where the calculated from this equation provides the basis for keeping the temperature of the distillation at the bubble point. This is accomplished by changing the temperature in proportion to the error in an analogous manner to a proportion controller action. Thus this can be represented by another differential equation dT ------- = Kc dx 2 PD-(11)

where a proper choice of the proportionality constant Kc will keep the error below a desired error tolerance. The calculation of Kc is a simple trial and error procedure for most problems. At the beginning Kc is set to a small value (say Kc = 1), and the system is integrated. If is too large, then Kc must be increased and the integration repeated. This trial and error procedure is continued until becomes smaller than a desired error tolerance throughout the entire integration interval. The temperature at the initial point is not specified in the problem, but it is necessary to start the problem solution at the bubble point of the initial mixture. This separate calculation can be carried out on Equation (D56) for x1 = 0.6 and x2 = 0.4 and the Antoine Equation (D55) using the POLYMATH Simultaneous Algebraic Equation Solver. The POLYMATH Report is as follows. Problem 11a1 - Initial Bubble Point... for Binary Batch Distillation NLE Solution
Variable Tbp xA PA PB xB yA yB Value 95.585087 0.6 1196.2189 485.67158 0.4 0.7869861 0.2130139 f(x) 2.654E-08 Ini Guess 90

NLE Report (safenewt) Nonlinear equations Explicit equations


[1] [2] [3] [4] [1] f(Tbp) = xA*PA+xB*PB-760*1.2 = 0 xA = 0.6 PA = 10^(6.90565-1211.033/(Tbp+220.79)) PB = 10^(6.95464-1344.8/(219.482+Tbp)) xB = 1-xA

Page PD-35

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

[5] yA = xA*PA/(760*1.2) [6] yB = xB*PB/(760*1.2)

where the resulting initial temperature is found to be T 0 = 95.5851 . The system of equations for the batch distillation can then be solved with the POLYMATH Simultaneous Differential Equation Solver for K c = 0.5 10 to yield the following POLYMATH Report. Problem 11 - DAE Equations for Binary Batch Distillation Calculated values of the DEQ variables
Variable x2 T L k2 Kc x1 k1 err initial value 0.4 95.5851 100 0.5325348 5.0E+05 0.6 1.311644 -3.646E-07 minimal value 0.4 95.5851 14.045555 0.5325348 5.0E+05 0.2 1.311644 -3.646E-07 maximal value 0.8 108.56926 100 0.7857526 5.0E+05 0.6 1.8566024 7.798E-05 final value 0.8 108.56926 14.045555 0.7857526 5.0E+05 0.2 1.8566024 7.747E-05
6

ODE Report (RKF45) Differential equations as entered by the user


[1] d(T)/d(x2) = Kc*err [2] d(L)/d(x2) = L/(k2*x2-x2) [1] [2] [3] [4] [5]

Explicit equations as entered by the user

k2 = 10^(6.95464-1344.8/(T+219.482))/(760*1.2) Kc = 0.5e6 x1 = 1-x2 k1 = 10^(6.90565-1211.033/(T+220.79))/(760*1.2) err = (1-k1*x1-k2*x2)

Independent variable
variable name : x2 initial value : 0.4 final value : 0.8

The final values from the Report indicate that 14.05 mol of liquid remain in the column when the concentration of the toluene reaches 80%. During the distillation the temperature increases from 95.6 C to 108.6 C . The error calculated from Equation PD-(10) increases from about 3.6 10 to 7.75 10 still small enough for the solution to be considered as accurate.
7 5

during the numerical solution, but it is

Page PD-36

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Approach 2 A different approach for solving this problem can be used because Equation (D56) can be differentiated with respect to x2 to yield ( k2 k1 ) dT ------- = -------------------------------------------------------------------------------------------------dx 2 B1 B2 ln ( 10 ) x 1 k 1 ----------------------+ x 2 k 2 ----------------------2 2 ( C1 + T ) ( C2 + T ) PD-(12)

Thus Equation PD-(12) can provide the bubble point temperature during the simultaneous integration with Equation (D54). The Report from the POLYMATH Simultaneous Differential Equation Solver is given below. Problem 11- Approach 2 - Differential Equations for Binary Batch Distillation Calculated values of the DEQ variables
Variable x2 L T k2 k1 x1 initial value 0.4 100 95.5851 0.5325348 1.311644 0.6 minimal value 0.4 14.041632 95.5851 0.5325348 1.311644 0.2001465 maximal value 0.8 100 108.57208 0.7856879 1.856466 0.6 final value 0.8 14.041632 108.57208 0.7856879 1.856466 0.2001465

ODE Report (RKF45)


Differential equations as entered by the user [1] d(L)/d(x2) = L/(k2*x2-x2) [2] d(T)/d(x2) = (k2-k1)/(ln(10)*(x1*k1*(-1211.033)/(220.79+T)^2+x2*k2*(1344.8)/(219.482+T)^2)) Explicit equations as entered by the user [1] k2 = 10^(6.95464-1344.8/(T+219.482))/(760*1.2) [2] k1 = 10^(6.90565-1211.033/(T+220.79))/(760*1.2) [3] x1 = 1-x2 Independent variable variable name : x2 initial value : 0.4 final value : 0.8

The above POLYMATH solution, Approach 2, gives essentially the same results as those determined in Approach 1. The POLYMATH problem solution files for this problem are D11a1.pol, D11a2.pol, and D11a3.pol.

Page PD-37

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Problem 12D Solution - Unsteady-state Heat Conduction in a Slab The appropriate equations for this solution are presented in the problem statement. (a) The problem then requires the solution of Equations (D-62), (D-63), and (D-65) which results in nine simultaneous ordinary differential equations and two explicit algebraic equation for the 11 temperature nodes. This set of equations can be entered into the POLYMATH Simultaneous Differential Equation Solver. The equation set as entered into POLYMATH is given below.

The plots of the temperatures in the first four sections, node points 2 5, are shown in Figure PD-(10). The transients in temperatures show an approach to steady state. The numerical results are compared to the hand calculations of a finite difference solution by Geankoplis1(pp. 4713) at the time of 6000 s in Table PD-(4). These results
Table PD-4 Results for Unsteady-State Heat Transfer in One-Dimensional Slab at t = 6000 s

Distance from Slab Surface in m 0 0.2 0.4 0.6 0.8 1.0

Geankoplis1 x = 0.20 m n 1 2 3 4 5 6 T in C 0.0 31.25 58.59 78.13 89.84 93.75

Method of Lines (a) x = 0.10 m n 1 3 5 7 9 11 T in C 0.0 31.71 58.49 77.46 88.22 91.66

Method of Lines (b) x = 0.05 m n 1 5 9 13 17 21 T in C 0.0 31.68 58.47 77.49 88.29 91.72

Page PD-38

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Figure PD-10 Temperature Profiles for Unsteady-State Heat Conduction in a One-Dimensional Slab

indicate that there is general agreement regarding the problem solution, but differences between the temperatures at the nodes increase as the nodes approach the insulated boundary of the slab. (b) The accuracy of the numerical solution can be investigated by doubling the number of sections for the numerical method of lines solution. This just involves adding an additional 10 equations given by the relationship in Equation (D62) and modifying Equation (D65) to calculate T21. The results for this change in the POLYMATH equation set are also summarized in Table PD-(4). Here the numerical solution is only slightly changed from the previous solution in part (a), which gives reassurance to the first choice of 10 sections for this problem. The temperature profiles are virtually unchanged. (c) The calculation of T1 at node 1 is required by the convection boundary condition for this case, and Equation (D68) can be entered into the equation set used in part (a) along with an equation for the ambient temperature T0. This equation set should indicate a somewhat slower response of the temperatures within the slab because of the additional resistance to heat transfer. A comparison with the approximate hand calculations by Geankoplis1 is summarized in Table PD-(5). In this case, the simplified hand calculations give results that have some error relative to the numerical method of lines solutions, which are in good agreement with each other.

Page PD-39

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

Table PD-5 Unsteady-State Heat Transfer with Convection in One-Dimensional Slab at t = 1500 s

Distance from Slab Surface in m 0 0.2 0.4 0.6 0.8 1.0

Geankoplis1 x = 0.20 m n 1 2 3 4 5 6 T in C 64.07 89.07 98.44 100.00 100.00 100.00

Method of Lines (a) x = 0.10 m n 1 3 5 7 9 11 T in C 64.40 88.13 97.38 99.61 99.96 100.00

Method of Lines (b) x = 0.05 m n 1 5 9 13 17 21 T in C 64.99 88.77 97.73 99.72 99.98 100.00

The POLYMATH problem solution files for this problem are D12a.pol, D12b.pol, D12c1.pol, and D12c2.pol.

Page PD-40

POLYMATH SOLUTIONS TO THE DEMONSTRATION PROBLEM SET

REFERENCES
1. Geankoplis, C. J. Transport Processes and Unit Operations, 3rd ed. Englewood Cliffs, NJ: Prentice Hall, 1993.

You might also like