It is a well known fact that the upper We limit encountered in continuum level viscoelastic flow simulations with typical constitutive equations for dilute polymeric solutions that predict bounded extensional viscosities in geometries with internal stagnation points on solid surfaces is a strong function of the strain hardening nature of the fluid. Specifically, the upper We limit decreases as the level of strain hardening of the fluid is increased. To provide further insight into the high We limitations in continuum level computations of this class of flows, we have performed extensive continuum and multiscale flow simulations in two benchmark flow problems, namely sedimentation of sphere in a tube and flow past a cylinder in a channel, utilizing the FENE-P (continuum and Brownian Configuration Fields (BCF)) and Giesekus (continuum) constitutive equations as well as the FENE (BCF) dumbbell micromechanical model. Extremely large stress gradients in the axial normal stress along the plane of symmetry in the wake of the cylinder and sphere are observed in FENE-P and Giesekus predictions using both multiscale and continuum numerical techniques at We of 0(1) where the numerical simulations begin to breakdown for significantly strain hardening fluids, i.e., b > 300, alpha < 0.005. The existence of very large localized polymeric stresses and stress gradients is shown to be the consequence of significant over prediction of macromolecular extension by the closed form constitutive equations in the extensionally dominated region of the flow. The inability of the aforementioned constitutive equations to accurately describe the flow microstructure coupling downstream of stagnation points in complex kinematics flows should motivate use of regularization techniques for the polymeric stress or utilization of more sophisticated constitutive equations in flows with strong straining components. Published by Elsevier B.V.