[Nek5000-users] Dynamo set up with weird results
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Jul 28 09:49:51 CDT 2015
Hi again,
I was able to solve the problem: I wasn't aware that the forcing in
userf is also added to the induction equation. To avoid this, one has to
use:
ffx=0.0d0
ffy=0.0d0
ffz=0.0d0
if (ifield .eq. 1) then
c----velocity forcing here
ffx=...
ffy=...
ffz=...
endif
as in the mhd example case. My example seems to be running fine (for now).
Cheers,
Jan
More information about the Nek5000-users
mailing list