Реализуйте типовую водородную модель стека топливного элемента
Simscape / Электрический / Специализированные Энергосистемы / Источники
Блок stack Топливного элемента реализует типовую модель, параметрированную, чтобы представлять большинство популярных типов стеков топливного элемента, питаемых водородом и воздухом.
Блок представляет две версии модели стека: упрощенная модель и подробная модель. Можно переключиться между этими двумя моделями путем выбора уровня в маске под Model detail level в диалоговом окне блока.
Эта модель основана на эквивалентной схеме стека топливного элемента, показанного ниже:
Упрощенная модель представляет конкретный стек топливного элемента, действующий при номинальных условиях температуры и давления. Параметры эквивалентной схемы могут быть изменены на основе кривой поляризации, полученной из таблицы данных производителя. Только необходимо ввести в маске значение напряжения в 0 и 1 А, номинала и максимальных рабочих точек, для параметров, которые будут вычислены. Диод используется, чтобы предотвратить поток отрицательного тока в стек. Типичная кривая поляризации состоит из трех областей:
Первая область представляет падение напряжения активации из-за замедления химических реакций, происходящих в поверхностях электрода. В зависимости от температурного и рабочего давления, типа электрода и используемого катализатора, эта область более или менее широка. Вторая область представляет резистивные потери из-за внутреннего сопротивления стека топливного элемента. Наконец, третья область представляет потери массового транспорта, следующие из изменения в концентрации реагентов, когда топливо используется.
Подробная модель представляет конкретный стек топливного элемента, когда параметры, такие как давления, температура, составы и скорости потока топлива и воздуха варьируются. Можно выбрать, который параметры варьироваться на Signal variation разделяют на области на диалоговом окне блока. Эти изменения влияют на напряжение разомкнутой цепи (Eoc), обменный ток (i0) и наклон Tafel, (э). Эок, i0 и A изменяется можно следующим образом:
где:
R = 8.3145 J / (молекулярная масса K)
F = 96 485 А s/mol
z = Количество движущихся электронов
En = Напряжение Нернста, которое является напряжением термодинамики ячеек и зависит от температур и парциальных давлений реагентов и продуктов в стеке (V)
α = Коэффициент передачи заряда, который зависит от типа электродов и используемых катализаторов
PH2 = Парциальное давление водорода в стеке (Па)
PO2 = Парциальное давление кислорода в стеке (Па)
k = Константа Больцманна = 1,38 × 10–23 J/K
h = Константа Планка = 6,626 × 10–34 J s
Δv = Фактор объема барьера активации (m3). Размер барьера активации (ΔG) вычисляется, принимая Δv = 1 м3.
ΔG = Размер барьера активации, который зависит от типа электрода и катализатора, использовал (J/mol)
T = Температура операции (K)
Kc = Напряжение, постоянное при номинальном условии операции
Эквивалентная схема эквивалентна для упрощенной модели, за исключением того, что параметры Eoc, i0 и Α должны быть обновлены онлайн как показано ниже:
Уровни преобразования (использования) водорода (UfH2) и кислород (UfO2) определяются в Блоке A можно следующим образом:
где
Pfuel = Абсолютное давление предоставления топлива (банкомат)
Пара = Абсолютное давление предоставления воздуха (банкомат)
Vlpm (топливо) = Топливная скорость потока (l/min)
Vlpm (воздух) = уровень Воздушного потока (l/min)
x = Процент водорода в топливе (%)
y = Процент кислорода в окислителе (%)
N = Количество ячеек
Постоянные 60000 прибывают из преобразования из скорости потока жидкости литра/min, используемой в модели к m3/s (1 литр/min = 1/60000 m3/s).
Парциальные давления и Напряжение Нернста определяются в Блоке B можно следующим образом:
и
где
PH2O = Парциальное давление водяного пара в стеке (банкомат)
w = Процент водяного пара в окислителе (%)
От парциальных давлений газов и Напряжения Нернста, могут быть вычислены новые значения напряжения разомкнутой цепи (Eoc) и обменного тока (i0).
Блок C вычисляет новое значение наклона Tafel (A).
Параметры α, ΔG, и Kc вычисляются на основе кривой поляризации при номинальных условиях операции наряду с некоторыми дополнительными параметрами, такими как КПД низкой теплоты сгорания (LHV) стека, состав топлива и воздуха, предоставляют давления и температуры. Они могут быть легко получены из таблицы данных производителя.
Номинальные уровни преобразования газов вычисляются можно следующим образом:
где:
ηnom = Номинальный КПД LHV стека (%).
Δh0(H2O (газ)) = 241,83 × 103 J/mol.
Vnom = Номинальное напряжение (V).
Inom = Номинальный ток (A).
Vlpm (воздух) имя = Номинальный уровень воздушного потока (l/min).
Pairnom = Номинальное абсолютное давление подачи воздуха (Па).
Tnom = Номинальная рабочая температура (K).
От этих уровней преобразования могут быть выведены номинальные парциальные давления газов и Напряжения Нернста. С Eoc, i0 и известный Α и принимая, что стек действует на постоянных уровнях преобразования или использований при номинальном условии, может быть определен α, ΔG, и Kc.
Если нет никакого топлива или воздуха во входе стека, он принят, что стек действует по фиксированной процентной ставке преобразования газов (номинальный уровень преобразования), то есть, предоставление газов настроено согласно току так, чтобы они всегда предоставлялись только необходимым, чем немного более стеком при любой загрузке.
Максимальный ток, который может поставить стек, ограничивается максимальными скоростями потока топлива и воздуха, который может быть достигнут. Кроме того максимальный ток, напряжение, выведенное стеком, уменьшается резко, когда более актуальный чертится.
Движущие силы топливного элемента представлены, если вы задаете время отклика и параметры для динамики потока (пиковое использование и соответствующее отклонение от номинала напряжения) на панели Fuel Cell Dynamics на диалоговом окне.
Время отклика (Td) 95% используются, чтобы смоделировать "заряд, двойной слой" явление из-за наращивания бросается на интерфейс электрода/электролита. Это влияет только на напряжение активации (NAln (ifc/i0)) как показано на эквивалентных схемах.
Пиковое использование (UfO2 (пик)) и соответствующее отклонение от номинала напряжения (Vu) используется, чтобы смоделировать эффект кислородного истощения (из-за воздушной задержки компрессора) на выходном напряжении ячейки. Напряжение Нернста изменяется должное с этой целью можно следующим образом:
где
K = постоянное отклонение от номинала напряжения
UfO2 (имя) = номинальное кислородное использование
K определяется можно следующим образом:
Текущий шаг и тесты прерывания должны быть сделаны на действительном стеке представлять точностью свою динамику. Рисунок ниже показывает ответ стека от этих тестов и обязательных параметров (Td, UfO2 (пик) и Vu).
Время отклика (Td) зависит от самого стека топливного элемента и обычно дается в таблице данных. Параметры для динамики потока (UfO2 (пик) и Vu) зависят от динамики внешнего оборудования (компрессор, регулятор и загрузки), и им не предоставляют производители, когда их значения меняются в зависимости от пользовательского приложения. Для симуляции примите значения UfO2 (пик) между 60% к 70% и Vu между 2-5% напряжения номинала стека.
Вот процедура, чтобы извлечь параметры из таблицы данных производителя стека топливного элемента. В данном примере таблица данных NetStack PS6 от NetStack используется:
Номинальная мощность стека составляет 6 кВт, и номинальное напряжение составляет 45 В. Следующие подробные параметры выведены из таблицы данных.
Напряжение на уровне 0 А и 1 А [Eoc, V1] = [65, 63]
Номинальная рабочая точка [Inom, Vnom] = [133.3, 45]
Максимальная рабочая точка [Iend, Продайте] = [225, 37]
Номинальный КПД стека (ηnom) = 55%
Рабочая температура = 65 0C
Номинальное давление предоставления [H2, Воздух] = [1.5 1]
Если данное давление относительно атмосферного давления, добавьте 1 панель, чтобы получить абсолютное давление.
Номинальный состав (%) [H2, O2, H2O (Воздух)] = [99.999, 21, 1]
Если воздух используется в качестве окислителя, примите 21% O2 и 1% H2O в случае, если их проценты не заданы.
Количество ячеек
Если не заданный, оцените его от формул ниже:
В этом случае,
Номинальный уровень воздушного потока
Если максимальный уровень воздушного потока дан, номинальная скорость потока жидкости может быть вычислена, приняв постоянное кислородное использование при любой загрузке. Ток, чертивший ячейкой, линейно зависим на уровне воздушного потока, и номинальной скоростью потока жидкости дают:
В этом случае,
В случае, если никакая информация не дана, примите уровень преобразования кислорода, чтобы быть 50% (когда это обычно имеет место для большинства стеков топливного элемента), и используйте формулы ниже, чтобы определить номинальный уровень воздушного потока.
Время отклика топливного элемента = 10 с
Примечание
Параметры [Eoc, V1], [Inom, Vnom], и [Iend, Продают], являются аппроксимированными и зависят от точности точек, полученных из кривой поляризации. Чем выше точность этих параметров, тем более закрыто симулированное напряжение стека к кривой таблицы данных. Инструмент, названный ScanIt (от amsterchem), может использоваться, чтобы извлечь точные значения из кривых таблицы данных.
Вышеупомянутыми параметрами кривая поляризации стека, действующего на фиксируемой номинальной ставке преобразования газов, закрывается для кривых таблицы данных как показано ниже: синяя пунктирная линия показывает, что симулированное напряжение стека и зеленая пунктирная линия показывают симулированную степень стека.
Выше максимального тока скорость потока газов, вводящих стек, максимальна, и напряжение стека уменьшается резко, когда более актуальный чертится.
Это предположения об этом блоке:
Газы идеальны
Стек питается водородом и воздухом
Стек оборудован системой охлаждения, которая обеспечивает температуру в катоде, и анод выходит устойчивый и равный температуре стека
Стек оборудован системой управления водными ресурсами, чтобы обеспечить влажность в ячейке на соответствующем уровне при любой загрузке
Падения напряжения ячейки происходят из-за кинетики реакции и транспорта заряда, когда большинство топливных элементов не действует в области массового транспорта
Перепады давления через каналы потока незначительны
Сопротивление ячейки является постоянным при любом условии операции
Это ограничения этого блока:
Движущие силы химической реакции, вызванные изменениями парциального давления химических разновидностей в ячейке, не рассматриваются
Выходная мощность стека ограничивается топливом и предоставленными уровнями воздушного потока
Воздействие температуры и влажность мембраны на внутреннем сопротивлении не рассматриваются
Поток газов или воды через мембрану не рассматривается
[1] Njoya, S. M. О. Тремблей и L.-A. Dessaint. Типовая модель топливного элемента для симуляции автомобилей на топливных элементах. Степень транспортного средства и Конференция по Движению, 2009, VPPC ’09, IEEE® . 7-10 сентября 2009, стр 1722–29.
[2] Motapon, S.N., О. Тремблей и L.-A. Dessaint. “Разработка типовой модели топливного элемента: приложение к симуляции автомобиля на топливных элементах”. Международный J. Силовой электроники. Издание 4, № 6, 2012, стр 505–22.