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.
Fluid pressure fringes
*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
$-----------------------------------------------------------------------------
$
$ 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
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