The characterization method of flow function is presented in Supplementary Note 11. Since the opening and closing time of solenoid valves cannot be ignored under high-frequency actions, we refer to the work of Xu et al. \cite{Xu.2017} and propose a time-based piecewise linear open degree function \(\varGamma \left( t \right) \)  to obtain the actual mass flow  through correction.