Basics : DEM coupling — Welcome to LS-DYNA Examples

Home » Knowledge Base » ICFD / Incompressible CFD » Basics Examples » Basics : DEM coupling

It is possible to couple the ICFD fluid solver with DEM particles. This feature can useful in cases involving mud, dust or snow deposition. The coupling is explicit and based on the particle drag coefficient and force.

Animated Result

Fluid velocity fringes

plate

 

 

Keywords

*TITLE
*CONTROL_DISCRETE_ELEMENT
*CONTROL_TERMINATION
*CONTROL_TIMESTEP
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_TITLE
*ICFD_BOUNDARY_FREESLIP
*ICFD_BOUNDARY_PRESCRIBED_VEL
*ICFD_BOUNDARY_PRESCRIBED_PRE
*ICFD_BOUNDARY_NONSLIP
*ICFD_CONTROL_DEM_COUPLING
*ICFD_CONTROL_TIME
*ICFD_INITIAL
*ICFD_MAT
*ICFD_PART
*ICFD_PART_VOL
*ICFD_SECTION
*INCLUDE
*KEYWORD
*MAT_RIGID_DISCRETE
*MESH_SIZE_SHAPE
*MESH_SURFACE_ELEMENT
*MESH_SURFACE_NODE
*MESH_VOLUME
*PARAMETER
*PART
*SECTION_SOLID
*END

Reduced Input

$-----------------------------------------------------------------------------
$
$ Example provided by 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 R9.0.0 (or higher) with double precision 
$X
$X------------------------------------------------------------------------------
$# UNITS: Dimensionless. 
$X------------------------------------------------------------------------------
$X
*keyword
*title
ICFD DEM Coupling
*include
mesh.k
*include
dem.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             PARAMETERS                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R    T_end     100.0
R  dt_plot      1.00
$
$--- Fluid
$
R  v_inlet       1.0
Rrho_fluid       2.0
R mu_fluid      0.01 
R dt_fluid     0.050
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                           ICFD CONTROL CARDS                                 $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_DEM_COUPLING
$#   ctype
         0
*ICFD_CONTROL_TIME
$#     ttm        dt
    &T_end &dt_fluid
$---+----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
$#     pid     secid       mid 
         1         1         1
*ICFD_PART
$#     pid     secid       mid 
         2         1         1
*ICFD_PART
$#     pid     secid       mid 
         3         1         1
*ICFD_PART_VOL
$#     pid     secid       mid 
        10         1         1
$#   spid1     spid2     spid3     
         1         2         3         
$---+----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         1         1         1           
*ICFD_BOUNDARY_PRESCRIBED_PRE
$#     pid      lcid        sf     death     birth     
         2         2     
*ICFD_BOUNDARY_FREESLIP
$#     pid   
         3
*ICFD_INITIAL
$#     pid        vx        
         0  &v_inlet 
*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           
         1         2         3
*MESH_SIZE_SHAPE
$#   sname
       box
$#   msize     pminx     pminy     pminz     pmaxx     pmaxy     pmaxz
     0.050        -1        -1        -1         1         1         1        
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             DATABASE (OUTPUT)                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
$#      dt
  &dt_plot
*END
 

Figures

figure_001

Download

Description

It is possible to couple the ICFD fluid solver with DEM particles. This feature can useful in cases involving mud, dust or snow deposition. The coupling is explicit and based on the particle drag coefficient and force.