24小时热门版块排行榜    

查看: 4152  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhs4276689

木虫 (正式写手)

[求助] 求利用matlab有限差分法解一维传热问题

这个问题困扰我一个多学期了,真心搞不懂了,主要是数学和软件功底都不好,因此只能求助大家了。谢谢先。

问题描述:熔融金属与金属基板相接触,只考虑一维方向,计算其凝固层厚度随时间的变化曲线。有方程(图1),其中,Ts为基板初始温度,Tl为熔体温度,除X(t)和t外其余都为已知量。

根据不同的熔体和基板初始温度,由上式X(t)的正负可出现两种情况:可能熔体凝固,也可能会使基板先熔化再凝固。其结果应该是如图2类似的曲线

文献中说利用有限差分法计算,原话直译为“该数值模型可利用有限差分法解决。中心有限差分用于计算区域,而前、后有限差分用于计算边界点。每个计算区域沿Z轴分为100个点,时间步长为10-3。Crank-Nicholson方程用于稳定。”。
边界条件如图3
固液界面定义使用Stefan 方程,如图4
初始条件如图5
详细可见图6

我目前只能搞懂文献中的其它问题,但是如何计算就不知道了,matlab应该可以实现有限差分的计算,但具体代码是怎么样的?

希望有人可以详细解惑,尽量详细,无以答谢,只有些许金币了。
求利用matlab有限差分法解一维传热问题
1.jpg


求利用matlab有限差分法解一维传热问题-1
2.jpg
求利用matlab有限差分法解一维传热问题-2
3.jpg


求利用matlab有限差分法解一维传热问题-3
4.jpg


求利用matlab有限差分法解一维传热问题-4
5.jpg
求利用matlab有限差分法解一维传热问题-5
6.jpg



[ Last edited by zhs4276689 on 2013-5-28 at 07:46 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuan1988hao

铜虫 (小有名气)

【答案】应助回帖

★ ★
zhs4276689: 金币+2, 就是想偷懒哦。。。呵呵 2013-06-11 07:46:45
把热传导方程离散成差分方程的形式,然后用迭代的方法求解差分方向组就行了。可你参考下西交大陶文铨编写的《数值传热学》,里面详细介绍整个过程如何实现,把例子看一下就懂了,很好理解的。
3楼2013-06-10 16:32:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

xiaoyueriyue

金虫 (正式写手)

【答案】应助回帖

★ ★ ★
dbb627: 金币+1, 感谢参与讨论 2013-06-07 19:45:24
zhs4276689: 金币+2, 谢谢支持,不过太简单了。。。呵呵 2013-06-11 07:46:12
有限差分用MATALB实现和VC基本上是一样的,先懂得差分格式,然后进行迭代就好了, 用循环语句就可以了
你找一下差分格式,带进去就可以实现了。
温度场的实现的算法应该是ADI,不过你的模型是一维模型,就用x方向的差分格式就行了
2楼2013-06-06 22:19:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们,祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见