Intermediate : SPHERIC Test 10
Fluid pressure fringes
*TITLE *KEYWORD *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *ICFD_BOUNDARY_NONSLIP *ICFD_CONTROL_IMPOSED_MOVE *ICFD_CONTROL_MESH *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_DATABASE_NODOUT *ICFD_DATABASE_TIMESTEP *ICFD_DEFINE_POINT *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE *LOAD_BODY_Y *MESH_BL *MESH_BL_SYM *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 R8.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: (kg/m/s)
$X------------------------------------------------------------------------------
$X
*KEYWORD
*TITLE
ICFD Sloshing 2
*INCLUDE
mesh3d_lat10_short.k
*INCLUDE
radial_vel_oil.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R T_end 6.0
R dt_plot 0.05
$
$--- Fluid
$
Rrho_fluid 900.
R mu_fluid 0.045
R dt_fluid 0.000
Rcfl_fluid 0.950
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 cfl lcidsf dtmin dtmax
&T_end &dt_fluid&cfl_fluid 0.05
*ICFD_CONTROL_OUTPUT
$# msgl
4
*ICFD_CONTROL_IMPOSED_MOVE
$# pid
0
$# alphal
100
$# ptid
1
*ICFD_CONTROL_MESH
$# mgsf
1.05
*ICFD_DEFINE_POINT
$# poid x y z
1 0 0 0
$---+----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 2
*ICFD_PART
$# pid secid mid
4 1 2
*ICFD_PART
$# pid secid mid
5 1 1
*ICFD_PART
$# pid secid mid
6 1 1
*ICFD_PART
$# pid secid mid
7 1 2
*ICFD_PART_VOL
$# pid secid mid
8 1 1
$# spid1 spid2 spid3 spid4
1 2 5 6
*ICFD_PART_VOL
$# pid secid mid
9 1 2
$# spid1 spid2 spid3 spid4
3 4 5 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
4
*ICFD_BOUNDARY_NONSLIP
$# pid
6
*ICFD_BOUNDARY_NONSLIP
$# pid
7
*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
$---+----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 4 6 7
*MESH_INTERF
$# volid
30
$# pid1
5
*MESH_BL
$# pid nelth blth blfe blst
6 1 0.002 1
*MESH_BL
$# pid nelth blth blfe blst
7 1 0.002 1
*MESH_BL_SYM
$# pid1
1
*MESH_BL_SYM
$# pid1
2
*MESH_BL_SYM
$# pid1
3
*MESH_BL_SYM
$# pid1
4
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_TIMESTEP
$# outlv
1
*ICFD_DATABASE_NODOUT
$# outlv dtout
1 0.01
$# nid1 nid2 nid3 nid4
1925 1923 1924 1926
*DATABASE_BINARY_D3PLOT
&dt_plot
*END
In this validation sloshing simulation, a partially filled tank is rolling with a prescribed angular motion. The impact pressure is compared with the experiment. Different fluids and different tank thicknesses are tried to capture the 3D effects. Warning : the mesh is fine, 8 CPUs or more are recommended for this case.
References :
[1] L. Delorme, A. Colagrossi, A. Souto-Iglesias, R. Zamora-Rodríguez, and E. Botia-Vera, “A set of canonical problems in sloshing. Part I: Pressure field in forced roll. Comparison between experimental results and SPH,” Ocean Engineering, vol. 36, no. 2, pp. 168–178, February 2009.
[2] A. Souto-Iglesias, E. Botia-Vera, A. Martin, and F. Pérez-Arribas, “A set of canonical problems in Sloshing. Part 0: Experimental setup and data processing,” Ocean Engineering, vol. 38, pp. 1823–1830, February 2011.
[3] Souto-Iglesias, A., E. Botia-Vera, and G. Bulian (2011, June). Repeatability and Two-Dimensionality of model scale sloshing impacts. In International Offshore and Polar Engineering Conference (ISOPE). The International Society of Offshore and Polar Engineers (ISOPE).
[4] 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.