I'm trying to run testcase 24 ('Volume flow boundary condition') from Elmer tutorials. I downloaded tutorial files, and without modifying anything, I proceed as follows:
1) generate mesh: ElmerGrid 1 2 mflow
2) compile the fortran file to impose prescribed volume flow: elmerf90 AddMassFlow.f90 -o AddMassFlow
then I run the testcase by typing 'ElmerSolver', but the computation crashes very quickly:
Code: Select all
ELMER SOLVER (v 5.5.0) STARTED AT: 2010/02/08 14:31:55
MAIN:
MAIN: ==========================================
MAIN: E L M E R S O L V E R S T A R T I N G
MAIN: Library version: 5.5.0 (Rev: 4322)
MAIN: ==========================================
MAIN:
MAIN:
MAIN: -----------------------
MAIN: Reading Model ...
MAIN: Done
MAIN: -----------------------
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: navier-stokes...done.
OptimizeBandwidth: Half bandwidth without optimization: 16
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 28
OptimizeBandwidth: Bandwidth optimization rejected, using original ordering.
OptimizeBandwidth: ---------------------------------------------------------
Loading user function library: [./AddMassFlow]...[AddMassFlow]
MAIN:
MAIN: -------------------------------------
MAIN: Time: 1/100 2.00000000000000004E-002
MAIN: -------------------------------------
MAIN:
SolveEquations: -------------------------------------
SolveEquations: Coupled system iteration: 1
SolveEquations: -------------------------------------
FlowSolve:
FlowSolve:
FlowSolve: -------------------------------------
FlowSolve: NAVIER-STOKES ITERATION 1
FlowSolve: -------------------------------------
FlowSolve:
FlowSolve: Starting Assembly...
FlowSolve: Assembly:
FlowSolve: Assembly done
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
FlowSolve: Dirichlet conditions done
AddMassFlow:
SolveWithLinearRestriction:
ERROR:: AddMassFlow: RestMatrix not associated
Thank you for any hints.
Best regards,
Martin Vymazal