CFD Thermocline-Analyse van een Thermische Buffertank

De Uitdaging

Een grootschalige thermische buffertank (warmteaccumulatortank) werd ontworpen als onderdeel van een stadsverwarmingsnetwerk in Nederland. De drukvaste cilindrische tank, met een interne diameter van 14 m en halfronde uiteinden, bevat circa 2.920 m³ water bij een werkdruk van 1,99 bar. Warm water van 120 °C en koud water van 70 °C worden in- en afgevoerd via twee cirkelvormige diffusors (Ø 4,5 m, hoogte 0,4 m) die nabij de boven- en onderkant van de tank zijn gepositioneerd.

De centrale engineeringvraag was of de thermocline — de dunne overgangszone die de warme en koude waterlagen scheidt — stabiel zou blijven tijdens het laden en ontladen bij een volumedebiet van 634 m³/h. Een bijzonder aandachtspunt was de nabijheid van de onderste diffusor tot de gebogen tankwand: de stroming zou een Coanda-effect kunnen veroorzaken, met grootschalige wervels die de thermische stratificatie zouden kunnen verstoren en de opslagefficiëntie van de tank zouden verminderen.

Onze Aanpak

We modelleerden het volledige fluïdedomein van de tank in Ansys Fluent met een halfsymmetrisch model om de transiënte rekentijden beheersbaar te houden. De mesh bestond uit polyedrische elementen met inflatielagen langs de tankwanden en diffusoroppervlakken, om het stromingsgedrag nabij de wand nauwkeurig op te lossen. Een dichtheidsgebaseerde solver werd geselecteerd om de door drijfkracht aangedreven stromingen correct vast te leggen, die ontstaan door de temperatuurafhankelijke dichtheidsverschillen over de thermocline.

Fase 1 — Studie van de diffusorpositie

In de eerste fase werden stationaire analyses uitgevoerd voor twee posities van de onderste diffusor: het oorspronkelijke ontwerp op 628 mm van de tankwand en een alternatief op 300 mm. Beide scenario's werden eerst doorgerekend met een uniforme watertemperatuur (zonder thermocline) en vervolgens met de werkelijke temperatuurverdeling, zodat we het effect van dichtheidsverschillen op het stromingspatroon konden isoleren. Snelheidscontourplots en stroomljnvisualisaties werden geëxtraheerd om wervelstructuren en het Coanda-effect te identificeren.

Fase 2 — Transiënt thermocline-gedrag

De tweede fase richtte zich op het dynamische gedrag van de thermocline zelf. Volledig transiënte geconjugeerde thermische analyses werden uitgevoerd voor twee bedrijfsscenario's: koud in (koud water via de onderste diffusor, warm water afgevoerd bovenaan) en warm in (warm water via de bovenste diffusor, koud water afgevoerd onderaan). De initiële thermocline-dikte was vastgesteld op 0,8 m, gedefinieerd als de zone waar de temperatuur meer dan 1% afwijkt van de bulk warm- of koudwatertemperatuur. 80 meetpunten verticaal verdeeld tussen de diffusors volgden het temperatuurprofiel over de tijd.

Simulatieduur

Elke transiënte analyse liep tot de thermocline de tegenoverliggende diffusor bereikte — circa 8.800 seconden voor het koud-in scenario en 7.800 seconden voor warm-in, overeenkomend met één volledige volumewissel tussen de diffusors.

Resultaten

Diffusorpositionering heeft minimale invloed

De stationaire analyses met uniforme watertemperatuur onthulden grootschalige wervels die de gehele tank besloegen voor beide diffusorposities, aangedreven door het Coanda-effect langs de gebogen bodemwand. Wanneer echter de werkelijke temperatuurverdeling werd toegepast (zonder de energievergelijking op te lossen, zodat de thermocline gefixeerd bleef), veranderde het stromingspatroon ingrijpend: de dichtheidsgradiënt over de thermocline fungeerde als een natuurlijke barrière, die de wervels beperkte tot het onderste gedeelte ruim onder de stratificatiezone. Het verschil tussen de 628 mm en 300 mm configuraties was verwaarloosbaar — beide leverden nagenoeg identieke stromingsbegrenzing op.

De thermocline blijft stabiel

De volledig transiënte analyses bevestigden dat de thermocline gedurende beide laadscenario's intact bleef. De grootschalige wervels die in de stationaire oplossingen werden waargenomen, verdwenen volledig in de transiënte simulaties — een belangrijke bevinding die de beperkingen van stationaire aannames voor dit type dynamisch thermisch opslagprobleem onderstreepte. In werkelijkheid verplaatst de thermocline zich continu, en een stationaire aanpak, die uitgaat van een stilstaand stromingsveld, overschat de sterkte van de recirculatie.

In het koud-in scenario groeide de thermocline-dikte geleidelijk van 0,8 m tot circa 2 m over de volledige cyclus van 8.800 s. In het warm-in scenario was de groei iets groter en bereikte circa 2,5 m over 7.800 s. Deze snellere verdikking trad op wanneer de thermocline de halfronde bodem van de tank binnentrad, waar de afnemende dwarsdoorsnede de overgangszone zijdelings samendrukte.

Opgeleverde Waarde

Dit project toonde het cruciale belang aan van het kiezen van de juiste simulatieaanpak. Een stationaire CFD-analyse alleen zou problematische grootschalige wervels door de gehele tank hebben voorspeld — wat mogelijk zou hebben geleid tot onnodige en kostbare ontwerpaanpassingen. De transiënte simulatie vertelde een fundamenteel ander en veel realistischer verhaal: de dichtheidsgestuurde stratificatie creëert een natuurlijk dempend effect dat menging onderdrukt en de thermocline intact houdt.

Het engineeringteam ontving duidelijk, kwantitatief bewijs dat de oorspronkelijke diffusorpositionering aanvaardbaar was. De analyse toonde aan dat geen geometrische herontwerp van de onderste diffusor nodig was, dat de thermische stratificatie stabiel blijft onder zowel laad- als ontlaadcondities, dat de halfronde tankgeometrie enige extra thermocline-verdikking veroorzaakt — een voorspelbaar effect dat meegenomen kan worden in de operationele regelstrategie, en dat de transiënte animatievisualisaties een effectief communicatiemiddel vormen voor projectbetrokkenen die niet vertrouwd zijn met CFD.

Door computational fluid dynamics te combineren met thermische analyse leverden we het vertrouwen dat nodig was om het buffertankontwerp zonder aanpassingen voort te zetten, wat het project zowel tijd als budget bespaarde.