pdepe solver boundary condition setup for discrete measured heat flux data

Consider 1D trainsient heat transfer into a slab(steel,0.3m thick) with uniform initial temperature T0, where the lower boundary is insulated and the upper boundary surface have a uniform along the surface but discrete in time heat flux term q(t)----supose q(t) includs the measured data point with sampling frequency 10Hz...How to set up the PDEPE solver to compute the temperature histry of a node inside the slab (x=0.15m）...

Thanks for all the help.

Answer 1

[Torsten]

By setting

pl = 0.0;
ql = 1.0;
qr = Area of the side face of the slab;
pr = q(t);

in the function where you define the boundary conditions for pdepe assuming that you defined

f = lambda*DuDx

in the function where you defined the differential equation.