[FLASH-USERS] Non-ideal MHD

[Owen Matthews]

Further to my previous question, I have managed to get the code to run for 
a few (~70) time steps by letting the time step get down to ~10^{-17}. It 
still hangs at this point though. There's nothing obviously wrong with the 
output just before hangning. (div B is very small for example.)

Any input gratefully received.

Owen.



On Tue, 17 Jul 2007, Owen Matthews wrote:

> 
> Hi,
> 
> I have a problem with non-ideal MHD hanging.
> 
> I'm trying a pretty simple test problem using the Orszag-Tang vortex as a 
> template but going to non-ideal MHD and to 3D cartesian (v and B not a 
> function of z). I'm using FLASH 2.4 with divb_diffuse, viscosity/constant 
> and magnetic_resistivity/const, but other than that it's pretty much the 
> Orszag-Tang test problem.
> 
> It works fine with very small values of resistivity and viscosity 
> (~10^{-12}) but for anything bigger it dies after just a couple of time 
> steps. It doesn't crash but gets stuck, while the time step remains at the 
> minimum length. 
> 
> Has anybody got any idea what might be wrong here?
> 
> Thanks in advance for your help,
> Owen.
> 
> 
> 

-- 
**********************************************************
Dr. Owen Matthews
Laboratory for Astrophysics    Tel (PSI): +41 56 310 42 56
Paul Scherrer Institut         Tel (ETH): +41 44 632 38 14 
Wurenlingen und Villigen       Fax:       +41 56 310 26 46  
CH-5232 Villigen PSI
Switzerland                    matthews at astro.phys.ethz.ch
----------------------------------------------------------
http://www.astro.phys.ethz.ch/staff/matthews/matthews.html
**********************************************************