🔬 physics

Development of a Quasi-Static MHD Solver for Turbulent Flows with Non-Uniform Applied Magnetic Fields in OpenFOAM

Para abordar los graves errores de cancelación discreta y la divergencia de la solución en campos magnéticos no uniformes que afectan a los resolvedores de MHD estándar de bajo número de Reynolds magnético, este artículo presenta **mhdturbFoamQS**, un nuevo resolvedor de OpenFOAM que elimina estos residuales mediante el cálculo directo de la fuerza de Lorentz a partir de la formulación del potencial eléctrico y la discretización implícita del amortiguamiento magnético, logrando así una alta precisión frente a referencias analíticas sin modificar los cierres de turbulencia estándar.

Autores originales: Wellington Fonseca

Publicado 2026-10-06

📖 5 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Wellington Fonseca

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina un mundo donde campos magnéticos invisibles pueden dirigir el flujo de metal líquido, de forma muy similar a como un conductor guía un río. Este es el reino de la magnetohidrodinámica, un campo de la física que estudia cómo se comportan los fluidos conductores de electricidad, como los metales líquidos o el agua salada, cuando se someten a fuerzas magnéticas. Estas interacciones son vitales para tecnologías que van desde sistemas de enfriamiento en reactores de fusión nuclear hasta métodos para agitar acero fundido sin tocarlo. En muchas de estas aplicaciones, el fluido se mueve a través de un campo magnético que no es perfectamente uniforme; en su lugar, la intensidad y la dirección del campo cambian de un punto a otro, creando un paisaje complejo e irregular que el fluido debe navegar. Para predecir cómo se moverá el fluido en estas condiciones, los científicos dependen de potentes simulaciones por computadora. Sin embargo, un fallo oculto en la forma estándar en que se construyen estas simulaciones las ha vuelto poco fiables cada vez que el campo magnético cambia a través del espacio, dando lugar a resultados que no son solo ligeramente erróneos, sino completamente equivocados.

Un investigador de la Universidad Federal de Pará en Brasil ha identificado y solucionado este problema específico mediante el desarrollo de un nuevo programa informático diseñado para manejar correctamente estos campos magnéticos irregulares. El problema surgió porque el método estándar utilizado para calcular la fuerza magnética sobre el fluido implica restar dos números muy grandes entre sí. Cuando el campo magnético es perfectamente uniforme, estos números se cancelan casi perfectamente, dejando la fuerza pequeña y correcta. Pero cuando el campo varía de una celda a otra en el modelo computacional, la cancelación ya no es exacta. La pequeña cantidad restante, que debería ser cero, actúa como una fuerza fantasma que se fortalece a medida que el campo magnético se vuelve más intenso. En pruebas con un fuerte campo de diez teslas, se descubrió que esta fuerza fantasma era casi cien veces mayor que la fuerza física real que el fluido debería sentir, lo que provocaba que la simulación fallara casi de inmediato.

Para resolver esto, el investigador, Wellington da Silva Fonseca, creó un nuevo solver llamado mhdturbFoamQS. En lugar de intentar restar dos números grandes para hallar la fuerza, este nuevo enfoque calcula primero la corriente eléctrica que fluye a través del fluido y luego utiliza esa corriente para determinar la fuerza directamente. Este método evita la peligrosa resta por completo, asegurando que no se creen fuerzas fantasma, incluso cuando el campo magnético es altamente irregular. El nuevo programa también maneja los límites del contenedor con mayor cuidado, asegurando que la corriente eléctrica no se escape a través de paredes aislantes, un detalle que intentos anteriores solían errar. Además, trata el efecto de amortiguación magnética —la forma en que el campo magnético frena al fluido— dentro del propio cálculo, lo que permite que la simulación se ejecute de forma estable sin necesidad de dar pasos de tiempo imposibles de reducir.

El nuevo solver fue sometido a una serie de pruebas rigurosas para demostrar que funciona. Primero, se comparó con soluciones matemáticas conocidas para un fluido fluyendo a través de un canal plano y una tubería circular. En estas pruebas, los resultados computacionales coincidieron con las respuestas teóricas exactas con un error de menos del medio por ciento, confirmando que la física subyacente se estaba resolviendo correctamente. Una prueba particularmente reveladora involucró un fluido en reposo en un campo magnético que no presentaba corrientes eléctricas fluyendo a través de él. La física dicta que tal fluido no debería moverse en absoluto. El método de simulación estándar causó que el fluido se acelerara a una velocidad significativa debido al error de cálculo, mientras que el nuevo solver mantuvo correctamente el fluido perfectamente quieto, devolviendo una velocidad de cero.

Finalmente, el investigador utilizó la nueva herramienta para simular un flujo turbulento de metal líquido dentro de un conducto rodeado por un anillo de imanes que creaban un campo magnético alterno de seis veces. Este es exactamente el tipo de entorno complejo y no uniforme donde los métodos antiguos fallan. La simulación se ejecutó con éxito, mostrando cómo el campo magnético aplanaba el flujo cerca de las paredes y aceleraba el núcleo del fluido. El estudio también exploró cómo diferentes modelos matemáticos de turbulencia respondían al campo magnético, encontrando que predecían diferentes niveles de aumento de presión y energía, lo que resalta que, si bien la fuerza magnética se calcula ahora correctamente, la forma en que la turbulencia interactúa con ella sigue siendo un área compleja para estudios futuros. Al solucionar el error de cálculo fundamental, este trabajo proporciona una base fiable para que ingenieros y científicos diseñen y analicen sistemas que involucren metales líquidos y campos magnéticos fuertes e irregulares, asegurando que sus modelos computacionales reflejen la realidad en lugar de artefactos numéricos.

Resumen Técnico: Desarrollo de un Solucionador MHD Cuasi-Estático para Flujos Turbulentos con Campos Magnéticos Aplicados No Uniformes en OpenFOAM

Planteamiento del Problema
Los solucionadores de magnetohidrodinámica (MHD) estándar de OpenFOAM suelen utilizar una formulación de inducción magnética donde la fuerza de Lorentz se evalúa mediante el tensor de tensiones de Maxwell. Aunque este enfoque es robusto para campos aplicados uniformes, sufre un severo error de cancelación discreta cuando el campo aplicado ( $B_0$ ) es no uniforme. En el límite continuo, la divergencia del tensor de tensiones de Maxwell y el gradiente de la presión magnética se cancelan exactamente. Sin embargo, en un entorno de volúmenes finitos discretos con campos espacialmente variables, esta cancelación es imperfecta. El residuo actúa como una fuerza de cuerpo espuria que escala con el cuadrado del campo aplicado ( $B_0^2$ ). Para configuraciones de ingeniería realistas (por ejemplo, un campo de 10 T en un conducto de 0.15 m), este residuo puede alcanzar aceleraciones de $\sim 10^3$ m s $^{-2}$ , dos órdenes de magnitud por encima de la aceleración de Lorentz física. Este error causa la divergencia de la solución en pocos pasos de tiempo y no se mitiga mediante el refinamiento de la malla o la reducción del paso de tiempo, lo que indica un fallo fundamental en la formulación en lugar de un problema de resolución numérica.

Metodología
Para abordar esto, el autor desarrolló mhdturbFoamQS, un solucionador de volúmenes finitos basado en la formulación de potencial eléctrico cuasi-estático (bajo número de Reynolds magnético, $R_m \ll 1$ ). La metodología central consiste en:

Cambio de Formulación: En lugar de calcular la fuerza de Lorentz a partir del tensor de tensiones de Maxwell, el solucionador resuelve una ecuación de Poisson para el potencial eléctrico ( $\phi$ ) derivada de la conservación de la carga ( $\nabla \cdot j = 0$ ). La densidad de corriente ( $j$ ) se reconstruye algebraicamente a partir de la ley de Ohm ( $j = \sigma(-\nabla\phi + u \times B_0)$ ), y la fuerza de Lorentz se evalúa directamente como $j \times B_0$ . Esto evita la sustracción de términos de gran magnitud $O(B_0^2)$ .
Manejo de Campos Generalizados: El campo aplicado $B_0$ se trata como un volVectorField general en lugar de una constante, lo que permite variaciones en la magnitud y dirección espacial.
Condiciones de Contorno Consistentes: Se aplica una implementación específica de la condición de pared aislante ( $j \cdot n = 0$ ) a ambos lados de la ecuación de potencial discreta. Esto asegura que el flujo de la fuerza electromotriz motriz se elimine correctamente del lado derecho, evitando la imposición inadvertida de una condición de electrodo.
Amortiguación Magnética Implícita: El término de amortiguación lineal estirado ( $-\sigma|B_0|^2 u$ ) derivado de la descomposición de la fuerza de Lorentz se trata implícitamente dentro de la matriz de momento usando fvm::Sp. Esto elimina la restricción de paso de tiempo asociada al tratamiento explícito de la amortiguación magnética, que de otro modo escalaría como $1/B_0^2$ .
Acoplamiento de la Turbulencia: El solucionador se integra con los cierres de turbulencia estándar de OpenFOAM (Spalart–Allmaras, $k$ – $\epsilon$ , $k$ – $\omega$ SST) sin modificar sus ecuaciones de transporte. Esto aisla la formulación electromagnética de los cambios en los modelos de turbulencia.

Contribuciones Clave

Identificación y Cuantificación: El artículo identifica y cuantifica el error de cancelación discreta en los solucionadores basados en inducción bajo campos no uniformes, demostrando que el error se satura en lugar de decaer con el refinamiento de la malla.
Desarrollo del Solucionador: La implementación de mhdturbFoamQS (344 líneas de C++), que incorpora la ruta del potencial cuasi-estático, el tratamiento consistente de paredes aislantes y el manejo de sumideros magnéticos implícitos.
Marco de Verificación: Un conjunto integral de verificación que incluye:
- Comparación con la solución exacta de Hartmann para canales planos (error en el factor de fricción < 0.17%, error en el perfil de velocidad < 0.5%).
- Comparación con la solución de Shercliff–Gold para tuberías circulares con paredes aislantes.
- Una "prueba nula" utilizando un campo aplicado no uniforme sin corriente donde el fluido está inicialmente en reposo. El solucionador cuasi-estático devuelve correctamente velocidad cero, mientras que el solucionador de inducción genera movimiento espurio.
- Pruebas de solución manufacturada que confirman la precisión de segundo orden tanto para el potencial como para la corriente reconstruida.
Reproducibilidad: El solucionador, los casos de verificación y los scripts de post-procesamiento se liberan bajo la licencia GPL v3 con un DOI de concepto en Zenodo, lo que garantiza que todos los números reportados sean independientemente reproducibles.

Resultados

Verificación: El solucionador reproduce soluciones analíticas para flujos laminares con alta precisión. En la prueba nula, el solucionador de inducción produjo velocidades espurias de hasta 0.36 m s $^{-1}$ (y 3.7 m s $^{-1}$ con valores exactos de pared) en un campo donde el fluido debería permanecer estacionario, mientras que el solucionador cuasi-estático mantuvo un cero de precisión de máquina.
Demostración Turbulenta: El solucionador se aplicó a un flujo turbulento en un conducto ( $Re \approx 10^4$ $R e \approx 1 0^{4}$ ) con un campo magnético alternante no uniforme y realista (estructura azimutal de seis pliegues). Las simulaciones se ejecutaron con tres cierres de turbulencia diferentes a intensidades de campo de hasta 10 T.
- El solucionador se mantuvo estable y no exhibió la divergencia observada en los enfoques basados en inducción.
- Los resultados mostraron un aumento monotónico de la presión promedio volumétrica y de la energía cinética turbulenta con la intensidad del campo.
- Se observó una sensibilidad significativa a la elección del cierre de turbulencia (por ejemplo, a 10 T, el aumento de presión varió desde +17.3% para $k$ – $\epsilon$ hasta +30.1% para Spalart–Allmaras), resaltando que la formulación electromagnética está desacoplada de la incertidumbre del modelado de la turbulencia.

Significancia y Reivindicaciones
El artículo sostiene que la formulación del potencial eléctrico cuasi-estático proporciona un método práctico y robusto para simular flujos MHD turbulentos con campos magnéticos espacialmente variables en OpenFOAM, eliminando eficazmente la patología de cancelación discreta inherente a las formulaciones de tensión de Maxwell. El trabajo enfatiza que evalúa el acoplamiento electromagnético y la formulación numérica en lugar de proponer nueva física de turbulencia; los cierres estándar se utilizan sin modificaciones para aislar los efectos numéricos. El estudio establece que el solucionador puede manejar campos magnéticos no uniformes de escala de ingeniería sin las fuerzas de cuerpo espurias que plagan las implementaciones existentes basadas en inducción, siempre que se cumpla la aproximación de bajo $R_m$ . El autor señala limitaciones, incluyendo la falta de términos magnéticos explícitos en los cierres de turbulencia, la restricción a paredes aislantes y el error de división de primer orden en el esquema de paso de tiempo, definiéndolos como direcciones para trabajos futuros.

Más como este