Problem with ArterySolve post file

Numerical methods and mathematical models of Elmer

Problem with ArterySolve post file

Postby shubh_shinde » 11 Oct 2017, 23:42

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,
Thank you in advance,
Shubham
Attachments
test_AS.sif
(4.35 KiB) Downloaded 4 times
shubh_shinde
 
Posts: 9
Joined: 13 Aug 2017, 23:58

Re: Problem with ArterySolve post file

Postby raback » 12 Oct 2017, 14:10

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
raback
Site Admin
 
Posts: 3139
Joined: 22 Aug 2009, 11:57
Location: Espoo, Finland

Re: Problem with ArterySolve post file

Postby shubh_shinde » 12 Oct 2017, 21:19

Hi Peter,
This worked great. Thank you.
-Shubham
shubh_shinde
 
Posts: 9
Joined: 13 Aug 2017, 23:58

Re: Problem with ArterySolve post file

Postby shubh_shinde » 13 Oct 2017, 00:49

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
(4.31 KiB) Downloaded 6 times
shubh_shinde
 
Posts: 9
Joined: 13 Aug 2017, 23:58


Return to ElmerSolver

Who is online

Users browsing this forum: No registered users and 5 guests