Non-equilibrium molecular dynamics are used to generate the flow of polymer solutions, specifically of Boger fluids, through a planar 2:1:2 contraction-expansion geometry. The solvent molecules are represented by Lennard-Jones particles, while linear molecules are described by spring-monomers with a finite extensible non-linear elastic spring potential. The equations for Poiseuille flow are solved using a multiple time-scale algorithm extended to non-equilibrium situations. Simulations are performed at constant temperature using Nose-Hoover dynamics. At simulation conditions, changes in concentration show no significant effect on molecular conformation, velocity profiles, and stress fields, while variations in the Deborah number have a strong influence on fluid response. Increasing the magnitude of the Deborah number (De), larger deformation rates are developed in the flow region. For a Deborah number of one, the non-dimensional pressure drop presents values lower than the correspondent Newtonian case. However, for large Deborah numbers, the pressure drop increases above the Newtonian reference. An effective excess pressure drop above the Newtonian value is predicted for Boger fluids along this geometry.