[Nek5000-users] Vorticity as initial condition

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

Hi Ismael,

There is no automatic way to do this in Nek.

However, a few years ago I coded up a little streamfunction-vorticity
solver inside userchk and tested it on a (one, only) 2D benchmark.

I've just posted this example to the repo as a demo on how to do
this.  It might shed some light on what you need to do.  Be aware,
however, that my demo had periodic boundary conditions, thus there
were no issues related to setting the boundary values, etc.  That
is a slightly deeper topic, but perhaps this will get you started.

Paul

PS - the new demo directory is examples/eddy_psi_omega




On Wed, 12 Feb 2014, nek5000-users at lists.mcs.anl.gov wrote:

> Hi,
> I have a vorticity field expression and I would like to use it as an
> initial condition in nek5000 (It a 2D case). Is there a way to do that
> directly in nek5000, or do I have to compute the inverse laplacian manually
> before with some other program?
>
> Thanks in advance!
> -- 
> Ismael
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>