Rough and textured surfaces are of paramount importance for lubrication, both in nature and in technology. While surface roughness relevantly influences both friction and wear, artificial surface texturing improves the performance of slider bearings as an energy efficiency action. The simulation of hydrodynamic lubrication by taking into account complex surfaces as boundaries requires the use of computational fluid dynamics (CFD) software able to predict the pressure and the velocity profile through the thickness of the fluid and at any point within the 3D domain. In the present study, a CFD-smoothed particle hydrodynamics (SPH) code is applied to simulate hydrodynamic lubrication for a linear slider bearing in the presence of a 3D rough surface, showing the capabilities of CFD-SPH in modelling such complex interaction phenomena. Numerical assessments involve the load capacity, the 3D fields of the velocity vector, and the pressure 3D field (both within the fluid domain and at the fluid-plate interface).Lubricants 2019, 7, 103 2 of 12 Based on the aforementioned theoretical framework, Williams and Symmons [3] developed a 1D computational fluid dynamics (CFD)-finite difference (FD) model to numerically reproduce the pressure longitudinal profiles within the fluid film of a linear slider. Dobrica and Fillon [4] developed a 2D CFD-finite volume method (FVM) code, alternatively using Navier-Stokes equations and Reynolds' equation for fluid films, and they validated it on the Rayleigh step bearings. They highlighted the importance of modelling the inertia terms, neglected by Reynolds' equation for fluid films. For step bearings, Vakilian et al. [5] found that neglecting the inertia terms in the momentum equations is responsible for underestimations at the leading edge and over-predictions at the trailing edge on the pressure field.Almqvist et al. [6] provided inter-comparisons between a FD model based on Reynolds' equation for fluid films and a commercial CFD-FVM code based on Navier-Stokes equations. Almqvist et al. [6,7] also provided report analytical solutions on pressure longitudinal profiles, velocity vertical profiles, friction force, and load-bearing capacity (the frictional coefficient is the ratio between these two forces) for both a linear slider and the Rayleigh step slider, with null Dirichlet's boundary conditions for pressure. They also derived the optimal geometric configuration for a linear slider to maximize the load capacity (e.g., [3]), also referred to as load-bearing capacity or load-carrying capacity. Further, they analysed the effects of the surface roughness by means of the homogenization technique [8].Rahmani et al. [9] presented an analytical approach based on Reynolds' equation for asymmetric partially textured slider bearings with surface discontinuities, to optimize the choice of the textures parameters with respect to the load capacity and the friction force.Papadopoulos et al.[10] used a 2D CFD-FVM code to optimize micro-thrust bearings with surface texturing by means of numerical inter...