Intermediate : Mold filling problem

This aim of this filling problem is to study the flow front and the filling times of the various compartments of the mold. Thermal effects (viscosity function of temperature, heat transfer) are not taken into account but can be added.

*KEYWORD
*TITLE
*ICFD_BOUNDARY_PRESCRIBED_VELOCITY
*ICFD_BOUNDARY_NONSLIP
*ICFD_CONTROL_FSI
*ICFD_CONTROL_OUTPUT
*ICFD_CONTROL_TIME
*ICFD_MAT
*ICFD_PART
*ICFD_PART_VOL
*ICFD_SECTION
*INCLUDE
*LOAD_BODY_Y
*MESH_SURFACE_ELEMENT
*MESH_SURFACE_NODE
*MESH_VOLUME
*PARAMETER
*END

$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LST-Ansys)
$
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X    Requires LS-DYNA MPP R11.0.0 (or higher) with double precision 
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s) 
*KEYWORD
*TITLE
Mold filling
*include
mesh_mold_scaled.k
*PARAMETER
R    T_end        4.
R  dt_plot      0.01
$
$--- Fluid
$
Rrho_fluid      7260
R mu_fluid    6.7e-3 
R dt_fluid    0.0001
R   dt_min      1e-8
R   dt_max      5e-4
R  grav         9.81
*DATABASE_BINARY_D3PLOT
$  dt/cycl      lcdt      beam     npltc
&dt_plot
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                           ICFD CONTROL CARDS                                 $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_OUTPUT
$#    msgl
         4
*ICFD_CONTROL_TURBULENCE
2
*ICFD_CONTROL_TIME
$#     ttm        dt
    &T_end &dt_fluid
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                    ICFD BOUNDARY/INITIAL/LOAD CONDITIONS                     $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_PRESCRIBED_VEL
3 , 1, 1, 1
*ICFD_BOUNDARY_PRESCRIBED_VEL
3 , 2, 1, 2
*ICFD_BOUNDARY_PRESCRIBED_VEL
3 , 3, 1, 3
*ICFD_BOUNDARY_NONSLIP
1
*ICFD_BOUNDARY_NONSLIP
2
*LOAD_BODY_z
100,1.
*DEFINE_CURVE
100
0., &grav
1.e5, &grav
*DEFINE_CURVE
 1 , 0, 1. , 1.130 , 0. , 0.
0.,0.601486
1e3,0.601486
*DEFINE_CURVE
 2 , 0, 1. , 1.130 , 0. , 0.
0.,0.056338
1e3,0.056338
*DEFINE_CURVE
 3 , 0, 1. , 1.130 , 0. , 0.
0.,-0.796905
1e3,-0.796905
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                       ICFD PARTS/ SECTION/ MATERIAL                          $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_PART
1,1,1
*ICFD_PART
2,1,4
*ICFD_PART
3,1,4
*ICFD_PART
4,1,4
*ICFD_PART_VOL
10,1,4
2,3,4
*ICFD_PART_VOL
11,1,1
1,4
*ICFD_MAT
 1, 0
*ICFD_MAT
$#     mid       flg        ro       vis          
         4         1&rho_fluid &mu_fluid 
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                            ICFD MESH KEYWORDS                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
20
1,2,3
*MESH_INTERF
20
4
*MESH_BL
1,1
*MESH_BL
2,1
*MESH_BL_SYM
3
*END