Preliminary 1 Results

The code appears to be running so it is now possible to obtain some prelimiary results for the base set of paramters to investigate, $E_{\text{APPLIED}}$=-0.16, -0.18, -0.20, -0.22, -0.24, -0.26, -0.28, -0.30; $c_{\theta}^{\infty}$=0.006, 0.012


In [12]:
from extremefill2D.fextreme.tools import get_treant_df, get_by_tags, get_by_uuid
import pandas
from extremefill2D.fextreme.plot import vega_plot_treant
from extremefill2D.fextreme import read_json
from toolz.curried import map, pipe
from pprint import pprint
from extremefill2D.fextreme.plot import vega_plot_treants

In [13]:
df = get_treant_df(['prelim1'], path='../../data')

In [14]:
treant = get_by_uuid('dbd3ed44', path='../../data')
print(treant)


/home/wd15/git/extremefill-project/extremefill2D/data/dbd3ed44-449f-4b48-be0a-730c3718f2eb

In [15]:
pprint(pipe(
    'dbd3ed44',
    get_by_uuid(path='../../data'),
    lambda x: x['params.json'].abspath,
    read_json
))


{'CFL': 0.1,
 'Nx': 100,
 'NxBase': 1000,
 'alpha0': 0.5,
 'alpha1': 0.3,
 'alpha_': 1,
 'appliedPotential': -0.3,
 'areaRatio': 0.093,
 'bulkCupric': 1000.0,
 'bulkSuppressor': 0.006,
 'capacitance': 0.3,
 'charge': 2,
 'current': 1e-07,
 'data_frequency': 10,
 'delta': 2.5e-05,
 'deltaRef': 0.03,
 'diffusionCupric': 2.65e-10,
 'diffusionSuppressor': 9.2e-11,
 'dt': 0.01,
 'dtMax': 1.0,
 'dtMin': 5e-08,
 'faradaysConstant': 96485.0,
 'featureDepth': 5.6e-05,
 'gamma': 2.5e-07,
 'gasConstant': 8.314,
 'i0': 50,
 'i1': 0.1,
 'kMinus': 1500000.0,
 'kPlus': 150.0,
 'kappa': 15.26,
 'levelset_update_ncell': 15,
 'omega': 7.1e-06,
 'rboundary': 2.8e-05,
 'rinner': 3.1e-06,
 'router': 8.5e-06,
 'shutdown_deposition_rate': 1e-13,
 'solver_tol': 1e-10,
 'spacing_ratio': 1.1,
 'sweeps': 4,
 'temperature': 298.0,
 'tol': 1e-10,
 'totalSteps': 100000000000000000000,
 'totalTime': 5000.0}

In [17]:
df_sort = df.sort(columns=['bulkSuppressor', 'appliedPotential'])
df_sort


/home/wd15/anaconda/envs/extreme/lib/python3.5/site-packages/ipykernel/__main__.py:1: FutureWarning: sort(columns=....) is deprecated, use sort_values(by=.....)
  if __name__ == '__main__':
Out[17]:
appliedPotential bulkSuppressor tags uuid
11 -0.30 0.006 [prelim1] dbd3ed44
3 -0.28 0.006 [prelim1] d7edd8f4
4 -0.26 0.006 [prelim1] 506cee39
7 -0.24 0.006 [prelim1] 1eb186bc
10 -0.22 0.006 [prelim1] 35a13610
12 -0.20 0.006 [prelim1] d3255296
6 -0.18 0.006 [prelim1] 63124df8
0 -0.16 0.006 [prelim1] 550e1938
13 -0.30 0.012 [prelim1] 17d5767e
2 -0.28 0.012 [prelim1] 95499a89
15 -0.26 0.012 [prelim1] f4e978e0
14 -0.24 0.012 [prelim1] 2431195c
5 -0.22 0.012 [prelim1] 5950f67d
9 -0.20 0.012 [prelim1] 068a2809
8 -0.18 0.012 [prelim1] 41727a23
1 -0.16 0.012 [prelim1] 766e328f

In [23]:
treants_sort = pipe(
    df_sort,
    lambda x: x.uuid,
    map(get_by_uuid(path='../../data')),
    list
)

In [26]:
out = vega_plot_treants(treants_sort[:8])
out.display()



In [27]:
out = vega_plot_treants(treants_sort[8:])
out.display()