## Problem with ArterySolve post file

Numerical methods and mathematical models of Elmer
shubh_shinde
Posts: 19
Joined: 13 Aug 2017, 23:58
Antispam: Yes

### Problem with ArterySolve post file

Dear all,
I am solving a fsi problem with ArteryOutlet boundary condition. I am writing the results in case.vtu file. The solver stops after one timestep and tells problem with the vtu file as below

Code: Select all

``````ComputeChange: SS (ITER=52) (NRM,RELC): ( 0.37410228E-08 0.53956913E-04 ) :: mesh update
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: artery outlet solver
OutletCompute: Starting
LumpedForceCompute: -------------------------------------
LumpedForceCompute: Computing Lumped Boundary:
LumpedForceCompute: -------------------------------------
OutletCompute: Lumping Boundary:  1
OutletCompute: Force: 0.2060E-06
OutletCompute: Area: 0.6283E-02
OutletCompute: Flux: 0.3912E-08
OutletCompute: Nonlinear iteration: 1
SolveSystem: Solving linear system
DirectSolver: Using direct method: banded
ComputeNorm: Computing norm of solution
ComputeChange: NS (ITER=1) (NRM,RELC): (  16.846774     0.16522544E-10 ) :: artery outlet solver
OutletCompute: Relative Change: 0.4130E-10
ComputeNorm: Computing norm of solution
ComputeChange: SS (ITER=52) (NRM,RELC): (  16.846774     0.16522544E-10 ) :: artery outlet solver
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: internalvtuoutputsolver
ResultOutputSolver: -------------------------------------
ResultOutputSolve: Saving with prefix: case
ResultOutputSolver: Working on mesh: .
ResultOutputSolver: Dimension of mesh is: 2
ResultOutputSolver: Creating list for saving - if not present
CreateListForSaving: Field Variables for Saving
CreateListForSaving: Scalar Field 1: pressure
CreateListForSaving: Vector Field 1: velocity
CreateListForSaving: Vector Field 2: displacement
CreateListForSaving: Vector Field 3: mesh velocity
CreateListForSaving: Vector Field 2 Complement: mesh update
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Saving results in VTK XML format with prefix: case
VtuOutputSolver: Saving number of partitions: 1
VtuOutputSolver: Number of active elements 1362 out of 1362
VtuOutputSolver: Number of geometry nodes 632 out of 632
VtuOutputSolver: Total number of geometry nodes to save:     632
VtuOutputSolver: Total number of dof nodes to save:     632
VtuOutputSolver: Total number of elements to save:    1362
VtuOutputSolver: Full filename base is: ././case
VtuOutputSolver: Setting offset for boundary entities: 100
VtuOutputSolver: Writing the vtu file: ././case0001.vtu
AscBinWriteInit: Initializing buffered ascii/binary writing
AscBinWriteInit: Writing in binary
AscBinWriteInit: Writing in double precision
AscBinWriteInit: Writing to unit number: 58
AscBinWriteInit: Size of buffer is: 1362
VtuOutputSolver: Writing nodal fields
VtuOutputSolver: Saving variable: pressure
VtuOutputSolver: Saving variable: velocity
VtuOutputSolver: Saving variable: displacement
VtuOutputSolver: Saving variable: mesh velocity
VtuOutputSolver: Number of nodal fields written: 4
VtuOutputSolver: Writing elemental fields
VtuOutputSolver: Number of elemental fields written: 0
VtuOutputSolver: Writing entity IDs for bodies and boundaries
VtuOutputSolver: Writing coordinates for each used node
VtuOutputSolver: Writing the elemental connectivity data
VtuOutputSolver: Writing nodal fields
VtuOutputSolver: Saving variable: pressure
VtuOutputSolver: Saving variable: velocity
VtuOutputSolver: Saving variable: displacement
VtuOutputSolver: Saving variable: mesh velocity
VtuOutputSolver: Writing elemental fields
AscBinWriteInit: Terminating buffered ascii/binary writing
VtuOutputSolver: All done for now
ResultOutputSolver: -------------------------------------
ResultOutputSolver: Working on mesh: 1d_mesh
ResultOutputSolver: Dimension of mesh is: 1
ResultOutputSolver: Recreating list for saving
CreateListForSaving: Field Variables for Saving
CreateListForSaving: Scalar Field 1: wnodal
ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format
VtuOutputSolver: Saving results in VTK XML format with prefix: case
VtuOutputSolver: Saving number of partitions: 1
VtuOutputSolver: Number of active elements 20 out of 20
VtuOutputSolver: Number of geometry nodes 21 out of 21
VtuOutputSolver: Total number of geometry nodes to save:      21
VtuOutputSolver: Total number of dof nodes to save:      21
VtuOutputSolver: Total number of elements to save:      20
VtuOutputSolver: Full filename base is: ./1d_mesh/case
VtuOutputSolver: Setting offset for boundary entities: 100
VtuOutputSolver: Writing the vtu file: ./1d_mesh/case0001.vtu
AscBinWriteInit: Initializing buffered ascii/binary writing
AscBinWriteInit: Writing in binary
AscBinWriteInit: Writing in double precision
AscBinWriteInit: Writing to unit number: 58
AscBinWriteInit: Size of buffer is: 21
At line 986 of file /build/elmerfem-csc-ULYpJ0/elmerfem-csc-8.2/fem/src/modules/ResultOutputSolve/VtuOutputSolver.F90 (unit = 58, file = '¸F7')
Fortran runtime error: File './1d_mesh/case0001.vtu' does not exist``````
I have made a 1D mesh in the solver itself as,

Code: Select all

``````Solver 4
Equation = Artery Outlet Solver
Time Derivative Order = 1
Stabilize = Logical False
Variable = Wnodal
Variable DOFs = 1
Procedure = "./ArteryOutlet" "OutletCompute"
Linear System Solver = Direct
Linear System Direct Method = banded
Linear System Iterative Method = CG
Linear System Convergence Tolerance = 1.0e-6
Linear System Max Iterations = 1000
Linear System Preconditioning = ILU0
Nonlinear System Max Iterations = Integer 1000
Nonlinear System Convergence Tolerance = Real 1.0e-8
Optimize Bandwidth = False

! These define the 1D mesh that is created on-the-fly for this solver
1D Mesh Create = Logical True
1D Element Order = Integer 1
1D Number Of Elements = Integer 20
1D Mesh Length = Real 0.005
1D Active Direction = Integer 1
1D Body Id = Integer 3
End``````
I have attached sif file for reference.
Kindly suggest,
Shubham
Attachments
test_AS.sif

raback
Posts: 3230
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

### Re: Problem with ArterySolve post file

Hi

The 1d mesh is created on-the-fly i.e. it is not on the disk. Hence when the code tries to save the vtu files to the mesh directory it fails. You have some options

1) Create the directory before running: mkdir 1d_mesh
2) Prescribe a different output directory in for the output solver: Output Directory = results

So the problem is that the code assumes the mesh directory to exist even though it doesn't.

-Peter

shubh_shinde
Posts: 19
Joined: 13 Aug 2017, 23:58
Antispam: Yes

### Re: Problem with ArterySolve post file

Hi Peter,
This worked great. Thank you.
-Shubham

shubh_shinde
Posts: 19
Joined: 13 Aug 2017, 23:58
Antispam: Yes

### Re: Problem with ArterySolve post file

Hi Peter,
I tried to solve the 2D blood flow in artery as shown in the Models manual, I get the convergence however the results are not proper. The results show only velocity at the inlet and not propagating through the geometry.
I have attached the sif file. Could you please suggest what am I missing?
The geometry I am working with is exactly same as that of shown in Models manual for artery flow.
Thanks,
Shubham
Attachments
test_AS.sif