Solve with linear restriction

Numerical methods and mathematical models of Elmer
Post Reply
Martin_
Posts: 11
Joined: 08 Feb 2010, 16:27

Solve with linear restriction

Post by Martin_ »

Dear Elmer users,

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
Does anyone have an idea how to fix this?

Thank you for any hints.

Best regards,

Martin Vymazal
Post Reply