<div dir="ltr"><div><div><div>Hey Juan!<br><br></div>It turns out I had been working on a similar problem for the past few days and came up with a small subroutine. It probably is not the most efficient way to do it, but it at least work when I use global direct and adjoint modes . I'll send it to you tomorrow, and if you wanna discuss, I can drop by ONERA next week.<br>

<br></div>Cheers.<br></div>JC<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014/1/21  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span><br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Neks,<br>
<br>
I would like to project the linearized NS equations onto a reduced basis formed by some POD modes (velocity vectors).<br>
We know that the perturbation equations can be read as du(t)/dx = A u(t), where<br>
<br>
// A = ( -U.grad(u) - u.grad(U) - grad(p) + 1/Re Lap(u)//<br>
<br>
To obtain the projected system, I need to compute the matrix Ar_i,j =  <U_i, AU_j> where U_j is the j mode (a vector containing the velocity  field).<br>
The step that  I  can't clearly see is the matrix-vector product A U_j.  Do  I  have to explicitily compute the matrix A  and then obtain the product? or is there a direct  way to  do  it  without this matrix (using the time-stepper)?<br>


<br>
Thanks in advance!<br>
Sincerely,<br>
Juan--<br>
<br>
<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr">Jean-Christophe Loiseau<br><a href="https://sites.google.com/site/loiseaujc/" target="_blank">Homepage</a><br></div>
</div>