O propósito desta tese é estudar um método numérico indireto para a solução de problemas de otimização de trajetórias espaciais em campo gravitacional central Newtoniano. Neste método, a resolução do problema de valor de contorno em dois pontos é baseada no conceito de variação segunda, motivo pelo qual ele é denominado de método da variação segunda.Inicialmente, este trabalho apresenta alguns fundamentos e resultados matemáticos da Teoria de Controle Ótimo. Estes conceitos são então utilizados no desenvolvimento matemático e formulação do método da variação segunda, incluindo sua extensão para problemas com vínculos nas variáveis de controle e a descrição do algoritmo final. Este algoritmo é então aplicado a problemas de transferências espaciais ótimas entre órbitas circulares coplanares próximas, considerando dois modelos clássicos de sistemas propulsivos: sistemas a velocidade de ejeção constante (VEC) e empuxo limitado; e sistemas com velocidade de ejeção modulável e potência limitada (PL). De modo geral, os resultados obtidos para os problemas são muito bons. Para o sistema PL, há grande concordância deles com resultados analíticos, provenientes de uma teoria linear, e numéricos, obtidos através de outros métodos. Para o sistema VEC, os resultados estão em acordo com a teoria de transferências orbitais ótimas e verifica-se que as duas formulações utilizadas, linearizada e não-linearizada, fornecem resultados bastante próximos, para problemas de transferência entre órbitas próximas.