24小时热门版块排行榜    

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

  • merry
    1 (33.3%)
  • me
    2 (66.7%)
查看: 1478  |  回复: 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能看出啥? +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +10 Tide man 2026-08-11 10/500 2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 固定端突然变了,今天 +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 +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见