use "generic steel" b(h) law - cf. tests done during ELEC0041 class :-)

Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro

@@ -460,8 +460,6 @@ Resolution {
   }
 }
 
-// Same PostProcessing for both static and dynamic formulations (both refer to
-// the same FunctionSpace from which the solution is obtained)
 PostProcessing {
   { Name Magnetodynamics2D_av; NameOfFormulation Magnetodynamics2D_av;
     PostQuantity {

Magnetodynamics/electromagnet.pro

@@ -55,7 +55,7 @@ Group {
 Function {
   DefineConstant[
     Current = {0.01, Name "Model parameters/Max. current [A]"},
-    murCore = {100, Visible !NonLinearCore,
+    murCore = {1000, Visible !NonLinearCore,
       Name "Model parameters/Mur core"}
     // by default, a linear ramp from 0 to 1 until 10 ms, then a constant value
     // of 1; any desired function can be specified
@@ -87,18 +87,24 @@ Function {
   nu[ Region[{Air, Ind, AirInf}] ]  = 1. / mu0;
 
   If(NonLinearCore)
-    Mat_core_b = {
-      0, 0.0902935, 0.194131, 0.293454, 0.397291,
-      0.496614, 0.600451, 0.690745, 0.794582, 0.893905,
-      0.997743, 1.09707, 1.19639, 1.2912, 1.39955,
-      1.49436, 1.54853, 1.59368, 1.63883, 1.693,
-      1.74266, 1.79684, 1.8465, 1.89616};
-    Mat_core_h = {
-      0, 4.336661847, 7.34801661, 9.72476491, 11.5026052,
-      13.1550314, 14.38338925, 15.72777105, 17.00435285, 17.9766304,
-      19.00397556, 20.54698273, 22.46623219, 25.98451136, 31.43230549,
-      42.54537896, 53.84673274, 66.63960707, 89.21748645, 138.2173146,
-      245.0317673, 454.3420034, 911.3848661, 1869.712483};
+    Mat_core_h() = {
+      0.0000e+00, 5.5023e+00, 1.1018e+01, 1.6562e+01, 2.2149e+01, 2.7798e+01, 3.3528e+01,
+      3.9363e+01, 4.5335e+01, 5.1479e+01, 5.7842e+01, 6.4481e+01, 7.1470e+01, 7.8906e+01,
+      8.6910e+01, 9.5644e+01, 1.0532e+02, 1.1620e+02, 1.2868e+02, 1.4322e+02, 1.6050e+02,
+      1.8139e+02, 2.0711e+02, 2.3932e+02, 2.8028e+02, 3.3314e+02, 4.0231e+02, 4.9395e+02,
+      6.1678e+02, 7.8320e+02, 1.0110e+03, 1.3257e+03, 1.7645e+03, 2.3819e+03, 3.2578e+03,
+      4.5110e+03, 6.3187e+03, 8.9478e+03, 1.2802e+04, 1.8500e+04, 2.6989e+04, 3.9739e+04,
+      5.9047e+04, 8.8520e+04, 1.3388e+05, 2.0425e+05, 3.1434e+05, 4.8796e+05, 7.6403e+05
+    } ;
+    Mat_core_b() = {
+      0.0000e+00, 5.0000e-02, 1.0000e-01, 1.5000e-01, 2.0000e-01, 2.5000e-01, 3.0000e-01,
+      3.5000e-01, 4.0000e-01, 4.5000e-01, 5.0000e-01, 5.5000e-01, 6.0000e-01, 6.5000e-01,
+      7.0000e-01, 7.5000e-01, 8.0000e-01, 8.5000e-01, 9.0000e-01, 9.5000e-01, 1.0000e+00,
+      1.0500e+00, 1.1000e+00, 1.1500e+00, 1.2000e+00, 1.2500e+00, 1.3000e+00, 1.3500e+00,
+      1.4000e+00, 1.4500e+00, 1.5000e+00, 1.5500e+00, 1.6000e+00, 1.6500e+00, 1.7000e+00,
+      1.7500e+00, 1.8000e+00, 1.8500e+00, 1.9000e+00, 1.9500e+00, 2.0000e+00, 2.0500e+00,
+      2.1000e+00, 2.1500e+00, 2.2000e+00, 2.2500e+00, 2.3000e+00, 2.3500e+00, 2.4000e+00
+    } ;
     Mat_core_b2 = Mat_core_b()^2;
     Mat_core_h2 = Mat_core_h()^2;
     Mat_core_nu = Mat_core_h() / Mat_core_b();
@@ -116,7 +122,7 @@ Function {
     nu[ Core ]  = 1. / (murCore * mu0);
   EndIf
 
-  sigma[ Core ] = 1e6 / 10;
+  sigma[ Core ] = 1e6 / 1000;
   sigma[ Ind ] = 5e7;
 
   Ns[ Ind ] = 1000 ; // number of turns in coil