tstart error with newline, debug shear experiment - sphere - GPU-based 3D discr… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit 31cb2c6665a9cca4fde1b0da0f1f45831540ac2f | |
parent b36c0fdbdc4167ec1b64d3da69799b9acbfc4229 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Thu, 6 Nov 2014 10:54:07 +0100 | |
start error with newline, debug shear experiment | |
Diffstat: | |
M python/halfshear-darcy-starter.py | 5 +++-- | |
M src/darcy.cpp | 3 ++- | |
2 files changed, 5 insertions(+), 3 deletions(-) | |
--- | |
diff --git a/python/halfshear-darcy-starter.py b/python/halfshear-darcy-starter… | |
t@@ -61,7 +61,8 @@ sim.setDampingTangential(0.0) | |
#sim.deleteAllParticles() | |
sim.fixvel[:] = -1.0 | |
-sim.initTemporal(total = 20.0, file_dt = 0.01, epsilon=0.07) | |
+#sim.initTemporal(total = 20.0, file_dt = 0.01, epsilon=0.07) | |
+sim.initTemporal(total = 1.0e-4, file_dt = 1.0e-5, epsilon=0.07) | |
# Fix lowermost particles | |
#dz = sim.L[2]/sim.num[2] | |
t@@ -70,6 +71,6 @@ sim.initTemporal(total = 20.0, file_dt = 0.01, epsilon=0.07) | |
sim.run(dry=True) | |
sim.run(device=device) | |
-#sim.writeVTKall() | |
+sim.writeVTKall() | |
#sim.visualize('walls') | |
#sim.visualize('fluid-pressure') | |
diff --git a/src/darcy.cpp b/src/darcy.cpp | |
t@@ -125,7 +125,8 @@ void DEM::checkDarcyStability() | |
//if (darcy.mu*time.dt/(dmin*dmin) > 0.5) { | |
if (time.dt >= 1.0/(2.0*alpha_max) * | |
1.0/(1.0/(dx*dx) + 1.0/(dy*dy) + 1.0/(dz*dz))) { | |
- std::cerr << "Error: The time step is too large to ensure stability in… | |
+ std::cerr | |
+ << "\nError: The time step is too large to ensure stability in " | |
"the diffusive term of the fluid momentum equation.\n" | |
"Increase the viscosity, decrease the time step, and/or increase " | |
"the fluid grid cell size." << std::endl; |