Hydrodynamic simulations of expanded warm dense foil heated by pulsed-power

Este artículo presenta un marco de modelado robusto que acopla simulaciones de circuitos eléctricos de potencia pulsada con códigos hidrodinámicos unidimensionales para diseñar y optimizar con precisión experimentos para la generación de materia densa caliente expandida en láminas metálicas delgadas confinadas dentro de celdas de zafiro.

Lo esencial

La visión general: ¿Qué es la "Materia Cálida Densa"?

Imagina un material que no es del todo un sólido, ni es del todo un líquido, ni es del todo un gas. Se encuentra en un punto intermedio extraño y desordenado llamado Materia Cálida Densa (WDM, por sus siglas en inglés). Piensa en ello como una pista de baile abarrotada donde todos se mueven rápido (calor) pero aún así chocan entre sí (densidad).

Los científicos necesitan comprender este estado de la materia para estudiar cosas como el interior de planetas gigantes o para construir mejores reactores de energía de fusión. Pero es difícil de estudiar porque es complicado de crear en un laboratorio y aún más difícil de predecir con matemáticas.

El experimento: El "Panqueque Eléctrico"

Los investigadores configuraron un experimento para crear este estado.

• La configuración: Tomaron una lámina de metal muy fina (como un papel de aluminio microscópico) y la colocaron entre dos placas gruesas y duras de zafiro (como el cristal de un reloj).
• La acción: Electrocutaron este sándwich de metal con un estallido de electricidad masivo y súper rápido (potencia pulsada).
• El resultado: La electricidad calienta el metal tan rápido (en menos de una millonésima de segundo) que se derrite, hierve y se convierte en un plasma caliente que se expande. Debido a que las placas de zafiro lo sujetan, el metal solo puede expandirse en una dirección, como un panqueque inflándose.

El problema: La "Caja Negra"

El desafío es que cuando se le da el choque eléctrico al metal, dos cosas suceden al mismo tiempo:

1. El circuito eléctrico: La electricidad fluye a través de los cables, el interruptor y el metal. A medida que el metal se calienta y cambia de forma, su capacidad para conducir electricidad cambia, lo que altera el flujo de la corriente.
2. El movimiento físico: El metal se calienta, se expande y se mueve. A medida que se mueve, cambia la forma del circuito, lo que cambia la electricidad nuevamente.

Es un bucle de retroalimentación. Si intentas calcular la electricidad sin saber cómo se mueve el metal, te equivocas. Si intentas calcular el movimiento del metal sin saber la electricidad, también te equivocas.

La solución: Un modelo de "Bicicleta Tándem"

Los autores construyeron un programa de computadora que actúa como una bicicleta tándem.

• Ciclista 1 (El modelo eléctrico): Esta parte simula el generador de energía, el interruptor y los cables. Calcula cuánta corriente está fluyendo.
• Ciclista 2 (El modelo hidrodinámico): Esta parte simula la lámina de metal. Calcula cómo el metal se calienta, se expande y cambia su densidad.

Estos dos ciclistas están unidos. Cada pequeña fracción de segundo, se comunican entre sí:

• "Oye, el metal acaba de calentarse y hacerse más delgado", dice el Ciclista 2.
• "Está bien, ajustaré el flujo de corriente porque el metal ahora es un peor conductor", dice el Ciclista 1.
• "De acuerdo, actualizaré el calor y la presión basándome en esa nueva corriente", dice el Ciclista 2.

Cómo lo probaron

Para asegurarse de que su "bicicleta tándem" funciona, probaron de tres maneras diferentes, como si revisaran el motor de un coche en distintos niveles:

1. La prueba de "Potencia Conocida": Le dieron a la computadora la medición real de la electricidad de un experimento real y preguntaron: "¿Pueden predecir cómo se mueve el metal?".

   • Resultado: Sí, muy bien. La computadora predijo la velocidad y la expansión del metal casi perfectamente. Esto les indicó qué reglas matemáticas ("Ecuaciones de Estado") describen mejor el comportamiento del metal.

2. La prueba de "Conductividad Conocida": Le dieron a la computadora la conductividad eléctrica real del metal (qué tan bien conduce) y preguntaron: "¿Pueden predecir la electricidad y el movimiento?".

   • Resultado: Sí. La computadora predijo con éxito el voltaje y la corriente, coincidiendo con el experimento real. Esto demostró que las dos partes del modelo se comunican correctamente.

3. La prueba de "Predicción Pura": Esta fue la más difícil. No le dieron a la computadora ningún dato del experimento real. Solo le dieron las leyes de la física y preguntaron: "¿Pueden predecir todo el experimento desde cero?".

   • Resultado: Fue muy cercano. La computadora predijo la velocidad, la corriente y el voltaje con buena precisión. Hubo pequeñas diferencias (como un error del 10% en el voltaje al final), pero la imagen general era correcta.

Por qué esto es importante

El artículo concluye que este modelo de computadora es una herramienta robusta y eficiente.

En lugar de simplemente adivinar cómo configurar futuros experimentos, los científicos ahora pueden usar este modelo de "bicicleta tándem" para diseñarlos. Pueden simular diferentes escenarios en la computadora para ver qué sucederá antes de encender la máquina real. Ayuda a comprender la física de la materia cálida densa sin tener que depender únicamente de experimentos costosos y difíciles.

En resumen: Construyeron un gemelo digital de una explosión eléctrica de alta velocidad. Demostraron que funciona comparándolo con explosiones reales, y ahora pueden usarlo para planificar futuros experimentos con confianza.

Resumen técnico

Resumen Técnico: Simulaciones Hidrodinámicas de Láminas Expandidas Calentadas por Potencia de Pulso

Planteamiento del Problema
La Materia de Densidad Cálida (WDM, por sus siglas en inglés), el régimen intermedio que une la materia condensada y el plasma, es crítica para campos que van desde las ciencias planetarias hasta la fusión por confinamiento inercial (ICF). Sin embargo, modelar este régimen es un desafío debido a la coexistencia de efectos de correlación iónica y degeneración electrónica, que ni la física de plasmas clásica ni las teorías convencentes de la materia condensada describen con precisión. Aunque existen diversas técnicas experimentales para producir WDM, los datos relativos al régimen expandido (densidad inferior a la nominal) son particularmente valiosos para validar las ecuaciones de estado (EOS) y los modelos de conductividad relevantes para las configuraciones de hohlraum de ICF.

El desafío específico abordado en este trabajo es el diseño y la simulación de un experimento de potencia de pulso donde láminas metálicas delgadas, confinadas dentro de una celda de zafiro, son calentadas por efecto Joule para alcanzar el régimen de WDM expandido. Diseñar estos experimentos es difícil porque requiere predecir simultáneamente la respuesta eléctrica del conductor de potencia de pulso y la evolución hidrodinámica del material calentado. Además, la interpretación precisa de la transición de sólido a plasma parcialmente ionizado requiere un marco numérico totalmente acoplado que vincule la deposición de energía eléctrica con la evolución termodinámica.

Metodología
Los autores desarrollaron un marco de modelado que acopla una descripción eléctrica de un sistema de potencia de pulso con un código hidrodinámico unidimensional (1D), ESTHER. El enfoque se basa en una geometría planar 1D, lo que simplifica el problema al asumir que la expansión ocurre en una única dirección preferencial debido al confinamiento por alta impedancia acústica (zafiro), permitiendo así despreciar los efectos de los campos magnéticos y reducir el costo computacional manteniendo la física esencial.

La metodología procede a través de tres etapas de modelado distintas:

1. Modelado de Circuito Eléctrico: El conductor de potencia de pulso, el interruptor de chispa (spark-gap) y la carga se modelan como un circuito RLC equivalente. La evolución de la corriente se resuelve numéricamente utilizando un esquema de diferencias finitas centradas de segundo orden para manejar la resistencia e inductancia dependientes del tiempo.

   • Conductor: Caracterizado por parámetros intrínsecos fijos ($R_f, L_f, C$).
   • Interruptor: Modelado mediante un modelo de arco de Braginskii modificado, donde la resistencia e inductancia de la chispa evolucionan según la integral de la corriente y el número de canales paralelos ($N$).
   • Validación: Este modelo eléctrico fue validado mediante mediciones de cortocircuito de cuatro generadores diferentes (EPP1, EPP2, EOLE, GEPI2) que cubren corrientes desde ~140 kA hasta ~5 MA.

2. Simulación Hidrodinámica: Se utilizó el código ESTHER, un resolvedor hidrodinámico Lagrangiano 1D, para simular la evolución termodinámica de la lámina. El código resuelve las ecuaciones de conservación de masa, momento y energía, incorporando una EOS y conducción térmica. El término fuente de calentamiento por efecto Joule se integra directamente en la ecuación de energía.

3. Estrategias de Acoplamiento: Se probaron tres enfoques para vincular los dominios eléctrico e hidrodinámico:

   • Enfoque A (Potencia Experimental): El término fuente es la potencia medida experimentalmente ($P = I_{exp} \times U_{exp}$). Esto aísla el rendimiento de diferentes modelos de EOS (SESAME, BLF, H´ebert et al.).
   • Enfoque B (Conductividad Experimental): El término fuente se deriva de la conductividad eléctrica experimental ($\sigma_{exp}$) combinada con el resolvedor del circuito. Esto valida el acoplamiento del resolvedor de corriente con la hidrodinámica.
   • Enfoque C (Totalmente Numérico): Un enfoque autosensible donde la conductividad eléctrica se calcula a partir de un modelo tabulado (Lee-More/Desjarlais modificado) basado en el estado termodinámico local (densidad y temperatura) de cada celda hidrodinámica. Esto elimina la dependencia de las señales de entrada experimentales.

Resultos Clave

• Validación Eléctrica: El modelo de circuito propuesto, incluyendo la dinámica de la chispa, mostró una excelente concordancia con las mediciones de corriente de cortocircuito a través de una amplia gama de generadores y condiciones operativas (50 kA a 5 MA).
• Selección de EOS: Al utilizar la potencia experimental como fuente, los modelos de EOS BLF y H´ebert et al. reprodujeron significativamente mejor las velocidades de expansión, densidades y energías internas que el modelo SESAME, particularmente respecto a la transición sólido-líquido y la subsiguiente meseta de velocidad.
• Validación Acoplada: Utilizando la conductividad experimental como término fuente, la simulación acoplada reprodujo con precisión las señales de corriente y voltaje experimentales. Los modelos BLF y H´ebert et al. proporcionaron la mejor concordancia para la evolución del voltaje, aunque se observaron ligeras discrepancias en los cálculos del espesor de la lámina.
• Simulación Totalmente Numérica: El enfoque totalmente autosensible (Enfoque C) reprodujo con éxito las tendencias generales de la velocidad de expansión, la corriente y el voltaje. Aunque se observaron ligeros desplazamientos temporales (atribuidos a la interacción entre la temperatura predicha por la EOS y el modelo de conductividad), las simulaciones enmarcaron eficazmente el comportamiento experimental. El estudio destacó que las discrepancias en el enfoque totalmente numérico probablemente surgen de la influencia combinada de la EOS (que afecta las predicciones de temperatura) y el modelo de conductividad, más que de una deficiencia única.

Significación y Reivindicaciones
El artículo afirma que este enfoque numérico constituye un "método robusto y eficiente" para diseñar y optimizar futuros experimentos de Materia de Densidad Cálida utilizando instalaciones de potencia de pulso. Al acoplar con éxito la respuesta eléctrica del conductor con la evolución hidrodinámica de la carga, el marco permite:

1. Optimización del Diseño: Predecir la trayectoria termodinámica de los objetivos antes de la experimentación.
2. Prueba de Modelos: Proporcionar una metodología directa para probar diferentes modelos de EOS y de conductividad eléctrica bajo condiciones relevantes.
3. Acceso a Cantidades No Medibles: Otorgar acceso a propiedades como la temperatura y la energía interna que son difíciles o imposibles de medir directamente en estos experimentos.

Los autores enfatizan que, si bien las simulaciones totalmente numéricas muestran ligeras discrepancias con respecto a los datos experimentales, el acuerdo general valida el marco como una herramienta poderosa para interpretar experimentos existentes y guiar diseños futuros, sin necesidad de inventar nuevas aplicaciones más allá del alcance de los experimentos de láminas de potencia de pulso descritos.