How can I solve for temperature distribution in a composite system (several pde with neumann boundary condtiions) using pdepe/any other PDE solver?

I have a system made of three layers of different materials, placed in parallel making a composite system. Heat is generated in the middle layer and is conducted to the first and bottom layer. I want to study temperature distribution in z direction with time. Hence I have three PDE's and each PDE works in different z intervals, namely z=0to a, z=a to b, and z=b to c. At the interfaces, my boundary condition is of constant heat flux, q1=q2 and temperature contuinuity, T1=T2. How can I solve three PDE's simultaneously, with these boundary conditions?

I am attaching my problem file. I would be grateful if anybody could help!

Answer 1

[Torsten]

Simply use pdepe, include the interface points within your mesh and make an if-statement to account for the different material properties.

Best wishes

Torsten.