“La turbulencia es el problema sin resolver más importante de la física clásica.”
Notas de campo de un investigador en CFD multifase compresible
START HERE
Elige tu primera ruta
Tres entradas claras para una primera visita: fundamentos, flujo compresible o un caso práctico.
FDM, FEM y FVM comparados
Descubre qué aproxima cada discretización y qué está diseñada para conservar.
Empezar esta rutaProblemas de Riemann y flujos
Desde el origen del cálculo de choques hasta la lógica de los solvers tipo Godunov.
Empezar esta rutaVerifícalo en OpenFOAM
Configura un escalón hacia atrás y analiza la región de recirculación resultante.
Empezar esta rutaMás leídos esta semana
…
La velocidad salió un sexto con τ=0.6: omitir la media fuerza F/2 en el forcing de Guo
Si se conserva la fuente de Guo pero se quita F/2 de la velocidad LBM, la fuerza efectiva se reduce por (1−1/2τ). En flujo de Kolmogorov con τ=0.6 la amplitud es 0.166.
Una pequeña fluctuación de la fracción de mezcla baja 964 K la temperatura media — β-PDF presunta y el infinito en Z=0
En Z̄ estequiométrica, la temperatura media con β-PDF queda 964 K por debajo de T(Z̄). Con a<1, una malla uniforme de 100 celdas además se pasa por 115 K. El porqué, con código.
Una interfaz de dos celdas pierde el 39% de su tensión superficial — la anisotropía del esfuerzo de Korteweg
En una interfaz plana líquido–vapor la presión normal es constante y solo la tangencial se hunde. Ese hundimiento es la tensión superficial; una red dz=1 pierde hasta el 39%.
[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.
Solo cambió τ y cambió la densidad de un gas en reposo — el esfuerzo (τ−½)² del forcing de Shan–Chen por desplazamiento de velocidad
Con la misma fuerza sobre una columna de gas D1Q3 en reposo, solo el forcing de Shan–Chen produce una densidad que depende de τ. La velocidad de equilibrio desplazada deja un esfuerzo (τ−½)²F²/ρ.
Dividir la constante inversa entre Kp disoció el 99.6% del AlO a 2500 K — el (p°/RT)^Δn oculto en la constante de equilibrio
Usar Kp en kr = kf/K desplaza el equilibrio de toda reacción con Δn≠0 por un factor (p°/RT)^Δn. Un reactor de disociación de AlO reproduce el error en SI y en CGS.
Al refinar la malla al doble, el error se triplicó — lo que oculta el criterio de parada por incremento en Jacobi
Si Jacobi se detiene cuando el incremento baja de una tolerancia, el error restante crece como 1/(1-ρ) al refinar. Una placa 2D compara Jacobi, SOR y ω óptimo.
Sumar un término en v y otro en v² acertó en los extremos y falló un 42 % en el medio — la ley de arrastre de dos términos de Prony y Re = 60
Bossut midió un arrastre proporcional a v² y Du Buat uno proporcional a v. La suma de Prony, a·v²+b·v, cambia de régimen en Re = 60 y allí sobrestima la velocidad terminal hasta un 42 %.
Darle su propio τ a la conducción dejó el calentamiento viscoso atrapado en Pr = 1 — el cuasi-equilibrio g* en LB de dos poblaciones
En Couette, un LB térmico de dos τ deja θ_max en 1.000 para cualquier Pr. Con el cuasi-equilibrio g* salen 0.498, 0.710, 1.000 y 1.997.
[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.