logo
Redian新闻
>未名空间
>GeoSpace - 地球与空间科学
>
谁微分方程厉害的?帮个忙
avatar
谁微分方程厉害的?帮个忙# GeoSpace - 地球与空间科学
k*n
1 楼
我在做个题,
就是INSTANTANIOUS HEATING OF A SEMI-INFINITE PLATE.
一维的.在SCHUBERT AND TURCOTTE的GEODYNAMICS书里有.
半无限的平面.
at time t=0, temperature is T0 everywhere.
at time t>0, temperature at the boundary is Ts at y=0.
at y-> infinite, temperature is still T0.
这个问题是有解了,就是那个什么ERROR FUNCTION.
但我现在的问题是要把Ts弄成TIME DEPENTENT.
假设那是地核的话,外面导热,温度要降低,
我自己算的好象是
ln(Ts) = sqrt (kappa*t/pi).
不过我改的程序好象是给出ln(Ts) proportional to t (not t^(1/2)).
有没有高手看一看啊?
avatar
a*r
2 楼
你没回答我上个帖子的问题啊，是否就是求解(*)和(**) ；
另外，“q可以求导得到”，是用的什么公式？

【在 k***n 的大作中提到】
: 我在做个题,
: 就是INSTANTANIOUS HEATING OF A SEMI-INFINITE PLATE.
: 一维的.在SCHUBERT AND TURCOTTE的GEODYNAMICS书里有.
: 半无限的平面.
: at time t=0, temperature is T0 everywhere.
: at time t>0, temperature at the boundary is Ts at y=0.
: at y-> infinite, temperature is still T0.
: 这个问题是有解了,就是那个什么ERROR FUNCTION.
: 但我现在的问题是要把Ts弄成TIME DEPENTENT.
: 假设那是地核的话,外面导热,温度要降低,
avatar
k*n
3 楼
http://wwwrses.anu.edu.au/~uli/Teaching/Heat/Time/InstHeat.html
http://www.geol.lsu.edu/Faculty/Nunn/gl7065/Chp4_notes.html
找到个连接.
就是这么问题,
唯一改变的是那个Ts现在不是固定的,而是改变的.
因为热量通过传导下去了.
所以这个Ts应该下降.
下降的应该是
rho*Cp*V*dTs=q*A*dt
rho-density
Cp-heat capacity
V-volume
q-heat flux at y=0.
A-area
那个q可以求导得出.
我最后算出来
ln(Ts) 正比于 t^(1/2)
可是模拟结果好象是
ln(Ts) 正比于 t
所以我不确定是不是我修改程序改的不对.
还是我的解有问题.

【在 a*******r 的大作中提到】
: 你没回答我上个帖子的问题啊，是否就是求解(*)和(**) ；
: 另外，“q可以求导得到”，是用的什么公式？
avatar
a*r
4 楼

：rho*Cp*V*dTs=q*A*dt （＊）
：那个q可以求导得出. （＊＊）
就是要求解（＊） （＊＊）两式？ （＊）应有个负号对吧。
要是已经解决了就算了。

【在 k***n 的大作中提到】
: http://wwwrses.anu.edu.au/~uli/Teaching/Heat/Time/InstHeat.html
: http://www.geol.lsu.edu/Faculty/Nunn/gl7065/Chp4_notes.html
: 找到个连接.
: 就是这么问题,
: 唯一改变的是那个Ts现在不是固定的,而是改变的.
: 因为热量通过传导下去了.
: 所以这个Ts应该下降.
: 下降的应该是
: rho*Cp*V*dTs=q*A*dt
: rho-density
avatar
k*n
5 楼
我现在要解TS
恩,应该是TS~ 1-EXP(T)吧,好象
我要仔细CHECK一下,看来.
avatar
k*n
6 楼
Ts～exp(-t)
是负t
这样就对了吧,
时间是0的时候,TS=1.

【在 a*******r 的大作中提到】
:
: ：rho*Cp*V*dTs=q*A*dt （＊）
: ：那个q可以求导得出. （＊＊）
: 就是要求解（＊） （＊＊）两式？ （＊）应有个负号对吧。
: 要是已经解决了就算了。
avatar
a*r
7 楼
你没回答我上个帖子的问题啊，是否就是求解(*)和(**) ；
另外，“q可以求导得到”，是用的什么公式？

【在 k***n 的大作中提到】
: Ts～exp(-t)
: 是负t
: 这样就对了吧,
: 时间是0的时候,TS=1.
avatar
k*n
8 楼
Ts～exp(-t)
是负t
这样就对了吧,
时间是0的时候,TS=1.

【在 a*******r 的大作中提到】
: 你没回答我上个帖子的问题啊，是否就是求解(*)和(**) ；
: 另外，“q可以求导得到”，是用的什么公式？
avatar
p*t
9 楼
把原题扫描或用word写下来发过来吧。BBS上看这种问题头疼。
顺便表扬alligator一下。

【在 k***n 的大作中提到】
: 我在做个题,
: 就是INSTANTANIOUS HEATING OF A SEMI-INFINITE PLATE.
: 一维的.在SCHUBERT AND TURCOTTE的GEODYNAMICS书里有.
: 半无限的平面.
: at time t=0, temperature is T0 everywhere.
: at time t>0, temperature at the boundary is Ts at y=0.
: at y-> infinite, temperature is still T0.
: 这个问题是有解了,就是那个什么ERROR FUNCTION.
: 但我现在的问题是要把Ts弄成TIME DEPENTENT.
: 假设那是地核的话,外面导热,温度要降低,
相关阅读
寻找在美国的有缘人，26岁征婚找朋友。unidentified_title菜鸟问题，关于沙漠地质工程Fresh Master求内推宇宙会不会只是某空间的一小部分？unidentified_titleSEG-Y display question伯克利的地球物理，地震学博士，比EECS的难念么？在NRCS下载soil survey maps中国地震局的局徽会不会改了如果把地球挖穿会怎样？前苏联挖了23年，得知结果后立马停手“地狱之门”科喇超深钻孔种无法解释的新发现海归啦！！即将回国，带不走的东西，回国前大清仓，低价大甩卖GAP一年记录贴：GIS篇Quartic Kernel Density Estimation哪个软件简单好用地球科学，如果想做教授，应该往哪个方向走。新发现！大洪水和大禹治水是真的有偿SCI文章引用人类有足以摧毁地球的杀伤力吗?ESRI在此领域是不是最好的公司？
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章，不代表同意其说法或描述，仅为提供更多信息，也不构成任何建议。文章信息的合法性及真实性由其作者负责，与Redian.news及其运营公司无关。欢迎投稿，如发现稿件侵权，或作者不愿在本网发表文章，请版权拥有者通知本网处理。