Parametric studies - Revision history

Jpye at 08:13, 29 March 2011

2011-03-29T08:13:06Z

← Older revision		Revision as of 08:13, 29 March 2011
Line 5:		Line 5:

	[[Image:Parameter-study.png]\|thumb\|400\|none\|Results of parametric study using [[ExtPy]] functionality]]		[[Image:Parameter-study.png]\|thumb\|400\|none\|Results of parametric study using [[ExtPy]] functionality]]
	<source lang="~~a4c~~">import extpy;		<source lang="py">import extpy;

	from pylab import *		from pylab import *

Jpye at 08:12, 29 March 2011

2011-03-29T08:12:43Z

← Older revision		Revision as of 08:12, 29 March 2011
Line 7:		Line 7:
	<source lang="a4c">import extpy;		<source lang="a4c">import extpy;

	from pylab import *		from pylab import *
	from solverreporter import *		from solverreporter import *

	def studyplot(self):		def studyplot(self):

	# following is an unfortunate necessity in the current system architecture:		# following is an unfortunate necessity in the current system architecture:
Line 76:		Line 76:
	browser.reporter.reportError("NO PLOTS CREATED")		browser.reporter.reportError("NO PLOTS CREATED")

	extpy.registermethod(studyplot)		extpy.registermethod(studyplot)
	#the above method can be called using "EXTERNAL studyplot(SELF)" in ASCEND.</source>		#the above method can be called using "EXTERNAL studyplot(SELF)" in ASCEND.</source>



	~~[[Category:Experimental]]~~
	[[Category:Documentation]]		[[Category:Documentation]]

Jpye at 08:12, 29 March 2011

2011-03-29T08:12:17Z

← Older revision		Revision as of 08:12, 29 March 2011
Line 4:		Line 4:
	The following code snippet implements something like the old <tt>[[STUDY]]</tt> method from the Tcl/Tk GUI. In this version, it is implemented as a plot, as shown above. At present the python script must be modified each time you want to change the parameters in your study. With some work, we hope to encapsulate this so that the [[STUDY]] can be created with only a METHOD call in the .a4l file.		The following code snippet implements something like the old <tt>[[STUDY]]</tt> method from the Tcl/Tk GUI. In this version, it is implemented as a plot, as shown above. At present the python script must be modified each time you want to change the parameters in your study. With some work, we hope to encapsulate this so that the [[STUDY]] can be created with only a METHOD call in the .a4l file.

	~~<div class="thumb tleft"><div class="thumbinner" style="width:402px;">~~[[Image:Parameter-study.png]~~] <div class="thumbcaption"><div class="magnify">[[File:Parameter-study.png~~\|~~<img src="/skins/common/images/magnify-clip.png" width="15" height="11" alt="" />]]</div>~~Results of parametric study using [[ExtPy]] functionality~~</div></div></div>~~		[[Image:Parameter-study.png]\|thumb\|400\|none\|Results of parametric study using [[ExtPy]] functionality]]
	<source lang="a4c">import extpy;		<source lang="a4c">import extpy;

Line 25:		Line 25:
	figure()		figure()

	#self.inlet.p.setRealValueWithUnits(20,~~"~~bar~~"~~);		#self.inlet.p.setRealValueWithUnits(20,"bar");
	#self.inlet.T.setRealValueWithUnits(100+273.15,~~"~~K~~"~~)		#self.inlet.T.setRealValueWithUnits(100+273.15,"K")

	betavals = [0.2177, 0.25, 0.275,0.3, 0.325,0.35,0.4,0.45,0.5, 0.7,0.9]		betavals = [0.2177, 0.25, 0.275,0.3, 0.325,0.35,0.4,0.45,0.5, 0.7,0.9]
Line 34:		Line 34:
	plots = 0		plots = 0
	for beta in betavals:		for beta in betavals:
	print ~~"~~BETA =~~ ~~%f~~"~~ % beta		print "BETA = %f" % beta
	self.OP.beta.setRealValue(beta)		self.OP.beta.setRealValue(beta)

Line 43:		Line 43:

	mdot = mdot_min + (mdot_max-mdot_min)*i/n		mdot = mdot_min + (mdot_max-mdot_min)*i/n
	self.mdot.setRealValueWithUnits(mdot,~~"~~kg/s~~"~~)		self.mdot.setRealValueWithUnits(mdot,"kg/s")

	try:		try:
	browser.sim.solve(browser.solver,SimpleSolverReporter(browser,~~"~~beta =~~ ~~%f, mdot =~~ ~~%f~~"~~%(beta,mdot)))		browser.sim.solve(browser.solver,SimpleSolverReporter(browser,"beta = %f, mdot = %f"%(beta,mdot)))

	except:		except:
	browser.reporter.reportError('Failed to solve for mdot =~~ ~~%d' % mdot)		browser.reporter.reportError('Failed to solve for mdot = %d' % mdot)

	continue		continue

	xdata.append(mdot)		xdata.append(mdot)
	ydata.append(-self.dp.as(~~"~~bar~~"~~))		ydata.append(-self.dp.as("bar"))

	if xdata:		if xdata:
Line 64:		Line 64:
	if plots:		if plots:

	xlabel(~~"~~Mass flow rate / [kg/s]~~"~~)		xlabel("Mass flow rate / [kg/s]")
	ylabel(~~"~~Pressure drop / [bar]~~"~~)		ylabel("Pressure drop / [bar]")
	title(~~"~~Pressure drop vs flow rate, various orifice ratios~~"~~)		title("Pressure drop vs flow rate, various orifice ratios")

	legend(loc='upper left')		legend(loc='upper left')
Line 74:		Line 74:
	show()		show()
	else:		else:
	browser.reporter.reportError(~~"~~NO PLOTS CREATED~~"~~)		browser.reporter.reportError("NO PLOTS CREATED")

	extpy.registermethod(studyplot)		extpy.registermethod(studyplot)
	#the above method can be called using ~~"~~EXTERNAL studyplot(SELF)~~"~~ in ASCEND.</source>		#the above method can be called using "EXTERNAL studyplot(SELF)" in ASCEND.</source>

Jpye at 09:37, 31 July 2010

2010-07-31T09:37:41Z

← Older revision		Revision as of 09:37, 31 July 2010
Line 1:		Line 1:
	~~<div class="pageexperimental">''This page documents an '''~~experimental''' ~~feature. You can help out by~~ '''~~testing it~~''~~' and <span class="plainlinks">~~[~~http://ascendwiki.cheme.cmu.edu/index.php?title=Parametric_studies&action=edit recording~~]~~</span> your experiences~~.''~~</div>~~		{{experimental}}
			''If you are interested in '''parametric models''', see the page on [[object-oriented modelling]] instead.''

	The following code snippet implements something like the old <tt>[[STUDY]]</tt> method from the Tcl/Tk GUI. In this version, it is implemented as a plot, as shown above. At present the python script must be modified each time you want to change the parameters in your study. With some work, we hope to encapsulate this so that the [[STUDY]] can be created with only a METHOD call in the .a4l file.		The following code snippet implements something like the old <tt>[[STUDY]]</tt> method from the Tcl/Tk GUI. In this version, it is implemented as a plot, as shown above. At present the python script must be modified each time you want to change the parameters in your study. With some work, we hope to encapsulate this so that the [[STUDY]] can be created with only a METHOD call in the .a4l file.

UploadBot: Restored page from Google Cache, uploaded by John Pye

2010-05-13T13:45:11Z

Restored page from Google Cache, uploaded by John Pye

New page

<div class="pageexperimental">''This page documents an '''experimental''' feature. You can help out by '''testing it''' and <span class="plainlinks">[http://ascendwiki.cheme.cmu.edu/index.php?title=Parametric_studies&action=edit recording]</span> your experiences.''</div>

The following code snippet implements something like the old <tt>[[STUDY]]</tt> method from the Tcl/Tk GUI. In this version, it is implemented as a plot, as shown above. At present the python script must be modified each time you want to change the parameters in your study. With some work, we hope to encapsulate this so that the [[STUDY]] can be created with only a METHOD call in the .a4l file.

<div class="thumb tleft"><div class="thumbinner" style="width:402px;">[[Image:Parameter-study.png]] <div class="thumbcaption"><div class="magnify">[[File:Parameter-study.png|<img src="/skins/common/images/magnify-clip.png" width="15" height="11" alt="" />]]</div>Results of parametric study using [[ExtPy]] functionality</div></div></div>
<source lang="a4c">import extpy;

from pylab import *
from solverreporter import *

def studyplot(self):

# following is an unfortunate necessity in the current system architecture:
browser = extpy.getbrowser()

# just check that all is ok
browser.do_solve()

ioff()
n = 20;

mdot_min = 1.7
mdot_max = 12
figure()

#self.inlet.p.setRealValueWithUnits(20,"bar");
#self.inlet.T.setRealValueWithUnits(100+273.15,"K")

betavals = [0.2177, 0.25, 0.275,0.3, 0.325,0.35,0.4,0.45,0.5, 0.7,0.9]

betalines = ['b.-','g.-','r.-' ,'c-','m-','y-','b--','g--','r--','c:','m:','y:']

plots = 0
for beta in betavals:
print "BETA = %f" % beta
self.OP.beta.setRealValue(beta)

xdata = []
ydata = []

for i in range(n):

mdot = mdot_min + (mdot_max-mdot_min)*i/n
self.mdot.setRealValueWithUnits(mdot,"kg/s")

try:
browser.sim.solve(browser.solver,SimpleSolverReporter(browser,"beta = %f, mdot = %f"%(beta,mdot)))

except:
browser.reporter.reportError('Failed to solve for mdot = %d' % mdot)

continue

xdata.append(mdot)
ydata.append(-self.dp.as("bar"))

if xdata:
plot(xdata,ydata,betalines[plots],label='%0.3f'%beta)

hold(1)
plots += 1

if plots:

xlabel("Mass flow rate / [kg/s]")
ylabel("Pressure drop / [bar]")
title("Pressure drop vs flow rate, various orifice ratios")

legend(loc='upper left')
ylim(ymax=4.0)
ion()

show()
else:
browser.reporter.reportError("NO PLOTS CREATED")

extpy.registermethod(studyplot)
#the above method can be called using "EXTERNAL studyplot(SELF)" in ASCEND.</source>

[[Category:Experimental]]
[[Category:Documentation]]