I’m trying to add a non-uniform heat flux to my simple model of a cylinder with a center hole. I know I need to create a function that describes the heat flux (Gaussian) but not sure how to do that in the context of my model. I believe the function can be applied as a thermalBC to Face 2 as shown below.
T_initial = 20;
t_end = 100;

model=createpde(‘thermal’,’transient’);

geo = multicylcinder([.3,2],.2,’void’,[1,0]);
modelGeometry = geo;
pdegplot(model,’FaceLabels’,on);
thermalProperties(model,’ThermalConductivity’,4,’MassDensity’,210,’SpecificHeat’,73);
thermalIC(model,T_initial);

% Heat Flux Function
% ——————
% function Qflux = @nonUniformHeatFlux
% ???
% ???
% ???
% end

thermalBC(model,’Face’,1,’Temperature’,20);
thermalBC(model,’Face’,2,’HeatFlux’,@nonUniformHeatFlux)

generateMesh(model)
pdemesh(model)

tlist = 0:10:t_end
modelresults = solve(model, tlist)

endI’m trying to add a non-uniform heat flux to my simple model of a cylinder with a center hole. I know I need to create a function that describes the heat flux (Gaussian) but not sure how to do that in the context of my model. I believe the function can be applied as a thermalBC to Face 2 as shown below.
T_initial = 20;
t_end = 100;

model=createpde(‘thermal’,’transient’);

geo = multicylcinder([.3,2],.2,’void’,[1,0]);
modelGeometry = geo;
pdegplot(model,’FaceLabels’,on);
thermalProperties(model,’ThermalConductivity’,4,’MassDensity’,210,’SpecificHeat’,73);
thermalIC(model,T_initial);

% Heat Flux Function
% ——————
% function Qflux = @nonUniformHeatFlux
% ???
% ???
% ???
% end

thermalBC(model,’Face’,1,’Temperature’,20);
thermalBC(model,’Face’,2,’HeatFlux’,@nonUniformHeatFlux)

generateMesh(model)
pdemesh(model)

tlist = 0:10:t_end
modelresults = solve(model, tlist)

end I’m trying to add a non-uniform heat flux to my simple model of a cylinder with a center hole. I know I need to create a function that describes the heat flux (Gaussian) but not sure how to do that in the context of my model. I believe the function can be applied as a thermalBC to Face 2 as shown below.
T_initial = 20;
t_end = 100;

model=createpde(‘thermal’,’transient’);

geo = multicylcinder([.3,2],.2,’void’,[1,0]);
modelGeometry = geo;
pdegplot(model,’FaceLabels’,on);
thermalProperties(model,’ThermalConductivity’,4,’MassDensity’,210,’SpecificHeat’,73);
thermalIC(model,T_initial);

% Heat Flux Function
% ——————
% function Qflux = @nonUniformHeatFlux
% ???
% ???
% ???
% end

thermalBC(model,’Face’,1,’Temperature’,20);
thermalBC(model,’Face’,2,’HeatFlux’,@nonUniformHeatFlux)

generateMesh(model)
pdemesh(model)

tlist = 0:10:t_end
modelresults = solve(model, tlist)

end thermal analysis, nonuniform-heat flux MATLAB Answers — New Questions

​