[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