Se construyó un modelo bidimensional de un reactor de fusión, a partir de las analogías que existen entre la teoría de matrices aleatorias y el log-gas de Coulomb. El modelo consiste en N0, cargas confinadas en un círculo, realizando movimientos angulares dictaminados por un movimiento browniano. En cada iteración, se revisa si un par de partículas están separadas menos de una distancia de fusión crítica. En caso positivo, hay un evento de fusión y se extraen las dos partículas del sistema. Se hicieron simulaciones, para N0 = 50 y N0 = 100. Se encontró que el sistema evoluciona a un número finito de partículas NF, independiente de la temperatura β. Para N = 50, la fracción final de partículas tiende a 0.65. Para N = 100 la fracción final tiende a 0.45. Este modelo se propone como una primera aproximación a los reactores de fusión toroidales.