[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