I try to utilize mesh refinement options of Elmer in a MagnetoDynamics2D model with mortar BC and rotating mesh part.
The attempt to increase mesh level by Mesh Levels = 2 in the Simulation section results in segmentation fault
Code: Select all
ELMER SOLVER (v 7.0) STARTED AT: 2014/07/12 20:24:33
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: Library version: 7.0 (Rev: 6814)
MAIN: HYPRE library linked in.
MAIN: Trilinos library linked in.
MAIN: MUMPS library linked in.
MAIN: =============================================================
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: case_PMSM.sif
Model Input: Unlisted keyword: [interpolation numeric epsilon] in section: [simulation]
Model Input: Unlisted keyword: [r outer] in section: [body 16]
Model Input: Unlisted keyword: [r inner] in section: [body 16]
Loading user function library: [RigidMeshMapper]...[RigidMeshMapper_Init0]
Loading user function library: [MagnetoDynamics2D]...[MagnetoDynamics2D_Init0]
Loading user function library: [MagnetoDynamics2D]...[BSolver_Init0]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init0]
Loading user function library: [SaveData]...[SaveScalars_Init0]
Loading user function library: [SaveData]...[SaveLine_Init0]
LoadMesh: Starting
ElmerAsciiMesh: Performing step: 1
LoadMesh: Base mesh name: PMSM/.
LoadMesh: Reading header info from file: PMSM/./mesh.header
LoadMesh: Number of nodes in mesh: 11097
LoadMesh: Number of bulk elements in mesh: 22120
LoadMesh: Number of boundary elements in mesh: 920
LoadMesh: Initial number of max element nodes: 3
ElmerAsciiMesh: Performing step: 2
LoadMesh: Reading nodes from file: PMSM/./mesh.nodes
LoadMesh: Dimension of mesh is: 2
LoadMesh: Scaling coordinates: 1.000E-03 1.000E-03
ElmerAsciiMesh: Performing step: 3
LoadMesh: Reading bulk elements from file: PMSM/./mesh.elements
ElmerAsciiMesh: Performing step: 4
LoadMesh: Reading boundary elements from file: PMSM/./mesh.boundary
LoadMesh: Performing node mapping
LoadMesh: Remapping bodies
LoadMesh: Minimum initial body index: 1
LoadMesh: Maximum initial body index: 16
LoadMesh: Remapping boundaries
LoadMesh: Minimum initial boundary index: 1
LoadMesh: Maximum initial boundary index: 2
CreateDiscontMesh: Number of discontinuous boundary elements: 848
CreateDiscontMesh: Number of discontinuous parent elements: 1696
CreateDiscontMesh: Number of candicate nodes: 848
CreateDiscontMesh: All candidate nodes used
CreateDiscontMesh: Number of discontinuous nodes: 848
CreateDiscontMesh: Number of primary staying nodes: 1696
CreateDiscontMesh: Number of primary moving nodes: 1688
CreateDiscontMesh: Family: 300 Decided: 1720, Undecided: 0
CreateDiscontMesh: Undecied boundary elements by family: 300 Decided: 1720, Undecided: 0
CreateDiscontMesh: Number of related elements moving: 840
CreateDiscontMesh: Number of related elements staying: 880
CreateDiscontMesh: All elements marked either moving or staying
EnlargeCoordinates: Increasing number of nodes from 11097 to 11945
EnlargeBoundaryElements: Increasing number of elements from 23040 to 23888
CreateDiscontMesh: Number of original bulk elements: 22120
CreateDiscontMesh: Number of original boundary elements: 1768
CreateDiscontMesh: Number of additional boundary elements: 848
EnlargeParallelInfo: Increasing global numbering size from 11097 to 11945
ElmerAsciiMesh: Performing step: 5
ElmerAsciiMesh: Performing step: 6
LoadMesh: Loading mesh done
LoadMesh: Elapsed time (CPU,REAL): 0.0600 0.0784 (s)
SplitMeshEqual: ******** Old mesh ********
SplitMeshEqual: Nodes : 11945
SplitMeshEqual: Bulk elements : 22120
SplitMeshEqual: Boundary elements : 1768
SplitMeshEqual: Edges : 34064
SplitMeshEqual: Faces : 0
SplitMeshEqual: ******** New mesh ********
SplitMeshEqual: Nodes : 46009
SplitMeshEqual: Bulk elements : 88480
SplitMeshEqual: Boundary elements : 3536
SplitMeshEqual: Elapsed time (CPU,REAL): 0.0200 0.0361 (s)
MAIN: -------------------------------------
Loading user function library: [RigidMeshMapper]...[RigidMeshMapper_Init]
Loading user function library: [RigidMeshMapper]...[RigidMeshMapper]
Loading user function library: [MagnetoDynamics2D]...[MagnetoDynamics2D_Init]
Loading user function library: [MagnetoDynamics2D]...[MagnetoDynamics2D]
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: mgdyn2d...done.
OptimizeBandwidth: Half bandwidth without optimization: 44429
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 1171
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [MagnetoDynamics2D]...[BSolver_Init]
Loading user function library: [MagnetoDynamics2D]...[BSolver]
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: mgdyn2dpost...done.
OptimizeBandwidth: Half bandwidth without optimization: 41935
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 1153
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init]
Loading user function library: [ResultOutputSolve]...[ResultOutputSolver]
Loading user function library: [SaveData]...[SaveScalars_Init]
Loading user function library: [SaveData]...[SaveScalars]
Loading user function library: [SaveData]...[SaveLine_Init]
Loading user function library: [SaveData]...[SaveLine]
MAIN:
MAIN: -------------------------------------
MAIN: Time: 1/31 5.55555555555999970E-005
MAIN: -------------------------------------
MAIN:
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: meshdeform
RigidMeshMapper: ---------------------------------------
RigidMeshMapper: Performing analytic mesh mapping
RigidMeshMapper: ---------------------------------------
RigidMeshMapper: Allocating new nodes of size: 46009
RigidMeshMapper: Number of nodes mapped: 20809
RigidMeshMapper: Coordinate mapping time: 0.68400000000000005
RigidMeshMapper: All done
WeightedProjectorDiscont: Creating Galerkin projector for discontinuous boundary
Segmentation fault (core dumped)
First order input mesh is used.
Where could be the problem?
BR, Pavel