Repository navigation
BUG: Reflective boudary doesn't reflect correctly BX3 in 2.5D #343
Description
Describe the issue:
Hi all,
By running MHD simulations in 2.5D (DIMENSIONS=2 and COMPONENTS=3), with the relfective boundary condition, I realised that the third component of the BX3 magnetic field is not reversed in the ghost cells as it should be (since it is tangential to the boundaries in IDIR and JDIR).
The reason seems to be that Idefix use only the cell-centered value of third component BX3 and not the face-centered one BX3s in 2.5D. So the reflective condition is applied in this loop
idefix/src/fluid/boundary/boundary.hpp
Lines 680 to 692 in ede2e1f
which does not reverse the magnetic field components. As I understand it, this is not a problem in 3D as they are overwitten later by the ReconstructVcField function, but in 2.5D, BX3 is left untouched. I fixed this by changing line 689 to
const int sign = (n == VX1+dir || (n >= BX1 && n != BX1+dir)) ? -1.0 : 1.0;which works well in my case, but I'm not sure if it is the best way to solve the problem.
Error message:
runtime information:
Idefix V2.2.01
Activity
fix reflective with MHD and DIMENSIONS < 3 (#343)
Fix axis and reflective boundary conditions (#345)
glesur commented on Jun 16, 2025
fixed and merged by #345