示例#1
0
    ':',
])
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',
示例#2
0
#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)):
示例#3
0
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
示例#4
0
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
示例#5
0
    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