我试图重建图如图6.2源(88页)。他们得到这些情节通过使用下面的速度和温度的解决方案(包括获得解决基本平行板坡问题)。
$ v_x = \压裂{2(\ρg sinα(\))^ {n}} {n + 1} (H ^ {n + 1} - (H-z) ^ {n + 1})美元
$ T = T_s + \压裂{问^{\补}_ {geo}} {\ kappa} (H-z) + \压裂{2啊^ {n + 3}(\ρg sinα(\))^ {n + 1}}{\卡帕(n + 2)}[1 - \压裂{z} {H} \压裂{1}{n + 3}(\压裂{H-z} {H}) ^ {n + 3}]美元
他们说为了创建图6.2,他们使用以下值。
- H = 100美元
- $ \α= 10 ^{\保监会}$
- T_s美元= -10 ^{\保监会}加元
- q $ ^{\补}_ {geo} = 50 W m ^ {2} $
- n = 3美元
- 一美元= 10 ^{-16}一个^ {1}Pa ^{3} = 10 ^{-16}一个^{1}公斤^ {3}m ^ {3} s ^ {6} $
- $ \ρ= 910公斤^ {3}$
- \ K = 2.1 W美元^ K ^ {1} {1} = 2.1 \ mbox{公斤}\ mbox {m} \ mbox{年代}^ {3}K ^ {1} $
- $ g = 9.81 m s ^ {2} $
我注意到,当我把这些值代入上面的公式我不得到相同的结果图。我看过其他来源与相似的图形,例如图2显示v_x美元在这里,所以我知道他们所做的是正确的,我觉得我失踪。
我仔细检查单位,他们的正确v_x美元方程。为元新台币方程我注意到我需要转换T_s美元先开尔文(即增加273)代入方程之前,然后我需要转换的最终结果元新台币方程回摄氏度(即减去273)。这还不给我一个匹配的图是什么图6.2中提供。
为什么这些数字不匹配时的值代入方程?我已经通过数学和编码一群,他们看起来好所以我觉得我的问题是源于一种误解的冰川动力学。有一些技巧,我失踪吗?谢谢!
这是原始图6.2
这就是我得到的
这是Matlab代码我曾经让我的结果,帮助。
H = 100;α= 10;= 10 e-16;ρ= 910;g = 9.81;k = 2.1;n = 3;q_geo = 50;Ts = -10;图()z_grid1 = linspace (0100401); v_x = (2.*A.*(1/(n+1)).*(rho.*g.*sin(alpha)).^n)*(H^(n+1) -(H-z_grid1).^(n+1)); %equation 6.13 plot(v_x,z_grid1) title('v_x graph') xlabel('velocity [m/a]') ylabel('Height above the base z [m]') figure() z_grid2 = linspace(0,100,401); T = (Ts+273) + (q_geo./kappa).*(H-z_grid2)+(2.*A.*H.^(n+3)).* (((rho.*g.*sin(alpha))^(n+1))./ (kappa.*(n+2))).*(1-(z_grid2 ./H) - (1/(n+3)).*((H-z_grid2)./H).^(n+3)); %equation 6.28 plot(T-273,z_grid2) title('T graph') xlabel('T [Degrees Celcius]') ylabel('Height above the base z [m]')
