Hi, apologies for the late reply. I have since been able to model the wave tank but now there appears to be another issue, which is the amount of particle mixing. With the piston (wave maker) deactivated, the only load which the particles are subjected to is a gravity loading, applied with the LOAD_BODY_Z keyword. In my mind, if there are no loads directly applied to the water within the tank then there should be no particle movement, but thats not the case. The particles appear to move spuradically and I do not know why. I have attached screenshots as well as my keyword input file (note that the piston is active in this file) to illustrate this. If you know how to correct this issue I would greatly appreciate a response. Also, please ignore part #6 contained in the file, it does not make contact with anything and is only there to serve as a visual aid on the wave height.
IMAGES WITH ONLY GRAVITY LOADING APPLIED
KEYWORD INPUT FILE
$# LS-DYNA Keyword file created by LS-PrePost(R) V4.9.12-16Dec2022
$# Created on Aug-8-2023 (13:38:55)
*KEYWORD
*TITLE
$# title
LS-DYNA keyword deck by LS-PrePost
*CONTROL_BULK_VISCOSITY
$# q1 q2 type btype tstype
0.0011.0000E-12 1 0 0
*CONTROL_SPH
$# ncbs boxid dt idim nmneigh form start maxv
1 11.00000E20 3 500 15 0.01.00000E15
$# cont deriv ini ishow ierod icont iavis isymp
0 0 0 0 0 0 1 100
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas nosol
15.0 0 0.0 0.01.000000E8 0
*CONTROL_UNITS
$# length time mass temp
m sec kg k
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
0.01 0 0 0 0
$# ioopt rate cutoff window type pset
0 0 0
*BOUNDARY_PRESCRIBED_MOTION_SET_ID
$# id heading
0Piston Motion
$# nsid dof vad lcid sf vid death birth
11 1 2 2 1.0 01.00000E28 0.0
*LOAD_BODY_Z
$# lcid sf lciddr xc yc zc cid
1 9.81 0 0.0 0.0 0.0 0
*PART
$# title
SphNode
$# pid secid mid eosid hgid grav adpopt tmid
1 1 1 1 0 0 0 0
*SECTION_SPH_TITLE
Water
$# secid cslh hmin hmax sphini death start sphkern
1 1.3 0.2 2.0 0.01.00000E20 0.0 0
*MAT_NULL_TITLE
Water
$# mid ro pc mu terod cerod ym pr
1 997.0 -10.0 0.001 0.0 0.0 0.0 0.0
*EOS_MURNAGHAN_TITLE
Water
$# eosid gamma k0 v0
1 7.0 65000.0 1.0
*PART
$# title
Piston
$# pid secid mid eosid hgid grav adpopt tmid
4 2 2 0 0 0 0 0
*SECTION_SHELL_TITLE
Piston
$# secid elform shrf nip propt qr/irid icomp setyp
2 16 1.0 2 1.0 0 0 1
$# t1 t2 t3 t4 nloc marea idof edgset
0.2 0.2 0.2 0.2 0.0 0.0 0.0 0
*MAT_RIGID_TITLE
Piston
$# mid ro e pr n couple m alias
2 100.0 20000.0 0.25 0.0 0.0 0.0
$# cmo con1 con2
1.0 5 7
$#lco or a1 a2 a3 v1 v2 v3
0.0
*PART
$# title
Container
$# pid secid mid eosid hgid grav adpopt tmid
5 2 4 0 0 0 0 0
*MAT_RIGID_TITLE
Container V2
$# mid ro e pr n couple m alias
4 100.0 20000.0 0.25 0.0 0.0 0.0
$# cmo con1 con2
1.0 7 7
$#lco or a1 a2 a3 v1 v2 v3
0.0
*PART
$# title
Wave Height
$# pid secid mid eosid hgid grav adpopt tmid
6 2 4 0 0 0 0 0
*HOURGLASS_TITLE
Water
$# hgid ihq qm ibq q1 q2 qb/vdc qw
1 11.00000E-6 0 1.5 0.06 0.1 0.1
*MAT_RIGID_TITLE
Container
$# mid ro e pr n couple m alias
3 100.0 20000.0 0.25 0.0 0.0 0.0
$# cmo con1 con2
1.0 5 7
$#lco or a1 a2 a3 v1 v2 v3
0.0
*DEFINE_BOX
$# boxid xmn xmx ymn ymx zmn zmx
1 -1.0 12.0 -1.0 1.0 -1.0 2.0
*DEFINE_CURVE_TITLE
Gravity
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 1.0
100.0 1.0
*DEFINE_CURVE_TITLE
Piston Motion
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 0.6 0.0795775 0.0 0.0 0 0
$# a1 o1
0.0 0.0
1.0 1.0
2.0 0.0
3.0 1.0
4.0 0.0
5.0 1.0
6.0 0.0
7.0 1.0
8.0 0.0
9.0 1.0
10.0 0.0
11.0 1.0
12.0 0.0
13.0 1.0
14.0 0.0
15.0 1.0
16.0 0.0
17.0 1.0
18.0 0.0
19.0 1.0
20.0 0.0
21.0 1.0
22.0 0.0
23.0 1.0
24.0 0.0
25.0 1.0
26.0 0.0
27.0 1.0
28.0 0.0
29.0 1.0
30.0 0.0