Wildfires Quasi-Implicit Alternative-Direction Simulations using Isogeometric Finite Element Method

Probleemstelling

Het modelleren van bosbranden is een complex probleem dat verschillende benaderingen vereist, variërend van empirische modellen tot volledig fysische simulaties. Bestaande methoden hebben echter significante beperkingen:

• Empirische modellen (zoals FARSITE) zijn rekenefficiënt maar missen een consistente fysische representatie van thermodynamische en transportprocessen, wat hun bruikbaarheid beperkt bij sterke koppeling tussen branddynamiek en atmosferische krachten.
• Volledig fysische CFD-modellen (zoals FIRETEC en WFDS) lossen de Navier-Stokes-vergelijkingen op en bieden hoge nauwkeurigheid, maar zijn extreem rekenintensief. Ze vereisen kleine tijdstappen en hoge-performance computing (HPC), waardoor ze ongeschikt zijn voor grootschalige voorspellingen of ensemble-simulaties.
• Data-gedreven en AI-benaderingen missen vaak robuustheid bij scherpe fronten en stijve reactietermen.

Er is een behoefte aan een model dat een goede balans biedt tussen fysische nauwkeurigheid en rekenkundige schaalbaarheid, specifiek voor het simuleren van de temperatuurverdeling en brandverspreiding op grote schaal.

Methodologie

De auteurs ontwikkelen een kwasi-impliciete brandverspreidingsmodel gebaseerd op de Isogeometrische Analyse (IGA) met B-spline elementen. De kern van de methode omvat:

1. Fysisch Model:

   • Het model lost een gereduceerde energiebalansvergelijking op voor het temperatuurscalarveld $T(x,t)$.
   • De vergelijking omvat warmtegeneratie (verbranding), transport (windgedreven advektie) en verliezen (geleiding, straling, convectie).
   • Niet-lineaire termen (zoals stralingsverlies $T^4$ en verbrandingsreacties) worden expliciet behandeld op tussentijdse tijdstippen om volledige niet-lineariteit in het lineaire systeem te vermijden.
   • Brandstofbeschikbaarheid wordt gemodelleerd via een ruimtelijk variërende coëfficiënt ($\eta(x)$) afgeleid van satellietbeelden (NDVI-data), waarbij groene kanalen worden geëxtraheerd om brandbaar materiaal te identificeren.

2. Ruimtelijke Discretisatie (Isogeometrie):

   • In plaats van standaard $C^0$ Lagrange-elementen, worden tensor-product B-splines gebruikt.
   • Dit biedt hogere continuïteit ($C^{p-1}$), wat leidt tot betere benaderingskwaliteit per vrijheidsgraad en verminderde numerieke dispersie, vooral belangrijk voor diffusie-gedomineerde operatoren en gladde temperatuurvelden.

3. Tijdsintegratie en Operator-Splitsing:

   • Om de hoge rekenkosten van impliciete methoden in 2D/3D te omzeilen, worden Alternating Direction Implicit (ADI) methoden gebruikt.
   • Twee specifieke schema's worden onderzocht:
     • Peaceman-Rachford: Een splitsing in twee richtingen (x en y) met een impliciete behandeling van de lineaire transport-diffusie-operatoren.
     • Strang-splitsing gecombineerd met Crank-Nicolson: Een symmetrische splitsing die tweede-orde nauwkeurigheid in de tijd biedt voor transport-gedomineerde regimes.
   • Door de tensor-product structuur van de B-splines, kunnen de multidimensionale systemen worden gereduceerd tot een reeks één-dimensionale lineaire systemen. De matrices hebben een Kroonek-product structuur, wat de oplossing zeer efficiënt maakt.

4. Software-implementatie:

   • De code is geschreven in C++ en is open-source beschikbaar (IGA-ADS bibliotheek).
   • De implementatie is parallelisabel voor shared-memory systemen (multithreading).

Belangrijkste Bijdragen

• Nieuwe Numerieke Schemata: De eerste toepassing van kwasi-impliciete, alternatieve-richtingsmethoden binnen een isogeometrische raamwerk specifiek voor bosbranden.
• Analyse van Stabiliteit en Nauwkeurigheid: Een gedetailleerde analyse toont aan dat het kwasi-impliciete schema 10 keer nauwkeuriger is dan volledig expliciete methoden voor kleine tijdstappen, terwijl het de stabiliteit behoudt.
• Open Source Software: De publicatie van een volledig werkende, open-source C++ code met installatie- en gebruikshandleidingen.
• Realistische Casestudies: Toepassing op twee recente, grote natuurrampen met vergelijking tegen satellietbeelden en meetgegevens.

Resultaten

De auteurs presenteren drie soorten numerieke experimenten:

1. Convergentie- en Stabiliteitstesten:

   • Met behulp van een "manufactured solution" (een bekende exacte oplossing) wordt de convergentie geanalyseerd.
   • De kwasi-impliciete methoden (Peaceman-Rachford en Strang/Crank-Nicolson) tonen een hogere nauwkeurigheid en betere stabiliteit dan expliciete methoden, zelfs bij grotere tijdstappen.
   • De methoden behouden een lineaire rekencomplexiteit van $O(N)$ per tijdstap, waarbij $N$ het totale aantal vrijheidsgraden is.

2. Vergelijking met FARSITE:

   • De simulaties tonen een vergelijkbare evolutie van het brandfront als het gevestigde FARSITE-model, maar met een hogere fysische consistentie.

3. Realistische Simulaties:

   • Chili (2024): Simulatie van de bosbranden in de Valparaíso-regio (Viña del Mar). Het model voorspelde de verspreiding van de brand nauwkeurig in overeenstemming met satellietbeelden en foto's van het verbrande gebied, hoewel het initiële brandpunt in de simulatie iets groter was dan in de realiteit.
   • Spanje (2019): Simulatie van de brand op Gran Canaria (Las Palmas). Het model reproduceerde de driehoekige vorm van het brandfront, veroorzaakt door de windrichting, en toonde hoe de brand zou zijn verspreid zonder brandbestrijding.

4. Schaalbaarheid:

   • Tests op een Intel Core Ultra 5 processor (14 kernen) tonen goede strong scaling (snelheidsverhoging) en weak scaling.
   • Voor kubieke B-splines werd een snelheidsverhoging (speedup) van ongeveer 7.5 bereikt bij 14 kernen.
   • Het is mogelijk om grote simulaties (bijv. $1500 \times 1500$ mesh) op een lokale workstation uit te voeren binnen enkele minuten tot een uur, afhankelijk van de polynoomgraad.

Betekenis en Conclusie

Dit onderzoek biedt een krachtig alternatief voor traditionele bosbrandmodellen. Door de combinatie van Isogeometrische Analyse (voor hoge nauwkeurigheid en gladheid) en Alternating Direction Implicit (ADI) methoden (voor rekenkundige efficiëntie), kunnen realistische, fysisch gebaseerde simulaties worden uitgevoerd op schaalbare hardware.

De methode overbrugt de kloof tussen de lage nauwkeurigheid van empirische modellen en de extreme rekenkosten van volledige CFD-simulaties. De beschikbaarheid van de open-source code maakt het mogelijk voor onderzoekers en beleidsmakers om gedetailleerde scenario's te analyseren, wat cruciaal is voor risicobeheer, evacuatieplanning en het begrijpen van de dynamiek van bosbranden in een veranderend klimaat. De auteurs wijzen er ook op dat toekomstig werk gericht zal zijn op de integratie van Variational Physics-Informed Neural Networks (VPINN).