clear all
T_ini=70; %initial temperature
T_melt=58; %melting temperature
C_liq=2.51*10^3; %specific heat of liquid [J/KgC]
C_sol=2.95*10^3; %specific heat of solid [J/KgC]
rho_liq=760; %density of liquid [Kg/m3]
rho_sol=818; %density of solid [Kg/m3]
k_sol=0.24; %thermal conductivity [W/mC]
k_liq=0.24; %thermal conductivity [W/mC]
lamda=226*10^3; %heat of fusion [J/Kg]
T_wall=33; %Wall temperature
dx=2/32;
dt=0.01;
t=20;
%% Initial Condition
h_ini=rho_sol*C_sol*(T_melt-T_ini); %sensible volumetric enthalpy [J/m3]
%% Boundary Condition
x=1;
t=1;
h_ff(x,t)=rho_liq*C_liq*(T_wall-T_ini);
%% Loop
H=1000;
n=1;
if H<0
T(n)=T_melt+H/(rho_sol*C_sol);
elseif 0<H<rho_liq
T(n)=T_melt;
else
T(n)=T_melt+(H-rho_liq*lamda)/(rho_liq*C_liq);
end
T(n)
没有合适的资源?快使用搜索试试~ 我知道了~
try1.zip_melting_zip
共1个文件
m:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 178 浏览量
2022-07-14
16:22:45
上传
评论
收藏 535B ZIP 举报
温馨提示
enthalpy change exlicit method
资源推荐
资源详情
资源评论
收起资源包目录
try1.zip (1个子文件)
try1.m 963B
共 1 条
- 1
资源评论
小波思基
- 粉丝: 72
- 资源: 1万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功