tdisplay values at a porosity of 0.3 - sphere - GPU-based 3D discrete element m… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit 48431f5c15a7f588b18bda59e3226a21515da310 | |
parent e00575a0c51ed0d98599781875be8f9e83f66bc1 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Wed, 12 Nov 2014 14:16:08 +0100 | |
display values at a porosity of 0.3 | |
Diffstat: | |
M python/sphere.py | 4 ++++ | |
1 file changed, 4 insertions(+), 0 deletions(-) | |
--- | |
diff --git a/python/sphere.py b/python/sphere.py | |
t@@ -3138,11 +3138,15 @@ class sim: | |
self.k_c[0] = k_c | |
phi = numpy.array([0.1, 0.9]) | |
k = self.k_c * phi**3/(1.0 - phi**2) | |
+ k_phi30 = self.k_c * 0.3**3/(1.0 - 0.3**2) | |
K = k * self.rho*numpy.abs(self.g[2])/self.mu | |
+ K_phi30 = k_phi30 = self.rho*numpy.abs(self.g[2])/self.mu | |
print('Hydraulic permeability limits for porosity phi = [0.1, 0.9]… | |
print('\tk = ' + str(k) + ' m*m') | |
+ print('\tk(phi = 0.3) = ' + str(k_phi30) + ' m*m') | |
print('Hydraulic conductivity limits for porosity phi = [0.1, 0.9]… | |
print('\tK = ' + str(K) + ' m/s') | |
+ print('\tK(phi = 0.3) = ' + str(K_phi30) + ' m/s') | |
else: | |
raise Exception('setPermeabilityPrefactor() only relevant for the ' | |
+ 'Darcy solver (cfd_solver = 1)') |