How to apply a time decaying internal heat source?

Hi
I’m trying to model the decay of nuclear fuel following an accident, I need to be able to apply an equation for the decay heat produced by the fuel to the internal heat source of the model.
This is my current code with P being the equation that should represent the decay heat.

thermalmodel= createpde(‘thermal’,"transient");
gm = multicylinder([4.655e-3 4.865e-3 5.59e-3],12.6e-3);
thermalmodel.Geometry = gm;
pdegplot(thermalmodel, "CellLabels",’on’,’FaceLabels’,’on’,"FaceAlpha",0.2);

thermalProperties(thermalmodel,’Cell’,1,’ThermalConductivity’,7,"MassDensity",10.97, "SpecificHeat",370); % Uranium
thermalProperties(thermalmodel,’Cell’,2,’ThermalConductivity’,252.48e-3 ,"MassDensity",0.178e-3,"SpecificHeat",5.193); % Helium
thermalProperties(thermalmodel,’Cell’,3,’ThermalConductivity’, 21.5,"MassDensity",6.55,"SpecificHeat",0.328); % Zircalloy

tlist = (10:1:100);

for a = 1 : length(tlist)
P(a) = 100e6*0.066*((tlist(a).^(-0.2))- ((10*24*60*60)+tlist(a)).^(-0.2));
end

internalHeatSource(thermalmodel,P,"Cell",1);

thermalIC(thermalmodel,873);

generateMesh(thermalmodel);
pdemesh(thermalmodel);
tic
thermalresults = solve(thermalmodel,tlist);
toc
Tcenter = interpolateTemperature(thermalresults,0,0,0,1:numel(tlist))
plot(tlist,Tcenter)Hi
I’m trying to model the decay of nuclear fuel following an accident, I need to be able to apply an equation for the decay heat produced by the fuel to the internal heat source of the model.
This is my current code with P being the equation that should represent the decay heat.

tlist = (10:1:100);

for a = 1 : length(tlist)
P(a) = 100e6*0.066*((tlist(a).^(-0.2))- ((10*24*60*60)+tlist(a)).^(-0.2));
end

internalHeatSource(thermalmodel,P,"Cell",1);

thermalIC(thermalmodel,873);

generateMesh(thermalmodel);
pdemesh(thermalmodel);
tic
thermalresults = solve(thermalmodel,tlist);
toc
Tcenter = interpolateTemperature(thermalresults,0,0,0,1:numel(tlist))
plot(tlist,Tcenter) Hi
I’m trying to model the decay of nuclear fuel following an accident, I need to be able to apply an equation for the decay heat produced by the fuel to the internal heat source of the model.
This is my current code with P being the equation that should represent the decay heat.

tlist = (10:1:100);

for a = 1 : length(tlist)
P(a) = 100e6*0.066*((tlist(a).^(-0.2))- ((10*24*60*60)+tlist(a)).^(-0.2));
end

internalHeatSource(thermalmodel,P,"Cell",1);

thermalIC(thermalmodel,873);

Cart

Cart