 |
tread noshear file as input when c_grad_p = 0.1 and sigma0=20e3 - sphere - GPU-… |
 |
git clone git://src.adamsgaard.dk/sphere |
 |
Log |
|
Files |
|
Refs |
|
LICENSE |
|
--- |
|
commit 53ebe460d5a88ab9eb12115197c0ca687da5c0f2 |
|
parent 66b4f94d6377edc0dbd6c2722bc3a7281e8a08e5 |
|
Author: Anders Damsgaard <[email protected]> |
|
Date: Wed, 24 Sep 2014 10:09:12 +0200 |
|
|
|
read noshear file as input when c_grad_p = 0.1 and sigma0=20e3 |
|
|
|
Diffstat: |
|
M python/shear-starter.py | 10 +++++++--- |
|
M python/sphere.py | 4 ++-- |
|
|
|
2 files changed, 9 insertions(+), 5 deletions(-) |
|
--- |
|
diff --git a/python/shear-starter.py b/python/shear-starter.py |
|
t@@ -25,9 +25,13 @@ sim = sphere.sim('diffusivity-sigma0=' + str(sigma0), |
|
fluid=True) |
|
sim.readlast() |
|
|
|
+if sigma0 == 20.0e3 and c_phi == 1.0 and c_grad_p == 0.1: |
|
+ sim.sid = 'shear-sigma0=20000.0-c_phi=1.0-c_grad_p=0.1-hi_mu-lo_visc-hw-no… |
|
+ sim.readlast() |
|
+ |
|
if fluid: |
|
sim.sid = 'shear-sigma0=' + str(sigma0) + '-c_phi=' + str(c_phi) + \ |
|
- '-c_grad_p=' + str(c_grad_p) + '-hi_mu-lo_visc-hw-noshear' |
|
+ '-c_grad_p=' + str(c_grad_p) + '-hi_mu-lo_visc-hw' |
|
else: |
|
sim.sid = 'shear-sigma0=' + str(sigma0) + '-hw' |
|
|
|
t@@ -39,8 +43,8 @@ sim.cleanup() |
|
sim.adjustUpperWall() |
|
sim.zeroKinematics() |
|
|
|
-#sim.shear(1.0/20.0) |
|
-sim.shear(0.0) |
|
+sim.shear(1.0/20.0) |
|
+#sim.shear(0.0) |
|
|
|
if fluid: |
|
#sim.num[2] *= 2 |
|
diff --git a/python/sphere.py b/python/sphere.py |
|
t@@ -3659,7 +3659,7 @@ class sim: |
|
:type disp: str |
|
''' |
|
|
|
- self.writebin() |
|
+ self.writebin(verbose=False) |
|
|
|
nd = '' |
|
if (disp == '2d'): |
|
t@@ -3668,7 +3668,7 @@ class sim: |
|
subprocess.call("cd .. && ./forcechains " + nd + "-f " + outformat \ |
|
+ " -lc " + str(lc) + " -uc " + str(uc) + " input/" + self.sid… |
|
+ ".bin > python/tmp.gp", shell=True) |
|
- subprocess.call("gnuplot tmp.gp && rm tmp.bin && rm tmp.gp", shell=Tru… |
|
+ subprocess.call("gnuplot tmp.gp && rm tmp.gp", shell=True) |
|
|
|
|
|
def forcechainsRose(self, lower_limit=0.25, graphics_format='pdf'): |
