← Back to team overview

dolfin team mailing list archive

Re: Strangeness in advection-diffusion demo

 

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?

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



Follow ups

References