Difference between revisions of "Assignment 7 - 2013"
Jump to navigation
Jump to search
Kevin Dunn (talk | contribs) |
Kevin Dunn (talk | contribs) |
||
| Line 12: | Line 12: | ||
A brute force approach to get the operating window for a very simple system: | A brute force approach to get the operating window for a very simple system: | ||
[[Image:operating-window-assignment-7-4N4-2013. | [[Image:operating-window-assignment-7-4N4-2013.jpg]] | ||
<syntaxhighlight lang="python"> | <syntaxhighlight lang="python"> | ||
| Line 30: | Line 30: | ||
pyplot.xlabel('Composition, A1') | pyplot.xlabel('Composition, A1') | ||
pyplot.ylabel('Total flow, F3') | pyplot.ylabel('Total flow, F3') | ||
pyplot.xlim([-0.05, 1.05]) | |||
pyplot.ylim([-5, 95]) | |||
pyplot.savefig('operating-window-assignment-7-4N4-2013.png') | pyplot.savefig('operating-window-assignment-7-4N4-2013.png') | ||
pyplot.show() | pyplot.show() | ||
</syntaxhighlight> | </syntaxhighlight> | ||
Revision as of 00:58, 19 November 2013
| Due date(s): | 15 November 2013, at class |
 (PDF)
|
Assignment questions |
| Other instructions | Tutorial date: 11 November 2013 |
A brute force approach to get the operating window for a very simple system:
import numpy as np
from matplotlib import pyplot
N=100
for Fs in np.arange(0, 60.00001,0.2): # The process can operate between 0 and 100% open for Fs
print(Fs)
for Fa in np.random.uniform(0,30.00001,N): # The process can operate between 0 and 100% open for Fa
flow = Fa + Fs
composition = Fa/(flow + 1E-9)
pyplot.scatter(composition, flow, s=2)
# then plot the data as a dot; every dot represents a valid operating point
pyplot.grid()
pyplot.xlabel('Composition, A1')
pyplot.ylabel('Total flow, F3')
pyplot.xlim([-0.05, 1.05])
pyplot.ylim([-5, 95])
pyplot.savefig('operating-window-assignment-7-4N4-2013.png')
pyplot.show()