ChE 201 - Spring 2005 - Individual Project #2 - 20 points

m-files Due at Start of Class on Thursday 27 January 2005

Write ONE function m-file that uses Matlab’s trapz function to integrate the heat capacity at constant pressure (C_p) between two temperatures. Use “Form 1” of the equation for Cp in Table B.2 of Felder and Rousseau. The function m-file should have the following input arguments: initial temperature, final temperature, and C_p coefficients A, B, C, and D. The function m-file should return the answer in kJ per mole. The function m-file should receive NO input from the keyboard via “input” statements, dialog boxes, etc. All input should be as arguments (the variable names that appear inside the parentheses in the first line of the function m-file). Also, the function m-file should not contain any sort of print statements. (HINT: there is a trapz tutorial on the ChE 201 homepage, and a discussion of numerical integration in Appendix A of Felder & Rousseau.)

Write an m-file that will find the outlet temperature of one stream in a heat exchanger, given the flow rates and inlet temperatures of both streams and the outlet temperature of the other stream. (HINT: this is a good time to use fsolve.) Assume that each stream will contain a single component. The flow rates, the temperatures, and any other required data can be written in assignment statement near the start of the m-file. Please use descriptive variable names. Don’t forget to check for phase changes! (If one of the streams solidifies, you should print an error message and stop.) Use your function m-file to do the integration for you.

Use the m-file you wrote to check whether the outlet temperature on the cold side of the heat exchanger in the unit operations laboratory is what you expect it to be. Be prepared to discuss the results in class.

Submit your m-files.

· Send mail to Dr. Young: valy@bobcat.ent.ohiou.edu.

· Return to ChE 201 home page.

(Last modified on 01/02/05)