<div dir="ltr">Hey NEKS<div>I am trying to run the conj HT example with variable properties where the density(param(1)), viscosity(param(2)), Cp(param(7)) and thermal conductivity(param(8)) are all functions of temperature.</div><div>I have made functions for all of these in the .usr file. However, the program blows up.</div><div>I am putting parts from the .usr file. please help</div><div><br></div><div><div> subroutine uservp (ix,iy,iz,ieg)</div><div> include 'SIZE'</div><div> include 'TOTAL'</div><div> include 'NEKUSE'</div><div><br></div><div> if (ifield.eq.1) then</div><div> param(1) = 298+(0.833/temp)</div><div> param(2) = (-15.837*temp*temp)+(167.47*temp)+159.626</div><div> utrans = param(1)</div><div> udiff = param(2)</div><div><br></div><div> else</div><div><br></div><div> param(7) = param(1)*((1.61*temp*temp)+(0.973*temp)+29.176)</div><div> param(8) = (-0.00177*temp*temp)+(0.0279*temp)+0.02479</div><div> utrans = param(7) ! thermal properties</div><div> udiff = param(8)</div><div><br></div><div> if (ieg .gt. nelgv) then ! properties in the solid</div><div> udiff = 0.1*param(8) ! conductivity</div><div> utrans = 1.0</div><div> endif</div><div><br></div><div> endif</div><div><br></div><div> return</div><div> end</div><div>c-----------------------------------------------------------------------</div><div> subroutine userf (ix,iy,iz,ieg)</div><div> include 'SIZE'</div><div> include 'TOTAL'</div><div> include 'NEKUSE'</div><div><br></div><div> ffx = 0.0</div><div> ffy = 0.0</div><div> ffz = 0.0</div><div><br></div><div> return</div><div> end</div></div><div>c-----------------------------------------------------------------------<br></div><div><div> subroutine userq (ix,iy,iz,ieg)</div><div> include 'SIZE'</div><div> include 'TOTAL'</div><div> include 'NEKUSE'</div><div><br></div><div> qvol = 0.0</div><div> if (ieg.gt.nelgv) qvol = 1.0</div><div><br></div><div> return</div><div> end</div></div><div>c-----------------------------------------------------------------------<br></div><div><div> subroutine userbc (ix,iy,iz,iside,ieg)</div><div>c NOTE ::: This subroutine MAY NOT be called by every process</div><div><br></div><div>C Set boundary conditions</div><div><br></div><div> include 'SIZE'</div><div> include 'TOTAL'</div><div> include 'NEKUSE'</div><div><br></div><div> ux = 4.0</div><div> uy = 0.0</div><div> uz = 0.0</div><div> TINF = 298.0</div><div> temp = 1.0</div><div><br></div><div> return</div><div> end</div><div>c-----------------------------------------------------------------------</div><div> subroutine useric (ix,iy,iz,ieg)</div><div><br></div><div>C Set initial conditions</div><div><br></div><div> include 'SIZE'</div><div> include 'TOTAL'</div><div> include 'NEKUSE'</div><div><br></div><div> ux = 0.0</div><div> uy = 0.0</div><div> uz = 0.0</div><div> temp = 0.0</div><div><br></div><div> return</div><div> end</div></div><div><br></div><div>Apoorva Vinayak Rudra<div>Graduate Research Assistant</div><div>Department of Chemical Engineering</div><div>City College of New York (CCNY)</div></div><div><br></div><div><br></div><div><div><div class="gmail_signature"><div dir="ltr"><div><br></div></div></div></div>
</div></div>