The aim of this thesis is to study numerically two phase flow using accurate interface tracking method. First, a two phase flow homogeneous isotropic turbulence is performed. This numerical configuration is used to study two phase flows turbulence and the ELSA model used in primary atomization modelling. Based on these results, modelling constants and source terms have been determined and validated. Then the evaporation process is studied and modelized by using a passive scalar and then by using a DNS formalism. Energy and species equations are added in the ARCHER code. Validation of the DNS performed by comparing the DNS results with theorical solution, such as the D2 law. Finally, limitations and interests of this approach are discussed and further improvements are proposed.