dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

On 11/05/10 19:22, Johan Hake wrote:
> On Tuesday May 11 2010 11:09:26 Garth N. Wells wrote:
> > On 11/05/10 18:11, Johan Hake wrote:
> > > On Tuesday May 11 2010 06:48:04 Anders Logg wrote:
> > > > The advection-diffusion demo (both C++ and Python versions) says in
> > > > the comment that a stabilized method is used, but the stabilization
> > > > seems to be missing.
> > >
> > > Yes, it looks like the stabilization part is missing.
> > >
> > > > I think I remember this being added at some point, by either Hake or
> > > > Garth. Any clue where this went?
> > >
> > > I have mentioned adding a SUPG stabilized demo a while ago. I do not know
> > > who added this one, eventhough my name comes up on bzr annotate :P. It
> > > looks like it is Kristian who is the original author?
> >
> > I added it a long time ago, but I don't know where it went.
> >
> > Kristian can't be the author since it pre-dates his arrival. Kristian
> > authored the DG advection-diffusion demo.
>
> Ok, maybe Kristian just added a Python version?

I forgot about that. Yes, he probably added the Python version.

Garth


> Johan
>
>
> > Garth
> >
> > > > There also seems to be a typo for the diffusive term which appears in
> > > > both LHS and RHS.
> > >
> > > The diffusive term looks correct to me. Note that the advective term is
> > > _not_ integrated by part. This implicitly sets a flux condition on the
> > > Neumann
> > >
> > > boundaries to:
> > >      dot(u*velocity,n)
> > >
> > > The diffusive flux across the boundaries are set to zero.
> > >
> > > Johan
> > >
> > > _______________________________________________
> > > Mailing list: https://launchpad.net/~dolfin
> > > Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> > > Unsubscribe : https://launchpad.net/~dolfin
> > > More help   : https://help.launchpad.net/ListHelp