Hello Peter,

I changed the definitions according to your suggestion, and the simulation starts - hoorray!

But right at the beginning the solver log window already fusses around:

*...*

MAIN: Reading Model ...

Model Input: Unlisted keyword: [parameter 1] in section: [solver 2]

Model Input: Unlisted keyword: [save materials] in section: [material 1]

Loading user function library: [SaveData]...[SaveMaterials_Init0]

MAIN: Done

MAIN: -----------------------

WARNING:: AddEquation: Time stepping method defaulted to IMPLICIT EULER

WARNING:: AddEquation: Time stepping method defaulted to IMPLICIT EULER

WARNING:: AddEquation: Time stepping method defaulted to IMPLICIT EULER

WARNING:: AddEquation: Time stepping method defaulted to IMPLICIT EULER

OptimizeBandwidth: ---------------------------------------------------------

OptimizeBandwidth: Computing matrix structure for: navier-stokes...done.

OptimizeBandwidth: Half bandwidth without optimization: 37955

OptimizeBandwidth:

OptimizeBandwidth: Bandwidth Optimization ...done.

OptimizeBandwidth: Half bandwidth after optimization: 1166

OptimizeBandwidth: ---------------------------------------------------------

WARNING:: AddEquation: Time stepping method defaulted to IMPLICIT EULER

WARNING:: ListFind:

WARNING:: ListFind: Requested property: [Equation], not found

WARNING:: ListFind:

Loading user function library: [SaveData]...[SaveMaterials_Init]

...
After that the first timestep seems to be calculated, but the simulation stops after writing the results of the first step with the following message:

*...*

SaveMaterials: Saving selected parameters as dofs

SaveMaterials: Parameter 1: density

SaveMaterials: All done

WriteToPost: Saving results in ElmerPost format to file ./case.ep.0

MAIN:

MAIN: -------------------------------------

MAIN: Time: 2/40 2.0000000000000000

MAIN: Estimated time left: 4 hours.

MAIN: -------------------------------------

MAIN:

SolveEquations: -------------------------------------

SolveEquations: Coupled system iteration: 1

SolveEquations: -------------------------------------

FlowSolve:

FlowSolve:

FlowSolve: -------------------------------------

FlowSolve: NAVIER-STOKES ITERATION 1

FlowSolve: -------------------------------------

FlowSolve:

FlowSolve: Starting Assembly...

[quad:16976] *** Process received signal ***

[quad:16976] Signal: Segmentation fault (11)

[quad:16976] Signal code: Address not mapped (1)

[quad:16976] Failing at address: (nil)

[quad:16976] [ 0] /lib/libpthread.so.0(+0xf8f0) [0x7f8f222c08f0]

[quad:16976] [ 1] /opt/elmer10/lib/libelmersolver-6.1.so(flowsolver_+0x743b) [0x7f8f23f11bc4]

[quad:16976] [ 2] /opt/elmer10/lib/libelmersolver-6.1.so(+0x14ab02) [0x7f8f23f6cb02]

[quad:16976] [ 3] /opt/elmer10/lib/libelmersolver-6.1.so(execsolver_+0x42) [0x7f8f23f6cb4b]

[quad:16976] [ 4] /opt/elmer10/lib/libelmersolver-6.1.so(__mainutils_MOD_singlesolver+0x67b) [0x7f8f23fc4dd8]

[quad:16976] [ 5] /opt/elmer10/lib/libelmersolver-6.1.so(__mainutils_MOD_solveractivate+0x43a) [0x7f8f23fc5234]

[quad:16976] [ 6] /opt/elmer10/lib/libelmersolver-6.1.so(__mainutils_MOD_solveequations+0x1983) [0x7f8f23fc6c0b]

[quad:16976] [ 7] /opt/elmer10/lib/libelmersolver-6.1.so(+0x27fabd) [0x7f8f240a1abd]

[quad:16976] [ 8] /opt/elmer10/lib/libelmersolver-6.1.so(elmersolver_+0x3079) [0x7f8f240a5a1f]

[quad:16976] [ 9] ElmerSolver_mpi(MAIN__+0xda) [0x40118e]

[quad:16976] [10] ElmerSolver_mpi(main+0x2a) [0x40153a]

[quad:16976] [11] /lib/libc.so.6(__libc_start_main+0xfd) [0x7f8f21f4cc4d]

[quad:16976] [12] ElmerSolver_mpi() [0x400ff9]

[quad:16976] *** End of error message ***

...
The density in fact appears in the list of displayable parameters in ElmerPost. In this respect the operation has succeeded indeed.

Just to give a complete picture of my setting, here is the complete sif file that is put to work:

Code: Select all

```
Header
CHECK KEYWORDS Warn
Mesh DB "." "."
Include Path ""
Results Directory ""
End
Simulation
Max Output Level = 4
Coordinate System = Cartesian
Coordinate Mapping(3) = 1 2 3
!Coordinate Scaling(3) = Real 2.5 2.5 2.5
Simulation Type = Transient
Output Intervals(2) = 1 1
Timestep intervals(2) = 20 20
Timestep Sizes(2) = 1 1
Steady State Max Iterations = 20
Solver Input File = case.sif
Post File = case.ep
End
Constants
Gravity(4) = 0 -1 0 9.82
Stefan Boltzmann = 5.67e-08
Permittivity of Vacuum = 8.8542e-12
Boltzmann Constant = 1.3807e-23
Unit Charge = 1.602e-19
End
Body 1
Name = "Body 1"
Equation = 1
Material = 1
Initial condition = 1
End
Solver 1
Equation = Navier-Stokes
!
Variable = String "Density"
Variable DOFs = 1
Variable = Flow Solution[Velocity:3 Pressure:1]
Bubbles = False
!
Stabilize = True
Steady State Convergence Tolerance = 1.0e-5
Nonlinear System Convergence Tolerance = 1.0e-4
Nonlinear System Max Iterations = 10
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1
Linear System Solver = Iterative
Linear System Iterative Method = BiCGStab
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0e-5
Linear System Preconditioning = ILU2
Linear System Residual Output = 1
End
Solver 2
Exec Solver = after timestep
Procedure = File "SaveData" "SaveMaterials"
Parameter 1 = String Density
End
Equation 1
Name = "NS"
NS Convect = False
Active Solvers(2) = 1 2
End
Initial Condition 1
Name = "Anfangsdruck"
Velocity 2 = 0
Velocity 1 = 0
Pressure = 100000
End
Material 1
Name = "Lubricant oil (coarse values)"
Heat Conductivity = 0.128
Reference Temperature = 293
Heat Capacity = 2093
!Density = Variable Pressure
! Real
! -500000 0.01
! 0 0.01
! 500 0.77
! 100000 0.86
! 1000000 0.86
! End
Density = Variable Pressure
Real MATC "max(0.0001 866/(1+0.09*(1/tx)))"
Compressibility Model = User Defined
Save Materials = Logical true
Viscosity = 20.0e-3
!Viscosity = 20.0e-6
End
Boundary Condition 1
Target Boundaries(1) = 2
Name = "Einlass"
!External Pressure = Variable Time; MATC "-100000 - 900000*tx"
External Pressure = Variable Time
Real
-1 -2000000
0 -3000000
10 -13000000
40 -43000000
End
Velocity 1 = 0
Velocity 2 = 0
End
Boundary Condition 2
Target Boundaries(5) = 3 4 5 7 8
Name = "ruhende Waende"
!Noslip wall BC = on
Velocity 1 = 0
Velocity 2 = 0
Velocity 3 = 0
End
Boundary Condition 3
Target Boundaries(1) = 9
Name = "Auslass"
!External Pressure = Variable Time; MATC "-1000 - 1000*tx"
External Pressure = -100000
Velocity 1 = 0
Velocity 2 = 0
End
```

What is wrong with it, please? Is the SaveMaterials solver compatible with parallel execution?

Kind regards,

Peter