<div dir="ltr"><br><div><br></div><div><div>Equation</div><div>Momentum Eq.1.1</div><div>Energy Eq.1.4</div><div>Passive scalar Eq.1.6</div><div>udiff</div><div>ρ</div><div>ρc p</div><div>(ρc p ) i</div><div>utrans</div><div>μ</div><div>k</div><div>k i</div><div>ifield</div><div>1</div><div>2</div><div>i-1</div></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>