Intermediate : Vena Cava — Welcome to LS-DYNA Examples

Home » Knowledge Base » ICFD / Incompressible CFD » Intermediate Examples » Intermediate : Vena Cava

This problem is part of a series of benchmarks proposed by the FDA and consists in studying the time averaged flow in a patientaveraged inferior vena cava.

Keywords

*KEYWORD
*TITLE
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_TITLE
*ICFD_BOUNDARY_NONSLIP
*ICFD_BOUNDARY_PRESCRIBED_VEL
*ICFD_BOUNDARY_PRESCRIBED_PRE
*ICFD_CONTROL_OUTPUT
*ICFD_CONTROL_MESH
*ICFD_CONTROL_TIME
*ICFD_DATABASE_FLUX
*ICFD_MAT
*ICFD_PART
*ICFD_PART_VOL
*ICFD_SECTION
*INCLUDE
*MESH_BL
*MESH_BL_SYM
*MESH_SURFACE_ELEMENT
*MESH_SURFACE_NODE
*MESH_VOLUME
*PARAMETER
*END

Reduced Input

$$$ Example provided by CJ Huang (Peggy) - LST
*KEYWORD
*TITLE
ICFD 3D patient-averaged vena cava flow at exercising condition
*INCLUDE
venacava_mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             PARAMETERS                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R    T_end      50.1
R  dt_plot      10.0
$
$--- Fluid
$
R  v_inlet 2.1386e-1
Rrho_fluid    1817.0
R mu_fluid5487.34e-6 
R dt_fluid     0.005
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                           ICFD CONTROL CARDS                                 $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_MPP_IO_NOD3DUMP
*ICFD_CONTROL_TIME
$#     ttm        dt
    &T_end &dt_fluid        1.   
*ICFD_CONTROL_MESH
$#    mgsf
        1.
*ICFD_CONTROL_OUTPUT
         0         7  &dt_plot
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                       ICFD PARTS/ SECTION/ MATERIAL                          $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_SECTION
$#     sid   
         1
*ICFD_MAT
$#     mid       flg        ro       vis          
         1         1&rho_fluid &mu_fluid       
*ICFD_PART_TITLE
iliac inlet left
$#     pid     secid       mid 
         1         1         1
*ICFD_PART_TITLE
iliac inlet right
$#     pid     secid       mid 
         2         1         1
*ICFD_PART_TITLE
IVC outlet
$#     pid     secid       mid 
         3         1         1
*ICFD_PART_TITLE
IVC wall
$#     pid     secid       mid
         4         1         1
*ICFD_PART_VOL
$#     pid     secid       mid 
        10         1         1
$#   spid1     spid2     spid3     
         1         2         3         4
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                    ICFD BOUNDARY/INITIAL CONDITIONS                          $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_PRESCRIBED_VEL
$#     pid       dof       vad      lcid       
         1         4         3         1   
*ICFD_BOUNDARY_PRESCRIBED_VEL
$#     pid       dof       vad      lcid
         2         4         3         1          
*ICFD_BOUNDARY_PRESCRIBED_PRE
$#     pid      lcid        sf     death     birth     
         3         2     
*ICFD_BOUNDARY_NONSLIP
$#     pid   
         4
*DEFINE_CURVE_TITLE
Velocity inlet
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         1                    &v_inlet
$#                a1                  o1
                 0.0                 1.0
             10000.0                 1.0
*DEFINE_CURVE_TITLE
Pressure outlet
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         2
$#                a1                  o1
                 0.0                 0.0
             10000.0                 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                            ICFD MESH KEYWORDS                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$#   volid     
         1
$#    pid1      pid2      pid3      pid4     
         1         2         3         4
*MESH_BL
$#     pid     nelth      blth      blfe      blst
         4         4       1.2    1.4e-4         3
*MESH_BL_SYM
$#    pid1
         1
*MESH_BL_SYM
$#    pid1
         2
*MESH_BL_SYM
$#    pid1
         3
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             DATABASE (OUTPUT)                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_FLUX
$#     pid
         1
*ICFD_DATABASE_FLUX
$#     pid
         2
*ICFD_DATABASE_FLUX
$#     pid
         3
*DATABASE_BINARY_D3PLOT
  &dt_plot
*END

Figures

venacava_rest_section.jpg

Download

– ICFD_example_venacava.zip

Description

Conference Paper