piece wise function for thermal properties in PDE solver

Hello all,

I need to basically create a piece wise function for heat capacity (C) and thermal conductivity (k) to accomodate different material properties across a boundary. For example, if x and y is within certain limits, then k=a+bT^3, else k = c + dT^3. Similary for C. I tried simply using an if then statement, but matlab will pass all of the geometrical points of the meshgrid to evaulate C and k and writing a for loop increased the solution time to be unreasonable. I would appreciate any thoughts or guidance on the issue.

Thanks

Answer 1

[Torsten]

When creating the geometry for your problem, I'd already take into account these two spatial regions where the different thermal properties apply. This way, you can separately apply different formulas to different RegionIDs for heat capacity and thermal conductivity.

[Brant]

That makes perfect sense! Thanks for the help