Recently, Faraday rotation measure (RM) at high Galactic latitude has been investigated, partly to explore the Galactic magnetic fields and partly to study the extragalactic magnetic fields. The Galactic contribution to RM comes from the global component as well as the turbulent component. So far the turbulent field was used to be analytically modeled with a Kolmogorov-type power spectrum. Here, we present the initial results of the work where the turbulent field is modeled using data of MHD turbulence simulations. Our work is intended to be applied to simulations of RM surveys with LOFAR, ASKAP, MeerKAT, and SKA.