Hi Mike,<br><br>Here are the .rea and .usr files of a case I did with 2 reacting scalars.  The files are modified versions of the rayleigh case from repo that simulates the Boussinesq equations, but I turned off the buoyancy term and added a second scalar with reaction.<br>

<br>In the userq routine in the .usr file, the qvol term specifies scalar production due to reaction.  The first scalar is called temp.  The next ones are ps(1), ps(2), etc.  The conditional statement (ifield.eq.2) specifies production of temp, the (ifield.eq.3) statement specifies production of ps(1).<br>

<br>The README of the shear4 example from repo, which has 4 passive scalars, describes how you can change the parameters relating to the multiple scalars.<br><br>hth,<br><br>David<br><br><br><div class="gmail_quote">On Thu, Jul 28, 2011 at 2:57 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">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;">All,<br>
<br>
I'm interested in using Nek for simulation of reacting flows.  Can anyone send me, or point me towards, a simple example that could get me going?  I guess I'm most interested to see the details of storage/operation with passive scalars.<br>


<br>
Thanks.<br>
--Mike<br>
<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>
</blockquote></div><br>