In this predominately predictive modelling finite volume/element study, a comparative analysis is performed for time-dependent and viscoelastoplastic flow in a circular contraction-expansion geometry of aspect-ratio 10:1:10. For this, a hybrid finite volume/element scheme is employed. A new and revised micellar model is investigated, under the denomination of BMP+_p, which reflects a bounded extensional viscosity response and an N1Shear-upturn at large deformation rates (lost in earlier model-variants), a versatile model capable of supporting plasticity, shear-thinning, strain softeninghardening and shear-banding. Many of these features are common to wormlike micellar and polymer solutions. Then, findings are contrasted against a De Souza model. Two flow regimes are addressed: plastic flow (low flow-rate Q≤1 units, solvent-fraction <10-2) and viscoelastic flow (larger-Q>1; minimised plasticity; =1/9); as quantified via flow-structure, yield-fronts and pressure-drops. Under the plastic regime, elasticity-increase causes asymmetry about the contraction-plane, whilst yieldstress and enhanced strain-hardening promote solid-like features, apparent through augmented unyielded-regions and rising pressure-drops. Concerning the viscoelastic regime and vortexstructures, extensional-deformation experienced correlates with hardening expectation in uniaxialextension, whilst streamline activity in vortex-cells correlates with normal-stress response in shear. Adjustment in strain-hardening/softening response with Q-rise, provides translation from weaker salient-corner vortex centres to stronger elastic corner-vortices; yet, when softening finally prevails, asymmetric upstream/downstream salient-corners vortex patterns are recovered. For strong-hardening and solvent-dominated ~0.8 fluids (as with Boger fluids), an intermediate lip-vortex-formation phase is noted, alongside coexistence of salient-corner vortices. Such a vortex-coexistence phase is distinctly absent in solute-concentrated fluids.