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 is to "operate" (simulate) the process at a variety of combinations of points. This is what the code does below.

More challenging simulations (e.g. in Aspen) would be done by picking only a few selected points to confirm your expectations on the process, and to simulate regions in the window where you are unable to predict the result.

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