Nuclear reactions in the plasma environment can be substantially different from those in conventional laboratory non-plasma cases, which have attracted considerable attention in the fields of fusion and astrophysics. To self-consistently model the nuclear reaction process during plasma dynamic evolution, an extended nuclear reaction calculation module is developed and included in two-dimensional particle-in-cell simulations. Through the self-consistent simulations, we systematically show that, apart from the plasma screening, the kinetic Weibel instability (WI) occurring in plasmas also results in significant enhancement of nuclear reactions, where the self-generated magnetic fields play a key role. Specifically, the self-generated magnetic fields in WI deflect ion motions, decreasing the relative velocity, and convert plasma kinetic energy to thermal energy, increasing the ion temperature. The simulation results show that, for the reaction with a sharp resonance peak in the cross section, the reaction product yield is enhanced four times due to the WI. For nuclear reactions that have more prominent resonance peaks in the cross section, like , it is expected that such enhancements can reach up to one or several orders of magnitude.