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
$$$ 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
Download
Description
