[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