[Nek5000-users] Nek5000-users Digest, Vol 94, Issue 1
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Fri Dec 2 01:01:20 CST 2016
Hi Sandeep,
Is your self-defined udiiff a non-dimensional one? If not, maybe you can try
to set dimensionless properties in uservp subroutine, then the non-positive
error may disappear.
Cheers,
Bolun
----------------------------------------------------------------
Bolun Xu
University of Science and Technology of China
Department of Modern Mechanics
Hefei, Anhui, China
-----邮件原件-----
发件人: nek5000-users-bounces at lists.mcs.anl.gov
[mailto:nek5000-users-bounces at lists.mcs.anl.gov] 代表
nek5000-users-request at lists.mcs.anl.gov
发送时间: 2016年12月2日 4:18
收件人: nek5000-users at lists.mcs.anl.gov
主题: Nek5000-users Digest, Vol 94, Issue 1
Send Nek5000-users mailing list submissions to
nek5000-users at lists.mcs.anl.gov
To subscribe or unsubscribe via the World Wide Web, visit
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
or, via email, send a message with subject or body 'help' to
nek5000-users-request at lists.mcs.anl.gov
You can reach the person managing the list at
nek5000-users-owner at lists.mcs.anl.gov
When replying, please edit your Subject line so it is more specific than
"Re: Contents of Nek5000-users digest..."
Today's Topics:
1. uservp MHD error (nek5000-users at lists.mcs.anl.gov)
----------------------------------------------------------------------
Message: 1
Date: Thu, 1 Dec 2016 18:31:45 +0100
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] uservp MHD error
Message-ID:
<mailman.7023.1480623493.3602.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="utf-8"
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/cdb6d
a49/attachment.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/cdb6d
a49/attachment.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/cdb6d
a49/attachment-0001.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/cdb6d
a49/attachment-0002.obj>
------------------------------
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
End of Nek5000-users Digest, Vol 94, Issue 1
********************************************
More information about the Nek5000-users
mailing list