应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 工科数学 » 差分法处理第三类非线性边界条件
101/1返回列表
查看: 1714  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

MSPX

木虫 (著名写手)

[求助] 差分法处理第三类非线性边界条件 已有2人参与

现在做一维传热模拟，要考虑辐射传热边界条件，温度T的梯度等于温度4次方的函数，请问在编程时如何处理？谢谢！

发自小木虫IOS客户端
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-10-30 00:35:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenyxtata

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
MSPX: 金币+2, ★有帮助 2015-10-30 09:57:04
feixiaolin: 金币+16, http://emuch.net/bbs/viewthread.php?tid=9565441&fpage=1&target=self&page=2 奖励 2015-11-01 21:04:04
把计算域往外延一个点，然后同时应用差分方程和边界条件处理，消去这个额外点的T值。
赞一下(1人) 回复此楼
2楼2015-10-30 07:30:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laosam280

禁虫 (正式写手)
★ ★
感谢参与，应助指数 +1
MSPX: 金币+2, ★有帮助 2015-10-30 09:57:16
本帖内容被屏蔽
3楼2015-10-30 08:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MSPX

木虫 (著名写手)
引用回帖:
2楼: Originally posted by shenyxtata at 2015-10-30 07:30:33
把计算域往外延一个点，然后同时应用差分方程和边界条件处理，消去这个额外点的T值。

但是这种做法只能提高一阶精度，并不能将边界点问题T^4表示成倒数第二个结点的函数，因此不好用内部结点的温度向量和其系数矩阵求解。
赞一下 回复此楼
4楼2015-10-30 09:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MSPX

木虫 (著名写手)
引用回帖:
2楼: Originally posted by shenyxtata at 2015-10-30 07:30:33
把计算域往外延一个点，然后同时应用差分方程和边界条件处理，消去这个额外点的T值。

但是这种做法只能提高一阶精度，并不能将边界点问题T^4表示成倒数第二个结点的函数，因此不好用内部结点的温度向量和其系数矩阵求解。
赞一下 回复此楼
5楼2015-10-30 09:44:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MSPX

木虫 (著名写手)
引用回帖:
3楼: Originally posted by laosam280 at 2015-10-30 08:32:27
如果求解区域固定，也可以做奇延拓或者偶延拓，获得周期边界条件，这样形成的系数矩阵就第一列和最后一列可以耦合起来。

周期性边界条件？傅里叶变换的条件？不是太懂啊？能具体说下吗？我传个图你帮我看看。

第一幅图是控制方程和边界条件的差分形式，i表示空间点（i=0和i=N分别表示边界结点），j表示时间点，T（0，j）和T（N，j）表示j时刻边界上的温度值。

第二幅图是写成矩阵的形式进行求解，如果只有对流条件，很容易将边界条件中T0/TN表示成T1/TN-1的表达式，再带入控制方程来消去TO和TN，这样就可以用内部结点计算温度值，然后再外插到边界温度值；  但是现在有T^4，温度值的4次方，不知该如何处理？谢谢！
差分法处理第三类非线性边界条件
1.png


差分法处理第三类非线性边界条件-1
2.png
赞一下 回复此楼
6楼2015-10-30 09:54:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MSPX

木虫 (著名写手)
还望大神们不吝赐教！！
赞一下 回复此楼
7楼2015-10-30 12:40:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MSPX

木虫 (著名写手)
顶一下

发自小木虫IOS客户端
8楼2015-10-30 18:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenyxtata

金虫 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by MSPX at 2015-10-30 09:44:40
但是这种做法只能提高一阶精度，并不能将边界点问题T^4表示成倒数第二个结点的函数，因此不好用内部结点的温度向量和其系数矩阵求解。...

既然边界是非线性的，整个问题就是非线性的。先尝试牛顿迭代吧。

发自小木虫Android客户端
赞一下 回复此楼
9楼2015-10-31 09:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MSPX

木虫 (著名写手)
引用回帖:
9楼: Originally posted by shenyxtata at 2015-10-31 09:43:44
既然边界是非线性的，整个问题就是非线性的。先尝试牛顿迭代吧。
...

谢谢，我去看看

发自小木虫IOS客户端
回复此楼
10楼2015-10-31 09:57:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 MSPX 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-02 00:15 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 19/950 2026-11-02 00:02 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 15:04 by BCdvBAEpt3zm
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
关闭
请填处理意见
关闭