Intermediate : Solitary wave — Welcome to LS-DYNA Examples

Home » Knowledge Base » ICFD / Incompressible CFD » Beta Examples » Intermediate : Solitary wave

This example features an incoming solitary wave impacting a floating box. It directly initializes the free surface level via an initial keyword rather than the classic way of meshing the initial free surface interface.

Animated Result

Fluid velocity fringes

dam

 

Keywords

*KEYWORD
*TITLE
*DATABASE_BINARY_D3PLOT
*DEFINE_CURVE_FUNCTION_TITLE
*DEFINE_CURVE_TITLE
*DEFINE_FUNCTION
*ICFD_BOUNDARY_FREESLIP
*ICFD_BOUNDARY_FSWAVE
*ICFD_BOUNDARY_NONSLIP
*ICFD_BOUNDARY_PRESCRIBED_MOVEMESH
*ICFD_CONTROL_IMPOSED_MOVE
*ICFD_CONTROL_OUTPUT
*ICFD_CONTROL_TIME
*ICFD_DEFINE_POINT
*ICFD_DEFINE_WAVE_DAMPING
*ICFD_INITIAL_LEVELSET
*ICFD_MAT
*ICFD_PART
*ICFD_PART_VOL
*ICFD_SECTION
*INCLUDE
*LOAD_BODY_Y
*MESH_SURFACE_ELEMENT
*MESH_SURFACE_NODE
*MESH_VOLUME
*PARAMETER
*END

Reduced Input

*KEYWORD
$—+—-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 &dt_sf
*ICFD_CONTROL_FSI
$# owc bt
0 &fsi_bt
*ICFD_CONTROL_OUTPUT
$# msglv
3
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
$ $
$ ICFD PARTS/ SECTION/ MATERIAL $
$ $
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
*ICFD_PART
$# pid secid mid
4 1 1
*ICFD_PART
$# pid secid mid
5 1 1
*ICFD_PART
$# pid secid mid
6 1 2
*ICFD_PART
$# pid secid mid
9 1 1
*ICFD_PART
$# pid secid mid
10 1 1
*ICFD_PART_VOL
$# pid secid mid
1 1 1
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
4 5 6 9 10 0 0 0
*ICFD_SECTION
$# sid
1
*ICFD_MAT
$# mid flg ro vis
1 0
*ICFD_MAT
$# mid flg ro vis
2 1&rho_fluid &mu_fluid
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
$ $
$ ICFD BOUNDARY/INITIAL CONDITIONS $
$ $
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
*ICFD_INITIAL_LEVELSET
1 0 1 0 0 0 0
*ICFD_BOUNDARY_NONSLIP
$# pid
10
*ICFD_BOUNDARY_FSI
10
*ICFD_BOUNDARY_NONSLIP
5
*ICFD_BOUNDARY_NONSLIP
6
*ICFD_BOUNDARY_PRESCRIBED_PRE
9,99
*DEFINE_FUNCTION
99
f(y)=-1000*9.8100*y
*ICFD_BOUNDARY_FSWAVE
4,4,&H_0,&Ampl
*ICFD_DEFINE_WAVE_DAMPING
1 2 5
*ICFD_DEFINE_POINT
1 5 0.3 0
*ICFD_DEFINE_POINT
2 1 0 0
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
$ $
$ ICFD MESH KEYWORDS $
$ $
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
*MESH_VOLUME
$# volid
20
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
4 5 6 9 10
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
$ $
$ DATABASE (OUTPUT) $
$ $
$—+—-1—-+—-2—-+—-3—-+—-4—-+—-5—-+—-6—-+—-7—-+—-8
*ICFD_DATABASE_DRAG
10
*END

Figures

figure_001

Download

Description

This example features an incoming solitary wave impacting a floating box. It directly initializes the free surface level via an initial keyword rather than the classic way of meshing the initial free surface interface.