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

• To: mathgroup at smc.vnet.net
• Subject: [mg87803] Re: heat equation through different media/problem with constant flux
• From: Oliver Ruebenkoenig <ruebenko at uni-freiburg.de>
• Date: Thu, 17 Apr 2008 06:58:22 -0400 (EDT)
• References: <fu6mi8\$nnm\$1@smc.vnet.net>

```Luigi,

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 + \(TavRInt - TavSInt\)\/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
Supplement.

http://www.imtek.uni-freiburg.de/simulation/mathematica/IMSweb/imsTOC/Application%20Examples/Finite%20Element%20Method/CompositeMaterialsDocu.html

There is also a 3D version of it.

http://www.imtek.uni-freiburg.de/simulation/mathematica/IMSweb/imsTOC/Application%20Examples/Finite%20Element%20Method/CompositeMaterials3DDocu.html

hth,

Oliver

Oliver Ruebenkoenig, <ruebenko AT uni-freiburg.de>

```

• 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