  In this work, we showed a calculation method of local stress tensor
applicable to non-equilibrium MD systems based on the Method of Plane (MoP).
From the relation between the macroscopic velocity distribution function and
the microscopic molecular passage across a fixed control plane, we derived a
method to calculate the basic properties of the macroscopic momentum
conservation law including the density, the velocity, the momentum flux, the
interaction and kinetic terms of the stress tensor defined on a surface with a
finite area. Any component of the streaming velocity can be obtained on a
control surface, which enables the separation of the kinetic momentum flux into
the advection and stress terms in the framework of MoP. We verified the present
method through the extraction of the density, velocity and stress distributions
in a quasi-1D steady-state Couette flow system and in a quasi-2D steady-state
system with a moving contact line. In our method, as opposed to volume average
method, the density, mass and momentum fluxes are defined on a surface, which
is essential to be consistent with the mass and momentum conservation laws in
dynamic systems.
