How to model a Phase Change Stefan problem with Simscape?

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]

[Mohamed El kholy]

Hello,

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

[Sreevatsan]

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

[Sooraj]

Hello Sreevatsan, I am working on PCM cooling with batteries and would love to communicate on it.

Answer 1

[Balavignesh]

Hi Elisa,

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

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

To answer your question.

1. As per Mass in thermal systems - MATLAB (mathworks.com), there are no options to change the mass parameter to density.
2. 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.

Hope this helps