[FLASH-USERS] ensuring div B = 0 when adjusting evolution

[Slavin, Jonathan]

Hi all,

I'm running MHD simulations where I set off a second supernova explosion
within a pre-existing remnant.  I'm using Simulation_adjustEvolution for
this.  However, I'm running into a problem with the magnetic field just
after initiating the second explosion.  I'm getting very large
pixel-to-pixel variations in Bx and By at the edge of the new expanding
blast wave.  I'm thinking that it could be because of non-zero div B in the
region in which the explosion is generated. So my question is, does anyone
have a suggestion for div B cleaning at the point that the explosion is
started, i.e. within Simulation_adjustEvolution, to prevent the problems
I'm having? I don't really expect the B field to be dynamically important
inside the remnant (beta >> 1), so accuracy of the B field is probably not
important there.  I should add that I'm doing these calculations in 2D
cylindrical symmetry (R-Z).  Thanks in advance for any help.

Regards,
Jon

-- 
________________________________________________________
Jonathan D. Slavin                 Harvard-Smithsonian CfA
jslavin at cfa.harvard.edu       60 Garden Street, MS 83
phone: (617) 496-7981       Cambridge, MA 02138-1516
cell: (781) 363-0035             USA
________________________________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170525/9925b9e6/attachment.htm>