Assignment 7 - 2013

From Engineering Economics and Problem Solving: 4N4
Revision as of 14:32, 22 November 2013 by Kevin Dunn (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
Due date(s): 15 November 2013, at class
Nuvola mimetypes pdf.png (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.

Operating-window-assignment-7-4N4-2013.jpg

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