To simulate the dynamics of two-dimensional dam-break flow on a dry horizontal bed, we use a smoothed particle hydrodynamics model implementing two advanced boundary treatment techniques: (i) a semi-analytical approach, based on the computation of volume integrals within the truncated portions of the kernel supports at boundaries and (ii) an extension of the ghost-particle boundary method for mobile boundaries, adapted to free-slip conditions. The trends of the free surface along the channel, and of the impact wave pressures on the downstream vertical wall, were first validated against an experimental case study and then compared with other numerical solutions. The two boundary treatment schemes accurately predicted the overall shape of the primary wave front advancing along the dry bed until its impact with the downstream vertical wall. Compared to data from numerical models in the literature, the present results showed a closer fit to an experimental secondary wave, reflected by the downstream wall and characterized by complex vortex structures. The results showed the reliability of both the proposed boundary condition schemes in resolving violent wave breaking and impact events of a practical dam-break application, producing smooth pressure fields and accurately predicting pressure and water level peaks. compared to 3D CFD (computational fluid dynamics) equations [1,2]. However, given their meshing features, assumptions regarding hydrostatic pressure, and neglect of vertical velocity and flow velocity uniformity along the vertical axis [3,4], traditional numerical solutions to SWEs cannot effectively reproduce the strong distortion of the free surface that occurs in dam-break flows.Mesh-free models, such as smoothed particle hydrodynamics (SPH), offer an alternative to solving SWEs with grid-based numerical methods and have, accordingly, been applied to several areas of computational fluid dynamics [5,6]. The SPH method presents different advantages: mesh deformation and cracking; calculation of the system's advection and transport (due to its Lagrangian nature); modeling of free surface and phase/fluid interface problems; the ability to manage very large deformations in high-energy phenomena (e.g., explosions, high-velocity impacts, and penetrations); applicability at multiple scales if coupled with molecular dynamics and dissipative particle dynamics; and greater suitability to 3D-modeling than mesh-based methods [7,8]. In addition, should a "weakly compressible" approach be adopted, the non-iterative SPH algorithms are simplified [8].The limits of SPH models are linked to slightly higher computational costs, complex procedures for the local refining of spatial resolution, and a relatively low accuracy for classical CFD applications where mesh-based models are well-validated (e.g., confined mono-phase flows [8]). However, they are effectively employed in several fields [9][10][11][12][13][14][15][16][17][18].Various flood propagation studies have focused on the use of the SPH method to investigate dam-break ...