thalfshear-darcy-fluid-pressures.py - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
 (HTM) git clone git://src.adamsgaard.dk/sphere
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
       thalfshear-darcy-fluid-pressures.py (1352B)
       ---
            1 #!/usr/bin/env python
            2 import sphere
            3 import shutil
            4 
            5 import seaborn as sns
            6 #sns.set(style='ticks', palette='Set2')
            7 #sns.set(style='ticks', palette='colorblind')
            8 sns.set(style='white', palette='Set2')
            9 sns.despine() # remove chartjunk
           10 
           11 
           12 sim = sphere.sim(fluid=True)
           13 
           14 def plotpressures(sid):
           15     sim.id(sid)
           16     outformat = 'pdf'
           17     sim.visualize('fluid-pressure',outformat=outformat, figsize=[3.74, 3.14])
           18     filename = sid + '-fluid-pressure.' + outformat
           19     shutil.copyfile(filename, '/home/adc/articles/own/2/graphics/' + filename)
           20 
           21 sids = [
           22 'halfshear-darcy-sigma0=20000.0-k_c=3.5e-15-mu=1.797e-06-velfac=1.0-shear',
           23 'halfshear-darcy-sigma0=20000.0-k_c=3.5e-14-mu=1.797e-06-velfac=1.0-shear',
           24 'halfshear-darcy-sigma0=20000.0-k_c=3.5e-13-mu=1.797e-06-velfac=1.0-shear',
           25 
           26 'halfshear-darcy-sigma0=20000.0-k_c=3.5e-15-mu=1.797e-07-velfac=1.0-shear',
           27 'halfshear-darcy-sigma0=20000.0-k_c=3.5e-15-mu=1.797e-08-velfac=1.0-shear',
           28 
           29 
           30 'halfshear-darcy-sigma0=80000.0-k_c=3.5e-15-mu=1.797e-06-velfac=1.0-shear',
           31 'halfshear-darcy-sigma0=80000.0-k_c=3.5e-14-mu=1.797e-06-velfac=1.0-shear',
           32 'halfshear-darcy-sigma0=80000.0-k_c=3.5e-13-mu=1.797e-06-velfac=1.0-shear',
           33 
           34 'halfshear-darcy-sigma0=80000.0-k_c=3.5e-15-mu=1.797e-07-velfac=1.0-shear',
           35 'halfshear-darcy-sigma0=80000.0-k_c=3.5e-15-mu=1.797e-08-velfac=1.0-shear'#,
           36 ]
           37 
           38 for sid in sids:
           39     plotpressures(sid)