MHD setup in Dedalus2

[Erin E Goldstraw]

Dear Dedalus Users,

I am trying to solve the MHD equations. 

The equations and boundary conditions are given in Equation_Summary.pdf. 

My current code setup is given in solar-fmhd-d2_minimal.py. 

Could someone please check this set up? 

I am not sure about the gauge conditions and keep encountering problems with numeric instabilities when running this code. 

Any advice would be greatly appreciate. I would like to understand the Dedalus methodology better.

Please and thank you,

Erin Goldstraw

[Sadokat Malikova]

Dear Erin Goldstraw,

From a quick glance, you are imposing perfectly conducting wall condition:
∂zBx = 0, ∂zBy = 0, Bz = 0.

In your code try with:

 bz=0  bc. 

I think this is the problem.

with best regards,

Sadokat

--
You received this message because you are subscribed to the Google Groups "Dedalus Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to dedalus-user...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/dedalus-users/056c8aa2-ee28-4e15-956d-34521e324933n%40googlegroups.com.

[Erin E Goldstraw]

Dear Sadokat,

Thank you for your suggestion. 

I think I have previously tried both d/dz(bz)=0 and bz=0 at one time or another without much affect. 

I will try again. 

Any other suggestions are welcome.

Thank you,

Erin