[Danath]

Can somebody explain me why OpenFoam uses

Su()=mass*3.0*constant::electromagnetic::mu0.value ()/p.rho()*magneticSusceptibility_/(magneticSusceptibility_ + 3)*HdotGradHInterp.interpolate(p.position(), p.currentTetIndices());

F = 3 * mass * μο / density_shere * χu / (χu + 3) Η * divH

in order to calculate paramagnetic force ?

i search bibliography but still i cant find something relevant

thanks in advance !