Correct geometry size for Inno's test

dG3D/benchmarks/Inno_EMT_test/InnoElecMagTherm.py

@@ -27,7 +27,7 @@ mag_mu_0 = 4.e-7 * math.pi # vacuum magnetic permeability
 mag_r_smp = 20.0 # relative mag permeability of SMP
 mag_mu_x = mag_mu_y = mag_mu_z = mag_r_smp * mag_mu_0 # mag permeability smp 
 magpot0_x = magpot0_y = magpot0_z = 0.0 # initial magnetic potential
-Irms = 10. # (Ampere)
+Irms = 2000. # (Ampere)
 freq = 1000. # (Hz)
 nTurnsCoil = 1000
 coilLx = coilLy = coilLz = 1. # m
@@ -73,14 +73,14 @@ mag_mu_x_vac = mag_mu_y_vac = mag_mu_z_vac = mag_r_vac * mag_mu_0 # mag permeabi
 magpot0_x_vac = magpot0_y_vac = magpot0_z_vac = 0.0 # initial magnetic potential
 
 # geometry
-geofile="smp.geo" # name of mesh file
-meshfile="smp.msh" # name of mesh file
+geofile="hollow_smp.geo" # name of mesh file
+meshfile="hollow_smp.msh" # name of mesh file
 
 # solver
 sol = 2  # Gmm=0 (default) Taucs=1 PETsc=2
 soltype = 1 # StaticLinear=0 (default) StaticNonLinear=1
-nstep = 50 # number of step (used only if soltype=1)
-ftime =1.0/freq # Final time (used only if soltype=1)
+nstep = 100 # number of step (used only if soltype=1)
+ftime =5.0/freq # Final time (used only if soltype=1)
 tol=1.e-3   # relative tolerance for NR scheme (used only if soltype=1)
 nstepArch=1 # Number of step between 2 archiving (used only if soltype=1)
 fullDg = 0 #O = CG, 1 = DG
@@ -234,6 +234,7 @@ mysolver.curlDisplacementBC("Edge",OutputPhysical,5,0.0)
 #mysolver.internalPointBuildView("sig_yz",IPField.SIG_YZ, 1, 1)
 #mysolver.internalPointBuildView("sig_xz",IPField.SIG_XZ, 1, 1)
 mysolver.internalPointBuildView("temperature",IPField.TEMPERATURE, 1, 1)
+mysolver.internalPointBuildView("w_T",IPField.THERMALSOURCE, 1, 1)
 mysolver.internalPointBuildView("jyx",IPField.THERMALFLUX_X, 1, 1)
 mysolver.internalPointBuildView("jyy",IPField.THERMALFLUX_Y, 1, 1)
 mysolver.internalPointBuildView("jyz",IPField.THERMALFLUX_Z, 1, 1)

dG3D/benchmarks/Inno_EMT_test/hollow_smp_data.geo

@@ -21,7 +21,7 @@ colorpp = "Ivory";
 
 
 // Smp
-RintSMP = 0.5*cm;
+RintSMP = 0.95*cm;
 RextSMP = 1*cm;
 H = 2*cm;
 
@@ -36,7 +36,7 @@ DefineConstant[
   wcoil = {3*cm, Name StrCat[pp, "1Coil width [m]"], Highlight Str[colorpp]}
   RintCoil = {4.5*cm, Name StrCat[pp, "2Coil inner radius [m]"], Highlight Str[colorpp]}
   RextCoil = {7.5*cm, Name StrCat[pp, "3Coil outer radius [m]"], Highlight Str[colorpp]}
-  RintSMP  = {0.5*cm, Name StrCat[pp, "4Sample inner radius [m]"], Highlight Str[colorpp]}
+  RintSMP  = {0.95*cm, Name StrCat[pp, "4Sample inner radius [m]"], Highlight Str[colorpp]}
   RextSMP  = {1*cm, Name StrCat[pp, "5Sample outer radius [m]"], Highlight Str[colorpp]}
   H     = {2*cm, Name StrCat[pp, "6Sample height [m]"], Highlight Str[colorpp]}
 ];