[Nek5000-users] Nek5000 pertubation mode
Paul Fischer
fischer at mcs.anl.gov
Wed Jun 3 04:49:38 CDT 2009
Dear Antonio,
Thanks for bringing this to my attention. I had similarly
found this as the issue last week and had started to investigate
when I got derailed by other projects. I'll try to sort this
out and hopefully repair the svn repo version of the code today.
Best regards,
Paul
On Wed, 3 Jun 2009, Antonios Monokrousos wrote:
> Dear Paul,
> I finally found what the problem with my case was. It was the
> de-aliasing. The problem appears first in the subroutine advabp (file
> perturb.f) where subroutine convop is called (file navier1.f). It seems
> that de-aliasing is not working for the pertubation/linear mode (?) When
> I deactivate the de-aliasing (from the rea file) the result looks correct.
> I tested it also with channel flow and it was the same problem: In
> pertubation mode, with de-aliasing a forcing appears in a few specific
> elements. I haven't been able to understand exactly why this happens but
> I suspect it has to do with the size of the arrays.
> Best regards
> Antonios
>
>
> --
> ---------------------------------------------------------------------------------------------
> Antonios Monokrousos
> Department of Mechanics email: antonios at mech.kth.se
> Royal Institute of Technology (KTH) phone: +46 8 790 68 76
> Osquars Backe 18 telefax: +46 8 796 98 50
> SE-100 44 Stockholm, SWEDEN
>
More information about the Nek5000-users
mailing list