Skip to content
← Series

Reseñas de artículos

Artículos de CFD y mecánica de fluidos leídos a través de su idea clave, sus ecuaciones y un código de reproducción.

43 partes · 2026.03.24 – 2026.11.01

  1. 01

    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.

    2026.03.24

  2. 02

    Esquema AP-IMEX-DG que libera el paso de tiempo de la velocidad del sonido a cualquier Mach y EOS, reproducido en NumPy.

    2026.05.03

  3. 03

    Cómo el modelo de seis ecuaciones con relajación rígida de presión esquiva las trampas del modelo difuso de cinco ecuaciones.

    2026.05.09

  4. 04

    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.

    2026.05.10

  5. 05

    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.

    2026.05.16

  6. 06

    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.

    2026.05.17

  7. 07

    Reconstrucción apropiada por onda: THINC en contactos, MP5 en choques, central en viscoso

    2026.05.23

  8. 08

    Primero resolver, luego dudar — reproducción del limitador MOOD de Ioriatti, Dumbser, Loubère (2020)

    2026.05.24

  9. 09

    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.

    2026.05.30

  10. 10

    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.

    2026.05.31

  11. 11

    Un precondicionador de complemento de Schur velocidad–presión mantiene plano el número de iteraciones cuando el Mach tiende a cero

    2026.06.06

  12. 12

    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.

    2026.06.07

  13. 13

    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.

    2026.06.13

  14. 14

    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.

    2026.06.14

  15. 15

    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.

    2026.06.20

  16. 16

    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/ε.

    2026.06.21

  17. 17

    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.

    2026.07.12

  18. 18

    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.

    2026.07.18

  19. 19

    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.

    2026.07.19

  20. 20

    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.

    2026.07.25

  21. 21

    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.

    2026.08.01

  22. 22

    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.

    2026.08.02

  23. 23

    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.

    2026.08.09

  24. 24

    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.

    2026.08.16

  25. 25

    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.

    2026.08.23

  26. 26

    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.

    2026.08.29

  27. 27

    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.

    2026.08.30

  28. 28

    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.

    2026.09.05

  29. 29

    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.

    2026.09.06

  30. 30

    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.

    2026.09.19

  31. 31

    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.

    2026.09.20

  32. 32

    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.

    2026.09.26

  33. 33

    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.

    2026.09.27

  34. 34

    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.

    2026.10.03

  35. 35

    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.

    2026.10.04

  36. 36

    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.

    2026.10.08

  37. 37

    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.

    2026.10.10

  38. 38

    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.

    2026.10.17

  39. 39

    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.

    2026.10.18

  40. 40

    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.

    2026.10.24

  41. 41

    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.

    2026.10.25

  42. 42

    Reproducción de la condición de no perturbación del modelo Baer–Nunziato basado en presión de Re y Abgrall (2022): si fracción volumétrica y masa no comparten flujo, la presión salta cientos de MPa.

    2026.10.31

  43. 43

    En flujos bifluidos a bajo Mach, la CFL acústica difumina la interfaz y multiplica los pasos. Reproducción 1D de cómo el IMEX de estado de referencia conserva una CFL material y una burbuja nítida.

    2026.11.01