python plot scripts to demo case: dt, iters, cp, kappa

Author: krebeljk

tutorials/demo/fill_pack/constant/plot_cp.py

@@ -0,0 +1,22 @@
+import numpy as np
+from io import StringIO
+import matplotlib.pyplot as plt
+
+# get data
+f = open('cpTable')
+bes = f.read()
+f.close()
+pod = bes.replace("("," ").replace(")", " ")
+tabela = np.loadtxt(StringIO(pod), usecols=(0,2))
+temp = np.add(tabela[:,0], -273.15)
+cp = tabela[:,1]
+
+# plot
+fig, ax = plt.subplots()
+ax.plot(temp, cp)
+ax.set(xlabel='T (degC)', ylabel='c_p (J/kgK)',
+       title='Specific heat at constant pressure')
+ax.grid()
+
+plt.show()
+# fig.savefig("cp.png")

tutorials/demo/fill_pack/constant/plot_kappa.py

@@ -0,0 +1,31 @@
+import numpy as np
+from io import StringIO
+import matplotlib.pyplot as plt
+
+# get data
+f = open('kappaTable')
+bes = f.read()
+f.close()
+pod = bes.replace("("," ").replace(")", " ")
+tabela = np.loadtxt(StringIO(pod), usecols=(0, 2, 4, 6))
+temp = np.add(tabela[:,0], -273.15)
+kappa1 = tabela[:,1]
+kappa2 = tabela[:,2]
+kappa3 = tabela[:,3]
+
+# plot
+fig, ax = plt.subplots()
+
+ax.set(xlabel='T (degC)'
+    ,ylabel='kappa (W/mK)'
+    ,title='Thermal conductivity')
+ax.grid()
+
+ax.plot(temp, kappa1, 'k', label='20 MPa')
+ax.plot(temp, kappa2, 'k--', label='80 MPa')
+ax.plot(temp, kappa3, 'k:', label='120 MPa')
+
+legend = ax.legend(loc='upper right', shadow=True, fontsize='x-large')
+
+plt.show()
+# fig.savefig("kappa.png")

tutorials/demo/fill_pack/plot_dt.py

@@ -0,0 +1,36 @@
+import numpy as np
+from io import StringIO
+import matplotlib.pyplot as plt
+import re
+import sys
+'''
+example usage:
+python plot_dt.py log.openInjMoldSim
+'''
+# get data
+cas = ''
+dcas = ''
+with open(str(sys.argv[1])) as origin_file:
+    data = origin_file.read()
+    for match in re.findall(r'(?m)^Time\s=.*', data):
+        cas = cas + match.split('=')[1] + '\n'
+    for match in re.findall(r'(?m)^deltaT\s=.*', data):
+        dcas = dcas + match.split('=')[1] + '\n'
+
+cas = np.loadtxt(StringIO(cas))
+dcas = np.loadtxt(StringIO(dcas))
+
+# plot
+fig, ax = plt.subplots()
+
+ax.set(xlabel='t [s]'
+      ,ylabel='dt [s]'
+      ,title='Casovni korak')
+ax.grid()
+
+ax.semilogy(cas, dcas, 'k', label='dt')
+
+# legend = ax.legend(loc='upper right', shadow=True, fontsize='x-large')
+
+plt.show()
+# fig.savefig("dt.png")

tutorials/demo/fill_pack/plot_iters.py

@@ -0,0 +1,39 @@
+import numpy as np
+from io import StringIO
+import matplotlib.pyplot as plt
+import re
+import sys
+'''
+example usage:
+python plot_iters.py log.openInjMoldSim
+'''
+# get data
+cas = ''
+iters = ''
+with open(str(sys.argv[1])) as origin_file:
+    data = origin_file.read()
+    for match in re.findall(r'(?m)^Time\s=.*', data):
+        cas = cas + match.split('=')[1] + '\n'
+    for match in re.findall(r'(?m)^PIMPLE:.*in.*iterations', data):
+        iters = iters + match.split(' ')[-2] + '\n'
+
+cas = np.loadtxt(StringIO(cas))
+iters = np.loadtxt(StringIO(iters))
+
+if cas.size == iters.size+1:
+    cas = cas[1:]
+
+# plot
+fig, ax = plt.subplots()
+
+ax.set(xlabel='t [s]'
+      ,ylabel='iters [1]'
+      ,title='Pimple iterations')
+ax.grid()
+
+ax.plot(cas, iters, 'k', label='iters')
+
+# legend = ax.legend(loc='upper right', shadow=True, fontsize='x-large')
+
+plt.show()
+# fig.savefig("iters.png")