The process of bioethanol production has been characterized with a structured and nonsegregated form of yeast growth dynamics. In this work, a geometric numerical method is applied to obtain the approximate solution of the oscillatory phenomena transpiring in the process of bioethanol production. This method is called group preserving scheme which is based on Lie group, proper for solving ordinary differential equations. In this regard, The Minkowski Cayley transformation is used to create this numerical method to get the approximate solutions of the problems. Moreover, figures are provided to show the reliability and accuracy of the proposed method.