Fluid pressure fringes
*KEYWORD *TITLE *BOUNDARY_PRESCRIBED_MOTION_RIGID_LOCAL *CONSTRAINED_EXTRA_NODES_NODE *CONSTRAINED_JOINT_CYLINDRICAL_ID *CONSTRAINED_JOINT_REVOLUTE_ID *CONSTRAINED_JOINT_STIFFNESS_GENERALIZED *CONTROL_ACCURACY *CONTROL_BULK_VISCOSITY *CONTROL_CONTACT *CONTROL_ENERGY *CONTROL_HOURGLASS *CONTROL_RIGID *CONTROL_SHELL *CONTROL_SOLID *CONTROL_TERMINATION *CONTROL_TIMESTEP *CONTROL_SOLUTION *DATABASE_BINARY_D3PLOT *DATABASE_BNDOUT *DATABASE_EXTENT_BINARY *DATABASE_FORMAT *DATABASE_GLSTAT *DATABASE_JNTFORC *DATABASE_MATSUM *DATABASE_RBDOUT *DEFINE_COORDINATE_NODES *DEFINE_CURVE *DEFINE_CURVE *DEFINE_CURVE_FUNCTION *DEFINE_VECTOR *ELEMENT_SHELL *ELEMENT_SOLID *INCLUDE *ICFD_BOUNDARY_FSI *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_MOVEMESH *ICFD_CONTROL_FSI *ICFD_CONTROL_MESH *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_DATABASE_TIMESTEP *ICFD_MAT *ICFD_SECTION *ICFD_PART *ICFD_PART *ICFD_PART *ICFD_PART *ICFD_PART *ICFD_PART *ICFD_PART *ICFD_PART_VOL *LOAD_BODY_Y *MAT_RIGID *MESH_BL *MESH_BL_SYM *MESH_INTERF *MESH_VOLUME *MESH_SURFACE_NODE *MESH_SURFACE_ELEMENT *NODE *PART_INERTIA *PART *SECTION_SHELL *SECTION_SOLID *END
$-----------------------------------------------------------------------------
$
$ Example provided by Iñaki (LSTC) and Nils Karajan (Dynamore)
$
$ 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
$===============================================================================
$ $
$ F L U I D P R O B L E M $
$ $
$===============================================================================
*KEYWORD
*TITLE
ICFD Tuned Liquid Damping Problem (SPHERIC Test 9)
*INCLUDE
mesh.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 cfl lcidsf dtmin dtmax
&T_end &dt_fluid &dt_min &dt_max
*ICFD_CONTROL_FSI
$# owc bt
0
*ICFD_CONTROL_OUTPUT
$# msglv
3
*ICFD_CONTROL_MESH
$# mgsf
1.05
$---+----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_oil &mu_oil
*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 1
*ICFD_PART
$# pid secid mid
5 1 2
*ICFD_PART
$# pid secid mid
6 1 2
*ICFD_PART
$# pid secid mid
7 1 2
*ICFD_PART
$# pid secid mid
8 1 2
*ICFD_PART_VOL
$# pid secid mid
10 1 1
$# spid1 spid2 spid3 spid4 spid5
1 2 3 4
*ICFD_PART_VOL
$# pid secid mid
11 1 2
$# spid1 spid2 spid3 spid4 spid5
4 5 6 7
$---+----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
2
*ICFD_BOUNDARY_NONSLIP
$# pid
3
*ICFD_BOUNDARY_NONSLIP
$# pid
5
*ICFD_BOUNDARY_NONSLIP
$# pid
6
*ICFD_BOUNDARY_NONSLIP
$# pid
7
*ICFD_BOUNDARY_FSI
$# pid
1
*ICFD_BOUNDARY_FSI
$# pid
2
*ICFD_BOUNDARY_FSI
$# pid
3
*ICFD_BOUNDARY_PRESCRIBED_MOVEMESH
$# pid dofx dofy dofz
5 0 0 1
*ICFD_BOUNDARY_FSI
$# pid
6
*ICFD_BOUNDARY_FSI
$# pid
7
$---+----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
1 2 3 5 6 7
*MESH_INTERF
$# volid
30
$# pid1
4
*MESH_BL
$# pid nelth blth blfe blst
6 1
*MESH_BL
$# pid nelth blth blfe blst
1 1
*MESH_BL_SYM
$# pid1
7
*MESH_BL_SYM
$# pid1
3
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_TIMESTEP
$# outlv
1
*END
The tank, initially in equilibrium, is free to rotate and its movement is induced by a weight that has a prescribed periodic movement along an initially horizontal rail attached to the rotation center of the tank which periodically displaces the center of gravity. The breaking waves and sloshing dynamics of the fluids placed inside the tank affect the damping characteristics of the system. It is recommended to use at least 16 CPUs for this problem.
References :
[1] G. Bulian, A. Souto-Iglesias, L. Delorme, and E. Botia-Vera, “SPH simulation of a tuned liquid damper with angular motion,” Journal of Hydraulic Research, vol. 48, no. Extra Issue, pp. 28–39, 2010.
[2] Botia-Vera, E., A. Souto-Iglesias, G. Bulian, and L. Lobovský (2010). Three SPH Novel Benchmark Test Cases for free surface flows. In 5th ERCOFTAC SPHERIC workshop on SPH applications.