<div dir="ltr">Hi again, NEK users.<div><br></div><div>I think I have found the problem. Following the table on page 17 of <a href="http://nek5000.github.io/NekDoc/Nek_users.pdf">http://nek5000.github.io/NekDoc/Nek_users.pdf</a>, I was mistakenly setting density in udiff and viscosity in utrans (field=1). Please correct this typo in the manual. NEK examples are however correct. </div><div><br></div><div>Regards</div><div>Antonio</div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-03-15 20:14 GMT+01:00 ANTONIO ANTORANZ PERALES <span dir="ltr"><<a href="mailto:anantora@ing.uc3m.es" target="_blank">anantora@ing.uc3m.es</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi,</div><div><br></div><div>After successfully running some DNS cases of a turbulent pipe with heat transfer in NEK (<a href="http://www.sciencedirect.com/science/article/pii/S0142727X15000533" target="_blank">http://www.sciencedirect.com/science/article/pii/S0142727X15000533</a>), I am now struggling to run a similar case but with variable fluid properties (Re and Pe being functions of temperature). </div><div><br></div><div>I am currently using Pn-Pn2 (lx2=lx1-2 in SIZE) method and stress formulation (IFSTRS=true), together with IFUSERVP=true (p030=1.0 in REA file) to be able to define the viscosity and conductivity functions. Initial flow conditions are taken from a previous constant fluid properties turbulent pipe case. </div><div><br></div><div>My problem is that, when setting IFUSERVP=true, the flow velocities go to zero in very few time steps, even though I was setting constant fluid properties in USERVP subroutine for trial purposes.</div><div><br></div><div>I have checked that the forcing in the momentum equation (subroutine USERF) is correctly set. </div><div><br></div><div>Could you please give me some indications as to what I can be doing wrong? Thanks in advance.</div><div><br></div><div>Regards</div><span class="HOEnZb"><font color="#888888"><div>Antonio</div><div><br></div><div><br></div><div><div><br></div>-- <br><div>ANTONIO ANTORANZ PERALES<br>Universidad Carlos III de Madrid</div>
</div></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">ANTONIO ANTORANZ PERALES<br>Universidad Carlos III de Madrid</div>
</div>