hazelsct wrote:This is a possible feature request though: it would be convenient if Elmer had an axis or symmetry boundary condition setting all of the normal fluxes, normal velocity, normal E and B field components, to zero.
This could indeed be implemented. Limiting to cases where x=0 is the symmetry axis would be rather straight-forward. However, the symmetry axis may also be some other axis. This would require some additional intelligence from the code.
The implementation could perhaps be put to the Solver_init code of each solver. If the symmetry flag would be present an appropriate Dirichlet condition would be added to the boundary list structure. And in GUI just a flag in BoundaryCondition / General as for the periodic stuff.