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
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
Kindly suggest,
Thank you in advance,
Shubham