Solvent design aiming at solution property induced surface stability: A case study using SiC solution growth

•A mesoscale simulation model for step motion is constructed.•Kinetic aspect of step bunching due to solute depletion is studied.•A non-dimensional number is proposed as a criterion of the tendency to step bunching.•A phase map for morphology stability is constructed as a pointer of solvent design.