[FLASH-USERS] Magnetic fields at refinement boundaries in USM

[Sean M. Couch]

I’ve also encountered this same issue at refinement jumps with USM. I could never get HLLD to behave. What I found was ok was to use HLLC with high-order electric field construction. Then, even better, I turned on magnetic resistivity with a small value. That smoothed things out nicely.

May “belief” at this point is that either: 1) staggered mesh CT simply does not behave well with AMR, or 2) the prolongation/restriction routines in PARAMESH for face variables is bugged somehow.

John Z., Is this for the mesh refinement versions of Athena and Athena++? Or uniform grids?

Sean

----------------------------------------------------------------------
Sean M. Couch
Assistant Professor
Department of Physics and Astronomy
Department of Computational Mathematics, Science, and Engineering
Facility for Rare Isotope Beams
Michigan State University
567 Wilson Rd, 3250 BPS
East Lansing, MI 48824
(517) 884-5035 --- couch at pa.msu.edu --- www.pa.msu.edu/~couch
On Mar 29, 2019, 7:08 AM -0400, John ZuHone <jzuhone at gmail.com>, wrote:
Hi Jonathan and all,

I also wonder if similar AMR + finite volume MHD codes exhibit this issue.

I can verify that this does not happen in Athena or Athena++. Not sure about Enzo or GAMER, which are the two codes which come to mind.

John
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20190329/5075afd7/attachment.htm>