Optimized tandem catalyst patterning for CO$_2$ reduction flow reactors

Deze studie toont aan dat het integreren van continuümtransportmodellering met op adjoint gebaseerde ontwerpoptimalisatie de prestaties van stroomreactoren voor CO₂-reductie aanzienlijk verbetert door strategische patronen van Ag- en Cu-katalysatoren aan te brengen om de ethyleenstroomdichtheid te maximaliseren, met name bij hoge spanningen en met meer patroongedeelten.

De kern

Stel je voor dat je probeert de perfecte taart te bakken, maar je keuken heeft een vreemde regel: je kunt niet alle ingrediënten tegelijk mengen. In plaats daarvan heb je twee aparte stations.

• Station A (Zilver): Dit station is uitstekend in het omzetten van rauwe bloem (Koolstofdioxide) in deeg (Koolmonoxide).
• Station B (Koper): Dit station is fantastisch in het omzetten van dat deeg in een heerlijke taart (Ethyleen, een waardevolle chemische stof).

Het probleem? Als je Station A ver van Station B plaatst, wordt het deeg weggeblazen door de wind (het stromende water in de reactor) voordat het Station B kan bereiken. Of, als je te veel Station A en te weinig Station B hebt, eindig je met een hoop deeg en geen taart.

Dit artikel gaat over het uitvinden van de perfecte lay-out voor deze twee stations om de meeste taart mogelijk te maken.

Het Grote Idee: "Tandem-katalyse"

De onderzoekers bestuderen een proces genaamd tandem-katalyse. Denk hierbij aan een assemblagelijn.

1. Zilver (Ag) fungeert als de eerste werknemer, die CO₂ omzet in CO.
2. Koper (Cu) fungeert als de tweede werknemer, die die CO oppakt en omzet in hoogwaardige producten zoals ethyleen (een bouwsteen voor plastics en brandstoffen).

In een traditionele opstelling kunnen deze werknemers gemengd zijn of geplaatst in grote, gescheiden blokken. De onderzoekers wilden weten: Als we de elektrode opsplitsen in vele kleine, afwisselende stroken van Zilver en Koper, en we de lengte van elke strook kunnen aanpassen, wat is dan het beste patroon om de meeste taart te krijgen?

Het Experiment: Een Digitale "Afregeling"-knop

In plaats van fysieke reactoren te bouwen en duizenden verschillende patronen te proberen (wat jaren zou duren), bouwde het team een computersimulatie.

Ze creëerden een digitale "stroomreactor" waar vloeistof over een vlak oppervlak stroomt. Ze gebruikten een slim computeralgoritme (zoals een super-geavanceerde GPS) om miljoenen verschillende patronen te testen. De computer zou:

1. Een patroon proberen (bijvoorbeeld een lange Zilver-strook, gevolgd door een korte Koper-strook).
2. Kijken hoeveel "taart" (ethyleen) er gemaakt werd.
3. De lengtes van de stroken lichtjes aanpassen.
4. Dit keer en keer herhalen totdat de absolute beste rangschikking gevonden was.

Wat Ze Vonden

De computer vond dat het "perfecte" patroon sterk afhankelijk is van hoe hard je het systeem duwt (de spanning) en hoe snel de vloeistof stroomt.

1. Het "Sterke Duw"-scenario (Hoge Spanning):
Toen ze het systeem hard duwden (met een sterke elektrische spanning), was het beste ontwerp om veel, heel veel kleine stroken te hebben (tot wel 12 secties) in plaats van slechts twee grote.

• Het Resultaat: Dit geoptimaliseerde patroon produceerde tot 65% meer ethyleen dan een eenvoudig, niet-geoptimaliseerd ontwerp.
• Waarom? Bij hoge snelheden beweegt de vloeistof snel. Als het Koper-segment te lang is, wordt het "deeg" (CO) helemaal aan het begin van de strook opgebruikt, en zit de rest van de Koper-strook inactief (een "dode zone"). Door de stroken korter en talrijker te maken, wordt vers deeg constant aangeleverd bij de Koper-werknemers, waardoor ze de hele tijd bezig blijven.

2. Het "Zachte Duw"-scenario (Lage Spanning):
Toen de duw zwakker was, zag het beste patroon er anders uit. Het bevoordeelde een zeer lange eerste Zilver-strook om een enorme hoop deeg te maken, gevolgd door een zeer lange laatste Koper-strook om alles op te eten, met kleine, snel wisselende stroken in het midden.

3. De Stroomsnelheid Is Belangrijk:

• Snelle Stroom: Als het water razendsnel voorbij stroomt, moet de reactie zeer sterk zijn (hoge spanning) om te voorkomen dat het deeg wegspoelt.
• Trage Stroom: Als het water traag is, heeft het deeg tijd om te bezinken, maar moet je oppassen dat je niet op verse ingrediënten zit.

De Geheime Ingrediënt: "Dode Zones" Vermijden

De belangrijkste reden waarom de geoptimaliseerde patronen zo goed werkten, is dat ze "dode zones" elimineerden.

Stel je een transportband voor waar de eerste paar werknemers druk bezig zijn, maar de laatste paar werknemers nietsdoen staan omdat de onderdelen op zijn. In de oude ontwerpen hadden de Koper-secties vaak deze dode zones aan het einde waar de CO op was.

De geoptimaliseerde ontwerpen van de computer herschikten de stroken zodat het "deeg" (CO) gelijkmatig werd verdeeld. Het zorgde ervoor dat elke centimeter van het Koper-oppervlak genoeg deeg had om aan te werken, waardoor de productie van het eindproduct werd gemaximaliseerd.

Samenvatting

Dit artikel is een "proof of concept". Het bouwde geen fysieke fabriek, maar het bewees dat het gebruik van wiskunde en computers om de lay-out van katalysatoren te ontwerpen aanzienlijk kan bijdragen aan hoe goed we CO₂ omzetten in bruikbare chemicaliën.

• Het Probleem: CO₂-reductie is lastig; tussenproducten gaan verloren of worden verspild.
• De Oplossing: Gebruik een computer om het perfecte patroon van afwisselende Zilver- en Koper-stroken te vinden.
• De Opbrengst: Door simpelweg de vorm van het katalysatoroppervlak te veranderen (niet de chemicaliën zelf), konden ze in hun simulatie de productie met tot 65% verhogen.

Het is als beseffen dat als je de meubels in een kamer herschikt, je veel sneller kunt bewegen, zelfs als je geen nieuwe meubels koopt.

Technische samenvatting

Technische Samenvatting: Geoptimaliseerde Tandemkatalysatorpatronen voor CO2-reductiestroomreactoren

Probleemstelling
Elektrochemische CO2-reductie (CO2R) biedt een route om overtollige hernieuwbare energie om te zetten in duurzame brandstoffen en chemicaliën, terwijl tegelijkertijd broeikasgasemissies worden verminderd. De commerciële schaalvergroting wordt echter belemmerd door prestatiebeperkingen in koper (Cu)-katalysatoren, specifiek hoge overpotentialen en een beperkte selectiviteit voor hoogwaardige meerkoolstofproducten (zoals ethyleen, ethanol). Deze beperkingen vloeien vaak voort uit een inefficiënt beheer van het transport van intermediaire species. Bij traditionele tandemkatalyse, waarbij een eerste katalysator (bijv. Ag) CO2 omzet in CO en een tweede katalysator (bijv. Cu) CO omzet in hoogwaardige producten, kunnen intermediaire species zoals CO verloren gaan door advectie of diffusie voordat ze de tweede katalysator bereiken. Hoewel eerdere studies diverse tandemkatalysatorconfiguraties hebben onderzocht (zoals verspreide deeltjes, afwisselende planaire secties), hebben deze grotendeels vertrouwd op voorwaartse analyse – het voorstellen en evalueren van specifieke ontwerpen – in plaats van het systematisch optimaliseren van de ruimtelijke rangschikking van katalysatoren om de prestaties te maximaliseren.

Methodologie
Deze studie presenteert een proof-of-concept raamwerk dat continuümtransportmodellering integreert met op adjoint gebaseerde ontwerpoptimalisatie binnen een vereenvoudigde tweedimensionale (2D) stroomreactoropstelling.

• Systeemopstelling: Het computationele domein modelleert een schuifstroom van een waterige bicarbonaat-elektrolyt (500 mM KHCO3, 34 mM CO2) over een actief elektrodeoppervlak. De elektrode bestaat uit afwisselende Ag- en Cu-secties. Ag-secties faciliteren de CO2 → CO-reactie, terwijl Cu-secties de omzetting van CO naar ethyleen (C2H4), ethanol (C2H6O), methaan (CH4) en waterstof (H2) faciliteren. CO is niet aanwezig in de inlaatstroom; het wordt gegenereerd door Ag en verbruikt door Cu.
• Besturende Vergelijkingen: Het model lost stationaire Nernst-Planck-vergelijkingen op voor species-transport (CO2, H+, OH-, HCO3-, CO3 2-, K+, CO, H2 en koolwaterstoffen), gekoppeld aan bulk-bicarbonaatbufferreacties. Oppervlaktekinetiek wordt gemodelleerd met Tafel-expressies voor zowel Ag- als Cu-katalysatoren. De vloeistofstroom wordt benaderd als een schuifstroom met een gedefinieerde schuifsnelheid, afgeleid van het Péclet-getal.
• Optimalisatieformulering: Het probleem wordt geformuleerd als een door partiële differentiaalvergelijkingen (PDE) beperkte optimalisatie. De ontwerpvariabelen zijn de lengtes van $N$ afwisselende Ag/Cu-secties ($l_j$). De doelstelling is het maximaliseren van de ruimtelijk gemiddelde ethyleenstroomdichtheid ($i_{C2H4}$).
• Algoritme: De optimalisatie maakt gebruik van de adjoint-methode om sensitiviteiten van de doelstelling ten opzichte van de katalysatorsectielengtes te berekenen. Een L-BFGS-B-kwasi-Newton-algoritme verfijnt de sectielengtes iteratief om de prestaties te maximaliseren. De studie onderzoekt de effecten van aangelegde spanning ($U_{app}$), elektrolytstroomsnelheid en het aantal patronensecties ($N$).

Belangrijkste Bijdragen

1. Ontwerpoptimalisatiekader: Het artikel introduceert een systematische optimalisatiebenadering voor tandemkatalysatorpatronen, die voorbijgaat aan voorwaartse analyse door ruimtelijke rangschikkingen van katalysatoren iteratief te verfijnen op basis van transport- en reactiebeperkingen.
2. Prestatieverbetering via Patronen: De studie toont aan dat het optimaliseren van de ruimtelijke verdeling van Ag- en Cu-secties de ethyleenproductie aanzienlijk verbetert, met name wanneer het aantal secties ($N$) wordt verhoogd en onder omstandigheden van sterke oppervlaktereacties (meer negatieve aangelegde spanningen).
3. Mechanistisch Inzicht: Door analyse van concentratievelden verduidelijkt het werk hoe geoptimaliseerde patronen "dode zones" (gebieden met lage CO-concentratie op Cu-oppervlakken) mitigeren en het benutten van de intermediaire CO-species verbeteren.

Resultaten

• Invloed van Aantal Secties ($N$): Het verhogen van het aantal secties van $N=2$ naar $N=12$ verhoogt de ethyleenstroomdichtheid aanzienlijk. Voor een aangelegde spanning van $-1,7$ V vs. SHE bereikt het geoptimaliseerde $N=12$-ontwerp tot een 65% stijging in ethyleenstroomdichtheid in vergelijking met een niet-geoptimaliseerd $N=2$-ontwerp met gelijke lengtes.
• Optimalisatie versus Aantal Secties: De prestatiewinst wordt ontleed in bijdragen van het verhogen van $N$ en het optimaliseren van de patronlengtes. Bij minder negatieve spanningen ($-1,35$ V) is de optimalisatie van het patroon de dominante factor. Bij meer negatieve spanningen ($-1,7$ V) is het verhogen van het aantal secties de primaire drijvende kracht van verbetering, hoewel optimalisatie nog steeds aanzienlijke winsten oplevert.
• Patroonkenmerken:
  • Bij $-1,35$ V (zwakkere reactie) vertonen geoptimaliseerde patronen voor hoge $N$ lange initiële Ag- en finale Cu-secties met korte, snel afwisselende intermediaire secties.
  • Bij $-1,7$ V (sterkere reactie) en lage stroomsnelheden geven geoptimaliseerde patronen de voorkeur aan langere Ag-secties en kortere Cu-secties om CO-productie te maximaliseren en CO-uitputting op Cu-oppervlakken te minimaliseren.
  • Bij $-1,7$ V en hoge stroomsnelheden blijven de geoptimaliseerde patronen dicht bij configuraties met gelijke lengtes, wat suggereert dat het initiële ontwerp met gelijke lengtes onder deze omstandigheden van hoog massatransport al bijna optimaal is.
• Afhankelijkheid van Stroomsnelheid: Het effect van de stroomsnelheid op de prestaties is spanningsafhankelijk. Bij $-1,35$ V verminderen hogere stroomsnelheden het CO-benuttingspercentage en de ethyleenproductie door het wegspoelen van zwak geproduceerd CO. Omgekeerd verhogen hogere stroomsnelheden bij $-1,7$ V de ethyleenproductie door reactanten effectiever aan te vullen, aangezien de sterke oppervlaktereacties hoge lokale CO-concentraties handhaven.

Betekenis en Claims
De auteurs claimen dat dit werk een proof of concept biedt voor het gebruik van continuümmodellering in combinatie met ontwerpoptimalisatie om de ruimtelijke rangschikking van katalysatoren in CO2R-elektrolyzers te sturen. De studie stelt vast dat:

• Het optimale reactorontwerp sterk afhankelijk is van bedrijfsomstandigheden (spanning, stroomsnelheid) en de mate van patronering ($N$).
• Het primaire mechanisme voor prestatieverbetering in geoptimaliseerde tandemkatalysatoren het minimaliseren is van zones met lage reactantconcentratie op het oppervlak van de secundaire katalysator, waardoor het benutten van het sleutelintermediaire product (CO) wordt gemaximaliseerd.
• Hoewel de huidige 2D-opstelling een vereenvoudiging is die niet alle complexiteiten van industriële gasdiffusie-elektroden (GDE's) vastlegt, biedt de methodologie een fundamentele stap naar het optimaliseren van realistischere, complexere elektrolyzergeometrieën.
• De integratie van high-performance computing met onderbouwde continuümmodellering dient als een krachtig hulpmiddel om de ontwerpruimte voor experimentele validatie te verkleinen en de efficiëntie van elektrochemische processen te verbeteren.

Het artikel concludeert dat toekomstig werk deze optimalisatiecapaciteiten moet uitbreiden naar 3D-omgevingen (bijv. GDE's), andere doelstellingsfuncties (selectiviteit, energie-efficiëntie) moet onderzoeken en bedrijfsomstandigheden naast katalysatorpatronering moet optimaliseren.