Result becames strange in an edited model

Numerical methods and mathematical models of Elmer
Post Reply
yuanxun
Posts: 17
Joined: 20 Feb 2019, 14:04
Antispam: Yes

Result becames strange in an edited model

Post by yuanxun »

Hi

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
And the new one:

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
Finally,here're the SIF file and VTU file.
Thanks for helping.
Last edited by yuanxun on 10 May 2019, 09:44, edited 1 time in total.
yuanxun
Posts: 17
Joined: 20 Feb 2019, 14:04
Antispam: Yes

Re: Didn't compute well

Post by yuanxun »

I‘m sorry that I’ve forgot the SIF file.
Thanks a lot.
Attachments
origin.zip
(1.51 KiB) Downloaded 229 times
new.zip
(1.5 KiB) Downloaded 229 times
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: Result becames strange in an edited model

Post by kevinarden »

The target bodies and target boundaries changed numbers, was that intentional? Difficult to review without mesh files.

Kevin
mzenker
Posts: 1999
Joined: 07 Dec 2009, 11:49
Location: Germany

Re: Result becames strange in an edited model

Post by mzenker »

Hi,

there is a kind of permutation in bodies and materials between the two cases. I don't have the time to analyze it ion depth. However, comparison of the two files shows also that in ine first one, there is no emissivity defined for iron. So if your iron body has a radiating boundary, this might cause a problem.
BTW you can assign the same material (e.g. material 2, which is air) to more than one body - no need to define it several times.

Matthias
yuanxun
Posts: 17
Joined: 20 Feb 2019, 14:04
Antispam: Yes

Re: Result becames strange in an edited model

Post by yuanxun »

The iron in first one doesn't have any radiation boundaries.
I directly use the matearials defined in original file and make them correspond to bodies correctly.(And air is divided to different heat conductivity)
Maybe there're something wrong with model?
Thank you.

Yuan
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Result becames strange in an edited model

Post by raback »

Hi

One thing I noticed is that you give "emissivity" for "air". Try to take it away.

Now the radiation BCs require 1) emissity and 2) normal direction. You can nicely do this simply by given emissity to the solid parts which truly form the surfaces and radiate outwards. If you also give air the emissity the code cannot find the direction of the normal using the material properties. This logic may not be apparent for the user, but if both sides have "emissivity" then the normal direction is arbitrary.

The other way (the original one) to give the needed information is to define the "emissivity" in the BC section, and define the direction of the normal by the "normal target boundary" keyword. This is more cumbersome, in my opinion.

-Peter
yuanxun
Posts: 17
Joined: 20 Feb 2019, 14:04
Antispam: Yes

Re: Result becames strange in an edited model

Post by yuanxun »

Hi Peter

If I take the emissivity away."WARNING:: GetParentMatProp: Property Emissivity not in either parents!"occurred.
Then I set the emissivity in BCs,the code in solver log becames like this:

Code: Select all

ELMER SOLVER (v 8.3) STARTED AT: 2019/05/15 16:19:16
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:     0.4820 (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:

HeatSolve: 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)    0.98    0.98
HeatSolve: iter:    1 Solve:    (s)    0.02    0.02
HeatSolve:  Result Norm   :    291.78228440658143
HeatSolve:  Relative Change :    2.0000000000000000
HeatSolve: 
HeatSolve: 
HeatSolve: -------------------------------------
HeatSolve:  TEMPERATURE ITERATION           2
HeatSolve: -------------------------------------
HeatSolve: 
HeatSolve: Starting Assembly...
HeatSolve: Assembly:

HeatSolve: 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.01    1.99
HeatSolve: iter:    2 Solve:    (s)    0.01    0.03
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:

HeatSolve: 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)    1.01    3.00
HeatSolve: iter:    3 Solve:    (s)    0.01    0.04
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):         3.93        3.93
ELMER SOLVER FINISHED AT: 2019/05/15 16:19:19
Post Reply