[Nek5000-users] Elsasser subroutines in the MHD solver

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Oct 29 09:23:50 CDT 2015


Hi,
I am simulating a flow in a rotating spheroidal shell (made from a cubed 
sphere grid) with a periodical forcing and no-slip boundaries for the 
velocity.  Without a magnetic field, this is working great. I recently 
tried my hand on dynamo simulations and tried to adapt quasi vaccum (n x 
B = 0) boundary conditions for the magnetic field. I have a local 
orthogonal coordinate system on the surface, where I set the two 
tangential components of B to 0 and use div B = 0 to find a value for 
the normal component . This vector is then transformed back to cartesian 
coordinates and used in userbc. This is of course a bit ad hoc and needs 
further testing.
However, I am not quite sure in how the MHD implementation in Nek5000 
actually works and wether it can be used for my problem. As far as I 
understand, the Elsässer formulation is mainly used for MHD turbulence. 
My main questions for further testing are
1) Is a kinematic dynamo possible with Nek5000? (we already have data 
from a finite volume code of kinematic dynamos for comparison)
2) Since the Elsasser variables combine u and B, are the forcings in 
userf actually seperate? I used "if (ifield .eq. 1)"  in userf, so that 
the forcing only adds to the velocity field, but the B field starts to 
grow before u.
As mentioned in the earlier mail, I already looked in the source code, 
but couldn't get any clear answers.

Cheers, Jan


More information about the Nek5000-users mailing list