Advanced : Flow through a funnel
In this coupled DEM particles-CFD free surface example, water is being poured in a funnel filled with Styrofoam spheres causing stirring and flow deviation.
https://www.dynaexamples.com/icfd/advanced-examples/funnel
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Advanced : Flow through a funnel
In this coupled DEM particles-CFD free surface example, water is being poured in a funnel filled with Styrofoam spheres causing stirring and flow deviation.
Fluid velocity fringes
*KEYWORD *TITLE *CONTROL_DISCRETE_ELEMENT *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *DEFINE_DE_TO_SURFACE_COUPLING *ELEMENT_DISCRETE_SPHERE_VOLUME *ELEMENT_SHELL *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_BOUNDARY_PRESCRIBED_VEL *ICFD_CONTROL_DEM_COUPLING *ICFD_CONTROL_MESH *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_MAT *ICFD_MODEL_POROUS *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE *LOAD_BODY *MAT_RIGID *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *PART *SECTION_SHELL *END
$-----------------------------------------------------------------------------
$
$ Example provided by Facundo and Iñaki (LSTC)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file as is.
$X Requires LS-DYNA MPP Dev 117500 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD Funnel flow
*INCLUDE
mesh_fluid.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_DEM_COUPLING
0
*ICFD_CONTROL_MESH
1.2
*ICFD_CONTROL_TIME
$# ttm dt cfl lcidsf dtmin dtmax
20.0 0.0 1.0 0 0.0 0.0
*ICFD_CONTROL_OUTPUT
$# msgl outl dtout lsppout
3
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_MAT
$# mid flg ro vis
1 1 1000.0 0.001567
*ICFD_MAT
$# mid flg ro vis
2 0 0.0 0.0
*ICFD_PART_TITLE
Wall1
$# pid secid mid
1 1 2
*ICFD_PART_TITLE
Outflow2
$# pid secid mid
2 1 2
*ICFD_PART_TITLE
Wall3
$# pid secid mid
3 1 2
*ICFD_PART_TITLE
Wall4
$# pid secid mid
4 1 1
*ICFD_PART_TITLE
Inflow5
$# pid secid mid
5 1 2
*ICFD_PART_VOL
$# pid secid mid
6 1 2
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
1 2 3 4 5 0 0 0
*ICFD_SECTION
$# sid
1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD BOUNDARY/INITIAL/LOAD CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_NONSLIP
$# pid
3
*ICFD_BOUNDARY_FREESLIP
1
*ICFD_BOUNDARY_FREESLIP
5
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
2 1 1.01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
4 4 1 2 1.0 01.00000E28 0.0
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
200.0 0.0
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 1.0
200.0 1.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$# volid
7
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 2 3 4 5 0 0 0
In this coupled DEM particles-CFD free surface example, water is being poured in a funnel filled with Styrofoam type spheres causing stirring and flow deviation.
