<div dir="ltr">Hello Ishmael<div>The code uses a kernel regularization technique to avoid the singularity in the integrand. This has a parameter delta; try to change it and see what happens. </div><div><br></div><div>If you use the P(n)-P(n-2) scheme, then you could integrate on the pressure mesh. This should avoid the singularity problem and you dont need the regularization. This may be a better approach.</div>
<div><br></div><div>Best</div><div>praveen</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Feb 12, 2014 at 11:04 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:23:13PM +0100) <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> :<br>

<div><div class="h5">> (Wed, Feb 12, 2014 at 06:39:58PM +0530) <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> :<br>
> > In a 2-d case involving a vortex dipole, I have used Biot-Savart law to<br>
> > calculate the velocity. This is very inefficient of course, and I could<br>
> > only run it in serial. If you want to see, it is here<br>
> ><br>
> > <a href="https://code.google.com/p/cfdlab/source/browse/trunk/nek5000/vortex_dipole/" target="_blank">https://code.google.com/p/cfdlab/source/browse/trunk/nek5000/vortex_dipole/</a><br>
> ><br>
> > See initcond.usr and the function usrchk.<br>
><br>
> Hi,<br>
><br>
> I had a try with your code, and it looks very strange with your example<br>
> case (A lot of noise). However it looks perfect with the one that I want,<br>
> which is only radial vorticity, so thank you very much for your help!<br>
<br>
</div></div>Sorry, I spoke too fast: I was looking twice at the same image...<br>
I have strange behavior of the initcond routine: you can see here the<br>
initial vorticity and the vorticity -> velocity -> vorticity.<br>
<a href="http://www.normalesup.org/~bouya/temp/vorticity.png" target="_blank">http://www.normalesup.org/~bouya/temp/vorticity.png</a><br>
<br>
They don't look the same at all, even if one can more or less "guess" the<br>
calculated value in the middle. Did you have the same kind of issues?<br>
Did I miss something? (I tried to follow the steps in the README file)<br>
<br>
Thanks!<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>