Resumen:
Una nueva aproximación euleriana-lagarangiana, en su forma de acople en dos vías,  para la simulación de flujo de burbujas, agua-aire  es presentada en la tesis, en la que se incluyen los efectos de las colisiones entre burbujas, así como las posibles roturas o coalescencia de burbujas. Esta aproximación utiliza el modelo Continuous Random Walk, CRW, para tener en cuenta las fluctuaciones de la velocidad. Esta aproximación se enmarca dentro de  un modelo de turbulencia k-epsilon para la fase continua del líquido.  En esta tesis se estudiarán los métodos para realizar el acople entre ambas aproximaciones, el efecto de la fuerza lift y de la dispersión turbulenta sobre la distribución de la fracción de huecos, así como los modelos  de coalescencia y rotura de burbujas que puedan ser empleados en este tipo de aproximación. 
Se ha partido de un código euleriano para simular la parte continua, y sobre él se ha acoplado la aproximación lagrangiana. Para que ese acople afecte a la fase continua sobre su solver ser han añadido  fuentes de momento y turbulencia. Además se ha modificado el volumen computacional de cada celda para que tenga en consideración el volumen ocupado por la fase dispersa. El acople en doble vía hace que los perfiles de velocidad y turbulencia de la fase continua se modifiquen notablemente y que se aproximen a los reales, lo que resulta básico para la correcta simulación de las fuerzas interfaciales.
La colisión entre burbujas, y burbujas y pared se ha incluido. Este efecto es necesario como paso previo a incluir los procesos de rotura o coalescencia de burbujas, aunque la colisión en sí tenga efectos limitados en al distribución de la fracción de huecos. El proceso de coalescencia se basa en el modelo de Chester ( 1991 ) , el modelo compara el tiempo de colisión con el tiempo de drenaje de la película entre burbujas para determinar si existe o no coalescencia. El modelo de rotura se basa en el modelo de Martínez-Bazán. 
Uno de los principales hitos de la tesis es el desarrollo realizado para incluir la fuente de turbulencia inducida en la fase continua debido a la presencia de burbujas y su inclusión el CRW. Esto permite evitar la utilización de los modelos de de fuerza d dispersión turbulenta y disponer así de un modelo mucho más realista.  Diferentes metodologías han sido probadas, y una expresión es propuesta que ofrece resultados muy acordes con los experimentales. En la tesis se hace uso de resultados experimentales para testear las soluciones aportadas en cada caso.