The solver just stops after calculating the gebhard factors in this case as well as for the crystal alone (radiating into the vacuum), but shows no error messages, except for some Matrix access errors way earlier, but these were often present in simulations that ran through.
Code: Select all
ELMER SOLVER (v 9.0) STARTED AT: 2024/02/27 20:56:43
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: 2024-02-22)
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: 0
LoadInputFile: Number of Initial Conditions: 1
LoadInputFile: Number of Materials: 2
LoadInputFile: Number of Equations: 1
LoadInputFile: Number of Solvers: 1
LoadInputFile: Number of Bodies: 2
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
ElmerAsciiMesh: Base mesh name: ./.
MapCoordinates: Scaling coordinates: 1.000E+00 1.000E+00 1.000E+00
LoadMesh: Elapsed REAL time: 0.1270 (s)
MAIN: -------------------------------------
AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: case
RadiationFactors: ----------------------------------------------------
RadiationFactors: Computing radiation factors for heat transfer
RadiationFactors: ----------------------------------------------------
RadiationFactors: Total number of Radiation Surfaces 4480 out of 9728
RadiationFactors: Temporarily updating the mesh.nodes file!
ViewFactors:
ViewFactors: ==================================================
ViewFactors: E L M E R V I E W F A C T O R S, W E L C O M E
ViewFactors: ==================================================
ViewFactors:
ViewFactors:
ViewFactors: Reading Model...
ViewFactors: Computing view factors as defined in file: 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: 0
LoadInputFile: Number of Initial Conditions: 1
LoadInputFile: Number of Materials: 2
LoadInputFile: Number of Equations: 1
LoadInputFile: Number of Solvers: 1
LoadInputFile: Number of Bodies: 2
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
ElmerAsciiMesh: Base mesh name: ./.
MapCoordinates: Scaling coordinates: 1.000E+00 1.000E+00 1.000E+00
LoadMesh: Elapsed REAL time: 0.1690 (s)
ViewFactors: Computing view factors for radiation body1
ViewFactors: Number of surfaces participating in radiation: 4480
ViewFactors: Computing viewfactors...
surfs: 4480, min(353)=0.00, max(29)=1.00, avg=0.99
ViewFactors: View factors computed in time (s): 10.40 10.40
ViewFactors:
ViewFactors: Viewfactors before manipulation:
ViewFactors: Minimum row sum: 0.000E+00
ViewFactors: Maximum row sum: 9.999E-01
ViewFactors: Symmetrizing Factors...
ViewFactors: View factors manipulated in time (s): 0.12
ViewFactors:
ViewFactors: Viewfactors after manipulation:
ViewFactors: Minimum row sum: 0.000E+00
ViewFactors: Maximum row sum: 9.999E-01
ViewFactors: Saving view factors in ascii mode
ViewFactors: *** ALL DONE ***
RadiationFactors: Loading view factors from ascii file: ././ViewFactors.dat
RadiationFactors: Computing factors...
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 932044 353
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 940739 365
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 940739 365
ERROR:: MakeMatrixIndex: Trying to access non-existent column: 966848 373
RadiationFactors: View factors filling (%) 9.2932E+01
RadiationFactors: Solution: 0 % done
RadiationFactors: Solution: 7 % done
RadiationFactors: Solution: 16 % done
RadiationFactors: Solution: 24 % done
RadiationFactors: Solution: 33 % done
RadiationFactors: Solution: 42 % done
RadiationFactors: Solution: 50 % done
RadiationFactors: Solution: 60 % done
RadiationFactors: Solution: 69 % done
RadiationFactors: Solution: 79 % done
RadiationFactors: Solution: 89 % done
RadiationFactors: Solution: 98 % done
RadiationFactors: Minimum Gebhart factors sum 0.000000E+00
RadiationFactors: Maximum Gebhart factors sum 8.910557E-04
RadiationFactors: Maximum share of omitted factors NaN
RadiationFactors: Gebhart factors filling (%) 1.992985E-05
RadiationFactors: Gebhart factors count 4
When running the radiation tutorial file as well as the radiation equilibrium between two crystals everything works just fine, as long as i dont define a surface radiating onto nothing.
I also ask myself if Elmer even can simulate heat loss through radiation without another body receiving the radiation, because i couldnt really find anything in the manuals on it.
The project folder is attached in a .zip file, but In have ommited the view factors file, because it was giant.
Thank you for your help in advance.