Flow simulations in pipelined channels and several kinds of parametrized configurations have a growing interest in many life sciences and industrial applications. Applications may be found in the analysis of the blood flow in specific compartments of the circulatory system that can be represented as a combination of few deformed vessels from reference ones, e.g. pipes. We propose a solution approach that is particularly suitable for the study of internal flows in hierarchical parametrized geometries. The main motivation is for applications requiring rapid and reliable numerical simulations of problems in domains involving parametrized complex geometries. The classical reduced basis (RB) method is very effective to address viscous flows equations in parametrized geometries (see, e.g., [10]). An interesting alternative foresees a combination of RB with a domain decomposition approach. In this respect, preliminary efforts to reduce the global parametrized problem to local ones have led to the introduction of the so-called reduced basis element method to solve the Stokes problem [6], and more recently to the reduced basis hybrid method [3] and to the static condensation method [7]. In general, we are interested in defining a method able to maintain the flexibility of dealing with arbitrary combinations of subdomains and several geometrical deformations of the latter. A further new contribution to this field is the computation of the reduced basis functions through an optimization greedy algorithm [11].