I'd like to do a simulation about thermal transfer in a wall.Once I changed the thickness of a kind of material,the computed result becames strange.
Condition and material setting is the same as original model,but there's something different in solver log.
Here're the code of the original one:
Code: Select all
ELMER SOLVER (v 8.3) STARTED AT: 2019/05/10 14:22:46
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: 8.3 (Rev: 432ee0dd, Compiled: 2018-12-18)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: Lua interpreted linked in.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: case.sif
LoadInputFile: Scanning input file: case.sif
LoadInputFile: Loading input file: case.sif
Loading user function library: [HeatSolve]...[HeatSolver_Init0]
LoadMesh: Base mesh name: ./.
LoadMesh: Elapsed REAL time: 1.2420 (s)
MAIN: -------------------------------------
AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: case
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: 21249
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 2041
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver]
MAIN:
MAIN: -------------------------------------
MAIN: Steady state iteration: 1
MAIN: -------------------------------------
MAIN:
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: heat equation
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
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 303945
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.02
10 0.2233E-03
20 0.8943E-04
30 0.4634E-04
40 0.2603E-04
50 0.5829E-02
60 0.1515E-04
70 0.1099E-04
80 0.9089E-05
90 0.4360E-04
100 0.7990E-05
110 0.6999E-05
120 0.7528E-05
130 0.3965E-04
140 0.5889E-05
150 0.5674E-05
160 0.4998E-05
170 0.4944E-05
180 0.4485E-05
190 0.6379E-05
200 0.4863E-05
210 0.1288E-04
220 0.4476E-05
230 0.4762E-05
240 0.4735E-05
250 0.4744E-05
260 0.4826E-05
270 0.5378E-05
280 0.5976E-05
290 0.6000E-05
300 0.5759E-05
310 0.7837E-05
320 0.6899E-05
330 0.1927E-03
340 0.2972E-04
350 0.2289E-04
360 0.6083E-05
370 0.1217E-04
380 0.1086E-04
390 0.1279E-04
400 0.5778E-05
410 0.1977E-04
420 0.6073E-04
430 0.6894E-05
440 0.1500E-04
450 0.9086E-05
460 0.8222E-05
470 0.1007E-04
480 0.1073E-04
490 0.9350E-05
500 0.1481E-04
501 0.1481E-04
NUMERICAL WARNING:: IterSolve: Too many iterations was needed.
ComputeChange: NS (ITER=1) (NRM,RELC): ( 277.84530 2.0000000 ) :: heat equation
HeatSolve: iter: 1 Assembly: (s) 1.83 1.83
HeatSolve: iter: 1 Solve: (s) 1.79 1.79
HeatSolve: Result Norm : 277.84529657907882
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
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 303945
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00
10 0.7313E-04
20 0.2223E-04
30 0.1214E-04
40 0.1534E-04
50 0.1832E-04
60 0.4477E-04
70 0.2298E-04
80 0.1770E-04
90 0.1166E-04
100 0.1072E-04
110 0.6969E-05
120 0.6870E-05
130 0.9831E-05
140 0.1056E-04
150 0.7407E-05
160 0.7001E-05
170 0.7582E-05
180 0.6518E-05
190 0.3861E-05
200 0.3510E-05
210 0.2995E-05
220 0.2987E-05
230 0.2888E-05
240 0.2649E-05
250 0.3034E-05
260 0.2289E-05
270 0.1250E-05
280 0.1264E-05
290 0.1342E-05
300 0.2417E-05
310 0.2528E-05
320 0.9870E-06
330 0.9702E-06
340 0.7179E-06
350 0.7172E-06
360 0.7085E-06
370 0.7335E-06
380 0.6136E-06
390 0.6935E-06
400 0.5606E-06
410 0.5159E-06
420 0.4920E-06
430 0.4035E-06
440 0.4116E-06
450 0.4737E-06
460 0.1206E-06
470 0.1495E-06
480 0.1955E-06
490 0.1421E-06
500 0.1588E-06
501 0.1588E-06
NUMERICAL WARNING:: IterSolve: Too many iterations was needed.
ComputeChange: NS (ITER=2) (NRM,RELC): ( 283.88470 0.21502868E-01 ) :: heat equation
HeatSolve: iter: 2 Assembly: (s) 1.89 3.72
HeatSolve: iter: 2 Solve: (s) 2.02 3.80
HeatSolve: Result Norm : 283.88469968026465
HeatSolve: Relative Change : 2.1502868429328144E-002
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 3
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:
: .Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 303945
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.01
10 0.3327E-06
20 0.1499E-06
30 0.3071E-06
40 0.1080E-06
50 0.9961E-07
60 0.1123E-06
70 0.1452E-06
80 0.3841E-06
90 0.4548E-06
100 0.7227E-06
110 0.1622E-06
120 0.1586E-06
130 0.4154E-06
140 0.1351E-06
150 0.1433E-06
160 0.1199E-06
170 0.8827E-07
180 0.2137E-05
190 0.1056E-06
200 0.2088E-06
210 0.9743E-07
220 0.1001E-06
230 0.1115E-06
240 0.1217E-06
250 0.1201E-06
260 0.7041E-07
270 0.7779E-07
280 0.8615E-07
290 0.7591E-07
300 0.8867E-07
310 0.1005E-06
320 0.7690E-07
330 0.6563E-07
340 0.6104E-07
350 0.6482E-07
360 0.6117E-07
370 0.6208E-07
380 0.5977E-07
390 0.8339E-07
400 0.6619E-07
410 0.5663E-07
420 0.5462E-07
430 0.3593E-07
440 0.8073E-06
450 0.3362E-07
460 0.3708E-07
470 0.3143E-07
480 0.2233E-07
490 0.1362E-07
500 0.2153E-07
501 0.2153E-07
NUMERICAL WARNING:: IterSolve: Too many iterations was needed.
ComputeChange: NS (ITER=3) (NRM,RELC): ( 284.01365 0.45414416E-03 ) :: heat equation
HeatSolve: iter: 3 Assembly: (s) 2.13 5.85
HeatSolve: iter: 3 Solve: (s) 2.16 5.97
HeatSolve: Result Norm : 284.01365354086988
HeatSolve: Relative Change : 4.5414416109434486E-004
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 4
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:
: .Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 303945
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00
10 0.2156E-06
20 0.5896E-08
30 0.1131E-07
40 0.1625E-07
50 0.1697E-07
60 0.1935E-07
70 0.5233E-07
80 0.5837E-07
90 0.1601E-06
100 0.4175E-06
110 0.6030E-07
120 0.5179E-07
130 0.6913E-07
140 0.8440E-07
150 0.1290E-06
160 0.1745E-06
170 0.6013E-08
180 0.6675E-08
190 0.1184E-07
200 0.4307E-08
210 0.3041E-08
220 0.3289E-08
230 0.2978E-08
240 0.2963E-08
250 0.2975E-08
260 0.2996E-08
270 0.1263E-08
280 0.2018E-08
290 0.1964E-08
300 0.1786E-08
310 0.1302E-08
320 0.1026E-08
330 0.1134E-08
340 0.1122E-08
350 0.9828E-09
360 0.8479E-09
370 0.8380E-09
380 0.8181E-09
390 0.8083E-09
400 0.8010E-09
410 0.7201E-09
420 0.6857E-09
430 0.7012E-09
440 0.6289E-09
450 0.6707E-09
460 0.3499E-09
470 0.7081E-09
480 0.1732E-09
490 0.4462E-10
490 0.4462E-10
ComputeChange: NS (ITER=4) (NRM,RELC): ( 284.02296 0.32751031E-04 ) :: heat equation
HeatSolve: iter: 4 Assembly: (s) 1.96 7.81
HeatSolve: iter: 4 Solve: (s) 1.95 7.92
HeatSolve: Result Norm : 284.02295543311391
HeatSolve: Relative Change : 3.2751030821153083E-005
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 5
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:
: .Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 303945
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.01
5 0.6560E-10
ComputeChange: NS (ITER=5) (NRM,RELC): ( 284.02296 0.37760813E-09 ) :: heat equation
HeatSolve: iter: 5 Assembly: (s) 1.98 9.79
HeatSolve: iter: 5 Solve: (s) 0.08 8.00
HeatSolve: Result Norm : 284.02295554036328
HeatSolve: Relative Change : 3.7760813098804609E-010
Loading user function library: [HeatSolve]...[HeatSolver_post]
ComputeChange: SS (ITER=1) (NRM,RELC): ( 284.02296 2.0000000 ) :: heat equation
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: internalvtuoutputsolver
ResultOutputSolver: -------------------------------------
ResultOutputSolve: Saving with prefix: case
ResultOutputSolver: Creating list for saving - if not present
CreateListForSaving: Field Variables for Saving
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Saving results in VTK XML format with prefix: case
VtuOutputSolver: Saving number of partitions: 1
ResultOutputSolver: -------------------------------------
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post]
ElmerSolver: *** Elmer Solver: ALL DONE ***
ElmerSolver: The end
SOLVER TOTAL TIME(CPU,REAL): 19.86 19.86
ELMER SOLVER FINISHED AT: 2019/05/10 14:23:06
Code: Select all
ELMER SOLVER (v 8.3) STARTED AT: 2019/05/10 14:25:39
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: 8.3 (Rev: 432ee0dd, Compiled: 2018-12-18)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: Lua interpreted linked in.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: case.sif
LoadInputFile: Scanning input file: case.sif
LoadInputFile: Loading input file: case.sif
Loading user function library: [HeatSolve]...[HeatSolver_Init0]
LoadMesh: Base mesh name: ./.
LoadMesh: Elapsed REAL time: 1.3600 (s)
MAIN: -------------------------------------
AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: case
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: 20431
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 1915
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver]
MAIN:
MAIN: -------------------------------------
MAIN: Steady state iteration: 1
MAIN: -------------------------------------
MAIN:
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: heat equation
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
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 287681
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.01
7 0.6821E-10
ComputeChange: NS (ITER=1) (NRM,RELC): ( 291.78228 2.0000000 ) :: heat equation
HeatSolve: iter: 1 Assembly: (s) 1.90 1.90
HeatSolve: iter: 1 Solve: (s) 0.04 0.04
HeatSolve: Result Norm : 291.78228440658143
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
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 287681
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00
2 0.1392E-10
ComputeChange: NS (ITER=2) (NRM,RELC): ( 291.78253 0.82785263E-06 ) :: heat equation
HeatSolve: iter: 2 Assembly: (s) 1.75 3.65
HeatSolve: iter: 2 Solve: (s) 0.00 0.04
HeatSolve: Result Norm : 291.78252595941331
HeatSolve: Relative Change : 8.2785263124128004E-007
HeatSolve:
HeatSolve:
HeatSolve: -------------------------------------
HeatSolve: TEMPERATURE ITERATION 3
HeatSolve: -------------------------------------
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly:
: .Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 287681
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.01
1 0.1178E-11
ComputeChange: NS (ITER=3) (NRM,RELC): ( 291.78253 0.21821426E-09 ) :: heat equation
HeatSolve: iter: 3 Assembly: (s) 2.04 5.68
HeatSolve: iter: 3 Solve: (s) 0.01 0.05
HeatSolve: Result Norm : 291.78252602308442
HeatSolve: Relative Change : 2.1821426439966022E-010
Loading user function library: [HeatSolve]...[HeatSolver_post]
ComputeChange: SS (ITER=1) (NRM,RELC): ( 291.78253 2.0000000 ) :: heat equation
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: internalvtuoutputsolver
ResultOutputSolver: -------------------------------------
ResultOutputSolve: Saving with prefix: case
ResultOutputSolver: Creating list for saving - if not present
CreateListForSaving: Field Variables for Saving
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Saving results in VTK XML format with prefix: case
VtuOutputSolver: Saving number of partitions: 1
ResultOutputSolver: -------------------------------------
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post]
ElmerSolver: *** Elmer Solver: ALL DONE ***
ElmerSolver: The end
SOLVER TOTAL TIME(CPU,REAL): 7.86 7.86
ELMER SOLVER FINISHED AT: 2019/05/10 14:25:47
Thanks for helping.