PDE-Toolbox: Time dependent applied Temperature

10 views (last 30 days)
Hello,
for the simulation of a brake disk, I want to apply the temperature for 7.2 seconds. Then there should only be convection, emissivity and Thermal Conductivity. How do I tell the toolbox to stop heating the disk after 7.2 seconds?
We've tried to put dT into a for-loop but the toolbox just takes the last known value.
Below there is a extract from our current program
%Variables
mm=200; %mass of the motorbike in kg
v0ms=250000/9; %velocity at the start in mm/s
c=460e6; %specific heat capacity
mb=0.43; %mass of the brake disk
thermalmodel = createpde('thermal','transient');
importGeometry(thermalmodel,'Model.stl');
generateMesh(thermalmodel);
thermalIC(thermalmodel,20);
%Equations
Ekin=0.5*mm*(v0ms^2); %kinetic energy of the bike
dT=Ekin/(c*mb); %Temperature change
tlist=0:1:20;
thermalBC(thermalmodel,'Face',95,'Temperature',dT); %Temperature, that should be applied for 7.2 seconds
%thermal losses
thermalmodel.StefanBoltzmannConstant=5.670373E-8;
thermalBC(thermalmodel,'Face',1:94,'ConvectionCoefficient',0.05,'AmbientTemperature',20);
thermalBC(thermalmodel,'Face',1:94,'Emissivity',0.4,'AmbientTemperature',20);%Konvektionskoeffizient
thermalBC(thermalmodel,'Face',96:100,'ConvectionCoefficient',0.05,'AmbientTemperature',20);
thermalBC(thermalmodel,'Face',96:100,'Emissivity',0.4,'AmbientTemperature',20);%Konvektionskoeffizient
thermalProperties(thermalmodel,'ThermalConductivity',lam,'MassDensity',rho,'SpecificHeat',c);
R = solve(thermalmodel,tlist);

Answers (1)

Kavya Vuriti
Kavya Vuriti on 8 Aug 2019
The temperature boundary condition can be specified as a function handle to specify temperature that depends on time and it can be passed to the function thermalBC.
dT= @(location,state) {code to set dT according to time}
thermalBC(thermalmodel,'Face',95,'Temperature',dT)
Where state is structure array and state.time field contains time at evaluation points.
For more information, refer the following link:

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!