24小时热门版块排行榜    

查看: 4344  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode,4个jtjc了 +4 ziyangfang 2026-08-19 5/250 2026-08-19 11:19 by ziyangfang
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +8 启萌科技 2026-08-12 30/1500 2026-08-19 10:22 by 启萌科技
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +52 医学老男孩 2026-08-13 114/5700 2026-08-19 09:38 by 医学老男孩
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见