2
\ begingroup美元

我试图重建图如图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]')
\ endgroup美元

    1回答1

    2
    \ begingroup美元

    在Matlab中,三角函数与角弧度,例如描述功能:

    输入角弧度,指定为一个标量、向量,矩阵,或多维数组。

    所以你要么保持代码但转换10°角弧度(改变α= 0.174533),或者您可以使用信德函数,它也计算正弦但需要角度在度作为输入。

    然后,你也有一个10倍误差的定义一个美元:10美元^ {-16}$应该注意的1 e-16。修复第一个图(速度剖面):

    在这里输入图像描述

    第二个看起来仍然奇怪……你也有一个错误的定义q $ ^{\补}_ {geo} $:单位是千瓦,所以它应该转换为q_geo = 0.05;。这个修正方程的第一行:如果你把其余的评论,它再现了线性,虚线,温度为-7.619°Cz美元m = 0。所以这个问题必须在第二部分,虽然我不能找出。

    \ endgroup美元
    5
    • 1
      \ begingroup美元 非常感谢您的帮助! !非常感激! \ endgroup美元
      - - - - - -k12345
      2月28日在13:26
    • 1
      \ begingroup美元 欢迎你!我找不到这个问题的第二部分方程6.28。我并不是说源是错误的,但它可以是错的。有一个完整的书勘误表纠正几个方程(这是完全正常的,我从来没见过一本书没有任何错误)。方程6.28不包括在勘误表,但它可能只是没有人发现/报告错误(如果有的话)。继续工作,如果你相信什么,你可以写信给作者。他们通常很高兴有人帮助他们解决他们的工作。 \ endgroup美元
      - - - - - -让-玛丽•Prival
      2月28日在14:50
    • 1
      \ begingroup美元 我非常欣赏这个! !谢谢你! \ endgroup美元
      - - - - - -k12345
      2月28日在阿联酋
    • 1
      \ begingroup美元 我发现这个问题在6.28 !只是更新了好玩。6.28中使用的值必须是美元Pa ^ {3} s ^{1} $为了单位匹配 \ endgroup美元
      - - - - - -k12345
      3月3日在十四11
    • 1
      \ begingroup美元 很有道理,好赶上!你可以写成一个答案,它是可以回答自己的问题,明白了www.hoelymoley.com/help/self-answer。它将有助于未来读者遇到同样的问题(虽然我承认它不太可能在你的情况下)。 \ endgroup美元
      - - - - - -让-玛丽•Prival
      在3月3日15:55

    你的答案

    通过点击“发布你的答案”,你同意我们服务条款,隐私政策饼干的政策

    不是你要找的答案?浏览其他问题标记问你自己的问题