Hmm...
>>>Why don't you upload a worksheet showing how far you have got.
I don;t have it. I have not started yet.
I try to recognize if this problem is solvable.
>>>Which section?
Through whole wall, for example three different properties layer.
Assumptions are simple:
- equation: diff(T(x,t),t)=a*diff(T(x,t),x,x)
- three material: (thermal conductivity (1,2,3), dinsity (1,2,3), heat capacity (1,2,3):
- initial condition T(x,0)=298 [K]
- boundary condition: T(0,t)=1473 [K], -lambda*dT/dx=h*((T(L,t)-T0)+emissivity*SBoltzman*(T(L,t)^2-T0^4),
I tried according to Mapleprimes suggestion using thermal diffisivity (lamba/cp*rho) as a piecewise function but got wrong results.
Problem lie for probably in interior boundary condition.
We have heat transfer betwen layer -lambda1*dT(x1,t)/dx=lambda2*dT(x2,t).
Last sugesstion from MaplePrimes was:
"Ah, I think I see the problem. It is not dT/dx that is continuous across the boundary, but the heat flux lambda*dT/dx. I think one way to tackle this is by a piecewise transformation of the position variable x, so that a displacement by dx in region j corresponds to a change of dX = dx/lambda in the new variable X. Then dT/dX is continuous across the boundary."
---
Generally speaking, I don;t have idea how to set up pdsolve!
wzel
