[Nek5000-users] Vorticity as initial condition

[nek5000-users at lists.mcs.anl.gov]

On Thu, Feb 13, 2014 at 4:10 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hi Praveen
>
> (Thu, Feb 13, 2014 at 01:02:38PM +0530) nek5000-users at lists.mcs.anl.gov :
> >                         This is very inefficient of course, and I could
> > only run it in serial.
>
> It turns out finally that I missed that in your first message, and now it
> works perfectly, thanks a lot!
>
> Best regards,
> --
> Ismael
>

Ok. But it is good to still vary delta and see how the vorticity changes.
You can compute error in recomputed vorticity and try to find a good value
of delta.

Also due to roundoff errors, it may be good to enforce continuity of
computed velocity. Looking at the eddy_psi_omega example, I think this is
done by

call dsavg(vx)
call dsavg(vy)

Best
praveen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140213/36d14e8f/attachment.html>