<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">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><div>Antonio</div><div><br></div><div><br></div><div><div><br></div>-- <br><div class="gmail_signature">ANTONIO ANTORANZ PERALES<br>Universidad Carlos III de Madrid</div>
</div></div>