tadded simulations for N=80kPa, decrease p_f axis lims - 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
---
(DIR) commit 01f8d706e9c77ddcd64baa72996b7758bf4180e5
(DIR) parent 3988d78163cbc3cafac1da217644cd7145716b86
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Wed, 15 Apr 2015 10:32:40 +0200
added simulations for N=80kPa, decrease p_f axis lims
Diffstat:
A python/halfshear-darcy-strength-di… | 49 +++++++++++++++++++++++++++++++
A python/halfshear-darcy-strength-di… | 51 +++++++++++++++++++++++++++++++
M python/halfshear-darcy-strength-di… | 2 +-
M python/halfshear-darcy-strength-di… | 4 ++--
A python/halfshear-darcy-velfac1-N80… | 49 +++++++++++++++++++++++++++++++
A python/halfshear-darcy-velfac1-N80… | 50 +++++++++++++++++++++++++++++++
6 files changed, 202 insertions(+), 3 deletions(-)
---
(DIR) diff --git a/python/halfshear-darcy-strength-dilation-rate-N80-cont.sh b/python/halfshear-darcy-strength-dilation-rate-N80-cont.sh
t@@ -0,0 +1,49 @@
+#!/bin/sh
+#PBS -N hs-d-r-N80
+#PBS -l nodes=1:ppn=3
+#PBS -l walltime=48:00:00
+#PBS -q qfermi
+#PBS -M adc@geo.au.dk
+#PBS -m abe
+
+# Grendel CUDA
+source /com/gcc/4.6.4/load.sh
+CUDAPATH=/com/cuda/5.5.22
+export PATH=$HOME/bin:$PATH
+export PATH=$CUDAPATH/bin:$PATH
+export LD_LIBRARY_PATH=$CUDAPATH/lib64:$CUDAPATH/lib:$LD_LIBRARY_PATH
+
+# Manually installed Python modules
+export PYTHONPATH=$HOME/.local/lib/python:$PYTHONPATH
+export PYTHONPATH=$HOME/.local/lib64/python:$PYTHONPATH
+
+# Manually installed Python
+#export PATH=/home/adc/.local/bin:$PATH
+
+# Shared Python2.7
+PYTHON=/com/python/2.7.6
+export PYTHONPATH=$PYTHON/lib:$PYTHONPATH
+export PATH=$PYTHON/bin:$PATH
+
+echo "`whoami`@`hostname`"
+echo "Start at `date`"
+
+ORIGDIR=/home/adc/code/sphere
+#WORKDIR=/scratch/$PBS_JOBID
+WORKDIR=$ORIGDIR
+
+#cp -r $ORIGDIR/* $WORKDIR
+
+cd $WORKDIR
+nvidia-smi
+rm CMakeCache.txt
+cmake . && make
+cd python
+python continue_sim.py halfshear-darcy-sigma0=80000.0-k_c=3.5e-13-mu=6.11e-07-velfac=1.0-shear 1 0 &
+python continue_sim.py halfshear-darcy-sigma0=80000.0-k_c=3.5e-14-mu=1.797e-07-velfac=1.0-shear 1 1 &
+python continue_sim.py halfshear-darcy-sigma0=80000.0-k_c=3.5e-15-mu=1.797e-08-velfac=1.0-shear 1 2 &
+wait
+
+#cp $WORKDIR/output/* $ORIGDIR/output/
+
+echo "End at `date`"
(DIR) diff --git a/python/halfshear-darcy-strength-dilation-rate-N80.sh b/python/halfshear-darcy-strength-dilation-rate-N80.sh
t@@ -0,0 +1,51 @@
+#!/bin/sh
+#PBS -N hs-d-r-N80
+#PBS -l nodes=1:ppn=3
+#PBS -l walltime=48:00:00
+#PBS -q qfermi
+#PBS -M adc@geo.au.dk
+#PBS -m abe
+
+# Grendel CUDA
+source /com/gcc/4.6.4/load.sh
+CUDAPATH=/com/cuda/5.5.22
+export PATH=$HOME/bin:$PATH
+export PATH=$CUDAPATH/bin:$PATH
+export LD_LIBRARY_PATH=$CUDAPATH/lib64:$CUDAPATH/lib:$LD_LIBRARY_PATH
+
+# Manually installed Python modules
+export PYTHONPATH=$HOME/.local/lib/python:$PYTHONPATH
+export PYTHONPATH=$HOME/.local/lib64/python:$PYTHONPATH
+
+# Manually installed Python
+#export PATH=/home/adc/.local/bin:$PATH
+
+# Shared Python2.7
+PYTHON=/com/python/2.7.6
+export PYTHONPATH=$PYTHON/lib:$PYTHONPATH
+export PATH=$PYTHON/bin:$PATH
+
+echo "`whoami`@`hostname`"
+echo "Start at `date`"
+
+ORIGDIR=/home/adc/code/sphere
+#WORKDIR=/scratch/$PBS_JOBID
+WORKDIR=$ORIGDIR
+
+#cp -r $ORIGDIR/* $WORKDIR
+
+cd $WORKDIR
+nvidia-smi
+rm CMakeCache.txt
+cmake . && make
+cd python
+# $ python shear-starter.py <DEVICE> <FLUID> <C_PHI> <C_GRAD_P> <SIGMA_0>
+#python halfshear-darcy-starter.py 0 1 1.0 3.5e-15 80000.0 1.797e-6 1.0 & from hs-d-v1-N80
+python halfshear-darcy-starter.py 0 1 1.0 3.5e-15 80000.0 6.110e-7 1.0 &
+python halfshear-darcy-starter.py 1 1 1.0 3.5e-15 80000.0 1.797e-7 1.0 &
+python halfshear-darcy-starter.py 2 1 1.0 3.5e-15 80000.0 1.797e-8 1.0 &
+wait
+
+#cp $WORKDIR/output/* $ORIGDIR/output/
+
+echo "End at `date`"
(DIR) diff --git a/python/halfshear-darcy-strength-dilation-rate.py b/python/halfshear-darcy-strength-dilation-rate.py
t@@ -429,7 +429,7 @@ for c in numpy.arange(0,len(mu_f_vals)):
ax1.set_xlim([0.0, 1.0])
ax1.set_ylim([0.0, 1.0])
ax2.set_ylim([0.0, 1.0])
- ax3.set_ylim([-200., 200.])
+ ax3.set_ylim([-150., 150.])
#ax1.set_ylim([0.0, 1.0])
#if pressures:
(DIR) diff --git a/python/halfshear-darcy-strength-dilation.py b/python/halfshear-darcy-strength-dilation.py
t@@ -233,7 +233,7 @@ for c in numpy.arange(0,len(k_c_vals)):
#alpha=alpha, color=color[c])
# smoothed
- ax1.plot(shear_strain[c][1:], smooth(friction[c], smooth_window)[1:], \
+ ax1.plot(shear_strain[c][1:-10], smooth(friction[c], smooth_window)[1:-10],\
label=label, linewidth=1,
alpha=alpha, color=color[c])
t@@ -440,7 +440,7 @@ for c in numpy.arange(0,len(k_c_vals)):
ax1.set_ylim([0.0, 1.0])
ax2.set_ylim([0.0, 1.0])
if fluid:
- ax3.set_ylim([-200., 200.])
+ ax3.set_ylim([-150., 150.])
#ax1.set_ylim([0.0, 1.0])
#if pressures:
(DIR) diff --git a/python/halfshear-darcy-velfac1-N80-cont.sh b/python/halfshear-darcy-velfac1-N80-cont.sh
t@@ -0,0 +1,49 @@
+#!/bin/sh
+#PBS -N hs-d-N80
+#PBS -l nodes=1:ppn=3
+#PBS -l walltime=48:00:00
+#PBS -q qfermi
+#PBS -M adc@geo.au.dk
+#PBS -m abe
+
+# Grendel CUDA
+source /com/gcc/4.6.4/load.sh
+CUDAPATH=/com/cuda/5.5.22
+export PATH=$HOME/bin:$PATH
+export PATH=$CUDAPATH/bin:$PATH
+export LD_LIBRARY_PATH=$CUDAPATH/lib64:$CUDAPATH/lib:$LD_LIBRARY_PATH
+
+# Manually installed Python modules
+export PYTHONPATH=$HOME/.local/lib/python:$PYTHONPATH
+export PYTHONPATH=$HOME/.local/lib64/python:$PYTHONPATH
+
+# Manually installed Python
+#export PATH=/home/adc/.local/bin:$PATH
+
+# Shared Python2.7
+PYTHON=/com/python/2.7.6
+export PYTHONPATH=$PYTHON/lib:$PYTHONPATH
+export PATH=$PYTHON/bin:$PATH
+
+echo "`whoami`@`hostname`"
+echo "Start at `date`"
+
+ORIGDIR=/home/adc/code/sphere
+#WORKDIR=/scratch/$PBS_JOBID
+WORKDIR=$ORIGDIR
+
+#cp -r $ORIGDIR/* $WORKDIR
+
+cd $WORKDIR
+nvidia-smi
+rm CMakeCache.txt
+cmake . && make
+cd python
+python continue_sim.py halfshear-darcy-sigma0=80000.0-k_c=3.5e-13-mu=1.797e-06-velfac=1.0-shear 1 0 &
+python continue_sim.py halfshear-darcy-sigma0=80000.0-k_c=3.5e-14-mu=1.797e-06-velfac=1.0-shear 1 1 &
+python continue_sim.py halfshear-darcy-sigma0=80000.0-k_c=3.5e-15-mu=1.797e-06-velfac=1.0-shear 1 2 &
+wait
+
+#cp $WORKDIR/output/* $ORIGDIR/output/
+
+echo "End at `date`"
(DIR) diff --git a/python/halfshear-darcy-velfac1-N80.sh b/python/halfshear-darcy-velfac1-N80.sh
t@@ -0,0 +1,50 @@
+#!/bin/sh
+#PBS -N hs-d-N80
+#PBS -l nodes=1:ppn=3
+#PBS -l walltime=48:00:00
+#PBS -q qfermi
+#PBS -M adc@geo.au.dk
+#PBS -m abe
+
+# Grendel CUDA
+source /com/gcc/4.6.4/load.sh
+CUDAPATH=/com/cuda/5.5.22
+export PATH=$HOME/bin:$PATH
+export PATH=$CUDAPATH/bin:$PATH
+export LD_LIBRARY_PATH=$CUDAPATH/lib64:$CUDAPATH/lib:$LD_LIBRARY_PATH
+
+# Manually installed Python modules
+export PYTHONPATH=$HOME/.local/lib/python:$PYTHONPATH
+export PYTHONPATH=$HOME/.local/lib64/python:$PYTHONPATH
+
+# Manually installed Python
+#export PATH=/home/adc/.local/bin:$PATH
+
+# Shared Python2.7
+PYTHON=/com/python/2.7.6
+export PYTHONPATH=$PYTHON/lib:$PYTHONPATH
+export PATH=$PYTHON/bin:$PATH
+
+echo "`whoami`@`hostname`"
+echo "Start at `date`"
+
+ORIGDIR=/home/adc/code/sphere
+#WORKDIR=/scratch/$PBS_JOBID
+WORKDIR=$ORIGDIR
+
+#cp -r $ORIGDIR/* $WORKDIR
+
+cd $WORKDIR
+nvidia-smi
+rm CMakeCache.txt
+cmake . && make
+cd python
+# $ python shear-starter.py <DEVICE> <FLUID> <C_PHI> <C_GRAD_P> <SIGMA_0>
+python halfshear-darcy-starter.py 0 1 1.0 3.5e-13 80000.0 1.797e-6 1.0 &
+python halfshear-darcy-starter.py 1 1 1.0 3.5e-14 80000.0 1.797e-6 1.0 &
+python halfshear-darcy-starter.py 2 1 1.0 3.5e-15 80000.0 1.797e-6 1.0 &
+wait
+
+#cp $WORKDIR/output/* $ORIGDIR/output/
+
+echo "End at `date`"