24小时热门版块排行榜    

查看: 4150  |  回复: 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的回帖

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的回帖

yuan1988hao

铜虫 (小有名气)

【答案】应助回帖

★ ★
zhs4276689: 金币+2, 就是想偷懒哦。。。呵呵 2013-06-11 07:46:45
把热传导方程离散成差分方程的形式,然后用迭代的方法求解差分方向组就行了。可你参考下西交大陶文铨编写的《数值传热学》,里面详细介绍整个过程如何实现,把例子看一下就懂了,很好理解的。
3楼2013-06-10 16:32:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhs4276689 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学336分求调剂 +3 zbzihdhd 2026-03-15 4/200 2026-03-16 23:52 by zbzihdhd
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 求调剂,一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见