Исследование формирования ближнего порядка в сплаве объемно-центрированного кубического Fe-Ga проведено с помощью Монте-Карло моделирования с использованием эффективных потенциалов взаимодействия, рассчитанных методами теории функционала плотности в ферромагнитном и парамагнитном состояниях. Установлено, что при концентрациях Ga, близких к границе двухфазной области, при T < T c формируется ярко выраженный ближний порядок типа D0 3, а при T < T c определенного ближнего порядка не наблюдается. Полученные результаты согласуются с экспериментальными наблюдениями методом рентгеноструктурного анализа. Обсуждается связь особенностей ближнего порядка в сплаве Fe-Ga с величиной магнитострикции.