Using pdepe for 1d transcient heat conduction through a composite wall

5 views (last 30 days)
I have four adjacent slabs each modelled by a pde. My problem is how to handle the boundary conditions at the interfaces. I know that the temperature and the first partial derivative are continuous across the perfect interface. I want to know how to define the bcfun function at the interfaces

Accepted Answer

Bill Greene
Bill Greene on 23 Apr 2016
Edited: Bill Greene on 23 Apr 2016
You want to model this multi-region wall with a single pde and then write your pde function (pdefun) so that it returns a different flux value (f) and capacitance (c) depending on the slab. You use the input x coordinate to test for the particular slab the coefficients are requested for. Also, as described in the documentation, you want to insure that you have mesh points (xmesh) at the intersection points of the slabs.
Then the BC function is defined to return boundary conditions only at the true boundaries of the wall.
Bill Greene
Bill Greene on 23 Apr 2016
My response should have said that you want to return both an f and a c depending on the x value. f is ki*DuDx and c is rhoi*ci where i=1,4 in your example. Your test on x is essentially correct but could be done better with an if/else block.

Sign in to comment.

More Answers (1)

Krishnaraj Sambath
Krishnaraj Sambath on 7 Oct 2016
@James: There's a subtlety to Bill's response. You need to update both c & f as rhoi*Cpi & ki, respectively, and not just f as ki / (rhoi * Cpi) and leaving c as 1 for all layers. Doing it the latter way will impose incorrect conditions at the interfaces.


Find more on Partial Differential Equation Toolbox in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!