This work proposes a numerical model that incorporates the effect of lubricant inertia on the hydrodynamic pressure distribution, fluid film reaction forces, and the fluid velocity component profiles for finite-length open-ended squeeze film dampers (SFDs). Firstly, the thin film flow equations for the SFD in presence of fluid inertia effects are introduced. Furthermore, a small first-order perturbation by means of the expressions for the fluid film velocity components and the lubricant pressure distribution that are expanded in power series of the squeeze film Reynolds number is applied to the flow equations. Subsequently the developed lubricant flow equations are solved to develop expressions for the velocity component profiles and the hydrodynamic pressure distribution in SFDs. The pressure expression is numerically solved by using Gauss-Seidel method with finite difference discretization. Moreover, the fluid film reaction forces are determined by numerically integrating the hydrodynamic pressure expressions over the journal surface. Additionally, the proposed pressure distribution expression and the numerical SFD forces are incorporated into a simulation model and the simulation results are compared with the existing models in the literature under different operating conditions, including eccentricity ratios and inertia effects (i.e., Reynolds numbers). The simulation results demonstrate the significant influence of both convective and temporal (i.e., unsteady) lubricant inertia terms on the SFD hydrodynamic pressure distribution and the fluid film reaction forces. Furthermore, the proposed SFD model is incorporated into a multi-mass flexible rotordynamic model to evaluate the effect of SFD fluid inertia on the mass unbalance induced steady-state vibrations of the rotor and the nodal transient orbits by implementing finite element method and transient modal integration with predictor-corrector solver. The results of the analysis demonstrate the significant effect of fluid inertia on the resonance frequencies of the rotor and the steady-state vibration amplitudes and the transient orbits at the resonance zone.