Este trabajo está dedicado a la investigación computacional de la deformación y ruptura de cuerpos de agua cilíndricos en el flujo de aire a alta velocidad detrás de ondas de choque incidentes. Se simularon configuraciones de columna simple y de columna en tándem en varias disposiciones reproduciendo el proceso de interacción choque/gota en un dispositivo de tubo de choque. Los cálculos se realizaron utilizando un solucionador de terceros recientemente desarrollado para flujos bifásicos compresibles en el marco de la caja de herramientas de volumen finito de código abierto OpenFOAM. El enfoque numérico se basa en el uso del método de volumen de fluidos para resolver la interfaz de fase, donde una técnica de discretización particular nos permite prevenir inestabilidades no físicas. El esquema numérico utiliza información más precisa de las velocidades de propagación locales para mantener una alta resolución y una baja viscosidad numérica. Comparaciones cualitativas y cuantitativas de los resultados con datos experimentales y numéricos de referencia demostraron un buen acuerdo para las principales características del proceso de interacción en términos de morfología, dinámica y ruptura de los cuerpos de agua en deformación.