<div dir="ltr">In a 2-d case involving a vortex dipole, I have used Biot-Savart law to calculate the velocity. This is very inefficient of course, and I could only run it in serial. If you want to see, it is here<div><br></div>
<div><a href="https://code.google.com/p/cfdlab/source/browse/trunk/nek5000/vortex_dipole/">https://code.google.com/p/cfdlab/source/browse/trunk/nek5000/vortex_dipole/</a><br></div><div><br></div><div>See initcond.usr and the function usrchk.</div>
<div><br></div><div>praveen</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Feb 12, 2014 at 6:26 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">(Wed, Feb 12, 2014 at 06:47:46AM -0600) <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> :<br>

<div class="">> I've just posted this example to the repo as a demo on how to do<br>
> this.  It might shed some light on what you need to do.  Be aware,<br>
> however, that my demo had periodic boundary conditions, thus there<br>
> were no issues related to setting the boundary values, etc.  That<br>
> is a slightly deeper topic, but perhaps this will get you started.<br>
><br>
> Paul<br>
><br>
> PS - the new demo directory is examples/eddy_psi_omega<br>
<br>
</div>Hi Paul,<br>
Thanks a lot for your help, I'll try to start from that!<br>
<br>
Best regards,<br>
<span class="HOEnZb"><font color="#888888">--<br>
Ismael<br>
</font></span><br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>