Tutorial de Simulación Electro-Térmica (Parte B): Malla Térmica, Condiciones de Contorno y Ajustes Experimentales

1. Introducción: alcance y enfoque de la Parte B

En Parte A ejecutó una simulación electro-térmica, identificó los principales mecanismos de generación de calor (calentamiento por transporte de portadores, calentamiento por recombinación y pérdidas parásitas), y aprendió a extraer salidas de temperatura de la red cristalina como función del sesgo aplicado. Esa primera parte se centró en qué calor se genera y dónde aparece en el dispositivo. En esta Parte B, nos centramos en la otra mitad del problema: cómo se transporta el calor fuera del dispositivo. Esto está gobernado por la malla térmica y, de forma crítica, por las condiciones de contorno térmicas. Las condiciones de contorno codifican cómo el dispositivo está montado, enfriado o conectado a un disipador en la práctica, y a menudo tienen un impacto de primer orden en el aumento de temperatura resultante. Dos simulaciones con comportamiento eléctrico idéntico y términos de generación de calor idénticos pueden producir perfiles de temperatura muy diferentes únicamente debido a distintas suposiciones de contorno térmico.

Por lo tanto, verá por qué la simulación electro-térmica normalmente utiliza un dominio térmico que es mayor que el dominio eléctrico, cómo se utilizan condiciones de contorno efectivas para representar disipadores de calor macroscópicos sin mallarlos explícitamente y por qué las cantidades de material dependientes de la temperatura se precalculan en una rejilla discreta de temperatura. Este tutorial continúa utilizando un ejemplo de diodo orgánico porque está fundamentado experimentalmente e incluye efectos medibles de auto-calentamiento, pero el flujo de trabajo es completamente general. Los mismos conceptos se aplican a semiconductores inorgánicos, dispositivos de película delgada, electrónica de potencia, fotodetectores y cualquier estructura donde la disipación de potencia eléctrica retroalimente el transporte a través de la temperatura.

Al final de la Parte B, debería comprender:

• por qué la malla térmica es distinta de la malla eléctrica,
• cómo las condiciones de contorno térmicas controlan la extracción de calor,
• cómo se utiliza la discretización de temperatura internamente para la eficiencia y
• por qué a veces debe incluirse el auto-calentamiento para coincidir con los datos experimentales.

2. Malla térmica y condiciones de contorno

Comience en la ventana principal de simulación y navegue a la cinta Thermal. Esta es la misma cinta introducida en la Parte A: contiene los conmutadores de fuentes de calentamiento y los controles principales de configuración (parámetros térmicos, malla térmica y condiciones de contorno). La cinta se muestra en ??.

2.1 La malla térmica y por qué no es la malla eléctrica

Haga clic en Thermal mesh para abrir el editor de malla térmica. A primera vista se parece a los editores de malla que ha visto en otros lugares en OghmaNano, pero está resolviendo un problema físico diferente: difusión de calor. El punto clave es que la simulación electro-térmica utiliza dos mallas, porque los problemas eléctricos y térmicos suelen vivir en escalas de longitud fundamentalmente diferentes:

• La malla eléctrica resuelve el transporte (drift–diffusion + Poisson + recombinación/trampas). Normalmente se concentra donde la acción eléctrica es más intensa: la capa activa y las interfaces cercanas, donde los campos y las densidades de portadores varían rápidamente. Estas variaciones pueden ocurrir en escalas de nanómetros a micrómetros.
• La malla térmica resuelve la ecuación de calor de la red cristalina. El calor no se detiene en la capa activa: se propaga por toda la pila, hacia contactos, sustratos, encapsulados y hacia cualquier elemento que elimine calor del dispositivo. Estas rutas suelen estar en escalas de milímetros a centímetros.

Esta discrepancia es la razón por la cual una sola “malla compartida” suele ser una abstracción incorrecta. El problema eléctrico necesita alta resolución en la región activa; el problema térmico necesita un dominio suficientemente grande para representar la trayectoria de flujo de calor. Por ejemplo:

• Un diodo de película delgada puede tener un espesor de capa activa de 100–300 nm, pero el sustrato puede tener 0.5–1 mm de espesor (vidrio) y el bloque de montaje o la etapa pueden tener centímetros.
• Un contacto metálico puede tener solo decenas de nanómetros de espesor eléctricamente, pero puede ser el principal dispersor lateral de calor térmicamente porque su conductividad térmica es grande en comparación con capas orgánicas u óxidos.

En este ejemplo, puede ver que la altura Y del dominio térmico abarca todo el espesor de la pila del dispositivo en lugar de estar restringida a la región eléctricamente activa. Esto es típico: el transporte térmico involucra todas las capas que conducen calor, no solo las capas donde fluye la corriente. A continuación configuramos cómo se permite que el calor salga de este dominio térmico utilizando condiciones de contorno. La malla define la región en la que el calor puede difundirse; las condiciones de contorno definen lo que ocurre en los bordes de esa región.

2.2 Condiciones de contorno: caras aislantes y un disipador efectivo

Abra el editor de condiciones de contorno mostrado en ??. En este ejemplo, la mayoría de las caras están configuradas como Neumann (a menudo mostrado como “Neumann (==0)”). Físicamente, una condición de contorno de Neumann especifica el flujo de calor normal en el límite. Cuando se establece en cero, impone:

\[ -k \nabla T \cdot \hat{n} = 0 \]

es decir, no fluye calor a través de ese límite. Estas caras se tratan como térmicamente aislantes. Esto no significa que el dispositivo esté térmicamente aislado en general; simplemente indica al solucionador térmico que esos límites no forman parte de la trayectoria prevista de eliminación de calor. La excepción aquí es la cara y_max, que está configurada como Heatsink. Este límite especifica una temperatura del dispositivo en y_max (aquí alrededor de 300 K), junto con una conductividad del disipador efectiva y una longitud del disipador (del orden de un milímetro en este ejemplo).

Conceptualmente, este es un modelo efectivo de eliminación de calor. Permite que el calor fluya fuera de la pila de dispositivos simulada hacia un disipador no mallado explícitamente, sin obligar a que la malla térmica incluya un objeto macroscópico. Esto es importante porque los disipadores de calor reales suelen ser órdenes de magnitud más grandes que los dispositivos de película delgada: mallar explícitamente un disipador a escala milimétrica o centimétrica con resolución submicrométrica sería computacionalmente ineficiente y físicamente innecesario. Cabe destacar que, en problemas térmicos, las condiciones de contorno determinan en gran medida cuán caliente se vuelve el dispositivo. Una ruta deficiente de extracción de calor puede hacer que la temperatura aumente muy rápidamente, mientras que una ruta eficiente puede mantener la temperatura cerca de la ambiental incluso bajo una disipación de potencia significativa.

Una analogía útil es pensar en el dispositivo como una bañera que se está llenando con agua. Los términos de generación de calor son el grifo: vierten energía en el sistema. Las condiciones de contorno térmicas son el desagüe:

• Si el desagüe está completamente abierto (excelente extracción de calor), el nivel de agua permanece bajo.
• Si el desagüe está parcialmente bloqueado (extracción de calor limitada), el nivel de agua sube.
• Si el desagüe está cerrado (sin extracción de calor), la bañera finalmente se desbordará.

En esta analogía, el nivel del agua es la temperatura de la red cristalina. Por lo tanto, el propósito de las condiciones de contorno térmicas no es cosmético ni secundario: codifican el entorno físico que en última instancia controla si el dispositivo funciona frío, templado o catastróficamente caliente.

2.3 “Puntos de malla” de temperatura (por qué se discretiza la temperatura)

El editor de malla térmica también muestra un rango de temperatura y un número de puntos de temperatura. Esto no es la malla térmica espacial: es la rejilla de temperatura utilizada para precálculos. Muchas partes del modelo eléctrico/térmico requieren que cantidades dependientes de la temperatura se evalúen repetidamente durante la solución acoplada, y es computacionalmente costoso recalcularlas desde cero en cada temperatura intermedia. Por ejemplo, dependiendo de las estadísticas elegidas y del modelo de densidad de estados, el solucionador puede necesitar relaciones dependientes de la temperatura para densidades de portadores, funciones cuasi-termodinámicas y cantidades de consulta relacionadas. En lugar de evaluarlas “en vivo” cada vez que el bucle acoplado actualiza \(T_L\), OghmaNano precalcula tablas de fondo en un número finito de puntos de temperatura y luego interpola entre ellas durante la ejecución. En este ejemplo las tablas se generan en 7 puntos entre 290 K y 350 K.

Esto convierte el rango de temperatura en una elección de modelado: debe abarcar cómodamente las temperaturas esperadas durante el auto-calentamiento. Si el dispositivo se calienta más allá de la temperatura superior, se corre el riesgo de artefactos de extrapolación (o limitación dependiendo de la configuración), lo cual rara vez es deseable. Como regla práctica: elija un rango que cubra las temperaturas de operación previstas con margen.

3. Parámetros térmicos (conductividad térmica y tiempos de relajación)

En la simulación electro-térmica, las condiciones de contorno definen cómo puede salir el calor del dispositivo, pero los parámetros del material definen cómo se mueve el calor dentro de la pila. En la práctica, una de las cantidades más importantes que controlan el aumento de temperatura y la formación de puntos calientes es la conductividad térmica \(k\) (a veces escrita \(\kappa\)).

Para ver o editar estos parámetros, abra la cinta Thermal y haga clic en Thermal parameters (el botón marcado con \(k\) / \(\kappa\)). Esto abre el editor de parámetros térmicos capa por capa mostrado en ??.

El editor expone tres parámetros clave para cada capa:

• Conductividad térmica \(k\): gobierna la eficiencia con la que el calor se difunde a través de esa capa. Las capas con bajo \(k\) actúan como cuellos de botella térmicos y pueden aumentar fuertemente el incremento de temperatura de red cristalina predicho.
• Tiempo de relajación de electrones y tiempo de relajación de huecos: solo son necesarios cuando se utiliza el modelo hidrodinámico / de transporte de energía, donde las temperaturas de los portadores pueden desviarse de la temperatura de la red cristalina. En el modelo estándar de calor de red cristalina no son necesarios porque las temperaturas de los portadores están fijadas a \(T_L\).

En la mayoría de los flujos de trabajo, las conductividades térmicas se toman de valores publicados (o medidos cuando están disponibles) y luego se refinan solo si existe una justificación experimental clara. Junto con las condiciones de contorno, \(k\) es uno de los principales determinantes de la temperatura absoluta del dispositivo y de la estructura espacial de los puntos calientes.

4. Cómo funciona el solucionador electro-térmico acoplado (y por qué es más lento)

Una simulación electro-térmica en OghmaNano consiste en dos solucionadores acoplados:

• un solucionador eléctrico completamente acoplado (drift–diffusion y Poisson, incluyendo recombinación y trampas cuando están habilitadas), y
• un solucionador térmico que resuelve la ecuación de difusión de calor para el campo de temperatura de red cristalina \(T_L\) utilizando términos de generación de calor calculados a partir de la solución eléctrica.

A un voltaje aplicado dado, el solucionador no ejecuta estas etapas una sola vez en secuencia. En su lugar, ejecuta un bucle externo autoconsistente:

1. resolver el problema de transporte eléctrico utilizando el campo de temperatura actual,
2. calcular los términos de generación de calor (transporte, recombinación, parásitos y opcionalmente ópticos),
3. resolver el problema de difusión térmica para actualizar \(T_L\), y
4. repetir hasta que tanto los residuos eléctricos como los térmicos hayan convergido.

Este acoplamiento iterativo es la razón por la cual las simulaciones electro-térmicas son más lentas que las puramente eléctricas: para cada punto de sesgo, la resolución de Newton eléctrica y la resolución térmica pueden ejecutarse múltiples veces antes de que el sistema acoplado se estabilice. El resultado es una solución multifísica genuinamente autoconsistente, en lugar de un cálculo eléctrico con la temperatura añadida como paso de posprocesamiento.

El acoplamiento es esencial porque los problemas eléctricos y térmicos suelen vivir en escalas de longitud físicas muy diferentes. El transporte eléctrico está dominado por la región activa, donde campos, densidades de portadores y tasas de recombinación varían rápidamente. El transporte térmico debe tener en cuenta toda la trayectoria de flujo de calor, incluidos contactos, sustratos, encapsulado y disipadores térmicos. Por esta razón OghmaNano proporciona controles independientes para la malla térmica y las condiciones de contorno térmicas.

También es importante reconocer que el marco electro-térmico implica naturalmente múltiples temperaturas. En general, el modelo distingue entre la temperatura de la red cristalina \(T_l\), que describe la temperatura de la red atómica, la temperatura de electrones \(T_e\), y la temperatura de huecos \(T_h\).

Los electrones y huecos son poblaciones de portadores móviles que pueden, en principio, tener sus propias energías locales, del mismo modo que tienen sus propios niveles cuasi-Fermi. En el modelo electro-térmico estándar utilizado aquí, las temperaturas de electrones y huecos están acopladas a la temperatura de la red cristalina, por lo que todos los subsistemas comparten un único campo de temperatura.

En regímenes extremos, OghmaNano también admite un modelo de transporte hidrodinámico en el que las temperaturas de electrones y huecos pueden desviarse de la temperatura de la red cristalina. Ese modelo avanzado se describe en otro lugar y solo es necesario en situaciones especializadas. Para la mayoría de las simulaciones de dispositivos, incluido este tutorial, el modelo electro-térmico basado en la red cristalina es apropiado y suficiente.

5. Auto-calentamiento y comparación con el experimento

Este dispositivo de ejemplo ha sido comparado con datos experimentales utilizando el modelo electro-térmico. No se realiza ningún ajuste en este tutorial. En su lugar, la comparación existente se utiliza para ilustrar cómo el auto-calentamiento altera la respuesta eléctrica del dispositivo, y por qué se requiere un tratamiento electro-térmico acoplado a corrientes de operación más altas.

El flujo de trabajo de comparación se accede mediante la barra de herramientas Automation, mostrada en ??. Al abrir la herramienta fit-to-experiment se muestran los datos experimentales de JV para este dispositivo, como se muestra en ??. Esta curva representa el comportamiento medido del dispositivo bajo sesgo.

Ejecute una iteración. Después de una única resolución electro-térmica, la curva JV simulada se superpone a los datos experimentales, como se muestra en ??. En este caso, la concordancia es cercana, lo que indica que el modelo captura la física dominante de transporte y térmica del dispositivo.

El papel del auto-calentamiento se vuelve claro repitiendo la misma comparación con el modelo térmico deshabilitado:

1. Regrese a la cinta Thermal.
2. Deshabilite el modelo térmico.
3. Ejecute una iteración de la comparación nuevamente.

La curva JV resultante ya no coincide con los datos experimentales, como se muestra en ??. La desviación es más pronunciada a sesgos más altos, donde las densidades de corriente son grandes y el auto-calentamiento por transporte de portadores se vuelve significativo.

Esta comparación demuestra un resultado central del modelado electro-térmico: una vez que los dispositivos operan en regímenes de alta densidad de corriente, la disipación de potencia eléctrica retroalimenta el transporte a través del aumento de temperatura. Las simulaciones eléctricas a temperatura fija no pueden capturar este efecto, mientras que el modelo electro-térmico acoplado sí puede.