---

In this study, the stochastic field method is developed for the direct numerical simulation of turbulent drag reduction by microfibers. For this purpose, the governing equations without any simplification are discretized on an Eulerian grid. A fifth-order upwind scheme is used for the discretization. A Monte-Carlo method is employed in the conformation space. Then, three-dimensional, time-dependent Navier-Stokes equations for the incompressible flow of a non-Newtonian fluid are numerically solved for a turbulent channel flow. Statistical quantities obtained by the proposed method are compared with those of a Lagrangian method and the high precision of the new method is demonstrated. The main advantage of the new method is its low computational cost.