Fluid pressure fringes
*KEYWORD *TITLE *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_SOLUTION *CONTROL_TERMINATION *CONTROL_TIMESTEP *DATABASE_BINARY_D3PLOT *DEFINE_CURVE *ELEMENT_SOLID *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_FSWAVE *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *ICFD_SOLVER_TOL_MOM *ICFD_SOLVER_TOL_PRE *INCLUDE *LOAD_BODY *MAT_ELASTIC *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *PART *SECTION_SOLID *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 R10.1 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD Partially submerged pillar
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 40.
R dt_plot 0.25
$
$--- Fluid
$
Rrho_fluid 1000
R mu_fluid 0.001
R dt_fluid 0.035
R grav 9.81
R fsi_bt 0.07
*INCLUDE
mesh.k
*INCLUDE
struc.k
$---+----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_FSI
$# owc bt
0 &fsi_bt
*ICFD_CONTROL_OUTPUT
$# msglv
3
*ICFD_CONTROL_TURBULENCE
2
0.13
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD SOLVER TOLERANCES $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_SOLVER_TOL_MOM
1.e-6 1.e-6
*ICFD_SOLVER_TOL_PRE
1.e-6 1.e-6
$---+----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
1 1 4
*ICFD_PART
$# pid secid mid
2 1 4
*ICFD_PART
$# pid secid mid
3 1 4
*ICFD_PART
$# pid secid mid
4 1 1
*ICFD_PART
$# pid secid mid
5 1 4
*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 1
*ICFD_PART
$# pid secid mid
10 1 4
*ICFD_PART
$# pid secid mid
11 1 4
*ICFD_PART_VOL
12 1 1
2, 4, 6, 7, 8, 9
*ICFD_PART_VOL
13 1 4
1, 2, 3, 5, 10, 11
*ICFD_MAT
$# mid flg ro vis
1 0
*ICFD_MAT
$# mid flg ro vis
4 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_BOUNDARY_FREESLIP
1
*ICFD_BOUNDARY_FREESLIP
3
*ICFD_BOUNDARY_NONSLIP
4
*ICFD_BOUNDARY_FSI
4
*ICFD_BOUNDARY_NONSLIP
5
*ICFD_BOUNDARY_FSI
5
*ICFD_BOUNDARY_FREESLIP
6
*ICFD_BOUNDARY_FSWAVE
7 2 2. 0.35 4. 20.
*ICFD_BOUNDARY_FREESLIP
8
*ICFD_BOUNDARY_FREESLIP
9
*ICFD_BOUNDARY_FREESLIP
10
*ICFD_BOUNDARY_FSWAVE
11 2 2. 0.35 4. 20.
*LOAD_BODY_Z
$# lcid sf
100 1
*DEFINE_CURVE_TITLE
Gravity force
$# lcid sidr sfa sfo offa offo dattyp
100 &grav
$# a1 o1
0.0 1.0
10000.0 1.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
20
1, 3, 4, 5, 6, 7, 8, 9
10, 11
*MESH_INTERF
20
2
*MESH_BL_SYM
1
*MESH_BL
4,1
*MESH_BL
5,1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_DRAG
5
*ICFD_DATABASE_DRAG
4
*DATABASE_BINARY_D3PLOT
&dt_plot
*DEFINE_CURVE
1 , 0, 1. , 1. , 0. , 0.
0.,5
1e3,5
*DEFINE_CURVE
2 , 0, 1. , 1. , 0. , 0.
0.,0.
1e3,0.
*END
In this example, a partially submerged pillar vibrates due to the periodic incoming waves. Strong FSI coupling is mandatory since the densities of the solid and the fluid are very close.