[Nek5000-users] [*] Re: MHD
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jun 28 22:56:56 CDT 2010
Hi Mani,
There are several things you have to do to get MHD working
SIZE:
lbx1=lx1, etc.
lbx2=lx2, etc.
lbelv=lelv,lbelt=lelt
ldimt=2 ! at least -- when there are no passive scalars
.rea/.re2:
p29 = -Rm ! magnetic Reynolds number
Boundary conditions for temperature (dummy field) and magnetic field
.usr:
in userchk...
ifheat = .false. ! unless you are also solving for heat/passive scalar advection-diffusion equation
if (istep.gt.0.and.mod(istep,iostep).eq.0)
$ call outpost(bx,by,bz,pm,t,' ') ! dumps B-field as velocity field into .f/.fld files with even numbers
in userbc/useric...
set B-field through something like
if (ifield.eq.1) then
ux = -sys ! Velocity components
uy = -cxc
uz = sxc + cys
elseif (ifield.eq.ifldmhd) then ! Magnetic field
ux = eps*sin(y)*sin(0.57*z)
uy = 0.
uz = 0.
endif
I am still working on 3D MHD example
Best,
Aleks
----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Sunday, June 27, 2010 11:10:55 AM GMT -06:00 US/Canada Central
Subject: Re: [Nek5000-users] [*] Re: MHD
Hi Aleks,
Any pointer to the switches that one needs to turn on to solve for
MHD using Nek5000?
Thanks,
Mani chandra
On 05/17/2010 11:29 PM, nek5000-users at lists.mcs.anl.gov wrote:
> Hi Mani,
>
> I'll set up an MHD example.
>
> Best,
> Aleks
>
>
> On Mon, 17 May 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Dear Paul and Stephan,
>>
>> I'd like to know how to setup nek5000 to solve for MHD equations.
>> Maybe an example?
>>
>> Regards,
>> Mani chandra
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
More information about the Nek5000-users
mailing list