 |
tfix lowermost particles so everything doesnt slide horizontally - sphere - GPU… |
 |
git clone git://src.adamsgaard.dk/sphere |
 |
Log |
|
Files |
|
Refs |
|
LICENSE |
|
--- |
|
commit 8443f37adc22d395277cf77eae9fc4ff30424426 |
|
parent b74bec6ccdc448acc13f019a5fdec30d8db66aea |
|
Author: Anders Damsgaard <[email protected]> |
|
Date: Wed, 13 Aug 2014 14:02:52 +0200 |
|
|
|
fix lowermost particles so everything doesnt slide horizontally |
|
|
|
Diffstat: |
|
M python/diffusivity-starter.py | 7 +++++++ |
|
M python/permeability-starter.py | 4 ++++ |
|
|
|
2 files changed, 11 insertions(+), 0 deletions(-) |
|
--- |
|
diff --git a/python/diffusivity-starter.py b/python/diffusivity-starter.py |
|
t@@ -39,13 +39,20 @@ for sigma0_str in sys.argv[4:]: |
|
sim.cleanup() |
|
sim.adjustUpperWall() |
|
sim.zeroKinematics() |
|
+ |
|
sim.consolidate(normal_stress = 10.0e3) |
|
+ |
|
sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic = True) |
|
sim.setFluidBottomNoFlow() |
|
sim.setFluidTopFixedPressure() |
|
sim.setDEMstepsPerCFDstep(10) |
|
sim.setMaxIterations(2e5) |
|
sim.initTemporal(total = 5.0, file_dt = 0.01, epsilon=0.07) |
|
+ |
|
+ # Fix lowermost particles |
|
+ I = numpy.nonzero(sim.x[:,2] < 1.5*dz) |
|
+ sim.fixvel[I] = 1 |
|
+ |
|
sim.run(dry=True) |
|
sim.run(device=0) |
|
#sim.writeVTKall() |
|
diff --git a/python/permeability-starter.py b/python/permeability-starter.py |
|
t@@ -37,6 +37,10 @@ for dp_str in sys.argv[4:]: |
|
z = dz*iz |
|
sim.p_f[:,:,iz] = p_bottom + dp/sim.L[2] * z |
|
|
|
+ # Fix lowermost particles |
|
+ I = numpy.nonzero(sim.x[:,2] < 1.5*dz) |
|
+ sim.fixvel[I] = 1 |
|
+ |
|
sim.setFluidTopFixedPressure() |
|
sim.setFluidBottomFixedPressure() |
|
sim.p_f[:,:,-1] = p_top |
