<div dir="ltr"><div>Dear FLASH users,</div><div><br></div>I have a question related to implementing a 'rigid body' (i.e. an object with reflective BC) for the unsplit staggered mesh (USM) MHD solver. It appears that the BC for B_x and B_z components of magnetic fields (i.e. in plane components for the 2D cylindrical geometry) are not assigned appropriately.<br>    <br>I'm running a modified 'laser-slab' simulation using the USM MHD solver, where I added a strong externally applied magnetic field along the Z-axis (2D cylindrical geometry) and defined a mesh using the 'rigid body' (take a look at the init.block file lines 166 - 191).  When the plasma flow reaches rigid bodies, the simulation stops due to a negative internal energy error. As shown in the plots below (for B_x and B_z), there are some abnormal values at the rigid body boundaries. <div><img src="cid:ii_l4rbt0853" alt="image.png" width="472" height="432"><br></div><div><img src="cid:ii_l4rbtxqf4" alt="image.png" width="472" height="434"><br></div><div><br>It seems that the issue is due to the boundary condition of rigid bodies for B_x and B_z (but not for Hydro), since the out of plane B_phi component works fine. So I'm wondering if there is a way to appropriately assign the B_x and B_z boundary conditions at the 'rigid body', because the FLASH manual (34.1.2) says "It is possible that the unsplit staggered mesh MHD solver can support the rigid body but we have not tested yet."  <br><br>Attached are my Simulation_initBlock.F90 and flash.par files. <div>Thank you!<br></div><div><br></div><div>Yigeng Tian, Graduate student</div><div>Department of Physics and Astronomy</div><div>University of Delaware</div></div></div>