Unable to get reasonable result for InductionHeating example in elmer-elmag

[atul]

Hello,

I ran the code from github repo: (below location) for the induction heating example.

https://github.com/ElmerCSC/elmer-elmag ... ionHeating.

The coil current e is of magnitude 350 to 50000 (units not known) and all quantities when viewed in ParaView are '0'

Below is the log file during the simulation run: Please advise on what could be wrong and how to correct it.

Thanks,
Atul

ELMER SOLVER (v 9.0) STARTED AT: 2023/09/06 09:24:08
ParCommInit: Initialize #PEs: 1

MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Version: 9.0 (Rev: Release, Compiled: 2022-08-02)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: Lua interpreter linked in.
MAIN: =============================================================

MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: case.sif

LoadInputFile: Scanning input file: case.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file

LoadInputFile: Loading input file: case.sif
LoadInputFile: Reading base load of sif file

LoadInputFile: Number of BCs: 4
LoadInputFile: Number of Body Forces: 2
LoadInputFile: Number of Initial Conditions: 0
LoadInputFile: Number of Materials: 3
LoadInputFile: Number of Equations: 2
LoadInputFile: Number of Solvers: 5
LoadInputFile: Number of Bodies: 3
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by area
ListTagKeywords: Setting weight for keywords!
ListTagKeywords: No parameters width suffix: normalize by volume

Loading user function library: [CoilSolver]...[CoilSolver_Init0]
Loading user function library: [MagnetoDynamics]...[WhitneyAVHarmonicSolver_Init0]
WhitneyHarmonicSolver_Init0: Setting element type to: "n:1 e:1"
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init0]
Loading user function library: [HeatSolve]...[HeatSolver_Init0]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init0]
ElmerAsciiMesh: Base mesh name: ./.

PermuteNodeNumbering: Performing node mapping

MapBodiesAndBCs: Minimum initial body index: 1
MapBodiesAndBCs: Maximum initial body index: 3

MapBodiesAndBCs: Minimum initial boundary index: 1
MapBodiesAndBCs: Maximum initial boundary index: 5

LoadMesh: Elapsed REAL time: 3.0200 (s)
MeshStabParams: Computing stabilization parameters

MeshStabParams: Elapsed REAL time: 0.1000 (s)
MAIN: -------------------------------------

Loading user function library: [CoilSolver]...[CoilSolver_Init]
Loading user function library: [CoilSolver]...[CoilSolver_bulk]
Loading user function library: [CoilSolver]...[CoilSolver]
OptimizeBandwidth: ---------------------------------------------------------

OptimizeBandwidth: Computing matrix structure for: coilsolver...done.
OptimizeBandwidth: Half bandwidth without optimization: 7197
OptimizeBandwidth:

OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 60
OptimizeBandwidth: ---------------------------------------------------------

Loading user function library: [MagnetoDynamics]...[WhitneyAVHarmonicSolver_Init]
Loading user function library: [MagnetoDynamics]...[WhitneyAVHarmonicSolver_bulk]
Loading user function library: [MagnetoDynamics]...[WhitneyAVHarmonicSolver]
OptimizeBandwidth: ---------------------------------------------------------

OptimizeBandwidth: Computing matrix structure for: mgdynamics...done.
OptimizeBandwidth: Half bandwidth without optimization: 309313
OptimizeBandwidth:

OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 25454
OptimizeBandwidth: ---------------------------------------------------------

Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init]
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_bulk]
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields]
OptimizeBandwidth: ---------------------------------------------------------

OptimizeBandwidth: Computing matrix structure for: mgdynamicscalc...done.
OptimizeBandwidth: Half bandwidth without optimization: 38961
OptimizeBandwidth:

OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 4074
OptimizeBandwidth: ---------------------------------------------------------

Loading user function library: [HeatSolve]...[HeatSolver_Init]
Loading user function library: [HeatSolve]...[HeatSolver_bulk]
Loading user function library: [HeatSolve]...[HeatSolver]
OptimizeBandwidth: ---------------------------------------------------------

OptimizeBandwidth: Computing matrix structure for: heat equation...done.
OptimizeBandwidth: Half bandwidth without optimization: 38961
OptimizeBandwidth:

OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 4074
OptimizeBandwidth: ---------------------------------------------------------

Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver]
ElmerSolver: Number of timesteps to be saved: 1

ListTagCount: Found number of normalized keywords: 0

CoilSolver: --------------------------------------
CoilSolver: Solving current distribution in a coil
CoilSolver: --------------------------------------
CoilSolver: Assuming that all coils are open!

CoilSolver: Coil system consists of 1 coils
CoilSolver: Set1 : 21 +nodes and 21 -nodes
CoilSolver: Computing the dummy potential field

DefaultFinishBulkAssembly: Saving bulk values for: coilsolver

DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 42
CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 115948
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 150
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.01

20 0.1218E-02

40 0.2221E-02

60 0.1193E-02

80 0.7458E-03

100 0.2065E-03

120 0.1909E-04

140 0.6552E-11
Idrs: 140 0.6552E-11

SolveLinearSystem: Calculating nodal loads
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.57815507 2.0000000 ) :: coilsolver
CoilSolver: Discrepancy of start and end coil currents: 1.3991E-09
CoilSolver: Initial coil current for coil 1: 5.0004E-05
CoilSolver: Coil potential multiplier: 1.9998E+06
CoilSolver: Computing current component: 1

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 115948
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 150
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.00

Idrs: 6 0.8842E-12
CoilSolver: Computing current component: 2

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 115948
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 150
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.00

Idrs: 7 0.1635E-12
CoilSolver: Computing current component: 3

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 115948
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 150
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.00

Idrs: 6 0.8070E-12
CoilSolver: Minimum current components: -1.4678E+06 -9.5086E+04 -1.5242E+06
CoilSolver: Maximum current components: 1.4777E+06 1.8718E+04 1.5009E+06
CoilSolver: Normalizing current density (as for stranded coil)

CoilSolver: Average current density: 1.3320E+06

CoilSolver: Finding fixing potential for elemental current density!

DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 1
CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 115948
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 150
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.01

20 0.1195E+01

40 0.1098E+01

60 0.8930E+00

80 0.7854E+00

100 0.8623E+01

120 0.1019E+00

140 0.3248E-06

Idrs: 157 0.4988E-11

SolveLinearSystem: Calculating nodal loads
CoilSolver: Fixing elemental current density to be divergence free!

CoilSolver: All done
CoilSolver: --------------------------------------
Loading user function library: [CoilSolver]...[CoilSolver_post]
MAIN:
MAIN: -------------------------------------
MAIN: Steady state iteration: 1
MAIN: -------------------------------------
MAIN:

DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 10458

SolveLinearSystem: Assuming complex valued linear system

SolveLinearSystem: Solution trivially zero!
Loading user function library: [MagnetoDynamics]...[WhitneyAVHarmonicSolver_post]
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.0000000 0.0000000 ) :: mgdynamics

MagnetoDynamicsCalcFields: ------------------------------
MagnetoDynamicsCalcFields: Computing postprocessed fields
MagnetoDynamicsCalcFields: Using potential variable: av

MagnetoDynamicsCalcFields: Eddy current power: 0.000000E+00
MagnetoDynamicsCalcFields: ElectroMagnetic Field Energy: 0.000000E+00
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_post]
ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.0000000 0.0000000 ) :: mgdynamicscalc

HeatSolver: -------------------------------------------
HeatSolver: Solving the energy equation for temperature

HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 1
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.31

20 0.6073E-06

Idrs: 25 0.4563E-08

ComputeChange: NS (ITER=1) (NRM,RELC): ( 225.46244 2.0000000 ) :: heat equation
HeatSolve: iter: 1 Assembly: (s) 3.05 3.05
HeatSolve: iter: 1 Solve: (s) 0.44 0.44
HeatSolve: Result Norm : 225.46243791746545
HeatSolve: Relative Change : 2.0000000000000000
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 2
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.05

20 0.5199E-07

Idrs: 24 0.4348E-08

ComputeChange: NS (ITER=2) (NRM,RELC): ( 272.22185 0.18790791 ) :: heat equation
HeatSolve: iter: 2 Assembly: (s) 3.01 6.06
HeatSolve: iter: 2 Solve: (s) 0.17 0.61
HeatSolve: Result Norm : 272.22184531911466
HeatSolve: Relative Change : 0.18790791261303133
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 3
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.04

20 0.1580E-07

Idrs: 21 0.1733E-08

ComputeChange: NS (ITER=3) (NRM,RELC): ( 286.83067 0.52262787E-01 ) :: heat equation
HeatSolve: iter: 3 Assembly: (s) 3.04 9.10
HeatSolve: iter: 3 Solve: (s) 0.16 0.77
HeatSolve: Result Norm : 286.83066644360690
HeatSolve: Relative Change : 5.2262786829916447E-002
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 4
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.05

Idrs: 19 0.5405E-08

ComputeChange: NS (ITER=4) (NRM,RELC): ( 291.23870 0.15250890E-01 ) :: heat equation
HeatSolve: iter: 4 Assembly: (s) 3.02 12.12
HeatSolve: iter: 4 Solve: (s) 0.15 0.92
HeatSolve: Result Norm : 291.23870248768640
HeatSolve: Relative Change : 1.5250889533305868E-002
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 5
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.04

Idrs: 16 0.7605E-08

ComputeChange: NS (ITER=5) (NRM,RELC): ( 292.58495 0.46118377E-02 ) :: heat equation
HeatSolve: iter: 5 Assembly: (s) 3.02 15.14
HeatSolve: iter: 5 Solve: (s) 0.14 1.06
HeatSolve: Result Norm : 292.58495246105491
HeatSolve: Relative Change : 4.6118377080377583E-003
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 6
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.04

Idrs: 19 0.8605E-09

ComputeChange: NS (ITER=6) (NRM,RELC): ( 292.97972 0.13483273E-02 ) :: heat equation
HeatSolve: iter: 6 Assembly: (s) 3.03 18.17
HeatSolve: iter: 6 Solve: (s) 0.15 1.21
HeatSolve: Result Norm : 292.97971887644917
HeatSolve: Relative Change : 1.3483272974531119E-003
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 7
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.04

Idrs: 6 0.7136E-08

ComputeChange: NS (ITER=7) (NRM,RELC): ( 293.06733 0.29899824E-03 ) :: heat equation
HeatSolve: iter: 7 Assembly: (s) 3.03 21.20
HeatSolve: iter: 7 Solve: (s) 0.08 1.29
HeatSolve: Result Norm : 293.06733239512988
HeatSolve: Relative Change : 2.9899824080886317E-004
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 8
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...

HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.05

Idrs: 5 0.6811E-08

ComputeChange: NS (ITER=8) (NRM,RELC): ( 293.09353 0.89376389E-04 ) :: heat equation
HeatSolve: iter: 8 Assembly: (s) 3.01 24.21
HeatSolve: iter: 8 Solve: (s) 0.08 1.37
HeatSolve: Result Norm : 293.09352686572930
HeatSolve: Relative Change : 8.9376389383801250E-005
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 9
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ...Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.05

Idrs: 6 0.6815E-08
ComputeChange: NS (ITER=9) (NRM,RELC): ( 293.10145 0.27031361E-04 ) :: heat equation
HeatSolve: iter: 9 Assembly: (s) 3.13 27.34
HeatSolve: iter: 9 Solve: (s) 0.09 1.46
HeatSolve: Result Norm : 293.10144968971338
HeatSolve: Relative Change : 2.7031360898484826E-005
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 10
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ....Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.05

Idrs: 4 0.8842E-08

ComputeChange: NS (ITER=10) (NRM,RELC): ( 293.10365 0.74924500E-05 ) :: heat equation
HeatSolve: iter: 10 Assembly: (s) 4.11 31.46
HeatSolve: iter: 10 Solve: (s) 0.07 1.53
HeatSolve: Result Norm : 293.10364574590756
HeatSolve: Relative Change : 7.4924500359304606E-006
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 11
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.05

Idrs: 5 0.7598E-08

ComputeChange: NS (ITER=11) (NRM,RELC): ( 293.10516 0.51773218E-05 ) :: heat equation
HeatSolve: iter: 11 Assembly: (s) 3.01 34.47
HeatSolve: iter: 11 Solve: (s) 0.08 1.61
HeatSolve: Result Norm : 293.10516324173227
HeatSolve: Relative Change : 5.1773218056379774E-006
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 12
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

: ..Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions

EnforceDirichletConditions: Dirichlet conditions enforced for dofs: 7849

CRS_IncompleteLU: ILU(1) (Real), Performing Factorization:

CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 1205263
CRS_IncompleteLU: ILU(1) (Real), filling (%) : 207
CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.05

Idrs: 2 0.8270E-08

ComputeChange: NS (ITER=12) (NRM,RELC): ( 293.10506 0.34622690E-06 ) :: heat equation
HeatSolve: iter: 12 Assembly: (s) 3.03 37.50
HeatSolve: iter: 12 Solve: (s) 0.07 1.68
HeatSolve: Result Norm : 293.10506176085812
HeatSolve: Relative Change : 3.4622689889019036E-007
Loading user function library: [HeatSolve]...[HeatSolver_post]
ComputeChange: SS (ITER=1) (NRM,RELC): ( 293.10506 2.0000000 ) :: heat equation

ResultOutputSolver: -------------------------------------
ResultOutputSolver: Saving with prefix: coilsolver
ResultOutputSolver: Working on mesh: .
ResultOutputSolver: Dimension of mesh is: 3
ResultOutputSolver: Creating list for saving - if not present
CreateListForSaving: Field Variables for Saving
CreateListForSaving: Scalar Field 1: coilpot
CreateListForSaving: Scalar Field 2: av re
CreateListForSaving: Scalar Field 3: av im
CreateListForSaving: Scalar Field 4: joule heating e
CreateListForSaving: Scalar Field 5: temperature
CreateListForSaving: Vector Field 1: coilcurrent e
CreateListForSaving: Vector Field 2: coilcurrent
CreateListForSaving: Vector Field 3: magnetic flux density re e
CreateListForSaving: Vector Field 4: magnetic flux density im e
CreateListForSaving: Vector Field 5: magnetic field strength re e
CreateListForSaving: Vector Field 6: magnetic field strength im e
CreateListForSaving: Vector Field 7: current density re e
CreateListForSaving: Vector Field 8: current density im e
CreateListForSaving: Vector Field 9: electric field re e
CreateListForSaving: Vector Field 10: electric field im e
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Saving results in VTK XML format with prefix: coilsolver
VtuOutputSolver: Saving number of partitions: 1

VtuOutputSolver: Writing the vtu file: results/coilsolver_inductor_t0001.vtu

VtuOutputSolver: Writing the vtu file: results/coilsolver_cylinder_t0001.vtu

VtuOutputSolver: Writing the vtu file: results/coilsolver_surrounding_t0001.vtu

ResultOutputSolver: -------------------------------------
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post]
CompareToReferenceSolution: Solver 1 PASSED: Norm = 5.78155071E-01 RefNorm = 5.78155074E-01
CompareToReferenceSolution: Relative Error to reference norm: 5.676369E-09
WARNING:: CompareToReferenceSolution: Solver 2 FAILED: Norm = 0.00000000E+00 RefNorm = 1.45863590E-02
CompareToReferenceSolution: Relative Error to reference norm: 1.000000E+00
WARNING:: CompareToReferenceSolution: Solver 4 FAILED: Norm = 2.93105062E+02 RefNorm = 6.58345530E+02
CompareToReferenceSolution: Relative Error to reference norm: 5.547854E-01
ReloadInputFile: Realoading input file
LoadInputFile: Loading input file:
WARNING:: CompareToReferenceSolution: FAILED 2 tests out of 3!

ElmerSolver: *** Elmer Solver: ALL DONE ***

ElmerSolver: The end
SOLVER TOTAL TIME(CPU,REAL): 75.97 75.97

[Rich_B]

Hello,

This line in the solver log:

MAIN: Version: 9.0 (Rev: Release, Compiled: 2022-08-02)

The version of Elmer that you have installed is over a year old. Try updating to the latest version and try running again.

Tip: save the solver log file to a text file and attach the text file to your post.

Rich.

[atul]

Hi Rich,

Thanks for your suggestions. I just updated Elmer to the latest and ran the Induction Heating example again and now the results seem okay. I have attached the log file (thanks for the tip:)
I have a hard time interpreting the actual values. For example coilcurrent is of the order of 10^6, which I assume has units of Amp/m^2. Is that correct?
The potential is of the order of 0.2, which I assume is in volts.
What are the units used for the magnetic flux?
Thanks and look forward to your response.

Atul

[raback]

Hi, If input parameters are in SI units and mesh as well (m) then all results will also be in SI units. -Peter