MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

2 domain PDE, NDSolve

  • To: mathgroup at
  • Subject: [mg88585] 2 domain PDE, NDSolve
  • From: "jeramyk" <jeramyk at>
  • Date: Fri, 9 May 2008 03:24:06 -0400 (EDT)
  • Organization: the University of Washington
  • Reply-to: "jeramyk" <jeramyk at>

I'm working on solving the basic 1D, transient heat transfer equation in a 
two layer slab configuration.  One side is insulated, the other has a 
constant heat flux, and there is not contact resistance between the two 

My formulation looks like this:
equation1:  T1(0,1)[z,t] == k1/rho1/Cp1*T1(2,0)[z,t]
equation2:  T2(0,1)[z,t] == k2/rho2/Cp2*T2(2,0)[z,t]
IC1:  T1[z,0] == T0
IC2:  T2[z,0] == T0
BC1:  T1(1,0)[0,t] == 0
BC2:  -k1*T2(1,0)[L2+L1,t] == qs
Match1:  k1*T1(1,0)[L1,t] == k2*T2(1,0)[L1,t]
Match2:  T1(1,0)[L1,t] == T2(1,0)[L1,t]

T1[z,t] is solve from 0 to L1, T2[z,t] from L1 to L1+L2.

I've tried NDSolve but it doesn't like the matching conditions (it 
interprets them as BC's and since they're not on the boundary it has issue 
with it).

There's an analytical solution but it's nearly unworkable...

Can this be solved in mathematica? 

  • Prev by Date: Re: substitution within a substitution list
  • Next by Date: Re: substitution within a substitution list
  • Previous by thread: Re: tablename[[x]]= xth row, how to get xth column
  • Next by thread: Solving Large Overdetermined Systems