CFD Stromings- en Thermische Analyse van een Verwarmde Buisreactor met Statische Mixer

De Uitdaging

Een chemisch verwerkingsbedrijf had een diepgaand inzicht nodig in het stromingsgedrag en de warmteoverdrachtsprestaties binnenin een compacte buisreactor. De reactor, met een lengte van 240 mm en een binnendiameter van 20 mm, bevat een spiraalvormige statische mixerschroef en een centraal geplaatst thermokoppel. De schroef induceert een wervelend stromingspad dat de menging van het procesfluïdum bevordert terwijl het door de verwarmde buis stroomt.

De opdrachtgever moest de interne stromingspatronen en snelheden valideren en bepalen hoe effectief de reactor warmte overdraagt van de extern verwarmde wand naar het fluïdum onder verschillende bedrijfscondities. Vier scenario's werden gedefinieerd, waarbij twee volumedebieten (6 ml/min en 19 ml/min) werden gecombineerd met twee viscositeitswaarden van het fluïdum (0,3 mPa·s en 1 mPa·s). Het begrijpen van het samenspel tussen debiet, viscositeit en warmteoverdracht is essentieel voor het optimaliseren van de prestaties en het opschalingspotentieel van de reactor.

Onze Aanpak

We bouwden een volledig 3D conjugaat warmteoverdrachtsmodel van de reactor in Ansys Fluent, inclusief zowel het fluïdumdomein als de vaste componenten: de buitenbuis uit Inconel 625 en de spiraalvormige schroef uit S316L roestvast staal. Een conjugaat warmteoverdrachtsaanpak was essentieel omdat de thermische interactie tussen de verwarmde buiswand, de metalen schroef en het stromende fluïdum de werkelijke temperatuurverdeling binnenin de reactor bepaalt.

Temperatuurafhankelijke fluïdumeigenschappen

Het procesfluïdum vertoont sterk temperatuurafhankelijke thermische eigenschappen. Tussen 60 °C en 600 °C daalt de dichtheid van 1.002 kg/m³ tot 504 kg/m³, terwijl de warmtecapaciteit bijna verdubbelt van 1.678 J/kg·K tot 3.497 J/kg·K. De thermische geleidbaarheid daalt met meer dan de helft over hetzelfde bereik. Elke eigenschap werd geïmplementeerd als een polynomiale benadering afgeleid uit gemeten datapunten, zodat de simulatie getrouw weergeeft hoe het gedrag van het fluïdum verandert naarmate het opwarmt over de reactorlengte.

Kalibratie van de warmte-invoer aan de hand van metingen

De reactor wordt verwarmd door een thermische straler van 2.300 W die de buis omgeeft. Door de lage emissiviteit van Inconel 625 (onder 0,1) en onbekende warmteverliezen naar de omgeving bereikt echter slechts een fractie van dit vermogen het fluïdum. De opdrachtgever verstrekte een gemeten fluïdumtemperatuur van 510 °C op 200 mm van de inlaat voor het scenario met 6 ml/min en een inlaattemperatuur van 200 °C. Via een energiebalans werd de werkelijke warmteopname door het fluïdum geschat op circa 73 W — slechts 3% van het nominale vermogen van de straler. We stemden de aangebrachte warmtebelasting op de buiswand af op 96,6 W om dit meetpunt te matchen, en pasten vervolgens dezelfde warmtebelasting toe op alle vier scenario's voor een consistente vergelijking.

Laminaire steady-state solver

Gezien de zeer lage debieten en de kleine buisdiameter liggen de Reynoldsgetallen in alle scenario's ruim onder de turbulente overgangsdrempel. Een steady-state laminaire solver was daarom passend, waarmee onnodige rekentijd wordt vermeden zonder in te boeten op nauwkeurigheid. Het polyhedrische mesh werd verfijnd rondom de oppervlakken van de spiraalvormige schroef om het complexe wervelende stromingspad en de thermische grenslagen op de schroef- en buiswanden nauwkeurig te capteren.

Resultaten

De stroomlijnenvisualisaties bevestigden dat de spiraalvormige schroef succesvol een goed gedefinieerde wervelstroming induceert doorheen de gehele reactor, zonder significante dode zones of stromingsstagnatie. De spiralerende beweging brengt het fluïdum herhaaldelijk in contact met de verwarmde buiswand, wat precies het menggedrag is dat de statische mixer beoogt te bereiken.

Effect van debiet op uitlaattemperatuur

Zoals verwacht heeft het debiet een doorslaggevende invloed op de uitlaattemperatuur. Bij het lagere debiet van 6 ml/min heeft het fluïdum meer verblijftijd in de verwarmde zone, waardoor het een gemiddelde uitlaattemperatuur van circa 531 °C bereikt — ruim boven de 510 °C die op 200 mm van de inlaat werd gemeten. Bij het hogere debiet van 19 ml/min beperkt de kortere verblijftijd de opwarming tot een gemiddelde uitlaattemperatuur van circa 323 °C. Ondanks dit temperatuurverschil is de totale warmteopname door het fluïdum marginaal hoger bij het hogere debiet (circa 81 W tegenover 80 W), wat wijst op een licht verbeterde convectieve warmteoverdracht bij de hogere snelheid.

Effect van viscositeit op warmteoverdracht

Een van de meest significante bevindingen was dat de viscositeit van het fluïdum vrijwel geen effect heeft op de warmteoverdrachtsefficiëntie. De temperatuurprofielen en opgenomen warmtewaarden waren nagenoeg identiek tussen de scenario's met 0,3 mPa·s en 1 mPa·s bij hetzelfde debiet. Viscositeit beïnvloedt voornamelijk het drukverlies over de reactor en daarmee het vereiste pompvermogen, maar beide waarden blijven zeer klein en verwaarloosbaar in de totale energiebalans.

De temperatuurplots op de dwarsdoorsneden bij 0, 120, 200 en 240 mm van de reactorinlaat onthulden hoe het thermische veld evolueert over de reactorlengte. Bij de inlaat bestaat er een sterke radiale temperatuurgradiënt tussen de hete buiswand en de koelere kernstroming. Ter hoogte van de uitlaat heeft de statische mixer de warmte herverdeeld over de doorsnede, waardoor de radiale niet-uniformiteit aanzienlijk is verminderd. Dit bevestigt dat de spiraalvormige schroef effectief is in het bevorderen van radiaal warmtetransport, een kritieke factor voor het bereiken van een consistente productkwaliteit in chemische verwerkingstoepassingen.

Geleverde Meerwaarde

Deze CFD-analyse verschafte de opdrachtgever gevalideerd, kwantitatief inzicht in de thermisch-hydraulische prestaties van de reactor over het volledige bedrijfsbereik. Het gekalibreerde model, verankerd aan een fysieke temperatuurmeting, geeft vertrouwen dat de voorspelde temperaturen en stromingspatronen representatief zijn voor de werkelijke omstandigheden.

De kernbevinding dat viscositeit geen invloed heeft op de warmteoverdrachtsefficiëntie vereenvoudigt het procesontwerp van de opdrachtgever: men kan de viscositeit van het fluïdum aanpassen aan andere procesvereisten zonder de thermische prestaties van de reactor te beïnvloeden. De gedetailleerde stroomlijnen- en contourvisualisaties bieden een helder inzicht in de mengkwaliteit die door de spiraalvormige schroef wordt bereikt, wat direct relevant is voor toekomstige ontwerpoptimalisatie — bijvoorbeeld het aanpassen van de schroefspoed, reactorlengte of verwarmingsvermogen voor verschillende doeltemperaturen en doorvoersnelheden.

Door CFD-stromingssimulatie te combineren met conjugaat warmteoverdrachtsanalyse leverden we een uitgebreide prestatiekarakterisatie in één simulatiecampagne — een veel efficiëntere en informatievere aanpak dan iteratief fysiek prototyping en testen.