24小时热门版块排行榜    

投票标题: 偏微分方程,常微分方程,应该还需软件协助,我已经搞不懂了。。。。算是编程离散求解 (单选)

  • merry
    1 (33.3%)
  • me
    2 (66.7%)
查看: 1473  |  回复: 8

ahui728

木虫 (文坛精英)

叫兽

[求助] 偏微分方程,常微分方程,应该还需软件协助,我已经搞不懂了。。。。算是编程离散求解

偏微分方程,常微分方程,应该还需软件协助,我已经搞不懂了。。。。算是编程离散求解

问题描述就是上面的情况,有人说有解析解?我不确定。

只要问题解决,无论解析的还是离散的。都行,都给金币。亦可另行约定其他报酬。

软件最好用 MATLAB、 maple 或 mathcad 吧,Fortran也行,但不太懂。。。
我有100双袜子,可以送你一双。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

1 首先问一个问题,eq 2 中为什么r和t相关?r不是空间变量吗?这和eq1 里T=T(r,t)作为r和t的函数的r有什么不同?
2 eq2 应该是状态方程,eq1 是能量方程(热传导方程), eq2 是什么方程?
2楼2014-12-15 23:11:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
3楼2014-12-16 08:10:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ahui728

木虫 (文坛精英)

叫兽

引用回帖:
2楼: Originally posted by pippi6 at 2014-12-15 23:11:38
1 首先问一个问题,eq 2 中为什么r和t相关?r不是空间变量吗?这和eq1 里T=T(r,t)作为r和t的函数的r有什么不同?
2 eq2 应该是状态方程,eq1 是能量方程(热传导方程), eq2 是什么方程?

我想还是先说一下物理现象。
球形水滴受热蒸发,主要考虑强烈的热辐射作用。开始时,水滴的温度和外界的温度相差很小。
eq2给出水滴半径随时间变化的关系。因为水滴受热蒸发的过程会使水滴温度变化,水滴温度的变化也会影响到eq2。

1.  “这和eq1 里T=T(r,t)作为r和t的函数的r有什么不同?”这个问题我有点懵,eq2中的r就是特指球的半径(暂时用R表示),应该是个时间变量,即球的半径R随时间t不断变化。这样的话,eq2中r用R代替更合理。

eq1中的r如果是空间变量的话,即r在[0,R]的区间变化,我也只关心r=R处的温度。

2. eq1确实是能量方程。eq2是个非稳态状态方程。

ps准确的说,我最后想要水滴半径R与时间t的曲线。
我的思路是不是有问题啊?被你一问,突然觉得好像哪里不对似得
我有100双袜子,可以送你一双。
4楼2014-12-16 12:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

引用回帖:
4楼: Originally posted by ahui728 at 2014-12-16 12:21:42
我想还是先说一下物理现象。
球形水滴受热蒸发,主要考虑强烈的热辐射作用。开始时,水滴的温度和外界的温度相差很小。
eq2给出水滴半径随时间变化的关系。因为水滴受热蒸发的过程会使水滴温度变化,水滴温度的变 ...

ok, 我觉得你把一般的作为空间变量的半径和水滴球体的外半径,姑且设为a,混起来了。你想得到的应该是a=a(t)。
5楼2014-12-16 14:05:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ahui728

木虫 (文坛精英)

叫兽

引用回帖:
5楼: Originally posted by pippi6 at 2014-12-16 14:05:43
ok, 我觉得你把一般的作为空间变量的半径和水滴球体的外半径,姑且设为a,混起来了。你想得到的应该是a=a(t)。...

嗯,是这个意思。
我看别人的文献最后也都是套用这几个公式的。
现在问题是 我要知道T和t的关系,才能搞定eq2。我的思路有问题吗?
暂时只能想到eq1解决水滴温度随时间变化关系
我有100双袜子,可以送你一双。
6楼2014-12-16 14:37:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ahui728

木虫 (文坛精英)

叫兽

引用回帖:
5楼: Originally posted by pippi6 at 2014-12-16 14:05:43
ok, 我觉得你把一般的作为空间变量的半径和水滴球体的外半径,姑且设为a,混起来了。你想得到的应该是a=a(t)。...

补充楼上的回复:

或者,一步一步来,eq1方程怎么解?能出数值解然后作图也行。
我有100双袜子,可以送你一双。
7楼2014-12-16 14:50:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ahui728: 金币+20, ★★★很有帮助, 亲,能帮忙搞一下吗?对代码命令流实在不感冒。看到MATLAB的界面就没心情了。我最近在用comsol试,进展不顺利。。 2014-12-18 09:03:46
引用回帖:
7楼: Originally posted by ahui728 at 2014-12-16 14:50:19
补充楼上的回复:

或者,一步一步来,eq1方程怎么解?能出数值解然后作图也行。...

看不出解析有可能,你先给出数值,我们可以试一试数值解
8楼2014-12-16 15:29:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ahui728

木虫 (文坛精英)

叫兽

引用回帖:
8楼: Originally posted by pippi6 at 2014-12-16 15:29:25
看不出解析有可能,你先给出数值,我们可以试一试数值解...

k1=4.75,k2=0.6,k3=2000
r=0.001
t=0时, T=300,
先看看T~t的关系曲线
我有100双袜子,可以送你一双。
9楼2014-12-16 15:47:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ahui728 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 21/1050 2026-08-12 16:08 by 云上清扬
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见