Implementação da distribuição de Maxwell

Com o intuito de aproximar melhor a simulação da teoria mais recente sobre difusão, o modo como o módulo das velocidades iniciais das moléculas é gerado foi alterado. Inicialmente, calculava-se um módulo de velocidade com base na equação $mv^2 = 3kT$, onde $m$ é a massa de uma molécula, $v$, sua velocidade, $k$, a constante de Boltzmann, e $T$, a temperatura ambiente. Porém, esta equação determina apenas a velocidade média de cada molécula num gás. O modelo cinético do gás prevê que as moléculas sigam uma distribuição de velocidade, conhecida como ``distribuição de velocidades de Maxwell''. Isso decorre da equipartição de energia quando uma colisão entre duas moléculas ocorre.

Aqui há dois pontos importantes a serem observados com relação à simulação. O primeiro é que está sendo assumido que a substância em difusão é um gás perfeito, que segue o modelo cinético. O segundo é que está se assumindo que as colisões que envolvem as moléculas são perfeitamente elásticas, ou seja, não há transferência de energia entre as moléculas e as paredes do meio de difusão, o que representaria um resfriamento no ambiente de difusão, e não há perda de energia na colisão entre duas moléculas.