[FLASH-USERS] Magnetic field boundary condition for NohCylindrical test problem

[Ernesto Zurbriggen]

Hi all!
Continuing with Greg's question, I have the following doubt.
Implementing boundary conditions for the magnetic field using the unsplit
staggered mesh,
in the boundary region:
Is essential to define the field either just in the cell center or just in
the cell faces, or in both places is necesary?

Thank you!




El vie., 24 ago. 2018 a las 11:59, Ernesto Zurbriggen (<
ezurbriggen at unc.edu.ar>) escribió:

> Hi Greg!
> I'm not in the details, but remember that in cylindrical geometry
> (IAXIS,JAXIS,KAXIS)=(r,z,theta).
>
> >>I am using the Unsplit Staggered Mesh MHD solver, and have included a
> modified Grid_bcApplyToRegion.F90 that alters the standard "reflect"
> condition on the HIGH face of IAXIS to:
> So, if you would like to modify the BC on the HIGH face of IAXIS,
>
> >> if(ivar == MAGZ_VAR) regionData(k-i,1:je,1:ke,ivar)= 1.0
> maybe you should put ivar == MAGX_VAR instead.
>
> Best!
>
> Ernesto
>
>
> El vie., 24 ago. 2018 a las 5:48, Boyle, Gregory (<gregory.boyle at desy.de>)
> escribió:
>
>> Dear FLASH community,
>>
>> I am attempting to modify the the 2D NohCylindrical test problem (32.2.4)
>> to include a permanent B_{\theta} condition at rightmost r boundary. This
>> is a first step towards including the effect of a longitudinal current.
>>
>> I am using the Unsplit Staggered Mesh MHD solver, and have included a
>> modified Grid_bcApplyToRegion.F90 that alters the standard "reflect"
>> condition on the HIGH face of IAXIS to:
>>
>> k = 2*guard+1
>> if(isFace)k=k+1
>> do i = 1,guard
>>    regionData(k-i,1:je,1:ke,ivar)= regionData(i,1:je,1:ke,ivar)*sign
>>    if(ivar == MAGZ_VAR) regionData(k-i,1:je,1:ke,ivar)= 1.0
>> end do
>>
>> However, when I look at the output it seems the condition of setting
>> MAGZ_VAR to 1.0 is being ignored.
>>
>> Any suggestions would be greatly appreciated!
>>
>> Cheers,
>> Greg.
>>
>
>
> --
> *Dr. Ernesto Zurbriggen*
>
> *Consejo Nacional de Investigaciones Científicas y Técnicas **(CONICET)**,
> Instituto de Astronomía Teórica y Experimental (IATE).*
>
> *Universidad Nacional de Córdoba **(UNC)**, **Observatorio Astronómico de
> Córdoba **(OAC)**.*
>
> *Teléfono: +54 0351 4331064-5, interno 222.*
>
> *Laprida 854, Bº Observatorio (CP 5000), Córdoba, Argentina.*
>
>

-- 
*Dr. Ernesto Zurbriggen*

*Consejo Nacional de Investigaciones Científicas y Técnicas **(CONICET)**,
Instituto de Astronomía Teórica y Experimental (IATE).*

*Universidad Nacional de Córdoba **(UNC)**, **Observatorio Astronómico de
Córdoba **(OAC)**.*

*Teléfono: +54 0351 4331064-5, interno 222.*

*Laprida 854, Bº Observatorio (CP 5000), Córdoba, Argentina.*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20180824/207a88a3/attachment.htm>