We deal with modelling of flows in channels or tubes with abrupt changes of the diameter. The goal of this work is to construct the FEM solution in the vicinity of these corners as precise as desired. We present two ways. The first approach makes use of a posteriori error estimates and the adaptive strategy. The second approach is based on the asymptotic behaviour of the exact solution in the vicinity of the corner and on the a priori error estimate of the FEM solution. Then we obtain the solution with desired precision also in the vicinity of the corner, though there is a singularity. Numerical results are demostrated on a 2D example.