manus_continuum_granular1_exp

commit 212f9f790c7f67593b0b641dc7be37dd9173817c
parent d082e478c0215f7acc848278f6d1a518896de824
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 30 Jun 2020 13:06:46 +0200

Include pore-skeleton compressibility in skindepth calculation

Diffstat:
M
fig-skin_depth/fig.gp
 | 
3
++-

1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/fig-skin_depth/fig.gp b/fig-skin_depth/fig.gp
@@ -12,6 +12,7 @@ set key bottom right #samplen 0.9
 
 # hydraulic parameters
 phi = 0.25
+alpha = 1e-8
 beta_f = 3.9e-10
 mu_f = 1.787e-3
 #k = 2e-17
@@ -23,7 +24,7 @@ f_min = 0.5*1.0/(24*365)/36000.0
 f_max = 2.0/3600.0
 
 # f: frequency, k: permeability
-skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*f))**(0.5)  # f: 1/s
+skindepth(f,k) = (k/((alpha + phi*mu_f)*beta_f*3.141592654*f))**(0.5)  # f: 1/s
 
 set xlabel "Forcing frequency {/:Italic f} [1/s]"
 set ylabel "Permeability {/:Italic k} [m^2]"