“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
Más leídos esta semana
…
Subir Δt cinco veces dio 1.9x de velocidad; diez veces no dio nada — la ventana que abre la tensión superficial implícita
Romper la restricción capilar de paso temporal deja otros techos en pie. Janodet et al. (2025) midieron el límite de estabilidad y el óptimo de rendimiento con razón de densidades 1000.
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.
Subir el número de Péclet de celda a 2 llevó la solución hasta -0.33 — donde desapareció la energía a minimizar
La advección rompe la simetría de la matriz de rigidez y con ella se va el principio de mínima energía potencial. Esta es la causa real de las oscilaciones de Galerkin y lo que compra la estabilización a cambio.
La fuerza de Coriolis hizo cero trabajo y el impulsor sumó 1800 J/kg — energía en un marco rotante
Lo que se conserva en un marco rotante es la rotalpía, no la entalpía total. Aquí se verifica numéricamente por qué Coriolis no hace trabajo y de dónde sale realmente la altura de una bomba.
Al quitar el −1/2 la viscosidad salió seis veces más grande — el Δt/2 que deja la discretización del LBM
Una integración trapezoidal y un cambio de variable dejan un Δt/2 en tres lugares distintos: la viscosidad, la velocidad y el esfuerzo. Una EDO escalar y una malla D2Q9 miden los tres.
La cuerda doblada aguantó, la cuerda cortada se pasó un 9% — el final numérico de la disputa de 1747
Las ondas viajeras de d'Alembert y la serie de senos de Bernoulli son la misma solución hasta que se trunca. El código muestra de dónde sale el sobrepico del 9%.
Tres en la pared, cinco en la esquina — contar las poblaciones que pierde un nodo de frontera en LBM
Tras el streaming, el número de poblaciones vacías en un nodo de lattice Boltzmann es igual al número de celdas sólidas a su alrededor: tres en una pared plana, cinco en una esquina cóncava, una en una esquina convexa.
Refinar la malla a la mitad adelantó el estallido al doble — autovalores complejos del modelo two-fluid de 6 ecuaciones
El modelo two-fluid de seis ecuaciones con presión única se vuelve complejo en cuanto las fases deslizan una respecto a la otra. Subir el coeficiente de presión interfacial por encima de 1 devuelve los autovalores al eje real.
Metí la presión tal cual y una interfaz en reposo empezó a temblar — dos formas del forcing en LBM no ideal
Al resolver un fluido de van der Waals con lattice Boltzmann, poner la presión p en el término de fuerza o el potencial químico mu coincide en el continuo y se separa en la red.
Bajar el humo 1 m costó 7 Pa de tiro — la estufa de Franklin y el presupuesto de flotabilidad
La boca de la chimenea no se mueve y aun así un tramo descendente rompe el tiro. Aquí se escribe el libro contable de la flotabilidad y se calcula cuándo el humo entra en la sala.
Un punto de cuadratura menos y la solución explotó — el piso de integración de DG y la base de Taylor
El número mínimo de puntos de Gauss en la integral de celda de DG lo fija el grado del integrando, no la malla. Pasar de tres puntos a dos no cambia nada; bajar a uno destruye el esquema.
La onda de choque nunca salió de la línea de partida — dónde se separan la forma conservativa y la primitiva
Escribir la misma ecuación como divergencia de flujo o como velocidad por gradiente cambia la velocidad de choque calculada. Refinar la malla ocho veces converge al valor equivocado.