[Revisión de artículo] MSTACS con pendiente 3 divergió en el paso 255 con Co=0.5 — esquemas compresivos y el límite de Courant
Se probó el esquema compresivo de la ec. (18) de MSTACS en una franja 1D. Sin redistribución diverge por encima de Co=0.44; la causa es el producto de la pendiente k y el número de Courant.
[Revisión de artículo] Con β=2 la interfaz se quedó en 14 celdas y una senoide quedó peor que con upwind de primer orden — el limitador compresivo de Tallois
Con un limitador compresivo β=2 en la fracción volumétrica, la interfaz de una placa sigue midiendo 14 celdas tras cinco periodos. El mismo limitador deja en una senoide un error L1 mayor que upwind de primer orden.
[Revisión de artículo] Mitad agua, mitad aire y el sonido se arrastra a 24 m/s — El término K∇·u del modelo de cinco ecuaciones de Murrone–Guillard
Con 50% de aire, la mezcla transmite el sonido a 23,6 m/s, más lento que cada fase. Sin el término K∇·u, la misma prueba 1D da 745 m/s.
[Revisión de artículo] El tetraedro que adelantó una esfera no se dispersó — un índice de achatamiento para cúmulos de fragmentos a Mach 20
Cuatro esferas a Mach 20 se separan según la postura del tetraedro. Los índices B y A en 40 000 orientaciones muestran qué número fija la separación.
El agua con sonido a 1625 m/s apareció a 18 m/s en el subsistema de presión — adónde lleva la rigidez el splitting de Toro–Vázquez
A 1 bar, el subsistema de presión del agua da A = 18.44 m/s y c = 1624.9 m/s. A/c queda en 0.011, y la rigidez se mueve entre los dos pasos IMEX.
Reducir 574 veces los cálculos de distancia costó 104 contactos — un solo grano fija el sector en DEM
La poda por sectores en la búsqueda de vecinos de DEM solo es segura si shell·s supera d_max, medido sobre 2400 granos en aglomeración.
Multipliqué los núcleos por 64 y la eficiencia cayó a 16.3% — la trampa de medir weak scaling refinando la malla
Por qué las iteraciones de Krylov crecieron 6.1 veces en un estudio de weak scaling que refina la malla, y cómo un nivel más de multigrid devuelve la curva a plana.
Al quitar K el agua se comprimió 1.86 veces más — un error de interfaz que refinar la malla no elimina
Derivación del factor de compresibilidad K del VOF compresible a partir del cierre isobárico, y medición de cuánto volumen de líquido se pierde al anularlo o al invertir su signo en una compresión 1D.
[Reseña de artículo] La presión pasó a implícita y Δt se redujo 100 veces — el segundo 1/M² oculto en la velocidad relativa
En un splitting all-Mach bifásico, el salto de entalpía esconde un segundo 1/M². El CFL sigue atado al Mach hasta integrar ese término en implícito.
[Reseña de artículo] 0.333 en tracción, 0 en compresión — dónde las partículas SPH se emparejan
El mismo SPH aguanta la compresión y empareja en tracción. Con h/dx = 1.2 la tasa vale 0.333; congelar el núcleo en la malla inicial la lleva a 0.
[Reseña de artículo] La función de color del mismo punto difirió en 0.544 según el camino — restricciones de rotacional y malla escalonada
Tras 400 pasos la malla escalonada deja el rotacional en 1e-12. El Rusanov colocado, con el mismo cizallamiento, abre una brecha de camino de 0.544.
La factura de un paso de tiempo 1001 veces mayor — por qué un esquema semiimplícito de presión no es L∞ estable
Un esquema semiimplícito de presión borra la velocidad del sonido del CFL y gana un paso 1001 veces mayor en M = 0.001. La oscilación no decae.
La presión saltó un 21% en un problema donde no debía pasar nada — el criterio de Abgrall y el término no conservativo
Que un flujo de presión y velocidad uniformes deje de serlo no es culpa del flujo numérico, sino de cómo se discretiza el término no conservativo.
[Reseña de artículo] Los dos fragmentos no se repelieron: volaron pegados — el ángulo de equilibrio de 144° de un par de esferas a Mach 20
En Mach 20 dos esferas en contacto no se repelen: se alinean en un ángulo y generan sustentación juntas. Ese ángulo sale de una aproximación newtoniana.
Subir Δt cinco veces dio 1.9x de velocidad; diez veces no dio nada — la ventana que abre la tensión superficial implícita
Con razón de densidades 1000, un paso capilar 5 veces mayor corrió 1.9 veces más rápido. Diez veces no aceleró. Janodet y otros (2025) midieron el techo.
Con CFL acústico 10 el precondicionador multiplicaba el error por 100 en cada sweep — la flecha que descarta un precondicionador por bloques
El Gauss-Seidel por bloques descarta una realimentación presión-velocidad cuya ganancia es exactamente el CFL acústico al cuadrado. Un modelo 1D lo mide.
Refinar la malla a la mitad adelantó el estallido al doble — autovalores complejos del modelo two-fluid de 6 ecuaciones
Con deslizamiento el two-fluid de seis ecuaciones se vuelve complejo. El coeficiente de presión interfacial debe pasar de 1 para volver al eje real.
[Reseña] Se soltó el grillete capilar y quedó CFL 0.05 — el techo real del transporte de interfaz en VOF
Tratar la tensión superficial de forma implícita rompe el límite capilar de paso temporal, pero el esquema compresivo que transporta la función de color lo frena en CFL 0.05.
[Reseña] Separar del todo advección y presión — la división TV del flux de Baer–Nunziato
Partir el flux de Baer–Nunziato deja el sonido solo en el sistema de presión. Toro Test 4 ganó a Rusanov por 1,5×; Noh llevó la presión a −0,12.
[Reseña de artículo] Cómo eliminar la subrelajación — linealización en un solver acoplado basado en presión
En el caso Sod, M/Co = 2.25 lleva la iteración con coeficiente fijo a P = 2.250, mientras Newton llega a la raíz física P = 1 en cinco iteraciones.
[Reseña de artículo] Una gota devuelve casi toda la onda sonora — impedancia acústica y un solver bifásico de interfaz sharp
La reflexión aire-agua devuelve el 99.889% de la energía acústica. Con solo 1% de gas, la velocidad de Wood cae de 1500 m/s a 120.5 m/s.
Devolver la nitidez a la interfaz en cada paso — Implementando el afilado anti-difusión (IST)
Tras 200 pasos, el borde de la burbuja de helio pasa de 6 a 20 celdas. El posprocesado IST vuelve a fijar ese grosor con el cociente ε/a.
Por qué una gota en reposo fluye sola — corrientes parásitas y tensión superficial well-balanced
Un error de curvatura del 3% lleva el residuo a 4.8e5, mientras la curvatura exacta lo deja en 4.1e-10. Esa brecha es la corriente parásita.
Menos ecuaciones, y dejó de reventar — el modelo de cuatro ecuaciones y la condición de equilibrio de interfaz
Transportar Γ de forma conservativa agita una presión uniforme. El modelo de cuatro ecuaciones lo mueve sin conservación y el error cae al redondeo.
[Reseña de artículo] Cuando una microburbuja arde como una estrella bajo el agua — Keller–Miksis y la fuerza secundaria de Bjerknes
Una burbuja de 5 μm tiene una frecuencia de Minnaert cerca de 650 kHz. A 0,1 R_E y n=1,4 la presión politrópica del gas ronda las 20 000 atm.
[Reseña] Cómo no matar el sonido a bajo Mach — esquema IMEX-AP de ondas
Con ε=0.01 en 64 celdas la amplitud upwind cae a 1.980e-7 tras 12800 pasos; la columna AP implícita queda en 1.00. Los pasos escalan como 1/ε.
[Reseña] Segundo orden en una sola etapa — un solver GRP para el modelo multifásico de Kapila
En N=320, un GRP de una sola etapa da orden 2.03 y un error L1 de 1.55e-5, frente a 2.31e-2 y orden 0.97 de Godunov en la misma onda.
Por qué un grito sobre el agua nunca llega abajo — Acústica a través de una interfaz difusa
Un salto de impedancia aire–agua cercano a 3600 refleja cerca del 99,9% del sonido. El ángulo crítico de 13,4° se comprueba en un tubo 1D.
[Reseña] Dibujar choques con tanh — Implementando la reconstrucción TENO-THINC
Con beta = 1.8, THINC-BVD mantiene un pico de 1.00 y un ancho de 30 celdas tras cinco vueltas; minmod TVD cae a 0.84 y 23 celdas.
[Reseña de artículo] Cómo una DQN domó los vórtices de un perfil con un chorro sintético — Control activo de flujo por aprendizaje por refuerzo
En 300 episodios la recompensa de juguete solo pasa de 3.12 a 3.25, y el chorro aprendido se detiene en 12 m/s, donde empieza la penalización.
[Reseña de artículo] GMRES corrió 1000 veces y no se detenía — precondicionamiento por bloques físico para un melt pool a toda velocidad
Un precondicionador de complemento de Schur velocidad–presión mantiene plano el número de iteraciones cuando el Mach tiende a cero
[Revisión de paper] Flujo bifásico sin Poisson — Acoplando la ecuación general de presión con MSTACS
Subir la velocidad del sonido artificial de 5 a 15 reduce el paso estable a un tercio. El error de volumen E_G de OS-MSTACS queda en 6,57e-8.
[Reseña] Que la velocidad del sonido no se coma el CFL — Peluchon (2017) y la división IMEX acústica-transporte
Sonido a 1500 m/s y flujo a 5 m/s: una celda de 1 mm exige un paso explícito de 0,67 μs. El split acústico deja el CFL exterior en la velocidad.
[Reseña] Sospechar de la celda después de resolverla — limitador MOOD subcell para DG
Primero resolver, luego dudar — reproducción del limitador MOOD de Ioriatti, Dumbser, Loubère (2020)
[Revisión de paper] Un solo esquema no basta para todas las ondas — Chamarthi (2025)
Reconstrucción apropiada por onda: THINC en contactos, MP5 en choques, central en viscoso
[Reseña] Por qué Roe se rompe a Mach bajo — Hope–Collins (2022) y las tres escalas de difusión artificial
A Mach 0.01 la escala convectiva deja una amplitud de 0.018, mientras la acústica conserva 0.974 y la mixta 0.971 tras un período.
[Revisión de artículo] Partir Newton en dos y la superficie libre se vuelve resoluble — Casulli–Zanolli (2012)
El día 1 de un pozo de 41 celdas muestra 39 celdas presurizadas y cero secas. Al día 10 quedan 0 y 12, y el Newton exterior baja de 5 a 1.
[Reseña] Cuando 30 K hacen mentir a Boussinesq — Demou (2018) NOB y la división de presión
En agua a ΔT=30 K el máximo impreso de |1-rho0/rho| es 0.0117, y el Poisson de presión sigue aún con coeficiente constante.
[Reseña de paper] Resolver con dos presiones, fundir en una sola — Saurel (2009) sobre relajación multifásica
Cómo el modelo de seis ecuaciones con relajación rígida de presión esquiva las trampas del modelo difuso de cinco ecuaciones.
[Reseña] Cuando la velocidad del sonido diverge y el paso de tiempo sobrevive — AP-IMEX para gases no ideales de Orlando–Bonaventura (2024)
Esquema AP-IMEX-DG que libera el paso de tiempo de la velocidad del sonido a cualquier Mach y EOS, reproducido en NumPy.
[Reseña] Elegir la reconstrucción con menor salto — Deng (2018) MUSCL-THINC-BVD
Deng (2018) construye MUSCL y THINC en cada celda y conserva la candidata con menor salto en las dos caras. Esa regla es BVD.
[Reseña de paper] Terminar la simulación antes de la tormenta — Casulli–Walters (2000), aguas someras semi-implícitas en mallas ortogonales no estructuradas
A 28 m de profundidad el paso explícito aguanta 0,2 s. Dejar implícito solo el gradiente de superficie libre es lo que sostiene un paso de cinco minutos.
[Reseña] Cuando WENO difumina los pequeños vórtices — Fu (2019) esquema TENO de baja disipación en volúmenes finitos
Un 10 por ciento más en un indicador de suavidad ya mueve los pesos WENO-JS. El corte 0/1 de TENO no cambia y el estencil de quinto orden sigue.
[Reseña de artículo] Physics-Informed Neural Networks para modelado de turbulencia
Sensores de presión reconstruyen un cilindro a Re=100 con error de velocidad inferior al 1%. El número se mueve si el peso físico ahoga el residuo viscoso.