Difference between revisions of "Assignment 7 - 2013"

From Engineering Economics and Problem Solving: 4N4
Jump to navigation Jump to search
 
Line 23: Line 23:
N=100
N=100
for Fs in np.arange(0, 60.00001,0.2): # The process can operate between 0 and 100% open for Fs
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
for Fa in np.random.uniform(0,30.00001,N): # The process can operate between 0 and 100% open for Fa
        flow = Fa + Fs
flow = Fa + Fs
        composition = Fa/(flow + 1E-9)
composition = Fa/(flow + 1E-9)
        pyplot.scatter(composition, flow, s=2)   
pyplot.scatter(composition, flow, s=2)   


# then plot the data as a dot; every dot represents a valid operating point
# then plot the data as a dot; every dot represents a valid operating point

Latest revision as of 14:32, 22 November 2013

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