Introducción

COMSOL Multiphysics ofrece una plataforma dinámica para modelar celdas de combustible proporcionando un conjunto versátil de herramientas y capacidades. Permite simulaciones exhaustivas, lo que permite a investigadores e ingenieros analizar diversos aspectos del rendimiento de las celdas de combustible. Desde comprender reacciones electroquímicas hasta estudiar el flujo de fluidos y transferencia de calor dentro de las celdas, COMSOL facilita investigaciones detalladas. Su enfoque de múltiples físicas permite la consideración simultánea de varios fenómenos físicos, brindando conocimientos sobre interacciones complejas dentro de las celdas de combustible.

Celdas de combustible de óxido sólido

Como ejemplo, revisamos el artículo escrito por Iliya K. Iliev et al. 2023 [1] con su trabajo titulado Numerical Simulation of Processes in an Electrochemical Cell Using COMSOL Multiphysics. Su estudio tiene como objetivo crear un modelo para Celdas de Combustible de Óxido Sólido (SOFC) que funcionen con combustible basado en carbono. Este modelo integral considera el calentamiento electroquímico, el flujo de fluidos no isotérmicos y múltiples fenómenos físicos como transporte de electrones y iones, difusión de gases, reacciones electroquímicas y transferencia de calor. Ofrece información sobre la concentración de reactantes, temperatura, distribución de corriente y cálculos de potencia de salida.

Los autores informan que las celdas alimentadas con hidrógeno superaron a las alimentadas con combustible basado en carbono, generando 1340 W/m2 a 0.25 A/cm2 en comparación con 1200 W/m2 para el metano a la misma densidad. El gas de síntesis alcanzó 1340 W/m2 a 0.3 A/cm2. Se observaron perfiles de temperatura distintos y cambios en la concentración de reactantes entre las celdas alimentadas con hidrógeno y carbono. El modelo se alineó estrechamente con la característica voltio-amperio de un prototipo industrial, mostrando una desviación menor al 5% en la mayoría de los rangos. Aunque es un paso significativo, los estudios futuros tienen como objetivo refinar suposiciones, incorporar geometría multicanales y considerar la radiación térmica para simulaciones más precisas, mejorando aún más el modelado de SOFC.

En la figura de la cabecera se muestra la fracción molar de O2 y H2 para un voltaje de 0.5 V. La imagen se obtuvo del modelo de la librería de aplicaciones de COMSOL: "Current Density Distribution in a Solid Oxide Fuel Cell, Application ID: 514", y que se estudia en el artículo citado en [1].

Comentarios finales

Las conclusiones de este estudio tienen una relevancia significativa en el campo de las Celdas de Combustible de Óxido Sólido (SOFC). Presentan un modelo integral, construido con COMSOL Multiphysics 6.1, que simula con precisión los procesos físicos complejos dentro de una SOFC. Al incorporar varios factores como el transporte de electrones y iones, la difusión de gases y las reacciones electroquímicas, este modelo sirve como una herramienta valiosa para comprender el comportamiento de las SOFC. La comparación del rendimiento de las SOFC cuando se alimentan con diferentes fuentes (hidrógeno, metano y gas de síntesis) brinda información crucial. Notablemente, los hallazgos ilustran salidas de potencia y distribuciones de temperatura variables entre diferentes tipos de combustible. Las reducciones de temperatura demostradas debido a reacciones endotérmicas en combustibles carbonáceos, en contraste con el hidrógeno, destacan la capacidad del modelo para simular comportamientos matizados.

Referencias

[1] Iliev, I.K.; Gizzatullin, A.R.; Filimonova, A.A.; Chichirova, N.D.; Beloev, I.H. Numerical Simulation of Processes in an Electrochemical Cell Using COMSOL Multiphysics. Energies 2023, 16, 7265. https:// doi.org/10.3390/en16217265

En nuestra última exploración de COMSOL Multiphysics v6.2, subrayamos una de sus características más destacadas: los modelos subrogados (surrogate models). Estos modelos, centrales en nuestra discusión de hoy, simplifican simulaciones complejas manteniendo la precisión. Ofrecen habilidades computacionales más rápidas, siendo inestimables para optimizar diseños, realizar estudios paramétricos y explorar extensos espacios de diseño eficientemente.

El Estudio de Entrenamiento de Modelos Subrogados en COMSOL Multiphysics introduce funcionalidades críticas para la creación de estos modelos eficientes. Énfasis en métodos estratégicos de muestreo DOE como el Latin hypercube sampling (LHS) para entrenamiento de modelos, equilibrando el tiempo de generación de datos con la precisión. Dentro de este estudio, se encuentran varios modelos subrogados como Deep Neural Network (DNN, ver fig. de la cabecera, donde se muestra una red neuronal de 5 capas. Imagen obtenida del ejemplo de COMSOL "A Tubular Reactor Surrogate Model"), Gaussian Process (GP) y Polynomial Chaos Expansion (PCE)*, cada uno con fortalezas únicas. Por ejemplo, DNN maneja conjuntos de datos más grandes, mientras que GP y PCE ofrecen estimaciones de incertidumbre. Estos modelos sobresalen en la interpolación de funciones multidimensionales, abordando intrincadas relaciones de datos no lineales. Las funciones de modelos subrogados extienden su utilidad más allá de aplicaciones y cuantificación de incertidumbre, encontrando aplicación en la representación de datos de materiales, optimización y más, permitiendo diferenciación en cuanto a parámetros de entrada.

La versión 6.2 ofrece ejemplos notables, como Tubular Reactor Surrogate Model application [1], Thermal Actuator Surrogate Model application [2] y Surrogate Model Training of a Battery Rate Capability model [3]. *El modelo subrogado DNN viene incluido en COMSOL Multiphysics, los modelos de GP y PCE están disponibles con el módulo de cuantificación de incertidumbre.

Referencias

[1] El Tubular Reactor Surrogate Model application muestra el importante aumento de velocidad en los cálculos al utilizar un modelo subrogado en lugar de un modelo de elemento finito exhaustivo
[2] El Thermal Actuator Surrogate Model application describe un método para agilizar el análisis multiphysics, utilizando un modelo subrogado diseñado para un modelo de geometría paramétrica completa.
[3] El Surrogate Model Training of a Battery Rate Capability Model demuestra el uso de una función de modelo subrogado, DNN, para entrenar y predecir la capacidad de carga de una celda de batería NMC/grafito.

Histograma por grupos realizado con Minitab

O cómo te ayudan un histograma y una serie temporal a llegar al fondo de la cuestión.

En anteriores artículos repasamos las principales características del histograma y cómo interpretarlo. En esta entrada vamos a explicar un caso real de cómo la visualización e interpretación de un histograma junto con una serie temporal ayudaron a encontrar un problema que existía en la fábrica desde hacía un año.

El problema

Desde el departamento de calidad recordaron que en las auditorías de revisión de registros de Excel siempre encontraban valores altos de fricción, pero dentro de tolerancias. Al hablarlo con los operarios, éstos afirmaron que se fiaban más de sus pruebas manuales de fricción que de los valores que medía la máquina. 

Analizando con pensamiento (estadístico) visual

Histograma de los datos de fricciónDesde calidad decidieron empezar por entender los datos históricos existentes de fricción que habían sido introducidos manualmente por dos operarios. Los datos de fricción medidos de forma automática estaban guardados en la máquina y a priori no era fácil acceder a ellos.

Como una primera forma de visualización de los datos entrados manualmente, se realizó el histograma al que se añadió una línea vertical con el límite de tolerancias en 4.2.

En el histograma se observó que había dos grupos de datos: unos datos alrededor del valor 1.9, muy lejos del límite de tolerancias y, por lo tanto, muy buenos, y otros que sospechosamente estaban muy cerca del límite de tolerancias, pero siempre por debajo.

Histograma de los datos por gruposAl indagar sobre el origen de los dos grupos se descubrió que correspondía cada uno a un operario, como se puede ver en el histograma por grupos de abajo. Los datos estaban siempre por debajo del límite de tolerancias (indicando que todo iba bien), pero siendo diferente dependiendo de cada operario: el operario 2 reportaba valores cerca tolerancias (pero siempre por debajo) mientras que el operario 1 reportaba valores más alejados de tolerancias (alrededor de 1.9).

Los operarios admitieron que los datos entrados manualmente eran valores “maquillados” porque los valores reales de fricción que daba la máquina eran muy malos, y ellos sabían que la fricción era buena, puesto que los dos elementos de la pieza se podían desplazar con las manos.

Finalmente, se logró obtener el registro de datos reales de la máquina y se pudieron comparar con los manuales. En este caso, para visualizar los datos se utilizó un gráfico de serie temporal. En el gráfico de serie temporal se grafica el tiempo en el eje horizontal (eje x) ya sea en formato fecha, hora, minutos desde el inicio… y en el eje vertical (eje y) se grafica el valor de la variable en cada instante de tiempo. Al tratarse de datos que tienen una relación de tiempo (unos van antes que los otros) se unen los puntos del gráfico con líneas para trasladarnos esta información de “orden”. En un solo gráfico se pueden representar varias series de temporales para facilitar la comparación.

Gráfico de dos series temporales: registro de máquina y registro de operarioPara comparar los datos manuales con los automáticos se representaron las dos series temporales en un solo gráfico al que se añadió una línea horizontal en el valor 4.2 que era el límite de tolerancias.

Se pudo ver que los valores de registro de máquina estaban siempre por encima de tolerancias (indicando que algo iba mal siempre). Por lo tanto, tenemos evidencias de que los dos operarios “maquillan los datos originales” pero de una forma diferente: uno los hacía “buenos” dejándolos cerca de tolerancias y el otro anotaba el valor que a criterio suyo se aproximaría al valor real.

Finalmente, se presentaron estos datos a ingeniería que lograron encontrar la causa raíz de la mala medición. ¡Había una fuga en la máquina que generaba mediciones aberrantes y que afectaba solo a esta medida y que llevaba 1 año sin haberse detectado! Los trabajadores, que sabían que los datos automáticos estaban mal, se habían acostumbrado a maquillar los datos, lo cual provoca registros de calidad falsos y síntomas de variación que conducen a hipótesis erróneas.

¿QUIERES APLICAR ESTAS HERRAMIENTAS A TUS DATOS?

Realizar este tipo de estudios no es una tarea difícil si tienes un buen conocimiento de las herramientas de visualización de datos que ofrece Minitab. Pero ¡ojo! No se trata de tener una noción profunda de cada herramienta (con todas sus parafernalias matemáticas) sino de ser capaz de imaginar qué herramienta sería la más adecuada en cada momento (y planificar una recogida de datos adecuada a la herramienta) y tener la habilidad de comprender e interpretar lo que dicen.

Precisamente esto es lo que te enseñaremos en el curso de Diagnóstico con visualización de datos con Minitab.

Por Stacey McDaniel

En los últimos años, un número cada vez mayor de agencias gubernamentales y policiales utilizan drones para obtener una vista aérea desde el cielo. Para ser efectivos, los drones deben poder maniobrar rápidamente hacia el objetivo previsto. Las maniobras rápidas les permiten navegar en entornos complejos y realizar sus tareas de manera eficiente. Los requisitos para convertirse en piloto certificado de drones varían según el país, pero todos requieren que se apruebe un examen escrito sobre conocimientos aeronáuticos administrado por la autoridad de aviación. Sorprendentemente, (en el momento de escribir esta publicación) no existe ningún requisito para demostrar el desempeño. Entonces, ¿cómo podemos comparar las habilidades de vuelo de dos pilotos certificados de drones?

GUIÓN:

Dos pilotos certificados de drones, Aleks y Dustin, son evaluados por la rapidez con la que posicionan el drone en un objetivo específico, en este caso a unos 150 metros de distancia de la posición de lanzamiento.

Queremos comparar para ver si la destreza de Aleks y Dustin es diferente entre sí. El tiempo se mide en segundos cuando comienzan (incluye elevación, movimiento hacia el objetivo). Una diferencia práctica es de 0,5 segundos, lo que significa que si un piloto puede colocar el dron en posición ½ segundo más rápido, será mejor. Aquí están los datos experimentales:

Con Minitab Statistical Software, podemos crear rápidamente estadísticas descriptivas.

Vemos que en diez vuelos, el tiempo promedio de Aleks para apuntar es de 15,86 segundos, mientras que el de Dustin es de 12,42 segundos.

La diferencia estimada es de 3,4 segundos. El intervalo de confianza del 95% nos dice que la verdadera diferencia promedio entre los pilotos es tan pequeña como 1,736 y tan grande como 5,131 segundos. Con base en este intervalo de confianza, podemos decir con confianza que los tiempos de vuelo no solo son diferentes de cero, sino que debido a que el extremo inferior del intervalo de confianza para la diferencia es mayor que nuestra diferencia práctica de 0,5 segundos, concluimos que Dustin es el mejor piloto de drones en esta prueba.

Un gráfico de caja del tiempo de vuelo proporciona una imagen del análisis:

En el futuro, podríamos realizar experimentos de seguimiento que incluyan múltiples distancias, obstáculos, diversas condiciones de vuelo y múltiples drones para que los utilicen los pilotos. Minitab Statistical Software proporciona análisis y visualizaciones, lo que facilita la medición del rendimiento.

La nueva actualización Maple Flow 2023.3 hace un uso efectivo de un motor de cálculo mejorado y añade nuevas características para que sea más eficiente crear y trabajar con hojas de trabajo de diseño.

  • Evaluación más rápida para los cálculos. Con el motor de Maple Flow optimizado, los cálculos se procesan más rápido, las operaciones matriciales son más fluidas y las actualizaciones de parámetros se muestran más rápidamente.
  • Nueva herramienta de dibujo de arco. Con esta adición a la barra de herramientas Dibujo, pueden añadirse líneas curvas a su hoja de trabajo y ajustar el grosor de las líneas y las puntas de flecha para mejorar diagramas y esquemas.
  • Mejora de la productividad al realizar múltiples tareas. La barra de progreso de ejecución ofrece una forma rápida de ver el estado de los cálculos de larga duración desde la barra de tareas de Windows o el Dock de Mac, lo que ayuda a gestionr el tiempo de forma eficaz.
  • Ayuda adicional sobre el producto. El nuevo contenido del sistema de ayuda facilita el inicio de funciones avanzadas, incluídas las técnicas de estadísticas y procesado de señales.
  • Actualizaciones de usabilidad: Varias mejoras hacen que Maple Flow sea más fácil de utilizar, incluyendo:
    • Reorganizar la hoja de trabajo. Ahora se puede mover una sección completa seleccionándola y arrastrándola.
    • Las actualizaciones del Selector de color facilitan la selección del color que se desea para el texto y los gráficos.

La nueva versión también cubre varias mejoras de la interfaz y correcciones de errores en respuesta a la solicitud de los usuarios.

En el mundo de las energías renovables, los módulos fotovoltaicos (PV) destacan como actores clave en el cambio hacia la sostenibilidad. Sin embargo, su rendimiento enfrenta un desafío significativo debido a la humedad en condiciones del mundo real. Comprender esta interacción es crucial para optimizar la longevidad y eficiencia de los paneles solares. En el presente estudio [1], se adopta un enfoque innovador al utilizar configuraciones idénticas en diferentes ubicaciones alrededor del mundo, ofreciendo una comparación única de diversos climas. Los módulos vienen equipados con sensores miniatura dentro ellos. Las configuraciones utilizadas proporcionan perspectivas del comportamiento real de la humedad. Ver figura de la cabecera, donde se muestran Mini módulos expuestos a condiciones exteriores [1].

Lo que distingue a este estudio es el uso estratégico de COMSOL para simulaciones exhaustivas de la entrada de humedad, utilizando el modelo de difusión de Fick y la segunda ley de difusión de Fick. Para mejorar la precisión, se crearon dos geometrías bidimensionales distintas.

La primera geometría, una sección transversal vertical del módulo, introdujo múltiples condiciones de contorno (BC). Una BC de Dirichlet en la interfaz aire-cubierta se calculó a partir de las condiciones del aire medidas a través de la solubilidad de la cubierta, utilizando la ley de Henry de sorción. Se asumió la constante presión parcial de vapor de agua entre el sensor de aire y la interfaz aire-cubierta. El valor de la humedad relativa (RH) en esa interfaz se recalculó utilizando la ecuación de Arden-Buck. Dos BC adicionales en la interfaz cubierta-encapsulante enfatizaron la continuidad del flujo y describieron el estado de equilibrio basado en la ley de Henry.

La segunda geometría simulada abordó la difusión horizontal frente a la celda, corrigiendo discrepancias observadas en simulaciones que dependen únicamente de la primera geometría. Esto también facilitó la integración de grietas en la celda. Las condiciones de contorno de Dirichlet en el borde de la celda incorporaron mediciones del sensor, con datos faltantes completados utilizando resultados de la simulación inicial. La temperatura del módulo, anclada a mediciones del sensor más cercano al centro de la celda, enriqueció la simulación con condiciones ambientales del mundo real.

Como conclusión, este trabajo presentó un novedoso modelo de simulación bidimensional en dos etapas, donde la investigación combinó la difusión de Fick en secciones transversales verticales y horizontales con condiciones de contorno cuidadosamente establecidas. Este enfoque innovador proporcionó una representación más precisa de la difusión de humedad frente a la celda fotovoltaica. Al comparar los resultados de la simulación con las mediciones, el estudio demostró que este enfoque logra un equilibrio entre la precisión de la simulación y la eficiencia computacional, destacando el papel clave de COMSOL Multiphysics en el avance de nuestra comprensión de la dinámica de la humedad en los módulos fotovoltaicos.

Referencias

[1] S. Mitterhofer et al., "Measurement and Simulation of Moisture Ingress in PV Modules in Various Climates," in IEEE Journal of Photovoltaics, doi: 10.1109/JPHOTOV.2023.3323808.
Web: https://ieeexplore.ieee.org/document/10295538.

La nueva versión contiene las siguientes implementaciones:

  • Modelos actualizados: Nuevos modelos US EPA AERMOD & AERMET (versión 23132) tanto en compilaciones de 32 como 64 bits.
  • Versión de AERMOD paralela actualizada: Nueva versión MPI (paralelo) para AERMOD 23132.
  • Elevaciones de terreno: Añadido soporte para fuentes líneas flotantes (buoyant line), RLINE, & RLINEXT cuando se importan elevaciones.
  • Notas adicionales: Revise las notas de la versión completas para una lista detallada de los cambios en el enlace inferior.