Uso de la programación Fortran para implementar métodos de simulación de dinámica molecular para estudiar las características del flujo de microfluidos

subrutina?init

sumv=0

sumv2=0

do?i=1,npart

x(i)=celosía-pos(i)

v(i)=(ranf(?)-0.5)

sumv=sumv+v(i)

sumv2=sumv2+v(i)**2

end?do

sum?=?sumv/npart

fs?=? sqrt(3*temp/sumv2)!?Factor de corrección de velocidad

do?i=1,npart

v(i)=(v(i)-sumv)* fs ?Corrección de velocidad

x(i)=x(i)-v(i)*dt?

end?do

regresar

end

Del código anterior, la fórmula sigue siendo muy simple. Puedes compararla una por una con las fórmulas del libro.