':', ]) colors = np.array([]) names = np.array(['01_dt-1', '02_dt-2', '03_dt-5', '04_dt-10']) #names = np.array(['01_dt1_p005_subgrid', '02_dt2_p005_subgrid', # '03_dt5_p005_subgrid', '04_dt10_p005_subgrid']) #names = np.array(['05_dt1_p005_Dborder', '06_dt2_p005_Dborder', # '07_dt5_p005_Dborder', '08_dt10_p005_Dborder']) ts = np.array([ 8.334e-06, 1.667e-05, 4.167e-05, 8.334e-05, 8.334e-06, 1.667e-05, 4.167e-05, 8.334e-05 ]) results = {} for nn in names: path = root_dir + '\\results\\output\\01_time_step\\' + nn + '\\' results[nn] = fn.load_obj(path + nn + '_results') for i in range(0, len(names)): temp = np.array(results[names[i]]['time']) temp *= scale results[names[i]]['time'] = temp.tolist() label = names #%% CH DISSOLUTION titles = [ 'Portlandite', 'Calcite', 'Calcium', 'Carbon', 'Average pH', 'Input C', 'Porosity' ] comp = ['portlandite', 'calcite', 'Ca', 'C', 'pH', 'C (1, 0)', 'avg_poros'] suffix = [ '_portlandite', '_calcite', '_calcium', '_carbon', '_average ph',
#names = np.array(['01_p005_ll0_100s', '02_p005_ll4_100s', # '03_p05_ll0_100s', '04_p05_ll4_100s']) names = np.array([ '01_p005_ll0_1000s', '01_p005_ll1_1000s', '02_p005_ll4_1000s', '03_p05_ll0_1000s', '03_p05_ll1_1000s', '04_p05_ll4_1000s' ]) porosity = np.array(['0.05', '0.05', '0.05', '0.5', '0.5', '0.5']) liqlayer = np.array(['0', '1', '4', '0', '1', '4']) label = np.array([ 'p0.05; l0', 'p0.05; l1', 'p0.05; l4', 'p0.5; l0', 'p0.5; l1', 'p0.5; l4' ]) results = {} for i in range(0, len(names)): path = paths[i] + names[i] + '\\' results[names[i]] = fn.load_obj(path + names[i] + '_results') for i in range(0, len(names)): temp = np.array(results[names[i]]['time']) temp *= scale results[names[i]]['time'] = temp.tolist() #%% PROPERTIES titles = ['Portlandite', 'Calcite', 'Calcium', 'Carbon', 'Input C', 'Porosity'] comp = ['portlandite', 'calcite', 'Ca', 'C', 'C (1, 0)', 'avg_poros'] suffix = [ '_portlandite', '_calcite', '_calcium', '_carbon', '_input_c', '_poros' ] for k in range(0, len(comp)): plt.figure(figsize=(8, 4)) for i in range(0, len(names)):
fname = 'dc' fpath = root_dir+'\\results\\output\\00_examples\\'#'\\results\\output\\00_default\\' fn.make_output_dir(fpath) name = '01_example_default' path = root_dir+'\\results\\output\\00_examples\\' + name + '\\' #name = '01_ie01_p05' #path = root_dir+'\\results\\output\\07_internal_energy\\' + name + '\\' #path = root_dir+'\\results\\output\\13_validation\\' + name + '\\' #name = '04_ccD13_Di' #name = '05_ccD13_Di_PSCS' linetype = '-' scale = 50 results = {} results = fn.load_obj(path + name +'_results') for n in ['time', 'portlandite', 'calcite']: temp = np.array(results[n]) temp *= scale results[n]= temp.tolist() concCA = np.load(path +'Ca.npy') concC = np.load(path +'C.npy') ph = np.load(path +'pH.npy') de = np.load(path +'De.npy') poros = np.load(path +'poros.npy') si = np.load(path +'SI.npy') phaseCC = np.load(path +'CC.npy') phaseCH = np.load(path +'CH.npy') #%% MAIN PROPERTIES
import sys,os root_dir = os.path.dirname(os.path.dirname(os.path.dirname(os.path.dirname(os.path.abspath(__file__))))) src_dir = os.path.dirname(os.path.dirname(os.path.dirname(os.path.abspath(__file__)))) sys.path.append(root_dir) sys.path.append(src_dir) import matplotlib.pylab as plt import numpy as np np.set_printoptions(precision=5, threshold=np.inf) import misc_func as fn #%% UPLOAD #name = '02_ccD12_long_prev' #name = '00_ccD11_prev' name = '04_ccD13_Di' path = root_dir+'\\results\\output\\13_validation\\' + name + '\\' results = fn.load_obj(path + '04_ccD13_Di' +'_results') scale = 100 #%% mm_CH = 74.093 mm_CC = 100.086 alpha = np.array([0.,0.3,0.55,1.05,1.25,1.55]) time = np.array([0.,2.,4.5,10.,12.,20.]) #days mass0 = 22.977 *1e-3 #[g] initial mass of CH crystal dmass = alpha*mass0/100 #delta mass increase dmm = mm_CC-mm_CH #[g/mol] molar mass difference xmol = dmass/dmm #[mol] mol of CH transformed to CC sa = 29.19e+6
os.path.dirname(os.path.dirname(os.path.dirname( os.path.abspath(__file__))))) src_dir = os.path.dirname( os.path.dirname(os.path.dirname(os.path.abspath(__file__)))) sys.path.append(root_dir) sys.path.append(src_dir) import matplotlib.pylab as plt import numpy as np np.set_printoptions(precision=5, threshold=np.inf) import misc_func as fn #%% UPLOAD #name = '02_ccD12_long_prev' name = '01_18days_D13' path = root_dir + '\\results\\output\\09_long\\' + name + '\\' results = fn.load_obj(path + '00_ccD11' + '_results') scale = 100 #%% mm_CH = 74.093 mm_CC = 100.086 alpha = np.array([0., 0.3, 0.55, 1.05, 1.25, 1.55]) time = np.array([0., 2., 4.5, 10., 12., 20.]) #days mass0 = 22.977 * 1e-3 #[g] initial mass of CH crystal dmass = alpha * mass0 / 100 #delta mass increase dmm = mm_CC - mm_CH #[g/mol] molar mass difference xmol = dmass / dmm #[mol] mol of CH transformed to CC sa = 29.19e+6 depth_sa = xmol * 0.0331 * 1e+15 / sa