[Nek5000-users] [*] Re: [*] Re: MHD

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Tue Jun 29 00:36:23 CDT 2010


  Thanks for the reply Alex. I'll await your example.

Regards,
Mani

On 06/29/2010 09:26 AM, nek5000-users at lists.mcs.anl.gov wrote:
> 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
> _______________________________________________
> 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