<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>