Shape optimization of a flow system is done to obtain the required effects, in the engineering fields. Most of these designs are accomplished by empirical or numerical analysis. In empirical analysis, it is difficult to obtain an optimal shape in the feasible design region. And, in numerical method, it usually needs much calculation expenses for shape optimization, because of design sensitivity analysis. In this study, we used the growth-strain method having only one distributed parameter such as a design variable. It optimizes a shape by making a distributed parameter such as dissipation energy uniform in a flow system, and then applied to two-flow systems. In order to overcome the stability occurred in numerical analysis performed by Azegami, the equation of volumic strain has been modified. Also, the shapes were compared with the known optimal shapes for the flow systems. Consequently, we confirm that the modified growth-strain method is very efficient and practical in shape optimization of the flow systems.