<div dir="ltr">Hi Guillaume,<div><br></div><div>Thank you for clear my misunderstanding! I am not sure which exact the numerical scheme adopted in nek5000 to solve the eigenfunctions of direct and adjoint problems. I myself did not implement new things, just run a steady cylinder flow supplied in cyl_restart and used the setups of lin_chan2D to continue my simulation. I enclosed the parameter files here, please check where I was wrong. Thanks!</div><div><br></div><div>Wei</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 26, 2018 at 9:38 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">
<div bgcolor="#FFFFFF" text="#000000">
<div class="m_-3367495712329327242moz-cite-prefix">Hi,<br>
<br>
If I understand the lin_channel2D example correctly, it is not
using any Krylov method but simply sets an initial perturbation
and lets it evolve. Are you using the Arnoldi algorithm?<br>
<br>
I can have a look if you upload your files somewhere.<br>
<br>
Guillaume<div><div class="h5"><br>
<br>
Le 24/02/2018 à 08:53, <a class="m_-3367495712329327242moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.<wbr>gov</a>
a écrit :<br>
</div></div></div>
<blockquote type="cite"><div><div class="h5">
<div id="m_-3367495712329327242contentDescription" style="line-height:1.5;text-align:justify;text-justify:inter-ideograph">
<div>Dear Neks,</div>
<div><br>
</div>
<div>Recently, I am trying to set up a linear stability analysis
of flow past a cylinder at very low Re, see Re=20, after I had
a converged base flow, I modified the example “lin_chan*” and
read the dns result as baseflow and updated the *usr file,
after testing groups of parameters in *.par, I obtained
different groups of growth rates, but none of them is closed
to the value documented by Luchini et tal, JFM or the value
given by other codes. Are there more details or rules to
choose the parameters when performing a linear stability, e.g.
Krylov space dimension, number of converged eigenvalues to be
dumped?</div>
<div><br>
</div>
<div>Thanks!</div>
<div>Wei</div>
<blockquote id="m_-3367495712329327242ntes-iosmail-quote" style="margin:0">
<div dir="ltr" class="m_-3367495712329327242gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr"> </div>
</div>
</blockquote>
</div>
<br>
<fieldset class="m_-3367495712329327242mimeAttachmentHeader"></fieldset>
<br>
</div></div><pre>______________________________<wbr>_________________
Nek5000-users mailing list
<a class="m_-3367495712329327242moz-txt-link-abbreviated" href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<wbr>gov</a>
<a class="m_-3367495712329327242moz-txt-link-freetext" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a>
</pre>
</blockquote>
<p><br>
</p>
</div>
<br>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><font size="4"><b><font face="arial narrow,sans-serif"><br></font></b></font><span><div><font face="Bookman Old Style"><font color="#7f7f7f">Wei HE 何偉<br></font></font></div><div><font face="Bookman Old Style" color="#7f7f7f">School of Aeronautics</font></div><div><font face="Bookman Old Style" color="#7f7f7f">Universidad Politécnica de Madrid</font></div></span><a href="http://www.cfm.upm.es/" target="_blank">http://www.cfm.upm.es/</a><br></div></div></div>
</div>