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