Red De Resistencia Térmica Y Simulación De CFD: Metodología De Ingeniería Para El Diseño Cuantitativo De Disipadores De Calor
Primero, de una red de resistencia térmica unidimensional a un campo de temperatura tridimensional
El punto de partida del diseño del disipador de calor es a menudo un diagrama de la red de resistencia térmica. La trayectoria del calor desde la unión de la viruta hasta el aire ambiente se descompone en: unión a la cáscara (Rhtjc, resistencia interna del paquete de viruta), cáscara al disipador de calor (Rhtcs, resistencia térmica TIM), disipador de calor al medio ambiente (Rhtsa, convección + radiación). Entre ellos, Rhtsa se puede descomponer en la resistencia térmica de difusión del substrato del disipador de calor (Rhtspread), la resistencia térmica de conductividad térmica unidimensional de las aletas (Rhtfin) y la resistencia térmica de convección (Rhtconv). Modelo de circuito de serie: resistencia térmica total = Rhtjc + Rhtcs + Rhtspread + Rhtfin + Rhtconv.
Este método de parámetros centralizados es rápido y efectivo en las estimaciones iniciales, pero el mayor inconveniente es que supone una distribución uniforme de la temperatura, cuando en realidad hay un violento efecto de difusión térmica 2D / 3D en el sustrato debajo del chip. Para chips avanzados con flujo de calor local de hasta 200 W / cm ², la resistencia térmica de difusión puede dominar e incluso llevar a la formación de "puntos calientes" en el sustrato, lo que hace que la temperatura local sea mucho más alta que la temperatura promedio. La simulación de CFD debe confiarse en este momento.
Las ecuaciones básicas de la simulación de dinámica de fluidos computacional
CFD resuelve tres ecuaciones diferenciales parciales acopladas:
Ecuación de continuidad
(Conservación de la masa): Haría / Haría + ∇ · (ρ u) = 0Ecuación de impulso
（Navier-Stokes）：∂(ρu)/∂t + ∇·(ρuu) = -∇p + ∇·(τ) + ρgEcuación de energía
: ○ (ρh) / ○ t + ∇ · (ρuh) = ∇ · (k ∇ T) + S _ h
Para la conducción de calor dentro de un disipador de calor, la ecuación de energía se reduce a la ecuación de conducción de calor sólido (con término de convección cero). Para el dominio del aire, es necesario resolver un modelo de turbulencia completo (el modelo k-ε más común o el modelo k-ω SST más avanzado) para capturar con precisión la velocidad y la temperatura dentro de la capa límite cercana a la pared, porque
El 80% del coeficiente de transferencia de calor convectivo depende de la subcapa viscosa en la capa límite con un espesor de solo unas pocas decenas de micras
。
III. Ansys Icepak: Frente a Superficies Complejas y Multifísica
Icepak se basa en el solver Fluent y utiliza una malla no estructurada (núcleo tetraédrico / hexagonal), que es altamente adaptable a geometrías curvas (por ejemplo, aletas circulares, conductos de aire de forma especial). La ventaja única de Icepak es que se puede acoplar sin problemas con Ansys Mechanical y Maxwell para electro-thermal-structural análisis de tres campos. Por ejemplo, los disipadores de calor en amplificadores de RF de alta potencia deben considerar tanto la fuente de calor distribuida espacialmente generada por pérdidas electromagnéticas (de Maxwell), el cambio en la resistencia térmica de contacto causado por la deformación térmica (de Mechanical) y la evaluación de la vida útil de fatiga bajo ciclos térmicos transitorios. Esta simulación acoplada es mucho más precisa que el análisis térmico aislado.
Como estrategia de malla, Icepak recomienda generar mallas prismáticas en la interfaz sólido-fluido, al menos 3-5 capas, para resolver la capa de temperatura de la capa límite. Para los disipadores de calor de CPU típicos, las mallas son típicamente de 5 millones a 20 millones, y el tiempo de solución es de aproximadamente 2-4 horas en una estación de trabajo de 16 núcleos.
FloTHERM: El rey de la eficiencia centrándose en la disipación electrónica del calor
Simcenter FloTHERM utiliza mallas cartesianas (mallas ortogonales), que se generan casi instantáneamente sin la intervención del usuario. Aunque la aproximación de la geometría curva produce un error de paso, este error puede controlarse dentro del rango aceptable de ingeniería para los disipadores de calor de aleta plana que se encuentran comúnmente en la electrónica de consumo (
FloTHERM
Centro de Comando
El módulo tiene un poderoso DOE (Diseño de Experimentos) y capacidades de optimización. Los ingenieros pueden definir funciones objetivas (resistencia térmica mínima o peso mínimo), establecer variables de diseño (altura de la aleta, espaciado, grosor, velocidad del ventilador) y dejar que el software repita automáticamente cientos de simulaciones para encontrar las leyes de las pocas fronteras vitales. Este proceso es casi imposible de hacer manualmente.
V. Trampas clave para simular condiciones límite
La precisión de la simulación depende en gran medida de la autenticidad de las condiciones de contorno de entrada. Aquí hay tres trampas comunes:
Error de suposición de la fuente de calor
: Simplifique el chip como una fuente de calor de superficie uniforme, ignorando la distribución de múltiples puntos de acceso en su interior. La práctica avanzada es utilizar el mapa de distribución de energía proporcionado por el fabricante del chip o calibrarlo a través de la medición del termopar.La convección natural no activa el término gravitacional
En el enfriamiento por convección natural, la flotabilidad es la única fuerza motriz. Sin activar el término de gravedad y ajustar la densidad del aire a la aproximación de Boussinesq, los resultados de la simulación predicen incorrectamente que hay poco flujo y la temperatura es inusualmente alta.La radiación es ignorada o sobreestimada
: Cuando la temperatura de la superficie es inferior a 100 ° C, la radiación generalmente representa solo el 5-15% de la disipación de calor total, lo que se puede simplificar. Pero si la superficie está ennegrecida con alta emisividad (emisividad> 0,9) y el caudal de aire es extremadamente bajo (
Seis, verificación de independencia de la red y criterios de convergencia
Cualquier simulación de CFD debe verificarse la independencia de la red antes del análisis formal. Método de operación: Genere tres conjuntos de rejillas gruesas, medianas y finas (el número de rejillas difiere al menos 2 veces) y calcule la temperatura de las ubicaciones clave (como la temperatura de unión de virutas). Diferencias entre los resultados de la rejilla y la rejilla fina
Los criterios de convergencia suelen fijarse de la siguiente manera: los residuos de energía caen por debajo de 1e-6, los residuos de impulso caen por debajo de 1e-4 y los cambios de temperatura en los puntos de control son inferiores a 0,01 ° C durante 100 iteraciones consecutivas.
VII. Calibración de circuito cerrado desde la simulación hasta el ensayo
La simulación nunca es lo mismo que la realidad física. El proceso de desarrollo más riguroso es: diseño de simulación térmica prueba térmica de muestra de producción de molde abierto (utilizando cámara de imagen térmica y termopar) prueba de comparación y parámetros de simulación de calibración de desviación de simulación (como correlación de convección del lado del aire, desviación de espesor TIM) diseño de corrección prueba secundaria. Después de dos rondas de bucle cerrado, la diferencia de temperatura entre simulación y prueba se puede controlar dentro de ±3 ° C. Esta base de datos de calibración es el activo de conocimiento central de la empresa.
VIII. Conclusión
La simulación térmica ha revolucionado el paradigma de diseño de los disipadores de calor, pasando de "experience-plus-test" a "diseño predictivo". Pero el software es solo una herramienta, y la verdadera experiencia radica en configurar el modelo físico correcto, interpretar con precisión los resultados de la simulación y calibrar continuamente el modelo a través de pruebas térmicas. En el futuro, con la proliferación de simulación asistida por IA y computación de alto rendimiento en la nube, la simulación térmica en tiempo real (Digital Twin) será posible: cada disipador de calor tendrá su gemelo digital, que refleja la temperatura de funcionamiento en tiempo real y predice la vida útil restante.
BQUQ es un fabricante profesional de disipadores de calor de metal, envíenos dibujos y nuestra empresa le cotizará en un plazo de 12 horas.
