dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #18269
Re: Strangeness in advection-diffusion demo
On 11 May 2010 20:23, Garth N. Wells <gnw20@xxxxxxxxx> wrote:
>
>
> 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.
Yes, I implemented a bunch of Python demos which was simply the C++
versions translated to Pydolfin syntax.
I also added a note on this at the top of all demo.py files to give
the original authors credit (and to avoid being hold responsible in
situations like this) but they must have been removed at some point.
Kristian
> 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
>
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin
> Post to : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help : https://help.launchpad.net/ListHelp
>
References