Axial piston machines find their use over a wide range of the power spectrum owing to their superior reliability, efficiency, and power density. They are also a key component in applications like reverse osmosis and firefighting wherein the working fluid is water. Utilizing low viscous fluid, such as water, as a working fluid poses challenges in designing the critical lubricating interfaces of the piston pumps. Specifically, low viscosity makes it difficult for the lubricating interfaces to provide sufficient bearing and sealing functions in challenging operating conditions. In order to maintain the lubricating interface performance in water hydraulic piston pumps, costly materials, and tight manufacturing tolerances are often utilized. To improve the efficiency and cost-effectiveness of these pumps, accurate numerical simulation tools that consider the fluid and structure interaction are needed to provide valuable insights into these lubricating interfaces. Although the Reynolds equation is a reliable method for determining the fluid pressure distribution in an oil-based piston pump, it assumes a laminar flow which may not be applicable to water piston machines. For example, in an inclined piston/cylinder interface of a water hydraulic pump, there may be regions in the film wherein the large gap height combined with the low viscosity of water induce turbulence effects. If the traditional Reynolds equation is used in such a scenario, it is likely to overestimate the leakage flow through the interface as it does not account for turbulence. Therefore, it is important to incorporate the effect of turbulence in the diffusive terms of the Reynolds equation to accurately describe the Poiseuille flow with high Reynolds numbers. The challenge is further compounded by the micromotion and deformation of the solid body, resulting in the unevenness of the gap height in the lubricating film. Therefore, the consideration of turbulence can only be applied regionally in such cases. The current study proposes a fluid-structure interaction model with the consideration of the localized turbulent effects. This modeling approach is applied to the piston/cylinder interface of an axial piston machine that uses water as the working fluid. The approach stems from the modification of the Poiseuille term to incorporate a function of the Reynolds number. The fluid dynamics considering the turbulence effect was validated against the solution of the Navier Stokes equation using commercial CFD software. The modified Reynold equation was implemented in an axial piston pump EHL model coupled with the multi-body dynamics. The simulation results from the novel pump model were compared to a measurement and the accuracy of the proposed model was found largely improved from the traditional laminar solution. The calculated flow rate was found to be 54.6% lower with the additional consideration of the turbulence effect in the studied case.