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]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] 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] simscape, matlab, phase change, heat transfer, stefan problem MATLAB Answers — New Questions

​