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