为了实现聚变点火条件，必须对等离子体进行辅助加热，中性束注入加热作为最有效的辅助加热手段之一，对其加热过程进行数值模拟是非常必要的。本文系统阐述了中性束加热的基本原理，介绍了中性束注入加热的五个过程：电离过程、衰减过程、漂移运动过程和快离子慢化过程。本文通过建立物理模型、采用合理的简化手段、利用Gill方法数值求解了一组中性束注入加热的基本方程组，模拟中性束注入加热过程中等离子体的电子和离子温度演化过程,计算结果与现有实验结果相符。进一步调整计算参数讨论不同注入中性束条件下的加热效率，计算结果表明：提高注入中性束的能量和等效流强能够提高注入后等离子体的电子、离子温度和电子、离子温度的爬升速度，使等离子体中的电子、离子以更快的速度加热到更高的温度。有限体积法是当前求解流动和传热问题最成功的数值计算方法，使用有限体积法数值求解Fokker-Planck方程，能够研究中性束注入加热等离子体过程中本底粒子的分布函数演化情况。本文叙述了求解Fokker-Planck方程的几种方法，分析了这几种方法的特点。介绍了求解偏微分方程的常见四种方法：有限差分法、有限元法、边界元法和有限体积法，分析了这四种方法的优缺点。利用有限体积法得到一维Fokker-Planck方程计算格式，计算中能够保证粒子数密度守恒。通过数值模拟一维超高斯分布电子自碰撞趋于麦克斯韦分布过程，证明有限体积法能够高效求解Fokker-Planck方程，求解过程能够保证粒子分布函数的非负性、粒子数密度守恒和能量守恒。进一步开发了求解二维Fokker-Planck方程的计算程序，用于数值模拟具体的中性束加热物理过程，目前程序还需要进一步调试。