Reduced-order turbulent flow solver to simulate streamwise periodic fins with iso-thermal walls

Dit artikel introduceert en verifieert een verminderde-orde turbulente stromingsoplosser in SU2 voor het efficiënt simuleren van warmtewisselaars met isotherme wanden en stroomwaartse periodieke vinnen, waarbij de resultaten nauwkeurig overeenkomen met die van volledige fin-array-simulaties.

De kern

De "Tijdmachine" voor Warmtewisselaars: Een Simpele Uitleg

Stel je voor dat je een gigantische, ingewikkelde machine bouwt om warmte te verplaatsen, zoals een radiator in een auto of een koelkast. Deze machines bestaan uit duizenden kleine, herhalende onderdelen (zoals vinnen of pinnen) waar lucht of water doorheen stroomt.

Om te weten of deze machine goed werkt, moeten ingenieurs simuleren hoe de lucht stroomt en hoe warmte wordt uitgewisseld. Maar hier zit het probleem: als je de hele machine in detail op de computer narekent, duurt het een hele dag (of langer) om één berekening te maken. Dat is te lang en te duur als je wilt experimenteren met verschillende ontwerpen.

De Oplossing: De "Oneindige Tunnel"

De auteurs van dit artikel (Nitish Anand en zijn collega's) hebben een slimme truc bedacht, een soort "tijdmachine" voor hun computer. In plaats van de hele machine te simuleren, simuleren ze slechts één klein stukje (één vin en de lucht eromheen).

Hoe werkt dat?
Stel je voor dat je door een tunnel loopt die aan beide kanten oneindig doorgaat. Als je in het midden van die tunnel loopt, maakt het niet uit of je 1 meter of 100 meter verder bent; de situatie is precies hetzelfde. De luchtstroom en de temperatuur herhalen zich.

In de echte wereld is dat niet zo: de lucht wordt warmer naarmate hij verder de machine in gaat. De auteurs hebben een wiskundige "magische formule" (een bronterm) bedacht die de computer vertelt: "Hé, doe alsof deze luchtstroom al 100 keer is geweest, maar pas de temperatuur en druk netjes aan alsof we in een oneindige tunnel zitten."

De Twee Werelden: Rustig en Turbulent

Deze truc was al bekend voor rustige luchtstromen (zoals een kalme rivier). Maar wat als de lucht wild gaat stromen, met wervelingen en chaos (zoals een storm)? Dat is turbulente stroming.

Tot nu toe wisten wetenschappers niet hoe ze die "magische formule" moesten aanpassen voor die wilde, turbulente lucht. Dit artikel is het eerste dat die formule heeft afgeleid en in de computercode (SU2) heeft gezet.

Het Experiment: De Proef op de Som

Om te bewijzen dat hun nieuwe formule werkt, hebben ze een test gedaan:

1. Ze hebben een hele rij van 11 vinnen op de computer gezet (de "oude, dure manier").
2. Ze hebben daarna slechts één vin gesimuleerd met hun nieuwe "tijdmachine-formule" (de "snelle, slimme manier").

Ze hebben dit gedaan voor twee situaties:

• Rustige lucht (Laminaire stroming): Hier bleek hun formule perfect te werken. De resultaten van de één vin waren identiek aan die van de hele rij.
• Wilde lucht (Turbulente stroming): Dit was de echte uitdaging. Ook hier bleek hun nieuwe formule perfect te kloppen! De computer zag geen verschil tussen de simpele één-ving-simulatie en de dure hele-rij-simulatie.

De Grootte van de Winst

Waarom is dit zo belangrijk?

• De dure manier (de hele rij simuleren) kostte ongeveer 1440 minuten (een dag).
• De slimme manier (één vin met hun formule) kostte slechts 30 minuten.

Dat is een tijdwinst van 48 keer!

Conclusie

Dit artikel is als het vinden van een snellere route naar je werk. Je komt op precies hetzelfde punt aan, maar je bent er 48 keer sneller. Hierdoor kunnen ingenieurs nu veel sneller nieuwe, betere warmtewisselaars ontwerpen en testen, zonder dat hun computer dagenlang hoeft te rekenen. Ze hebben de sleutel gevonden om de wiskunde van de "wilde lucht" te temmen in een simpele, snelle simulatie.

Technische samenvatting

Titel: Verminderde orde turbulent stromingsoplosser voor het simuleren van stroomwaartse periodieke vinnen met isotherme wanden

Auteurs: Nitish Anand, Praharsh Pai Raikar, Carlo De Servi
Publicatie: Preprint (Elsevier), maart 2026

---

1. Probleemstelling

Het thermohydraulische ontwerp van warmtewisselaars (HEX) is traditioneel gebaseerd op vereenvoudigde methoden (zoals $\epsilon$-NTU en LMTD) en empirische correlaties. Hoewel deze methoden standaard zijn, leiden ze vaak tot conservatieve en suboptimale ontwerpen vanwege onzekerheden in de schatting van warmteoverdrachtscoëfficiënten en wrijvingsfactoren.

Computational Fluid Dynamics (CFD) biedt een nauwkeuriger alternatief, maar het simuleren van complexe HEX-geometrieën met vele parallelle kanalen is computatierijk en kostbaar. Een veelgebruikte aanpak om deze kosten te drukken is het gebruik van stroomwaartse periodieke (SP) modellen. Deze modellen gaan ervan uit dat de stroming en temperatuurverdeling na een bepaalde afstand van de inlaat zelfgelijkend (self-similar) worden.

Het specifieke gat in de literatuur:
Hoewel SP-modellen goed gevestigd zijn voor laminaire stroming en wanden met een voorgeschreven warmtestroom (heat-flux), ontbreekt er tot nu toe een wiskundige formulering en implementatie voor turbulente stromingen in kanalen met isotherme wanden (constante wandtemperatuur). Dit is een kritieke beperking, aangezien veel warmtewisselaars (zoals condensatoren) werken met isotherme wanden.

2. Methodologie

De auteurs hebben een Streamwise Periodic Reynolds-Averaged Navier-Stokes (SP-RANS) oplossing ontwikkeld en geïmplementeerd in de open-source CFD-software SU2.

Kernpunten van de methodologie:

• Governing Equations: De incompressibele Navier-Stokes-vergelijkingen worden aangevuld met brontermen (source terms) om de periodieke randvoorwaarden te simuleren.
  • Momentum: Een constante drukgradiënt ($\Delta p/L$) wordt toegevoegd als bronterm.
  • Energie: Voor isotherme wanden wordt een gereduceerde temperatuur ($\theta$) ingevoerd. De auteurs hebben een nieuwe bronterm voor de energie-equatie afgeleid die specifiek geldt voor turbulente stroming.
• Afleiding van de Turbulente Bronterm:
  De nieuwe bronterm voor de energie-equatie ($S_{\theta, turb}$) houdt rekening met het effect van turbulentie op de warmtediffusie. De formule combineert de laminaire term met een correctie voor de turbulentieviscositeit ($\mu_{turb}$) en het turbulent Prandtl-getal ($Pr_{turb}$):
  $$S_{\theta, turb} = S_{\theta, lam} - \theta \lambda_L \frac{\nabla \mu_{turb} \cdot \hat{i}_L}{\rho Pr_{turb}}$$
  Hierbij is $\lambda_L$ de exponentiële vervalconstante van het temperatuurverschil tussen het bulkfluid en de wand, die iteratief wordt opgelost.
• Validatie: De nieuwe solver werd gevalideerd door de resultaten te vergelijken met een volledige simulatie van een fin-array (11 vinnen) met een standaard RANS-solver.
• Testgeval: Een offset pin-fin array (gestaggerde opstelling) met een straal $r$.
  • Laminaire stroming: Reynoldsgetal ($Re$) = 100.
  • Turbulente stroming: $Re$ = 10.000.
  • Randvoorwaarde: Isotherme wanden ($T_w = 300$ K).

3. Belangrijkste Bijdragen

1. Wiskundige Formulering: De eerste openbare afleiding van de brontermen voor de energie-equatie van een SP-solver voor turbulente stroming met isotherme wanden.
2. Implementatie: De succesvolle integratie van deze termen in de open-source CFD-suite SU2.
3. Validatie: Een uitgebreide verificatie die aantoont dat de SP-solver zowel voor laminaire als turbulente regimes nauwkeurige resultaten levert die overeenkomen met volledige array-simulaties.
4. Inzicht in Ontwikkeling: Het identificeren van de exacte locatie waar de stroming hydraulisch en thermisch ontwikkeld is voor zowel laminaire als turbulente regimes in een pin-fin array.

4. Resultaten

• Nauwkeurigheid: De resultaten van de SP-solver (gebaseerd op één eenheidscel) kwamen perfect overeen met de resultaten van de volledige array-simulatie (11 vinnen) voor zowel snelheid, druk als temperatuurverdeling.
  • Voor de temperatuur werd de gereduceerde temperatuur ($\theta$) uit de SP-solver omgerekend naar de werkelijke temperatuur en vergeleken met de array-simulatie; de overeenkomst was exact.
  • Ook de turbulentie-eigenschappen (zoals de wervelviscositeit) werden correct voorspeld.
• Stromingsontwikkeling:
  • Laminaire stroming: De stroming werd hydraulisch en thermisch ontwikkeld na vinnen #5-#8.
  • Turbulente stroming: De stroming vereiste meer afstand om zich te ontwikkelen. De stroming werd hydraulisch ontwikkeld na vinnen #5, maar thermisch ontwikkeld pas na vinnen #7-#9. Dit suggereert dat inlaat-effecten bij turbulente stroming langer aanhouden dan bij laminaire stroming.
• Prestatieverschillen:
  • De temperatuurvervalconstante ($\lambda_L$) was ongeveer een orde van grootte kleiner bij turbulente stroming (0.0337 m⁻¹) dan bij laminaire stroming (0.2258 m⁻¹), wat wijst op een andere warmteoverdrachtsdynamiek.
  • De drukval ($\Delta p$) was bij laminaire stroming (6.41 Pa) ongeveer twee keer zo hoog als bij turbulente stroming (3.24 Pa) voor de specifieke testcondities.
• Rekenkosten:
  • Volledige Array: ~1440 minuten (ongeveer 1 dag).
  • SP-Solver (één cel): ~30 minuten.
  • Dit resulteert in een reductie van de rekentijd met een factor ~48.

5. Betekenis en Conclusie

Dit onderzoek opent de deur voor efficiënter ontwerp en optimalisatie van complexe warmtewisselaars. Door het gebruik van de ontwikkelde SP-RANS-solver kunnen ingenieurs:

• De enorme rekenkosten van het simuleren van volledige HEX-arrays vermijden.
• Betrouwbare thermohydraulische prestaties voorspellen voor geometrieën met isotherme wanden, zelfs in het turbulente regime.
• Optimalisatiealgoritmen (zoals vormoptimalisatie of topologieoptimalisatie) veel sneller en praktischer toepassen op complexe HEX-ontwerpen.

De studie bevestigt dat de afgeleide brontermen correct zijn en dat de SP-methode een krachtig alternatief is voor standaard CFD-simulaties, mits de stroming voldoende ver van de inlaat en uitlaat is om periodieke aannames te rechtvaardigen.