Difference between revisions of "Dynamic models - 2014"
Jump to navigation
Jump to search
Kevin Dunn (talk | contribs) |
Kevin Dunn (talk | contribs) |
||
| Line 9: | Line 9: | ||
| vimeoID6 = 84972209 | | vimeoID6 = 84972209 | ||
| vimeoID7 = 85162288 | | vimeoID7 = 85162288 | ||
| vimeoID8 = | | vimeoID8 = 85354975 | ||
| vimeoID9 = | | vimeoID9 = | ||
| vimeoID10 = | | vimeoID10 = | ||
| Line 39: | Line 39: | ||
| video_download_link7_MP4_size = 294 M | | video_download_link7_MP4_size = 294 M | ||
| video_notes7 = | | video_notes7 = | ||
| video_download_link8_MP4 = | | video_download_link8_MP4 =http://learnche.mcmaster.ca/media/2014-3P4-Class-04A.mp4 | ||
| video_download_link8_MP4_size = M | | video_download_link8_MP4_size = 271 M | ||
| video_notes8 = | | video_notes8 = | ||
| video_download_link9_MP4 = | | video_download_link9_MP4 = | ||
Revision as of 15:52, 29 January 2014
| Class date(s): | 13 to 31 January 2014 | ||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
Readings and preparation for class; video and audio files
- Main reference: Marlin textbook, Chapter 3, Chapter 4, Chapter 5, and Appendix B
- Alternative reference: Seborg textbook, Chapters 2, 3, 4 and 5
| Date | Class number | Video and audio files | Main concepts | Reading (other) | Handout | |
|---|---|---|---|---|---|---|
| 13 January | 02A | Video (361M) | Audio (42M) |
|
|
None |
| 15 January | 02B | Video (364M) | Audio (43M) |
|
None | |
| 17 January | 02C | Video (227M) | Audio (42M) |
|
|
Table of Laplace transforms (from Seborg) |
| 20 January | 03A | Video (367M) | Audio (43M) |
| ||
| 22 January | 03B | Video (305M) | Audio (42M) |
|
None | None |
| 24 January | 03C | Video (317M) | Audio (42M) |
|
None | None |
| 27 January | 04A | Video (294M) | Audio (41M) |
|
None | None |
| 29 and 31 January | 04A |
If you wish to read ahead, we will cover
|
None | None | ||
* The book by Seborg et al. is easily available new or second hand, as it was the prescribed textbook in 2013 (Marlin's book was prescribed in 2012). I will make reference to the chapters from Seborg on the website as well.
Test your understanding before and after class with these resources from Dr. Thomas Marlin. This website also contains full Powerpoint slides for each chapter from his textbook. Use this as a resource to get a different teaching perspective on the same topic. It's quite OK if someone else's approach is a "better fit" for you than my approach.
Computer code: class 02B, 15 January
In a file called cstr_height.m:
function d_depnt__d_indep = cstr_height(indep, depnt)
% Dynamic balance for the CSTR height
% indep: the independent ODE variable, such as time or length or the reactor
% depnt: a VECTOR of dependent variables
%
% Returns:
%
% d(depnt)
% ---------- = a vector of ODEs
% d(indep)
% Assign some variables for convenience of notation: one row per DEPENDENT variable
h = depnt(1);
% Constant and other equations
A = 0.5; % m^2
F_i = 0.8; % m^3/min
R = 15; % min/m^2
F_o = h/R; % m^3/min
% Output from this ODE function must be a COLUMN vector, with n rows
% n = how many ODEs in this system?
n = numel(depnt); d_depnt__d_indep = zeros(n,1);
% Specify every element in the vector below: 1, 2, ... n
d_depnt__d_indep(1) = 1/A * (F_i - F_o);
The call the above model file from the "driver"; you can call this file anything, e.g. ODE_driver.m:
% The independent variable always requires an initial and final value:
indep_start = 0.0; % s
indep_final = 50.0; % s
% Set initial condition(s): for integrating variables (dependent variables)
h_depnt_zero = 4.0; % i.e. h(t=0) = 3.0
IC = [h_depnt_zero];
% Integrate the ODE(s):
[indep, depnt] = ode45(@c02B_linear, [indep_start, indep_final], IC);
% Plot the results:
clf;
plot(indep, depnt(:,1))
grid('on')
hold('on')
xlabel('Time [min]')
ylabel('Tank height')
legend('h')
title('Tank height with time')
% Does it match the analytical equation?
height = 12 - 8.*exp(-indep/7.5);
hold on
plot(indep, height, 'r.')