24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1609  |  回复: 9

MSPX

木虫 (著名写手)

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

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

发自小木虫IOS客户端
回复此楼
已阅   回复此楼   关注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值。
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

木虫 (著名写手)

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 上海大学实验技术岗位非升即走 +9 嘻嘻哈哈乐呵呵 2026-05-15 9/450 2026-05-19 01:47 by jurkat.1640
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 评审有感 +6 popular289 2026-05-18 11/550 2026-05-18 22:40 by iaeyuan
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 7/350 2026-05-17 19:42 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:06 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 8/400 2026-05-17 07:26 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见