Un nuevo método de identificación de superficie libre para el método MPS 3D

Por qué importa rastrear las superficies de agua en movimiento

Desde olas que chocan contra muros de contención hasta el vaivén de combustible dentro de tanques espaciales, los ingenieros confían cada vez más en simulaciones por ordenador para predecir cómo se mueven los líquidos y qué fuerzas ejercen. Muchas de estas simulaciones usan un enjambre de partículas virtuales en lugar de una malla rígida para representar el agua. Pero hay un problema sutil: el ordenador debe decidir constantemente qué partículas están en la superficie exterior del líquido, donde actúa la presión del aire circundante. Si esa decisión es incluso ligeramente incorrecta, las presiones predichas sobre paredes y estructuras pueden volverse ruidosas o engañosas. Este estudio presenta una nueva forma de identificar esas partículas de superficie en simulaciones tridimensionales, haciendo que el agua virtual se comporte más como la real.

Una forma diferente de trazar la línea entre agua y aire

En el método Moving Particle Semi-Implicit (MPS), el agua se modela como innumerables partículas diminutas que se mueven e interactúan. Los enfoques tradicionales deciden si una partícula está en la superficie libre simplemente contando cuántos vecinos tiene. Menos vecinos suele significar que está en o cerca de la superficie. Sin embargo, esta regla empírica puede fallar cuando las partículas están desigualmente espaciadas, como ocurre con frecuencia en flujos turbulentos o cerca de esquinas afiladas. Los autores, en cambio, se centran en cómo se disponen espacialmente las partículas vecinas, usando una cantidad llamada divergencia de posición relativa, que mide cuánto "se extienden" los vecinos alrededor de una partícula dada. Si los vecinos tienden a situarse mayormente en un lado, el valor de la divergencia cambia de una manera que indica una probable superficie.

Añadir dirección a la prueba de superficie

Aunque la medida de divergencia mejorada coincide muy bien con las expectativas teóricas, los autores muestran que en tres dimensiones un umbral simple sobre este valor no puede separar perfectamente partículas interiores y de superficie, especialmente cerca de aristas pronunciadas y dentro de cavidades. Para solucionarlo, introducen un paso adicional que estima la dirección local de la superficie, o vector normal, a partir de las partículas circundantes. Usando esta dirección, el método define una región en forma de cono que apunta hacia el exterior desde cada partícula candidata. Si ese cono no contiene otras partículas dentro de cierta distancia, se considera que la partícula yace en la superficie libre. Este enfoque combinado, llamado RPD+NV, emplea tanto cómo se distribuyen los vecinos como dónde hay espacio vacío, ofreciendo una imagen más fiable de la verdadera interfaz entre líquido y aire.

Probar el método en formas complejas y agua en calma

Los investigadores prueban primero su método en formas geométricas estáticas y puras: un cubo rematado por una cúpula, un cubo con un hueco en forma de cuenco y un tubo en forma de S que contiene una cavidad interna. Estos casos desafían al algoritmo con superficies tanto convexas como cóncavas. Al comparar las partículas de superficie detectadas con las superficies verdaderas conocidas, muestran que el nuevo método captura correctamente regiones convexas y cóncavas y resuelve con precisión cavidades ocultas. Luego aplican el enfoque a un contenedor hidrostático simple y verifican que, con elecciones razonables del espaciamiento de partículas y del paso temporal, la presión predicha en el fondo coincide muy de cerca con el valor teórico. En estas condiciones de calma, el nuevo detector de superficie produce menos oscilaciones de presión que los métodos antiguos, manteniendo casi el mismo coste computacional.

Capturar olas, impactos y vaivenes

A continuación, el equipo somete el método a situaciones altamente dinámicas. En una simulación tridimensional de rotura de presa, una columna de agua colapsa y golpea una pared, creando superficies libres que cambian rápidamente, salpicaduras y bolsas de aire. En comparación con criterios tradicionales basados solo en el recuento de vecinos o en la divergencia, el método RPD+NV rastrea superficies más continuas y detalladas y produce campos de presión más suaves y realistas. De forma crucial, la presión máxima de impacto en la pared se aproxima mucho más a las mediciones de laboratorio. En un caso de vaivén amortiguado, donde las ondas dentro de un tanque se van asentando gradualmente, el nuevo método conduce a la evolución más suave de las fuerzas sobre las paredes, lo que indica una reducción del ruido numérico. Finalmente, en una prueba que imita el vaivén de combustible dentro de un tanque anular que se balancea, las formas de la superficie y los picos de presión de la simulación concuerdan bien con imágenes de cámara de alta velocidad y datos de sensores de presión, nuevamente sin picos espurios de presión en la superficie.

Qué significa esto para las simulaciones del mundo real

Para lectores fuera de la dinámica de fluidos computacional, el mensaje práctico es que este trabajo mejora la forma en que los ordenadores deciden dónde termina el agua y comienza el aire en modelos basados en partículas. Al combinar una medida refinada de cómo rodean los vecinos a cada partícula con una prueba direccional de "cono vacío", el nuevo método RPD+NV reduce drásticamente las clasificaciones erróneas de partículas de superficie. El resultado son predicciones de presión más limpias y estables en una variedad de flujos, desde agua en reposo hasta impactos violentos y complejos vaivenes tridimensionales. Dado que el método añade poco coste computacional, ofrece una mejora práctica para ingenieros que simulan estructuras costeras, el movimiento de buques, tanques de almacenamiento y otros sistemas donde capturar con precisión superficies de agua en movimiento es esencial.

Cita: Geng, C., Wang, Wh., Heng, My. et al. A new free surface identification method for 3D MPS method. Sci Rep 16, 13829 (2026). https://doi.org/10.1038/s41598-026-44218-9

Palabras clave: simulación de fluidos basada en partículas, detección de superficie libre, impacto numérico de olas, vaivenes en tanques, método de partículas en movimiento