[Nek5000-users] Primer for reacting flows

[nek5000-users at lists.mcs.anl.gov]

Hi Mike,

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.

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).

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.

hth,

David


On Thu, Jul 28, 2011 at 2:57 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> All,
>
> 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.
>
> Thanks.
> --Mike
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110728/981a21a9/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dd.rea
Type: application/octet-stream
Size: 4507 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110728/981a21a9/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dd.usr
Type: application/octet-stream
Size: 17269 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110728/981a21a9/attachment-0001.obj>