MathGroup Archive 2008

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

Search the Archive

Re: heat equation through different media/problem with constant flux


On Thu, 17 Apr 2008, Luigi B wrote:

> Dear All,
> I am trying to solve the heat conduction problem in a sequence of
> three media with different properties. For that I am using NDSolve
> with my own grid. The code (without the 'tedious' definition of the
> time dependent boundary conditions) is:
> \!\(NDSolve[{=A2=DF\_t u[x,
>    t] == alfa[x]*=A2=DF\_{x, 2}u[x, t], u[
>        x, 0] == TavSInt[0] + \(TavRInt[0] - TavSInt[0]\)\/L*
>      x, u[0, t] == TavSInt[
>        t], \([L, t]\) == TavRInt[
>          t]}, u, {x, 0, L}, {t, 0, tmax}, MaxSteps -> 50000, Method -
>> \
> {"\<MethodOfLines\>", \ "\<SpatialDiscretization\>" -> {\ \
> "\<TensorProductGrid\>", "\<Coordinates\>" -> {mygrid}}}]\)
> However, i still do not get a satisfactory result. Probably because I
> am not including the condition that at the interface between two media
> the heat flux is constant. How can I do this?

The thing is here, that you have a jump condition between the two
materials. A crude workaround would be to refine the mesh at the
interface. Another option would be to smooth out the interface.

Jose LLuis G=C3=B3mez-Mu=C3=B1osome time a ago did a finite element version of a
similar problem. You can find it as part of the IMTEK Mathematica

There is also a 3D version of it.



Oliver Ruebenkoenig, <ruebenko AT>

  • Prev by Date: Re: A kernel, multiple notebooks, and Global?
  • Next by Date: Re: Re: Package to 'graft' Lists onto Matrices
  • Previous by thread: Basic algebraic identities
  • Next by thread: Re: heat equation through different media/problem with constant flux