[Nek5000-users] uservp MHD error

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Dec 1 11:31:45 CST 2016


Hello All,

I want to run MHD case with different magnetic diffusivities in different
region of the domain.
Can we run MHD run with different magnetic diffusivities?
Assuming that udiff=magnetic diffusivity (eta), in the uservp file I added
the following

      1.0000000E+00     p30 > 0 ==> properties set in uservp()

      subroutine uservp (ix,iy,iz,ieg)
      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'

      if (ifield.eq.1) then ! velocity
utrans = param(1)
udiff  = param(2)
      elseif (ifield.eq.ifldmhd) then ! B-field
utrans = 1.0
        R1 = sqrt(x*x+y*y);
if(R1 .le. 2.0) then
udiff= param(29)
else
udiff= 1000.0D0*param(29)
endif
 endif

      return
      end

I get an error

 done :: set initial conditions

  ERROR:  Non-positive diffusivity (    0.00    ) specified for field 2,
group 0 element
  ERROR:  Non-positive diffusivity (    0.00    ) specified for field 2,
group 0 element

Attached are the .usr, .rea and the logfile

When I set the p30=0.000 and just utrans=0.0 udiff=0.0 in the uservp, the
simulation is running fine.
Thank you in advance.

Cheers,
Sandeep
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161201/cdb6da49/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gpf.usr
Type: application/octet-stream
Size: 12369 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161201/cdb6da49/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gpf.rea
Type: application/octet-stream
Size: 14433846 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161201/cdb6da49/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gpf.log.16
Type: application/octet-stream
Size: 18825 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161201/cdb6da49/attachment-0005.obj>


More information about the Nek5000-users mailing list