based on your example, changed some naming to fit my .sif, did some flux measurement tests. Complete project is attached below, with 'make' (defaults to 16 CORES for fine mesh). Thing to explore is which per-conditioners would in this case improve the convergence speed. With BiCGStabL polynomial degrees higher than 5 internal error arises. But this is not the major issue now.
The two flux measurement areas of the Torus appear to be symmetrical (also reported by Elmer and is ofc slightly mesh dependent), however calculating the average should, according to the orientation (also a good question how to define orientation of a surface), provide either double value so around 0.07 T or ->0 T, but it's just something.
It is also unclear to me how to report average flux for each boundary separately; now I've done more runs to enable one, or the other, or both, and the last I've changed mesh grid to finest (12 from curvatures, in .geo.opt currently in the zip project).
Expected flux density in this core with u_r = 3000 and wire current of 1.256 A should be ~0.037 T
Code: Select all
Variables in columns of matrix: results/coreWire.dat
1: simulation: cpu time (s)
2: simulation: real time (s)
3: res: eddy current power
4: res: electromagnetic field energy
5: res: magnetic flux average
6: res: magnetic flux density average
7: res: magnetic flux area
2.052951300000E+001 1.744762611389E+001 5.964347131448E-004 7.683346874482E-007 -1.574121712043E-007 -2.915040207464E-003 5.400000000045E-005
2.155061100000E+001 1.811320304871E+001 5.964347146686E-004 7.692430960868E-007 5.177143995973E-007 1.917460739233E-002 2.700000000023E-005
2.167175200000E+001 1.824741506577E+001 5.964347146686E-004 7.692430960868E-007 5.177143995973E-007 1.917460739233E-002 2.700000000023E-005
2.172054300000E+001 1.825896120071E+001 5.964347146686E-004 7.692430960868E-007 -1.575314886442E-007 -2.917249789682E-003 5.400000000045E-005
2.181652600000E+001 1.827390098572E+001 5.964347146686E-004 7.692430960868E-007 -6.752458882415E-007 -2.500910697170E-002 2.700000000023E-005
2.193741600000E+001 1.868077683449E+001 5.964347146686E-004 7.692430960868E-007 -1.575314886442E-007 -2.917249789682E-003 5.400000000045E-005
5.879800900000E+001 4.814381885529E+001 5.964347154959E-004 7.697399605983E-007 -1.575966774580E-007 -2.918456989939E-003 5.400000000045E-005
1.719310000000E+001 1.474894309044E+001 5.964347148393E-004 7.693471671225E-007 -1.575451448709E-007 -2.917502682770E-003 5.400000000045E-005
1.730557200000E+001 1.490097379684E+001 5.964347148393E-004 7.693471671225E-007 -1.575451448709E-007 -2.917502682770E-003 5.400000000045E-005
8.483373020000E+002 7.520740270615E+002 6.236439989278E-004 8.431919786786E-007 -2.571142490298E-007 -4.599135063742E-003 5.590491374276E-005
Uros