How to model a Phase Change Stefan problem with Simscape?

Question

Elisa Revello 2023-1-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1903145-how-to-model-a-phase-change-stefan-problem-with-simscape

评论： Sreevatsan 2024-6-25

I am trying to model a Phase Change Material with Stefan problem (attached in "PCM_1D Stefan problem.png" file) using Simscape library.

I assumed two thermal mass blocks (for the liquid and solid phase respectively) and conductive heat transfer to describe the phase change problem, as depicted in the scheme below. Particularly, I don't understand how to change from mass parameter to density parameter in the thermal mass block (because from my input data I only have the material density) and how to specify the initial and boundary conditions in "PCM_model.slx"?

Thanks in advance if someone could help me.

Input data for the PCM:

rho = 1370; % density [kg/m^3]

k_s = 0.830; % thermal conductivity solid phase [W/(m K)]

k_l = 0.660; % thermal conductivity liquid phase[W/(m K)]

cp_s = 1.69; % specific heat capacity solid phase [kJ/(kg K)]

cp_l = 1.96; % specific heat capacity liquid phase [kJ/(kg K)]

L = 227; % latent heat of the phase change [kJ/kg]

T_melt = 115+273; % melting temperature of the PCM [K]

alpha_s = k_s/(rho*cp_s); % thermal diffusivity solid phase [m^2/s]

alpha_l = k_l/(rho*cp_l); % thermal diffusivity liquid phase [m^2/s]

R_pcm = 0.05; % PCM radius [m] (supposed)

T_0 = -10+273; % initial temperature of the PCM [K]

Input data for the IMD (geometrical properties):

D_imd = 254*(10^-3); % diamiter [m]

h_imd = 186*(10^-3); % height [m]

T_wall = 135+273; % Temperature of the stator external surface [K]

2 个评论
显示无隐藏无

Mohamed El kholy 2023-4-15

Hello,

Could you reach out to me since I m working in the samet thing and we can discuss about it

Sreevatsan 2024-6-25

Hello Mohamed , iam working on modelling the PCM for battery cooling, can you help me with modelling.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Balavignesh 2023-2-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1903145-how-to-model-a-phase-change-stefan-problem-with-simscape#answer_1171765

Hi Elisa,

I understand that you are trying to simulate phase change model using simscape and seek clarifications on

How to change the mass parameter to density in the thermal mass block
How to specify initial and boundary conditions

To answer your question.

As per Mass in thermal systems - MATLAB (mathworks.com), there are no options to change the mass parameter to density.
You can set the initial conditions and boundary conditions by specifying the initial target (Set Priority and Initial Target for Block Variables) and nominal values (System Scaling by Nominal Values) respectively.