presented at event SIAM-AMS Conference on Numerical Solution of Field problems in Continuum Physics Conference