Thermohydraulics have played an important role in a variety of industrial devices. In particular, constitutive equations have been developed to describe the interactions of gas and liquid at the interfaces of various flow patterns. Such behavior used to be solved numerically by a set of finite difference equations. However, numerical simulation of interfacial motion in a flow is still a incomplete in many cases. The cellular automaton method is expected to be an effective technique for simulating fluid motion numerically. Because of discretized space, time and variables, the discrete dynamic system can be defined, and the simple structure of dynamics allows it to be dealt with easily with a large scale of freedom. These enable us to realize discrete molecular dynamics on the computers. We have demonstrated that an automaton solution of two-component counter-flow results in turbulent behavior which scatters droplets at its interface. After evaluation of the turbulence intensity, velocity distribution of the averaged solution and an inertial energy, it can be concluded that this method is useful to simulate two-component flow.