[dune-pdelab] Unphysical solution

[Oliver Sander]

Dear conf86,
this is a mailing list for technical specific questions regarding the 
development and use of
the dune-pdelab module.  Vague questions in the style of 'my code ain't 
working -- help me'
are out of place here.  Please go ask your thesis supervisor about stuff 
like that.
To be quite honest, we are even less inclined to help you because you 
consistently refuse
to give your full name and affiliation.  In this community this is 
interpreted as a sign of
disrespect.  I doubt you will get any more answers here until you share 
your identity.
Regards,
Oliver

On 21.10.2015 11:45, conf86 at web.de wrote:
> I have the following problem:
> A sperical body has a radially dependent initial temperature u_0 from a shock wave as given in u_0.png. There is no further heat source and the body simply cools down.
> So I have a*\partial_t u - \Delta u = 0 and Neumann boundary conditions for the outflow (black body radiation).
>
> However as seen in u_t.png the solution looks very unphysical, because the heat should not distribute on the surface in this fashion.
> I did a model calculation with a constant initial temperature beforehand and it looked ok, so I really don't know what was going wrong here.
> However the gradient of the initial condition is rather high, could this be the reason?
> I double-checked possible sign errors and similar, but I couldn't find any.
>
> Do you have any idea what could be reason for this or where to look at?
>
>
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20151021/61cbd972/attachment.htm>