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.
https://www.dynaexamples.com/icfd/intermediate-examples/vena-cava
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
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.
*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
$$$ 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
