Solving two-temperature model using pdepe

Hi,

Am not so good in programming but am trying to do a simuation of laser ablation using two-temperature model. However, am not able to define boundary conditions correctly which is giving me a false output due to very limited temperature change. Can someone please help me code the boudary conditions for the following equation;

The code am having for the boudary condition is;

function [pl,ql,pr,qr] = pdebc(xl,ul,xr,ur,t) % Boundary Conditions

pl = [0;ul(1)];

ql = [300;300];

pr = [ur(1)-300;ur(2)];

qr = [300;300];

end

and the one for my initial condition is;

function u0 = pdeic(x) % Initial Conditions

u0 = [300; 300];

end

[Luis]

Hey, Do you upload the code to Github or somewhere, I'm a new user of Matlab and I'll like to do the same simulation, nevertheless i can't find more information than this post to move forward.

Thx in advance.

Answer 1

[Saurav Chaudhary]

You can refer following MATLAB answer to get an idea about how to define boundary conditions for the equations.

https://www.mathworks.com/matlabcentral/answers/17349-pdepe-and-bcfun?s_tid=answers_rc1-1_p1_Topic

Also the following documentation link contains multiple example that may help you with coding boundary conditions and getting more insight into pdepe.

https://www.mathworks.com/help/matlab/ref/pdepe.html