MathGroup Archive 2008

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

Search the Archive

Re: 2 domain PDE, NDSolve

  • To: mathgroup at smc.vnet.net
  • Subject: [mg88606] Re: 2 domain PDE, NDSolve
  • From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
  • Date: Fri, 9 May 2008 07:16:56 -0400 (EDT)
  • Organization: The Open University, Milton Keynes, UK
  • References: <g00v09$gkr$1@smc.vnet.net>

jeramyk wrote:

> 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 
> layers.
> 
> 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).

<snip>

Please, could you post genuine Mathematica expressions as well as what 
you tried with NDSolve? (I believe that by writing T1(0,1)[z,t] you mean 
actually Derivative[0, 1][T1][z, t], but I am too lazy today to correct 
all the expressions :-)

Regards,
-- Jean-Marc


  • Prev by Date: Re: tablename[[x]]= xth row, how to get xth column
  • Next by Date: Re: tablename[[x]]= xth row, how to get xth column
  • Previous by thread: Re: Display commands in a package
  • Next by thread: Re: 2 domain PDE, NDSolve