A high ion temperature plasma in the Large Helical Device is examined in the case in which the ion temperature gradient mode is unstable. The nonlinear gyro-kinetic simulation is performed to evaluate the turbulent ion heat diffusivity with the kinetic electron response. It is clarified that the decay time of zonal flows [S. FerrandoMargalet et al., Phys. Plasmas 14, 122505 (2007)] decreases radially outward due to the trapped electron and the ion energy transport increases outward. To reduce the computational cost for applying to the dynamical transport simulation, an extended transport model for the ion heat diffusivity in terms of the mixing length estimate and the characteristic quantity for the linear response of zonal flows is proposed. Turbulent transport is one of the most critical issues for plasma confinement in magnetic fusion devices. Theoretical expressions for the turbulent transport due to the various instabilities were reviewed [1]. Recently, a large number of the gyro-kinetic simulations which are applied to the turbulent transport have been performed in toroidal plasmas, e.g. [2][3][4][5][6][7][8]. The gyro-kinetic analysis results in tokamak [9][10][11] and helical [12][13][14] plasmas have been studied with the experimental observations. Gyro-kinetic simulations of helical plasmas require a large number of mesh points along the field line in order to capture the helical ripple structure. Therefore, the gyro-kinetic analysis in helical plasmas consumes the larger computer resources than for tokamaks. In helical plasmas, it is not practical to perform the nonlinear gyro-kinetic simulation at each time step of the dynamical transport code. The reduced model, which reproduces the nonlinear gyro-kinetic analysis results, is needed for the dynamical transport simulation in helical plasmas.The GKV code [15] has been used to examine the ion temperature gradient (ITG) mode and zonal flows in the Large Helical Device (LHD) for studying the turbulent transport [13]. The gyro-kinetic simulation with the adiabatic electron is performed for the high ion temperature LHD discharge (shot number 88343 [16]). The ion energy flux by the ITG mode instability agrees with the experimental results [13,14]. The reduced model for the ion heat diffusivity is proposed [14] by the simulation with the adiabatic electron. This reduced model is the function of the author's e-mail: toda@nifs.ac.jp linear growth rate for the ITG mode and the zonal flow decay time [17,18]. In tokamak plasmas, the gyro-kinetic analysis at each time step is globally performed in the dynamical transport simulation [19,20]. How to apply the reduced model of the turbulent ion heat diffusivity from the gyro-kinetic simulation with the adiabatic electron to the transport code has been shown in helical plasmas [21]. The simulation with the kinetic electron shows the larger ion energy flux than the experimental results in the high-T i LHD discharge [22]. On the other hand, the electron and ion energy fluxes of the simulation results with the kinetic el...