See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/311273939

Estudio del Efecto Altura en un Sistema de Empaque de Baterías de Ion Litio

Conference Paper · October 2016

CITATIONS READS
0 10,507

5 authors, including:

Marcelo Cortés Abdiel Mallco

University of Antofagasta University of Antofagasta
40 PUBLICATIONS   405 CITATIONS    18 PUBLICATIONS   19 CITATIONS   

SEE PROFILE SEE PROFILE

Jorge Reyes-Marambio
University of Chile
9 PUBLICATIONS   108 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Subdaily hydrological alteration in Chile View project

Ayllu solar View project

All content following this page was uploaded by Marcelo Cortés on 01 December 2016.

The user has requested enhancement of the downloaded file.

 Estudio del Efecto Altura en un Sistema de Empaque
de Baterías de Ion Litio

Marcelo Cortés-Carmona, Abdiel Mallco C. Williams Calderón, Jorge Reyes, Rodrigo Palma
Universidad de Antofagasta Universidad de Chile
Centro de Desarrollo Energético Antofagasta Centro de Energía
Antofagasta, Chile Santiago, Chile
[email protected], [email protected] [email protected], [email protected]

Resumen—Uno de los efectos más notorios en la pérdida de Las baterías de ion-litio pueden ser consideradas en la
rendimiento y capacidad de sistemas térmicos se produce por el actualidad como uno de los dispositivos ideales para el
efecto de la altura, la que provoca que la densidad de un fluido almacenamiento de energía en el rango de potencias medias,
compresible y la presión atmosférica se vea disminuida las que se utilizan principalmente en autos eléctricos y
considerablemente, lo que acarrea una reducción de potencia aplicaciones estacionarias para apoyo de la red eléctrica
para sistemas eléctricos y térmicos. Ante esto, los empaques de (Battery Energy Storage Systems, BESS). Esto se debe gracias
baterías de ion Litio que cuenten con refrigeración forzada a sus favorables ventajas, como su alta densidad de potencia,
mediante un fluido compresible, se ven afectados directamente
su larga vida útil y su moderada tasa de descarga [4]. Sin
con las condiciones ambientales incidiendo en la pérdida de
embargo, la temperatura de las baterías de ion-litio varía
capacidad de refrigeración del sistema afectando el rendimiento
del empaque cuando trabaja a gran altura. El objetivo de esta
dependiendo del contexto en que se encuentren, ya sea
investigación es cuantificar el impacto que tiene la altitud sobre aumentando o disminuyendo, afectando así las propiedades de
el nivel del mar en el diseño de un pack de baterías de ion litio la batería. Por ejemplo, la temperatura aumenta durante los
que cuentan con refrigeración forzada. procesos de carga y descarga de la batería.
Las futuras tecnologías vehiculares, tienen como objetivo
Palabras claves— Li-ion batteries; Thermal modeling; Thermal introducir vehículos de alto rendimiento eléctrico, de esa
behavior; Simulation
manera se requiere de sistemas inteligentes para controlar y
mantener la temperatura de su empaque de baterías dentro de
I. INTRODUCCIÓN rangos deseables de funcionamiento. Durante el
Las baterías de iones de litio son dispositivos que se usan funcionamiento con altas tasas de carga y descarga, las baterías
principalmente en el almacenamiento de energía para de iones de litio pueden calentarse hasta 55°C o más, sobre
Vehículos eléctricos híbridos (HEV), Vehículos eléctricos todo cuando las celdas se apilan en módulos y paquetes. Por lo
híbridos enchufables (PHEV) y Vehículos Eléctricos (EV), tanto, se requiere de un sistema que pueda modelar la
debido a su larga vida útil, velocidad de carga y densidad generación térmica de la batería, de manera tal que la batería
energética. Sin embargo, uno de los factores que hace que estos pueda mantener su temperatura dentro del rango deseable y
tipos de vehículos no tengan un precio competitivo en el garantizar un funcionamiento óptimo del vehículo en términos
mercado automotriz, es el costo de la batería, la que representa de seguridad, rendimiento y vida útil de la batería.
aproximadamente el 50 % del costo total del vehículo [1]. Ante Uno de los efectos más notorios en la pérdida de
esta eventualidad es necesario desarrollar tecnologías que rendimiento y capacidad de sistemas térmicos se produce por el
puedan reducir los costos, mejorar la densidad energética y que efecto de la altura, la que provoca que la densidad de un fluido
tenga particularmente el control de los sistemas de gestión de compresible y la presión atmosférica se vea disminuida
baterías como también la gestión térmica del sistema [2]. considerablemente acarreando una reducción de potencia para
Uno de los requerimientos importantes para lograr la sistemas eléctricos y térmicos [5]. Ante esto, los empaques de
gestión de la vida útil de la batería, se relaciona con el control baterías de ion litio que usan refrigeración forzada mediante un
de la temperatura de la batería y de esta forma limitar el rango fluido compresible, se ven afectados directamente perdiendo
de temperatura de funcionamiento. En la práctica, sólo la capacidad de refrigeración, y disminución del rendimiento.
temperatura de la superficie de la celda puede ser medida, y se El objetivo de esta investigación es cuantificar el impacto
utiliza para el monitoreo de seguridad. Además, la degradación que tiene la altura sobre el nivel del mar en un diseño de un
de las baterías de iones de litio se ve influenciada por la pack de baterías de ion litio que cuenta con refrigeración
temperatura del núcleo, el cual está más cerca del electrodo [3], forzada. La hipótesis que se plantea es que el efecto de la
de esta forma, éste parámetro proporcionará una referencia más altitud es suficientemente significativo como para justificar que
precisa para la estimación de la vida útil de la batería. las baterías deberían diseñarse orientadas a mercados
específicos en relación a la altitud de trabajo, lográndose con

EnerSOL 2016 1
esto diseños más económicos cuando se conoce que el uso de segundo circuito RC representan la dinámica de difusión del
estas estará orientado principalmente a ciudades cercanas a las litio en el electrolito. Las variables de estado del modelo
costa. eléctrico son SOC, V1 y V2. La corriente I es la variable de
entrada, mientras que la variable de salida VT es el voltaje en
II. METODOLOGIA terminales el que se define por (4).
En esta sección se presenta un modelo eléctrico y térmico (4)
de una celda de ion litio que permite estimar el desempeño de VOCV representa el voltaje de circuito abierto, mientras que RS
esta. La dinámica eléctrica se representa por un modelo de un representa la resistencia interna de la celda. VOCV sólo depende
circuito eléctrico consistente en una fuente de voltaje, una del estado de carga de la batería, sin embargo, los parámetros
resistencia y dos circuitos resistivos capacitivos en paralelo del circuito equivalente dependerán del estado de carga, la
(OCV-R-RC-RC). Para el fenómeno térmico se utiliza un temperatura, y la dirección de la corriente. En [6] se muestran
modelo de dos estados que captura la dinámica de la parámetros en función del SOC, temperatura y la dirección de
temperatura en el núcleo y en la superficie de la batería. La la corriente. Los valores de los coeficientes corresponden a la
dependencia entre ambos modelos se establece mediante la celda 26650 fabricada por la empresa A123 Systems.
generación de calor que se produce al circular corriente
eléctrica por la celda, obteniéndose así un modelo electro- La Temperatura de la celda de la batería es conducida por
térmico de la batería. la generación de calor Q definido en (5).

A. Modelo eléctrico (5)

El estado de carga de una celda “state of charge”, (SOC) se El calor generado en la celda está definido por el
define como la integración de la corriente que circula por la calentamiento por efecto Joule y la energía disipada en el
celda durante un periodo de tiempo, lo que queda definido por electrodo de mayor potencial.
(1):

(1) B. Modelo térmico

La distribución de temperatura radial dentro de una celda
La capacidad nominal de la celda, Cn [Ah] se puede obtener cilíndrica puede ser descrita por ecuaciones diferenciales
realizando ciclos de carga y descarga de acuerdo a las ordinarias, las cuales consideran la generación de calor y de
recomendaciones que entregue el fabricante. transferencia de esta. Un modelo térmico de dos estados
El comportamiento eléctrico de una batería se puede térmicos simplificado es el que define en (6) y (7),
representar por un circuito eléctrico equivalente como el que se
muestra en la Fig. 1. Este tipo de estructura es una buena (6)
elección para representar el modelo de la batería. Los dos
circuitos RC representan la constante de tiempo rápida y lenta (7)
para el establecimiento de la tensión en terminales de la celda
[6].
Donde [°C] y [°C] representan los estados de la
temperatura del centro y de superficie, respectivamente. Los
parámetros que caracterizan este modelo son los siguientes;
[°C] es la temperatura de entrada del aire y la generación de
calor del modelo eléctrico. [J/K] es la capacidad de calor en
el núcleo, [J/K] es la capacidad de calor en la superficie,
[K/W] es la resistencia de conducción entre el núcleo y la
superficie, [K/W] es resistencia de convección en torno a la
celda. La resistencia de convección depende de las
Fig. 1. Modelo del circuito equivalente eléctrico de la batería de ion-Litio condiciones de flujo, y puede ser modelado para diferentes
Las ecuaciones (2) y (3) representan el modelo eléctrico de tipos de fluidos. La Fig. 2 muestra una representación gráfica
la batería. del modelo de transferencia de calor.
Para efecto de interacción con las variables eléctricas de la
(2) celda, la temperatura utilizada por el modelo equivalente del
circuito es la media de las temperaturas del núcleo y de la
(3) superficie definido como Tm [°C] en (8).

Los estados V1 y V2 son los voltajes de los condensadores, (8)

Rs es la resistencia interna, los parámetros R1 [Ω], C1 [F]
corresponden al primer circuito RC y representan la dinámica
de difusión del litio en la fase sólida. R2 [Ω], C2 [F] en el

EnerSOL 2016 2
 D. Resistencia convectiva entre celdas
Cuando se tiene un arreglo de celdas un tema que debe ser
resuelto es la determinación de la resistencia de convección
entre celdas. Una estimación aproximada de este valor se
puede obtener usando el procedimiento establecido por
Zhukauskas [7]. Esta metodología fue diseñada para
intercambiadores de calor, pero debido a que la configuración
geométrica utilizada en intercambiadores de calor es similar a
la que se usa en arreglo de celdas circulares, también se puede
aplicar en baterías.
Fig. 1. Modelo Térmico de una Celda Zhukauskas (1972, 1987) proporcionó una amplia teoría
sobre el ordenamiento de tubos alineados y escalonados como
la transferencia de calor y caída de presión. La transferencia
C. Representación de un arreglo de celdas
de calor desde una celda depende de su ubicación dentro del
La transferencia de calor hacia o desde un banco de tubos empaquetamiento y dentro de rangos bajos del número de
en flujo cruzado es relevante para numerosas aplicaciones Reynolds. Los tubos en la primera fila tienen un
industriales, como la generación de vapor en una caldera o el comportamiento similar a las de las filas inmediatamente
enfriamiento en un serpentín de un acondicionador de aire. En superiores. Para números de Reynolds mayores, la turbulencia
la Fig. 4 se muestra un arreglo de celdas en forma lineal y del flujo conduce a una mayor transferencia de calor desde los
escalonada. Normalmente un fluido se mueve sobre los tubos. tubos interiores que de los de la primera fila. La transferencia
de calor se hace invariante con la ubicación del tubo después
de la tercera o cuarta fila en el régimen de flujo mixto, esto se
produce por encima del número de Reynolds ReDmax.
E. Efecto de la altitud en el diseño de baterías
Uno de los efectos más notorios en la pérdida de
rendimiento y capacidad de sistemas térmicos se produce por
el efecto de la altura, la que provoca que la densidad de un
fluido compresible y la presión atmosférica se vea disminuida
considerablemente acarreando una reducción de potencia para
Fig. 2. Arreglo geométrico en forma lineal y escalonada sistemas eléctricos y térmicos. Ante esto, los empaques de
baterías de ion litio que usan refrigeración forzada mediante
Las filas de los tubos de un banco están escalonadas o un fluido compresible, se ven afectados directamente
alineadas en la dirección de la velocidad del fluido. La perdiendo capacidad de refrigeración, y disminución del
configuración se caracteriza por el diámetro del tubo D y por rendimiento.
la separación transversal ST y la separación longitudinal SL Desde un punto de vista mecánico, la diferencia
medidas entre los centros de los tubos. Las condiciones del fundamental entre líquidos y gases consiste en que estos
flujo dentro del banco están dominadas por los efectos de últimos, dada su naturaleza de fluido compresible, pueden ser
separación de la capa limite y por las interacciones de estelas, comprimidos. Su volumen, por tanto, no es constante y
que a su vez influyen en la transferencia de calor por consiguientemente tampoco su densidad. Teniendo en cuenta
convección. el papel fundamental de esta magnitud física en la estática de
A diferencia del modelo de una celda individual, en este se fluidos, se comprende que el equilibrio de los gases haya de
requiere modelar la transferencia de calor por convección considerarse separadamente de los líquidos.
entre celdas vecinas. De (6) y (7) se aprecia que la temperatura Una forma simple de expresar el cambio de presión con la
interna no se vincula con la temperatura superficial de las altura esta expresada por la ecuación general de la
otras celdas, por lo tanto, se puede deducir que la trasferencia hidrostática, que generalmente se aplica en fluidos no
de calor entre celdas depende de las temperaturas superficiales compresibles, pero que también es aplicable para fluidos
entre celdas y el coeficiente de convección entre estas. De esta compresibles.
manera se tiene que para las celdas i, j la transferencia de calor (10)
entre celdas vecinas queda establecida por (9).
Dónde es la diferencia de Presión entre dos puntos, es
(9) la diferencia de alturas, la densidad del fluido y la
aceleración de la gravedad.
En general en los fluidos compresibles la densidad varía de
Donde Rci,j es la resistencia de convección entre celdas i y j, Ωi
manera casi proporcional a la presión. Por ello si la
es en conjunto de celdas vecinas a la celda i.
integración de (10) supone grandes cambios de presión, la
densidad debe ser considerada como variable. Para analizar

EnerSOL 2016 3
este tema, se consideran dos situaciones específicas en los ambiente . De estos resultados se puede determinar que el
procesos de integración. coeficiente de convección h es en promedio 19 % mayor a
En la primera, se considera un gas estático a temperatura nivel del mar que a 3.000 m.s.n.m. Por otra parte, la Fig. 8
constante. En este cuyo caso, (10) se combina con la ecuación muestra los resultados obtenidos cuando se considera el efecto
de estado del gas y se integra para obtener la presión y la de la temperatura ambiente. En este caso, la temperatura
densidad del gas en función de la altura sobre el nivel del mar. promedio a nivel del mar es 20 °C, mientras que a 3.000
Un caso más interesante, pero menos habitual consiste en m.s.n.m es 12 °C. Este análisis en equivalente a comparar las
obtener la presión y densidad del gas considerando que la condiciones de las ciudades de Antofagasta y Calama. Cuando
presión varía linealmente con la altura. En esta situación, la la temperatura es considerada por el modelo, el coeficiente de
presión y densidad quedarán expresadas en función de la convección h es en promedio 21% mayor a nivel del mar con
altura y de la temperatura ambiente. En (11) y (12) se respecto a 3.000 m.s.n.m. Esto implica que considerar la
muestran las ecuaciones de presión y densidad del gas en temperatura en el modelo produce una diferencia de 2 % con
función de estas dos variables. respecto a modelos que no la consideran.
(11)

(12)

Donde, es la presión a una altura H y temperatura en

[Pa], es la presión a nivel del mar, 101325 [Pa], es la
densidad del aire a nivel del mar, 1.2 [kg/m3]. es la
aceleración de gravedad, 9.8 [m/s2], el gradiente térmico
0.00650 [°K/m], es la constante de gases, que para el caso
del aire es 287 [ ], es la temperatura a nivel del mar [°K]
Fig. 6. Coeficiente de convección a nivel del mar
y H es la altura sobre el nivel del mar en [m]. La Fig. 5 muestra
la densidad en función de la altitud y la temperatura del
entorno,

Fig. 7. Coeficiente de convección a 3000 m.s.n.m. sin considerar la

temperatura

Fig. 5. Coeficiente de convección a nivel del mar

En la Fig. 5 se aprecia que la densidad disminuye cuando la

altitud crece. Asimismo, la densidad disminuye cuando la
temperatura aumenta.

III. RESULTADOS
En esta sección se realizan diferentes análisis en los cuales Fig. 8. Coeficiente de convección a 3000 m.s.n.m. considerando la
se explora el comportamiento de parámetros y variables temperatura
relevantes en el diseño de baterías cuando cambia la
temperatura del entorno y la altitud. La Fig. 6 muestra el Este efecto es aún más cuantificable cuando se compara la
coeficiente de convección en función de la separación vertical diferencia de presión entre la entrada al módulo de celdas y la
ST [m] parametrizado con distintas velocidades de salida. La diferencia de presión tiene un mayor impacto en el
refrigeración a nivel del mar. De la misma manera, en la Fig. 7 diseño del arreglo de celdas dado que define la capacidad del
se muestran los resultados para una altura de 3000 m.s.n.m. En ventilador afectando el rendimiento energético de la batería.
estas simulaciones no se ha considerado la temperatura La Fig. 9 presenta la caída de presión ∆P a nivel de mar en

EnerSOL 2016 4
función de la separación vertical ST, parametrizada a distintas
velocidades de refrigeración. Asimismo, la Fig. 10 muestra la
caída de presión a 4000 m.s.n.m. La caída de presión
disminuye conforme aumenta la altitud, de esta manera se
tiene, que a 4000 m.s.n.m esta disminuye en un 27% cuando la
velocidad de refrigeración es 5 m/s.

Fig. 11. Temperatura de pack de 40 celdas a nivel del mar

Fig. 9. Caída de presión a nivel del mar en función de la separación vertical

ST.

Fig. 12. Temperatura de pack de 40 celdas a 4000 m.s.n.m

Otro ejercicio que se realizó para cuantificar el efecto de la

altitud fue calcular la potencia del ventilador. Para esto se
comparó la potencia que se requiere para el ventilador a nivel
del mar con la requerida a 4000 metros de manera que la
temperatura sea la misma. De esta manera se tiene que la
temperatura máxima a nivel del mar con una velocidad de 0.1
m/s es 30.7 °C. Para lograr una temperatura similar a 4000 m
se requiere una velocidad cercana a los 0.5 m/s. La potencia
Fig. 10. Caída de presión a 4000 m.s.n.m. función de la separación vertical requerida para ventilación en este caso es 87 veces la requerida
ST. a nivel del mar. Esto permite concluir que el efecto de la altitud
influye considerablemente en el consumo de energía utilizado
Otra variable que se puede analizar en el diseño de un pack en refrigeración.
de batería es la temperatura. Para cuantificar esta variable, se
realizaron diferentes simulaciones con distintas IV. CONCLUSIONES
configuraciones. La Fig. 11 muestra los resultados obtenidos a
En esta investigación se ha analizado el efecto que tiene la
nivel del mar para una configuración escalonada consistente altitud y temperatura ambiente en el diseño de un banco de
de 40 celdas conectadas eléctricamente en serie y organizadas baterías de ion litio. Los resultados obtenidos permiten concluir
en 10 columnas y 4 filas. Se utilizó una celda cilíndrica que la altitud tiene un mayor efecto que la temperatura
LiFePO4, tipo ANR 26650 fabricada por A123 System. Las ambiente en el diseño. En lo que se refiere al coeficiente de
características de la celda y parámetros de simulación son: convección h, este es en promedio 21% mayor a nivel del mar
Voc = 3.35 (V), 2.5 Ah, D = 2.6 cm, SL = 3.3 cm, ST = 3.4 cm, que con respecto a 3.000 m.s.n.m. Por otro lado, la temperatura
T∞ = 24 °C, V = 0.1 m/s. Como corriente de carga se aplica un ambiente sólo tiene un efecto de un 2% en la variación de la
pulso de 3C en descarga durante 15 min, y 0.5C en carga densidad. En lo que se relaciona con la temperatura en el
durante 15 min. El tiempo de simulación es 60 min. núcleo de la celda, está experimenta un aumento del orden de 1
Asimismo, la Fig. 12 muestra los resultados de la simulación °C entre nivel del mar y 4000 m.s.n.m. Finalmente, la energía
para una altitud de 4000 m.s.n.m. requerida para mantener la misma temperatura de la celda que
En las Fig. 11 y 12 se aprecia que el efecto de la altitud se produce a nivel del mar requiere un consumo en ventilación
produce un incremento de temperatura cercano a 1 °C. De esta 87 veces mayor. Este efecto, permite justificar la hipótesis de
manera se puede concluir que la altura produce un efecto no que es conveniente realizar diseños personalizados en base de
significativo en la temperatura. altitud, como una forma de lograr diseños más económicos.

EnerSOL 2016 5
 AGRADECIMIENTOS [3] P. Antonucci, and V. Antonucci, “Electrochemical energy storage” in
Electrochemical energy storage in the emerging era of smart grids,
Los autores agradecen el generoso apoyo financiero Rosario Carbone, Eds. 2011, pp. 2.
entregado por los proyectos CONICYT/ FONDAP/15110019 [4] N. Baba, H. Yoshida, M. Nagaoka, C. Okuda, S. Kawauchi, “Numerical
“Solar Energy Research Center” SERC-Chile e INNOVA simulation of thermal behavior of lithium-ion secondary batteries using
the enhanced single particle model,” Journal of Power Sources, vol.
CORFO 12IDL2-16296. 252, pp. 214-228, April 2014.
[5] M. Lapuerta, O. Armas, J. R. Agudelo, y C. A. Sánchez, “Study of the
Altitude Effect on Internal Combustion Engine Operation. Part 1:
REFERENCES Performance”, Información Tecnológica, vol. 17 N° 5, pp. 21-30, 2006.
[1] S. Al Hallaj, H. Maleki, J. S. Hong, and J. R. Selman, “Thermal [6] H.Perez and J. Sieguel, "Parameterization and validation of an integrated
modeling and design considerations of lithium-ion batteries,” Journal of electro-thermal cylindrical LPF battery model," ASME 2012 5th Annual
Power Sources, vol. 83, pp. 1-8, October 1999. Dynamic Systems and Control Conference, 2012.
[2] S. Al-Hallaj, J. R. Selman, “Thermal modeling of secondary lithium [7] F. Incropera, and D. De Witt, Fundamentals of heat and mass transfer,
batteries for electric vehicle/hybrid electric vehicle applications,” John Wiley and Sons Inc.,New York, NY, 1985.
Journal of Power Sources, vol. 110, pp. 341-3488, August 2002.

EnerSOL 2016 6

View publication stats