GPU-native extension of a high-order compressible solver to multicomponent reacting flows

De moderne wetenschap heeft een punt bereikt waarop de krachtigste instrumenten voor ontdekking niet alleen telescopen zijn die naar buiten kijken of microscopen die naar binnen kijken, maar enorme digitale laboratoria die draaien op supercomputers. Binnen het vakgebied van de fluïdummechanica, dat bestudeert hoe gassen en vloeistoffen bewegen, bestaat er een specifieke uitdaging die lang als een van de moeilijkste problemen wordt beschouwd: het simuleren van een vuur. In tegen tegenstelling tot een eenvoudige wind die over een vleugel blaast, omvat een vuur een chaotische mix van kolkend turbulentie, razendsnelle chemische explosies en de beweging van vele verschillende soorten gasmoleculen tegelijkertijd. Om dit te begrijpen, moet men zich voorstellen dat men probeert het pad van elke individuele molecuul zuurstof, waterstof en stikstof te volgen terwijl ze botsen, verbranden en wervelen, terwijl de lucht zelf wordt samengedrukt en verhit. Dit met perfecte nauwkeurigheid doen vereist een niveau van rekenkracht dat ooit onmogelijk was, maar de opkomst van gespecialiseerde computerchips genaamd GPU's heeft een nieuwe deur geopend. Deze chips, oorspronkelijk ontworpen voor het renderen van videogames, zijn nu de motoren die de volgende generatie wetenschappelijke simulaties aandrijven, waardoor onderzoekers complexe fysische verschijnselen met ongekende details kunnen modelleren.

Een team onderzoekers van de Universiteit van Rome en de Europese Ruimtevaartorganisatie heeft een belangrijke stap voorwaarts gezet in dit veld door een geavanceerde computercode te upgraden om deze vurige, meerkomponentige stromingen aan te kunnen. Hun werk richt zich op een solver genaamd STREAmS-2.0, een programma ontworpen om te simuleren hoe samendrukbare gassen zich gedragen wanneer ze met hoge snelheden bewegen. Hoewel de vorige versie van deze software uitstekend was in het afhandelen van enkelvoudige soorten gas, kon deze nog niet omgaan met de complexiteit van een echte chemische reactie waarbij verschillende gassen mengen, verbranden en veranderen in nieuwe stoffen. De onderzoekers hebben dit instrument nu uitgebreid om mengsels van gassen te kunnen verwerken die met elkaar reageren, een capaciteit die essentieel is voor het ontwerpen van schonere en efficiëntere raketmotoren en straalmotoren. Ze hebben niet alleen een nieuwe functie toegevoegd; ze hebben de wiskundige motor van de software herbouwd om ervoor te zorgen dat deze efficiënt kan draaien op de enorme arrays van grafische processoren die in de krachtigste supercomputers van vandaag te vinden zijn.

De kern van hun prestatie ligt in de manier waarop ze de computer hebben geleerd om met de chaos van een reactief mengsel om te gaan. In een niet-reactieve stroming volgt de computer simpelweg hoe lucht beweegt. In een reactieve stroming moet de computer ook bijhouden hoe snel verschillende chemische reacties plaatsvinden, wat miljoenen keren sneller kan gebeuren dan de lucht zelf beweegt. Dit creëert een moeilijk wiskundig probleem waarbij de computer gelijktijdig moet oplossen voor twee zeer verschillende snelheden. Het team loste dit op door het probleem in twee delen te splitsen: één deel dat de beweging van het gas afhandelt en een ander dat de chemie afhandelt. Vervolgens gebruikten ze een specifieke, zeer efficiënte methode om het chemische deel op te lossen, wat een complexe reeks vergelijkingen omvat die beschrijven hoe moleculen uiteenvallen en weer samenklonteren. Door bepaalde chemische eigenschappen vooraf te berekenen en op een manier op te slaan die de grafische processoren direct kunnen ontsluiten, hebben ze ervoor gezorgd dat de computer geen tijd verspilt aan het wachten op gegevens. Deze aanpak zorgt ervoor dat de simulatie stabiel en accuraat blijft, zelfs wanneer de chemische reacties extreem snel en onvoorspelbaar zijn.

Om te bewijzen dat hun nieuwe instrument werkt, hebben de onderzoekers het onderworpen aan een reeks strenge tests, waarbij ze de resultaten vergeleken met bekende wiskundige oplossingen en gegevens van andere betrouwbare simulaties. Ze testten alles van eenvoudige eendimensionale problemen, zoals een schokgolf die door een buis met gas reist, tot complexe driedimensionale scenario's met kolkende vortexen van brandstof en lucht. In één test simuleerden ze een vlamfront dat zich door een mengsel van waterstof en zuurstof beweegt, om te controleren of de computer de temperatuur en de snelheid van de vlam nauwkeurig kon voorspellen. In een andere test modelleerden ze een turbulente wolk van gas waarin verschillende soorten stoffen mengden en reageerden. In elk geval kwamen de resultaten van hun nieuwe software met opmerkelijke precisie overeen met de referentiedata. De simulaties lieten zien dat de software in staat was om de vorming van schokgolven, de diffusie van verschillende gassen en de vrijgekomen hitte door chemische reacties correct te vangen, wat bevestigt dat de onderliggende fysica correct in code is vertaald.

Naast nauwkeurigheid maakte het team zich evenveel zorgen over hoe snel de software op echte hardware kon draaien. Ze testten de solver op de Leonardo supercomputer, een van de krachtigste machines ter wereld, met gebruik van maximaal 512 high-performance grafische processoren. Ze ontdekten dat de software bijna perfect schaalt: naarmate ze meer processoren toevoegden, nam de tijd die nodig was om een simulatie te voltooien recht evenredig af, wat betekent dat de machine efficiënter werd naarmط de groei ervan. Dit is een zeldzame en waardevolle eigenschap in high-performance computing, waar het toevoegen van meer processoren vaak leidt tot verminderde rendementen door communicatievertragingen tussen hen. De onderzoekers ontdekten dat de software deze hoge efficiëntie behoudt zolang elke processor genoeg werk heeft, specifiek vereist dat elke chip ten minste 10% van zijn beschikbare geheugen gebruikt om productief te blijven. Wanneer de werklast te klein is, besteden de processoren meer tijd aan het wachten op gegevens dan aan het uitvoeren van berekeningen, maar wanneer het probleem groot genoeg is, werkt het systeem met bijna perfecte efficiëntie.

De studie onthulde ook waar de computer de meeste tijd doorbrengt. In simulaties zonder chemische reacties is de processor bezig met het berekenen van hoe het gas beweegt en hoe drukgolven zich voortplanten. Echter, wanneer de chemische reacties worden ingeschakeld, wordt een aanzienlijk deel van de rekenkracht — ongeveer 40% van de totale tijd — besteed aan het berekenen van de chemische brontermen, oftewel de snelheden waarmee moleculen worden gecreëerd of vernietigd. Dit benadrukt de enorme computationele kosten van het modelleren van verbranding, maar de onderzoekers toonden aan dat hun methode deze last effectief afhandelt zonder de gehele simulatie te vertragen. Het werk bevestigt dat het mogelijk is om hoogwaardige numerieke methoden, die een hoge nauwkeurigheid bieden, te combineren met complexe chemische kinetica en de massale parallelle kracht van moderne grafische processoren.

Vooruitblikkend willen de onderzoekers deze nieuwe capaciteit gebruiken om een laboratoriumschaal jetflame in volledige detail te simuleren, waarbij ze hun resultaten direct vergelijken met fysieke experimenten om het model verder te valideren. Deze volgende stap zal dienen als de definitieve test voor het vermogen van de software om real-world verbrandingsverschijnselen te voorspellen. Het succes van dit project toont aan dat de combinatie van geavanceerde wiskunde, chemische fysica en moderne hardware een volwassenheid bereikt waarbij zij de moeilijkste problemen in de luchtvaarttechniek kunnen aanpakken. Door een instrument te bieden dat in staat is om reactieve stromingen op een massale schaal accuraat en efficiënt te simuleren, legt dit werk de basis voor het ontwerpen van voortstuwingssystemen die niet alleen krachtiger, maar ook schoner en efficiënter zijn, wat bijdraagt aan de voortdurende inspanningen om de milieu-impact van de luchtvaart en ruimtevaart te verminderen.