Neste trabalho apresentamos uma revisão de métodos de simulação de energia eletrostática de sistemas de cargas e uma proposta de adaptação de algoritmo ultilizado na literatura de sistemas gravitacionais para estudo das propriedades estatísticas de sistemas coulombianos. Na primeira parte do estudo, revisamos os fundamentos teóricos do método de Ewald e suas condições de aplicabilidade, procurando esclarecer as referências mais importantes no assunto, que são de difícil compreensão, gerando equívocos na utilização do termo de dipolo. Detalhamos o estudo sobre a análise da convergência da série em que a técnica se baseia, bem como sua interpretação física mostrando a equivalência entre as duas abordagens . Na segunda parte do trabalho analisamos os fundamentos do Fast Multipole Method desenvolvido para interação gravitacional, para o qual construímos programas em linguagem C para uma versão na rede. Criamos um algoritmo que denominamos Fast Multipole Monte Carlo (FMMC) e desenvolvemos um programa para cálculo das propriedades termodinâmicas de sistemas coulombianos. Os programas são testados comparando resultados para a energia e propriedades térmicas do modelo LRPM com resultados de simulação através de cálculo direto.