Intermediate : Dam break on shallow wet bed — Welcome to LS-DYNA Examples

Home » Knowledge Base » ICFD / Incompressible CFD » Intermediate Examples » Intermediate : Dam break on shallow wet bed

At t=0, the gate is released with a prescribed velocity allowing the collapse of a water column and propagation of a dam break wave. The depth of the wet bed (shallow/deep) influences on the shape of the propagating flow. Warning : the mesh is fine, 8 CPUs or more are recommended for this case.

Animated Result

Fluid pressure fringes

dam

 

 

Keywords

*KEYWORD
*TITLE
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_TITLE
*DEFINE_FUNCTION
*ICFD_BOUNDARY_FREESLIP
*ICFD_BOUNDARY_NONSLIP
*ICFD_BOUNDARY_PRESCRIBED_PRE
*ICFD_CONTROL_ADAPT_SIZE
*ICFD_CONTROL_IMPOSED_MOVE
*ICFD_CONTROL_OUTPUT
*ICFD_CONTROL_TIME
*ICFD_DATABASE_NODOUT
*ICFD_MAT
*ICFD_PART
*ICFD_PART_VOL
*ICFD_SECTION
*INCLUDE
*LOAD_BODY_Y
*MESH_INTERF
*MESH_SURFACE_ELEMENT
*MESH_SURFACE_NODE
*MESH_VOLUME
*PARAMETER
*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: (kg/m/s) 
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD Dam break on shallow wet bed
*INCLUDE
mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             PARAMETERS                                       $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R    T_end       0.6
R  dt_plot      0.01
$
$--- Fluid
$
Rrho_fluid      1000
R mu_fluid     0.001 
R dt_fluid     0.000
R   dt_min      1e-8
R   dt_max      5e-4
R  grav         9.81
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                           ICFD CONTROL CARDS                                 $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_TIME
$#     ttm        dt
    &T_end &dt_fluid
*ICFD_CONTROL_OUTPUT
$#    msgl
         4
*ICFD_CONTROL_IMPOSED_MOVE
$#     pid      lcvx      lcvy     
         6         2         4
*ICFD_CONTROL_IMPOSED_MOVE
$#     pid      lcvx      lcvy     
         7         2         4
*ICFD_CONTROL_IMPOSED_MOVE
$#     pid      lcvx      lcvy     
         8         2         4
*ICFD_CONTROL_IMPOSED_MOVE
$#     pid      lcvx      lcvy     
         9         2         4
*ICFD_CONTROL_IMPOSED_MOVE
$#     pid      lcvx      lcvy     
        10         2         4
*ICFD_CONTROL_ADAPT_SIZE
$#   asize       nit
         1        -5
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         2                    
$#                a1                  o1
                 0.0                 0.0
             10000.0                 0.0
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         4                    
$#                a1                  o1
                 0.0                 0.0
               0.005                 1.5
                0.11                 1.5
                0.12                 0.0
                   1                 0.0
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         5                    
$#                a1                  o1
                 0.0                   2
               0.010                   2
               0.020                   5
               0.030                  10
               0.100                  10
               0.110                 100
                   1                 100
$---+----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_MAT
$#     mid       flg            
         2         0      
*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
$#     pid     secid       mid 
         4         1         2
*ICFD_PART
$#     pid     secid       mid 
         5         1         2
*ICFD_PART
$#     pid     secid       mid 
         6         1         1
*ICFD_PART
$#     pid     secid       mid 
         7         1         1
*ICFD_PART
$#     pid     secid       mid 
         8         1         1
*ICFD_PART
$#     pid     secid       mid 
         9         1         2
*ICFD_PART
$#     pid     secid       mid 
        10         1         2
*ICFD_PART_VOL
$#     pid     secid       mid 
        12         1         1
$#   spid1     spid2     spid3     spid4     spid5     spid6     spid7
         1         2         3         6         7         8       
*ICFD_PART_VOL
$#     pid     secid       mid 
        13         1         2
$#   spid1     spid2     spid3     spid4     spid5     spid6     spid7
         3         4         5         9        10         
$---+----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   
         1
*ICFD_BOUNDARY_NONSLIP
$#     pid   
         4
*ICFD_BOUNDARY_FREESLIP
$#     pid   
         5
*ICFD_BOUNDARY_FREESLIP
$#     pid   
         6
*ICFD_BOUNDARY_FREESLIP
$#     pid   
         7
*ICFD_BOUNDARY_FREESLIP
$#     pid   
         8
*ICFD_BOUNDARY_FREESLIP
$#     pid   
         9
*ICFD_BOUNDARY_FREESLIP
$#     pid   
        10
*ICFD_BOUNDARY_PRESCRIBED_PRE
$#     pid      lcid
         2         3
*LOAD_BODY_Y
$#    lcid        sf       
         1        1          
*DEFINE_CURVE_TITLE
Gravity force
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp
         1                    &grav
$#                a1                  o1
                 0.0                 1.0
             10000.0                 1.0
*DEFINE_FUNCTION
$#    fcic
         3
$# function
f(y)=-rho_fluid*grav*y+grav*18
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                            ICFD MESH KEYWORDS                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$#   volid     
        30
$#    pid1      pid2      pid3      pid4      pid5      pid6      pid7      pid8        
         1         2         4         5         6         7         8         9
$#   pid10
        10
*MESH_INTERF
$#   volid     
        30
$#    pid1                
         3  
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$                                                                              $
$                             DATABASE (OUTPUT)                                $
$                                                                              $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
&dt_plot
*END
 

Figures

figure_001

Download

Description

At t=0, the gate is released with a prescribed velocity allowing the collapse of a water column and propagation of a dam break wave. The depth of the wet bed (shallow/deep) influences on the shape of the propagating flow. Warning : the mesh is fine, 8 CPUs or more are recommended for this case.

References :

[1] I. M. Jánosi, D. Jan, K. G. Szabó, and T. Tél, “Turbulent drag
reduction in dam-break flows,” Experiments in Fluids, vol. 37, no. 2,
pp. 219–229, 2004.

[2] M. Gomez-Gesteira, A. J. C. Crespo, B. D. Rogers, R. A. Dalrymple,
J. M. Dominguez, and A. Barreiro, “Sphysics – development of a
free-surface fluid solver – part 2: Efficiency and test cases,” Comput.
Geosci., vol. 48, pp. 300–307, Nov. 2012