<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto"><div></div><div>You can provide your own implementation of the thermal divergence in userqtl in the .rea file. It's no longer hardcoded.</div><div><br></div><div>The MV boundary condition is currently not supported with PnPn.</div><div><br>On 19 May 2017, at 23:55, "<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>" <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br><br></div><blockquote type="cite"><div>


<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Kopano HTML builder 1.0">
<title></title>


<!-- Converted from text/plain format -->

<p><font style="font-family: courier" size="2">
Marco, <br>
<br>
>    - the low Mach does not seem to work for PnPn-2 formulation. Following this remark, what would be the correct set up with the low Mach approach to >compute the full stress tensor with variable viscosity?<br>
<br>
       You will have to use the Pn-Pn formulation with IFSTRS set to true.<br>
<br>
>    - the user manual reads 'The implementation of the equation if state for the Low Mach formulation is for the moment hard-coded to be the ideal gas >equation of state p = ρRT. This allows for both variable density and variable viscosity.’ Does that mean the density is always computed using the ideal gas >equation of state? What does happen when the density is specified in uservp?<br>
<br>
       For variable density, the RHS of the velocity constraint is non-zero, i.e. div u = -1/rho (D rho/D t) .  This was implemented in Nek5000 assuming an ideal gas equation of state regardless of what you prescribe in uservp.  You will need to modify qthermal.f to account for whatever equation of state you need.<br>
<br>
-Dillon<br>
<br>
<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">mailto:nek5000-users-bounces@lists.mcs.anl.gov</a>] On Behalf Of <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
Sent: Friday, May 19, 2017 2:51 PM<br>
To: <a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
Subject: [Nek5000-users] Low Mach configuration in Nek5000<br>
<br>
Hello,<br>
<br>
I am trying to simulate two membranes that produce a synthetic jet with Nek5000. I have been using the low Mach configuration in Nek5000 (the flow conditions are low Mach). The low Mach set up in the low_mach example was helpful but I still have a few questions:<br>
<br>
    - the low Mach does not seem to work for PnPn-2 formulation. Following this remark, what would be the correct set up with the low Mach approach to compute the full stress tensor with variable viscosity?<br>
<br>
    - the user manual reads 'The implementation of the equation if state for the Low Mach formulation is for the moment hard-coded to be the ideal gas equation of state p = ρRT. This allows for both variable density and variable viscosity.’ Does that mean the density is always computed using the ideal gas equation of state? What does happen when the density is specified in uservp?<br>
<br>
    - since I have a moving boundary condition ‘mv’ and I solve the heat equation, what would be a good boundary condition for the temperature at the membranes? I do not think ‘mv’ boundary type is supported for the heat equation. I have tried with ’t’ and ‘I’ without noticing any significant differences in the numerical solution.<br>
<br>
Thanks,<br>
<br>
Marco<br>
<br>
    <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">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><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">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</font>
</p>

</div></blockquote></body></html>