Preconditioning for a Cahn-Hilliard-Navier-Stokes model for morphology formation in organic solar cells

Technische Samenvatting: Preconditionering voor een Cahn-Hilliard-Navier-Stokes-model voor morfologievorming in organische zonnecellen

1. Probleemstelling

De fabricage van organische fotovoltaïsche cellen (OPV's) is afhankelijk van het drogen van een oplosmiddelgebaseerd mengsel dat een polymeer (elektrondonor), een non-fullereen acceptor (NFA) en een oplosmiddel bevat. De prestaties van het resulterende apparaat zijn kritisch afhankelijk van de morfologie van de dunne film, die ontstaat door een complex samenspel van spontane fase-scheiding en oplosmiddelverdamping. Het experimenteel optimaliseren van de productiecondities om een gewenste morfologie te bereiken is tijdrovend en vaak onhaalbaar vanwege de gevoeligheid van het proces voor materiaaleigenschappen en procesparameters.

De wiskundige uitdaging ligt in het modelleren van deze morfologie-evolutie. Het systeem omvat:

• Macroscopische vloeistofdynamica: De beweging van het vloeistofmengsel.
• Microscopische fase-scheiding: De scheiding van donor, acceptor en oplosmiddelcomponenten.
• Verdamping: De overgang van oplosmiddel van vloeistof naar gas en de diffusie hiervan in de omringende lucht.

Deze fysische processen zijn gekoppeld, wat leidt tot een systeem van partiële differentiaalvergelijkingen (PDE's) die, wanneer gediscretiseerd, resulteren in grootschalige, slecht geconditioneerde lineaire en niet-lineaire systemen. Het efficiënt oplossen van deze systemen, met name met betrekking tot veranderingen in discretisatieparameters (roostergrootte, tijdstap), is de primaire numerieke hindernis die in dit werk wordt aangepakt.

2. Methodologie

Wiskundige Modellering
De auteurs stellen een continuümmodel voor dat drie hoofdonderdelen koppelt:

1. Navier-Stokes vergelijkingen: Beschrijven de macroscopische snelheid ($v$) en druk ($p$) van het vloeistofmengsel.
2. Cahn-Hilliard vergelijkingen: Beschrijven de evolutie van volumefracties ($\phi_i$) voor het polymeer, de NFA, het oplosmiddel en de lucht, gedreven door chemische potentialen ($\mu_i$) afgeleid van een vrije energiefunctionaal.
3. Allen-Cahn vergelijking: Beschrijft de evolutie van een damporde-parameter ($\Phi_v$) om de faseverandering van het oplosmiddel van vloeistof naar gas te modelleren. Deze vergelijking is niet-massa-behoudend, wat de afvoer van oplosmiddel naar de omgeving reflecteert.

Het model incorporeert een vrije energiefunctionaal gebaseerd op de Flory-Huggins theorie, die interpoleert tussen vloeistof- en gasfasen met behulp van een gladde functie van de damporde-parameter. Randvoorwaarden omvatten no-slip voor de snelheid en specifieke flux-condities aan de bovenste rand om de uitstroom van oplosmiddel en de instroom van lucht te modelleren.

Discretisatie

• Ruimtelijk: De vergelijkingen worden gediscretiseerd met behulp van de Eindige Elementen Methode (FEM). De auteurs gebruiken een gesplitste formulering om aan de continuïteitseisen van de basisfuncties te voldoen. Taylor-Hood elementen worden gebruikt voor de Navier-Stokes vergelijkingen (snelheid van graad $\ell+1$, druk van graad $\ell$), terwijl continue Lagrange-elementen worden gebruikt voor de fasevelden en chemische potentialen.
• Temporeel: Er wordt een semi-impliciet tijdstap-schema gebruikt. Lineaire termen worden impliciet behandeld, terwijl niet-lineaire termen voortvloeiend uit het vrije energiepotentiaal expliciet worden behandeld.

Preconditioneringsstrategie
De kernbijdrage van het artikel is de ontwikkeling van een robuuste gepreconditioneerde iteratieve oplosser voor de resulterende saddle-point systemen.

• Block Preconditioning: De auteurs construeren blok-diagonaal preconditioners die de saddle-point structuur van de gediscretiseerde vergelijkingen benutten.
• Schur Complement Benadering: Voor de Cahn-Hilliard en Navier-Stokes blokken benadert de preconditioner de Schur complement.
  • Voor Cahn-Hilliard wordt de Schur complement benaderd door een product van matrices (matching approach) om directe inversie van een som van grote matrices te vermijden.
  • Voor Navier-Stokes wordt een pressure-convection-diffusion (PCD) preconditioner gebruikt, gebaseerd op de commutator van de operatoren.
• Lineaire Solvers: De diagonale blokken en Schur complement benaderingen worden geïnverteerd met behulp van Algebraic Multigrid (AMG) methoden (specifiek Ruge-Stüben en Smoothed Aggregation). Het totale systeem wordt opgelost met de Generalized Minimal Residual (GMRES) methode.

3. Belangrijkste Bijdragen

• Nieuwe Modelformulering: Het artikel presenteert een gekoppeld model dat vloeistofdynamica, fase-scheiding en oplosmiddelverdamping koppelt, specifiek voor de context van morfologievorming in organische zonnecellen, waarbij eerdere Cahn-Hilliard-Navier-Stokes modellen worden uitgebreid met verdampingseffecten via een Allen-Cahn vergelijking.
• Robuuste Preconditionering: De auteurs ontwerpen en implementeren een preconditioneringstechniek die ervoor zorgt dat het aantal GMRES-iteraties klein en onafhankelijk blijft van de roostergrootte ($h$) en de tijdstap ($\tau$). Dit "parameter-robuuste" gedrag is cruciand voor de schaalbaarheid van de solver in grootschalige 3D-simulaties.
• Verbeteringen in Numerieke Stabiliteit: Om numerieke instabiliteiten aan te pakken die voortkomen uit de logaritmische vrije energiefunctionaal, stellen de auteurs een polynomiale benadering voor (geïnspireerd door double-well potentialen) die de gewenste fysische kenmerken (fase-scheidingsminima) behoudt terwijl de numerieke stabiliteit wordt verbeterd.

4. Numerieke Resultaten

De auteurs valideren het model en de solver via uitgebreide numerieke experimenten in 1D, 2D en 3D domeinen:

• Solver Prestaties: In alle geteste scenario's (variërende roostergroottes van $50^2$ tot $400^2$ in 2D, en $10^3$ tot $50^3$ in 3D, en variërende tijdstappen), bleven de GMRES-iteratieaantallen laag en stabiel. De runtime schaalde lineair met de vrijheidsgraden, wat de efficiëntie van de preconditioner bevestigt.
• Fysisch Gedrag:
  • Verdamping: Simulaties van pure oplosmiddelverdamping toonden de verwachte uitputting van het oplosmiddel en de vorming van een damplaag.
  • Fase-scheiding: In het volledige model (polymeer, NFA, oplosmiddel) legden de simulaties succesvol de fase-scheiding van de actieve lagen vast. Het gebruik van de polynomiale vrije energie-benadering maakte stabiele simulaties van complexe morfologieën in 3D mogelijk, die voorheen gevoelig waren voor numerieke instabiliteiten met de logaritmische potentiaal.
• 3D Morfologie: Het artikel demonstreert succesvol de vorming van 3D-morfologieën, waarbij de evolutie van polymeer- en NFA-domeinen tijdens het droogproces wordt getoond.

5. Betekenis en Claims

Het artikel beweert dat nauwkeurige modellering van morfologievorming een onmisbaar instrument is voor het optimaliseren van OPV-productie, aangezien experimentele trial-and-error onpraktisch is. De betekenis van dit werk ligt in:

1. Mogelijk maken van 3D-simulaties: Door een robuuste preconditioneringstechniek te bieden, stellen de auteurs de simulatie van het morfologievormingsproces in drie dimensies mogelijk, wat een noodzakelijke stap is naar realistische device modellering.
2. Parameter Robuustheid: De voorgestelde solver is robuust ten opzichte van discretisatieparameters, waardoor deze geschikt is voor adaptieve roosterverfijning en variërende tijdstap-strategieën die vereist zijn voor complexe fysische simulaties.
3. Kwalitatieve Vastlegging: Het model legt kwalitatief het gedrag van filmmorfologie tijdens het drogen vast, inclusief de interactie tussen fase-scheiding en verdamping.

De auteurs concluderen dat hoewel het model de essentiële fysica vangt, de mogelijkheid om deze simulaties in 3D uit te voeren een krachtig hulpmiddel biedt voor het onderzoeken van geoptimaliseerde productiecondities. Ze beweren niet het inverse probleem van het vinden van optimale condities te hebben opgelost, maar bieden het numerieke kader dat nodig is om het forward probleem van morfologie-evolutie te bestuderen.