Artificial Viscosity vs. Laminar+SPS issues
From what I understood, DualSPHyiscs has two ways to model turbulent flows including dissipation: 1) Artificial viscosity which is dependent on a factor of alpha which is set in the input files, and 2) Laminar+SPS which implements Smagorinsky's eddy viscosity model and is dependent on Cs (Smagorinsky coefficient of 0.12 set in JSph.cpp) and CI (Blin constant of 0.0066 also set in JSph.cpp) (see 3.2.1 and 3.2.2 - https://github.com/DualSPHysics/DualSPHysics/wiki/3.-SPH-formulation).
I'm trying to implement the Laminar+SPS viscosity in the input file of the Dambreak 2D example. But in the input file when I switch the viscotreatment from 1 (artificial viscosity - default setting) to 2 (Laminar + SPS), my model is still sensitive to alpha and the viscobound factor. The alpha term is only present in the artificial viscosity equations (see Eq. 8 - https://github.com/DualSPHysics/DualSPHysics/wiki/3.-SPH-formulation) so theoretically if I switch to Laminar+SPS, alpha should not matter.
Additionally, it is sensitive to the viscobound factor, but I'm not sure exactly what the viscobound factor refers to, since I cannot find it in the DualSPHyiscs formulation nor which viscosity treatment it is associated with.