El artículo considera la modelización numérica de los procesos de condensación y evaporación homogénea y heterogénea en flujos multifásicos utilizando el método de momentos. Los procesos no estacionarios de dinámica de gases y transiciones de fase en el plano bidimensional y en regiones axisimétricas se describen mediante un sistema general de ecuaciones. El sistema de ecuaciones se amplía añadiendo dos ecuaciones. Una describe la evolución de la fracción de masa total de la sustancia que se condensa; la otra describe la evolución de la fracción de masa de las partículas sólidas. Se utiliza un modelo de humectación instantánea para modelar la nucleación heterogénea. Se emplea el modelo de Gyarmathy para la aproximación de la tasa de crecimiento promedio de las gotas. La condensación heterogénea se modela en función de la función de distribución de impurezas extranjeras. Se propone un enfoque para calcular la evaporación en el caso heterogéneo. Se presenta una comparación de los modelos propuestos con un experimento numérico. Se lleva a cabo una simulación numérica de la condensación homogénea-heterogénea en un eyección de dinámica de gases.