A rigours mathematical model based on the bivariate population balance frame work (the base of LLECMOD ''Liquid-Liquid Extraction Column Module'') for the steady state and dynamic simulation of pulsed liquid-liquid extraction columns is developed. The model simulates the coupled hydrodynamic and mass transfer for pulsed (packed and sieve plate) extraction columns. It is implemented using visual digital FORTRAN and then integrated into the LLECMOD program. Experimental validated correlations are used for the estimation of the droplet terminal velocity in extraction columns based on single and swarm droplet experiments in laboratory scale devices. Additionally, recently published correlations for turbulent energy dissipation, droplet breakage and coalescence frequencies are discussed as being used in this version of LLECMOD. In a case study, LLECMOD is used here to simulate the steady state performance of pulsed extraction columns with two chemical test systems recommended by the European Federation of Chemical Engineering (water-acetone-n-butyl acetate and water-acetone-toluene) and an industrial test system. Model predictions are successfully validated against steady state and transient experimental data, where good agreements are achieved. The simulated results (holdup, mean droplet diameter and mass transfer profiles) compared to the experimental data show that LLECMOD is a powerful simulation tool, which can efficiently predict the dynamic and steady state performance of pulsed extraction columns.