“…A constraint function is applied to ensure that all the nodes at the top of the pile have the same vertical displacement. The size of the mesh elements and the domain are selected to satisfy classical rules
24 : per Rayleigh wavelength,
, there should be a minimum of 10 elements to have a high spatial resolution of propagating waves, and the minimum domain size should be greater than
to prevent reflections at the boundaries. The mesh widths are 120 m for
, 190 m for
and 380 m for
.…”