Assignment 7 - 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()