timprove stress plot w different permeabilities - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
HTML git clone git://src.adamsgaard.dk/sphere
DIR Log
DIR Files
DIR Refs
DIR LICENSE
---
DIR commit 669749b24e19b7624a788903e16ffff634579f6d
DIR parent 22f5eb64dee87919b6ab90d4b1e8d5419e3892e6
HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 14 Apr 2015 19:47:43 +0200
improve stress plot w different permeabilities
Diffstat:
M python/halfshear-darcy-strength-di… | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
---
DIR diff --git a/python/halfshear-darcy-strength-dilation-rate.py b/python/halfshear-darcy-strength-dilation-rate.py
t@@ -35,8 +35,7 @@ k_c = 3.5e-15
# 5.0e-8 results present
#mu_f_vals = [1.797e-06, 1.204e-06, 5.0e-8, 1.797e-08]
#mu_f_vals = [1.797e-06, 1.204e-06, 3.594e-07, 1.797e-08]
-#mu_f_vals = [1.797e-06, 1.797-07, 1.797e-08]
-mu_f_vals = [1.797e-06, 1.797-07, 1.797e-08]
+mu_f_vals = [1.797e-06, 1.797e-07, 1.797e-08]
#mu_f_vals = [1.797e-06, 1.204e-06, 1.797e-08]
#mu_f_vals = [1.797e-06, 3.594e-07, 1.797e-08]
velfac = 1.0
t@@ -283,7 +282,7 @@ for c in numpy.arange(0,len(mu_f_vals)):
if zflow:
ax3.set_ylabel('$\\boldsymbol{v}_\\text{f}^z h$ [ms$^{-1}$]')
if pressures:
- ax3.set_ylabel('Mean fluid pressure $\\bar{p}_\\text{f}$ [kPa]')
+ ax3.set_ylabel('Fluid pressure $\\bar{p}_\\text{f}$ [kPa]')
#ax4.set_ylabel('Particle contact force $||\\boldsymbol{f}_\\text{p}||$ [N]')
#ax1.set_xlim([200,300])