As it is known, the Womersley function models velocity as a function of radius and time. It has been widely used to simulate the pulsatile blood flow through circular ducts. In this context, the present study is focused on the introduction of a simple function as an approximation of the Womersley function in order to evaluate its accuracy. This approximation consists of a simple quadratic function, suitable to be implemented in most commercial and non-commercial computational fluid dynamics codes, without the aid of external mathematical libraries. The Womersley function and the new function have been implemented here as boundary conditions in OpenFOAM ESI software (v.1906). The discrepancy between the obtained results proved to be within 0.7%, which fully validates the calculation approach implemented here. This approach is valid when a simplified analysis of the system is pointed out, in which flow reversals are not contemplated.