Error Creating Data Files for Coupled Momentum Method FSI

Provides a system for patient-specific cardiovascular modeling and simulation.
POST REPLY
User avatar
Ryan Pewowaruk
Posts: 14
Joined: Mon Nov 13, 2017 1:04 pm

Error Creating Data Files for Coupled Momentum Method FSI

Post by Ryan Pewowaruk » Tue Jul 17, 2018 9:23 am

I have been running into problems trying to perform FSI simulations with the coupled momentum method, even with the cylinder example from the tutorial.

When I create the Data files, two errors appear

Code: Select all

warning: illegal maxChild (410994/)

Code: Select all

************************************************************
 ****************** f a t a l    e r r o r ******************
 ************************************************************

 routine cg
Simvascular still says that data files have been created but the displacement.vtp file has only zero values when opened in paraview.

Thanks
Top
User avatar
Justin Tran
Posts: 109
Joined: Mon Sep 30, 2013 4:10 pm

Re: Error Creating Data Files for Coupled Momentum Method FSI

Post by Justin Tran » Thu Aug 02, 2018 8:26 am

Hi Ryan,

Thank you for your question, and apologies for the late reply. The two errors and warnings (about the illegal max child and the fatal error) are actually expected and usually do not cause any problems for users. We are using a CG linear solver from like the 1970s to solve for the initial displacement field. The function thinks it is expecting some kind of data, but it does not need it.

Though the fact that your displacement field is zero is curious. Can you send me the entire output of your presolver so that I may investigate?
Top
User avatar
Ryan Pewowaruk
Posts: 14
Joined: Mon Nov 13, 2017 1:04 pm

Re: Error Creating Data Files for Coupled Momentum Method FSI

Post by Ryan Pewowaruk » Wed Aug 15, 2018 7:52 am

Here is the complete output from the pre-solver. Thanks!

Code: Select all

LINE 0000: mesh_and_adjncy_vtu mesh-complete/mesh-complete.mesh.vtu

LINE 0001: set_surface_id_vtp mesh-complete/mesh-complete.exterior.vtp 
LINE 0002: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA.vtp 2
LINE 0003: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_1.vtp
LINE 0004: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_1a.vt
LINE 0005: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_1b.vt
LINE 0006: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_1c.vt
LINE 0007: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_2.vtp
LINE 0008: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_2a.vt
LINE 0009: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_3.vtp
LINE 0010: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_3a.vt
LINE 0011: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_3b.vt
LINE 0012: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_3c.vt
LINE 0013: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_4.vtp
LINE 0014: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_4a.vt
LINE 0015: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_4b.vt
LINE 0016: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_LPA_4c.vt
LINE 0017: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_MPA.vtp 1
LINE 0018: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_MPA_2.vtp
LINE 0019: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_1.vtp
LINE 0020: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_1a.vt
LINE 0021: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_1b.vt
LINE 0022: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_2.vtp
LINE 0023: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_2a.vt
LINE 0024: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_3.vtp
LINE 0025: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_3a.vt
LINE 0026: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_3b.vt
LINE 0027: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_3c.vt
LINE 0028: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_3d.vt
LINE 0029: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_3e.vt
LINE 0030: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_4.vtp
LINE 0031: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_5.vtp
LINE 0032: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_5a.vt
LINE 0033: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_5b.vt
LINE 0034: set_surface_id_vtp mesh-complete/mesh-surfaces/cap_RPA_5c.vt
LINE 0035: fluid_density 1.06

LINE 0036: fluid_viscosity 0.04

LINE 0037: initial_pressure 0

LINE 0038: initial_velocity 0.0001 0.0001 0.0001

LINE 0039: prescribed_velocities_vtp mesh-complete/mesh-surfaces/cap_MP
LINE 0040: bct_analytical_shape womersley

LINE 0041: bct_period 0.7

LINE 0042: bct_point_number 71

LINE 0043: bct_fourier_mode_number 15

LINE 0044: bct_flip

LINE 0045: bct_create mesh-complete/mesh-surfaces/cap_MPA.vtp cap_MPA.f
LINE 0046: bct_write_dat bct.dat

LINE 0047: bct_write_vtp bct.vtp

LINE 0048: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA.vtp 0

LINE 0049: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_1.vtp 0

LINE 0050: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_1a.vtp 0

LINE 0051: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_1b.vtp 0

LINE 0052: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_1c.vtp 0

LINE 0053: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_2.vtp 0

LINE 0054: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_2a.vtp 0

LINE 0055: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_3.vtp 0

LINE 0056: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_3a.vtp 0

LINE 0057: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_3b.vtp 0

LINE 0058: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_3c.vtp 0

LINE 0059: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_4.vtp 0

LINE 0060: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_4a.vtp 0

LINE 0061: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_4b.vtp 0

LINE 0062: pressure_vtp mesh-complete/mesh-surfaces/cap_LPA_4c.vtp 0

LINE 0063: pressure_vtp mesh-complete/mesh-surfaces/cap_MPA_2.vtp 0

LINE 0064: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_1.vtp 0

LINE 0065: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_1a.vtp 0

LINE 0066: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_1b.vtp 0

LINE 0067: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_2.vtp 0

LINE 0068: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_2a.vtp 0

LINE 0069: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_3.vtp 0

LINE 0070: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_3a.vtp 0

LINE 0071: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_3b.vtp 0

LINE 0072: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_3c.vtp 0

LINE 0073: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_3d.vtp 0

LINE 0074: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_3e.vtp 0

LINE 0075: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_4.vtp 0

LINE 0076: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_5.vtp 0

LINE 0077: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_5a.vtp 0

LINE 0078: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_5b.vtp 0

LINE 0079: pressure_vtp mesh-complete/mesh-surfaces/cap_RPA_5c.vtp 0

LINE 0080: deformable_wall_vtp mesh-complete/walls_combined.vtp

warning: illegal maxChild (251316)
warning: illegal maxChild (251316)
LINE 0081: deformable_thickness .5

LINE 0082: deformable_E 1e6

LINE 0083: deformable_nu 0.5

LINE 0084: deformable_kcons 0.833333

LINE 0085: deformable_pressure 0

LINE 0086: deformable_solve_displacements

 N=      126471
 ndim=     1386582
 maxnz=     1386582
 before call nw=      652353
 before call inw=      379413

     initial iterative parameters
     preprocessor and preconditioner parameters
          iparm(12) =              4    (nstore)
          iparm(13) =              0    (iscale)
          iparm(14) =              0    (iperm )
          iparm(15) =              1    (ifact )
          iparm(16) =              0    (lvfill)
          iparm(17) =              0    (ltrunc)
          iparm(18) =              2    (ipropa)
          iparm(19) =             -1    (kblsz )
          iparm(20) =             -1    (nbl2d )
          iparm(21) =              1    (ifctv )
          iparm(22) =              1    (iqlr  )
          iparm(23) =              2    (isymm )
          iparm(24) =              0    (ielim )
          iparm(25) =              1    (ndeg  )
          rparm(13) = 0.00000000D+00    (timfac)
          rparm(14) = 0.00000000D+00    (timtot)
          rparm(15) = 0.35500000D-11    (tol   )
          rparm(16) = 0.00000000D+00    (ainf  )


 ************************************************************
 ****************** f a t a l    e r r o r ******************
 ************************************************************

 routine cg
Top
User avatar
Justin Tran
Posts: 109
Joined: Mon Sep 30, 2013 4:10 pm

Re: Error Creating Data Files for Coupled Momentum Method FSI

Post by Justin Tran » Fri Aug 17, 2018 8:59 am

Hi Ryan,

Thank you for posting your complete output! One thing that stood out to me was this command:

LINE 0085: deformable_pressure 0

This command sets the scalar pressure that will be used for the displacement solve. CMM requires that the deformable wall simulation starts out with an initial displacement configuration that is in equilibrium with the displacement field. The pre-solver solves a static problem for a displacement field that is in equilibrium with a spatially uniform pressure field. This property sets the value for that pressure field. We recommend setting this parameter to the spatially averaged pressure in the last step of your rigid wall simulation. Setting this parameter to zero is likely the cause of your error.

Hope that helps! Let us know how it goes.
Top
POST REPLY