Resum: Una nova aproximació euleriana-lagarangiana, en la seva forma d'acoblament en dues vies, per a la simulació de flux de bombolles, aigua-aire és presentada a la tesi. També s'inclouen els efectes de les col.lisions entre bombolles, així com les possibles trencaments o coalescència de bombolles. Aquesta aproximació utilitza el model Continuous Random Walk, CRW, per tenir en compte les fluctuacions de la velocitat. Aquesta aproximació s'emmarca dins d'un model de turbulència k-epsilon per a la fase contínua del líquid. En aquesta tesi s'estudiaran els mètodes per realitzar l'acoblament entre les dues aproximacions, l'efecte de la força lift i de la dispersió turbulenta sobre la distribució de la fracció de buits, així com els models de coalescència i trencament de bombolles que puguin ser emprats en aquest tipus d'aproximació. S'ha partit d'un codi eulerià per simular la part contínua, i sobre ell s'ha acoblat l'aproximació lagrangiana. Perquè aquest acoblament afecti al solver de la fase contínua s'ha de afegir fonts de moment i turbulència. A més s'ha modificat el volum computacional de cada cel.la perquè tingui en consideració el volum ocupat per la fase dispersa. El acoblament en doble via fa que els perfils de velocitat i turbulència de la fase contínua es modifiquin notablement i que s'aproximin als reals, el que resulta bàsic per a la correcta simulació de les forces interfacials. La col.lisió entre bombolles, i bombolles i paret s'ha inclòs. Aquest efecte és necessari com a pas previ a incloure els processos de ruptura o coalescència, encara que la col.lisió en si tingui efectes limitats en la distribució de la fracció de buits. El procés de coalescència es basa en el model de Chester (1991), el model compara el temps de col.lisió amb el temps de drenatge de la pel.lícula entre bombolles per determinar si existeix o no coalescència. El model de ruptura es basa en el model de Martínez-Bazán. Un dels principals fites de la tesi és el desenvolupament realitzat per incloure la font de turbulència induïda en la fase contínua a causa de la presència de bombolles i la seva inclusió el CRW. Això permet evitar la utilització dels models de força de dispersió turbulenta i disposar així d'un model molt més realista. Diferents metodologies han estat provades, i una expressió és proposta que ofereix resultats molt d'acord amb els experimentals. En la tesi es fa ús de resultats experimentals per testejar les solucions aportades en cada cas.