Impact of perturbed eddy-viscosity modeling on stability and shape sensitivity of the hydro-turbine vortex rope using linearized Reynolds-averaged Navier-Stokes equations

Deze studie toont aan dat het consistent lineariseren van het turbulentiemodel, in plaats van het 'bevriezen' ervan, essentieel is om de vormgevoeligheid van de vortex-touw-instabiliteit in een waterturbine correct te voorspellen en overeen te laten komen met experimentele trends.

De kern

Stel je voor dat je een enorme, krachtige watermolen hebt die elektriciteit maakt. Soms, als de waterstroom niet precies op het perfecte tempo loopt (bijvoorbeeld als er te weinig water is), begint er een gevaarlijk fenomeen op te treden: er vormt zich een grote, draaiende waterwervel in de uitlaatpijp. In de vakwereld noemen ze dit een "vortex rope" (wervelkabel).

Deze wervel is als een onrustige, dansende slang die tegen de wanden van de pijp slaat. Het veroorzaakt trillingen, lawaai en kan zelfs de molen beschadigen. De ingenieurs willen deze wervel stoppen of verzwakken, maar hoe doe je dat zonder de hele molen te slopen?

Dit onderzoek is een zoektocht naar het perfecte antwoord op de vraag: "Waar moeten we de pijp een beetje verbreden of versmallen om die wervel te kalmeren?"

Hier is hoe de onderzoekers dit aanpakken, vertaald in simpele taal:

1. De Simulatie: Een Digitale Proefpersoon

De onderzoekers bouwen geen nieuwe molen, maar maken een extreem gedetailleerde digitale kopie. Ze gebruiken wiskundige formules (de Navier-Stokes vergelijkingen) om te simuleren hoe het water stroomt.

Maar water is chaotisch en turbulent. Om dit hanteerbaar te maken, gebruiken ze een trucje: ze kijken niet naar elk individueel waterdruppeltje, maar naar het gemiddelde gedrag van de stroom. Dit noemen ze "RANS" (Reynolds-averaged Navier-Stokes).

2. Het Grote Geheim: De "Viskeuze" Wolk

Bij turbulente stroming is er een lastig probleem: hoe beschrijf je die kleine, chaotische werveltjes die de grote stroom beïnvloeden? De onderzoekers gebruiken een model dat dit simuleert met een soort "virtuele dikte" of wrijving van het water. Dit noemen ze "turbulentie" of "wervelviscositeit".

Hier komen twee manieren om dit te berekenen:

• De "Bevroren" Methode (De oude manier): Je doet alsof die wervelviscositeit statisch is. Het is als een bevroren beeld: je kijkt naar de wrijving en zegt: "Oké, dit is de wrijving, en die verandert niet, hoe de stroom ook beweegt." Het is makkelijk, maar misschien niet helemaal eerlijk.
• De "Gestoorde" Methode (De nieuwe manier): Je neemt aan dat de wrijving wel meebeweegt met de stroom. Als de stroom verandert, verandert ook de wrijving. Het is als kijken naar een levend, trillend weefsel in plaats van een stenen muur.

3. Het Experiment: De "Shape Sensitivity"

De kern van dit onderzoek is het testen van een specifieke vormverandering. Stel je voor dat je de binnenkant van de pijp een heel klein beetje kunt vervormen (zoals een zachte rubberen slang). De onderzoekers willen weten: "Als ik hier een beetje dikker maak, wordt de wervel dan rustiger of juist wilder?"

Dit noemen ze vormgevoeligheid (shape sensitivity).

4. De Verrassende Bevinding

Wat ze ontdekten, is fascinerend en een beetje tegenintuïtief:

• Voor de "Wervel" zelf (de frequentie en groei): Het maakt bijna geen verschil welke methode je kiest. Of je nu de "bevroren" of de "gestoorde" methode gebruikt, de computer voorspelt bijna hetzelfde gedrag voor de wervel zelf. Het is alsof je een auto rijdt; of je de bandenspanning exact meet of schat, de auto rijdt ongeveer hetzelfde.
• Voor de "Oplossing" (waar je moet ingrijpen): Hier wordt het cruciaal! De twee methoden geven volledig tegenovergestelde adviezen.
  • De "bevroren" methode zegt: "Maak de pijp op dit punt dunner om de wervel te stoppen."
  • De "gestoorde" methode zegt: "Nee, maak het juist dikker op dat punt!"

5. De Realiteitstest

De onderzoekers hebben hun digitale modellen vergeleken met echte experimenten in een waterlaboratorium.

• Het resultaat? De "gestoorde" methode had gelijk. Als je de pijp inderdaad dikker maakte op het voorgestelde punt, werd de wervel rustiger.
• De "bevroren" methode gaf een verkeerd advies. Als je die had gevolgd, had je de situatie waarschijnlijk verergerd.

De Metafoor: De Dansende Danseres

Stel je de waterwervel voor als een danseres die op een podium (de pijp) draait.

• De bevroren methode kijkt alleen naar de danseres en denkt: "Als ik de vloer hier een beetje schuiner maak, stopt ze met draaien."
• De gestoorde methode kijkt ook naar de danseres, maar ziet ook hoe de vloer zelf reageert op haar bewegingen. Ze ziet dat als ze de vloer schuin maakt, de vloer zelf ook gaat trillen en de danseres juist meer aanmoedigt om te draaien.
• De gestoorde methode zegt: "Nee, maak de vloer juist steviger en dikker op die plek, dan dempt de trilling van de vloer haar dans."

Conclusie

De boodschap van dit onderzoek is simpel maar belangrijk voor de toekomst:
Als je wilt voorspellen hoe je een complexe, turbulente stroming (zoals in een watermolen, een vliegtuigvleugel of een auto) kunt verbeteren door de vorm te veranderen, mag je de kleine details van de turbulentie niet "bevroren" laten. Je moet meedenken met de turbulentie.

Alleen door die dynamische wrijving mee te nemen in de berekeningen, krijg je het juiste advies over waar je de vorm moet aanpassen. Het is het verschil tussen een gok doen en een wetenschappelijk onderbouwde oplossing vinden.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Impact of perturbed eddy-viscosity modeling on stability and shape sensitivity of the hydro-turbine vortex rope using linearized Reynolds-averaged Navier–Stokes equations", geschreven in het Nederlands.

Probleemstelling

In hydro-turbines, specifiek bij deelslast-operatie (part-load), ontstaat er een fenomeen genaamd de "vortex rope" (wervelkoord). Dit is een niet-lineair verzadigde globale instabiliteit die leidt tot grote drukpulsaties en structurele trillingen, wat een significant risico vormt voor de veilige werking van de turbine. Het controleren en mitigeren van deze instabiliteit is cruciaal.

Linear Stability Analysis (LSA) en adjoint-gebaseerde gevoeligheidsanalyse (sensitivity analysis) bieden een fysisch onderbouwd kader om deze instabiliteiten te begrijpen en flow control strategieën te ontwikkelen (bijv. via vormoptimalisatie). Echter, bij toepassing op volledig turbulente stromingen (zoals in een Francis-turbine) rijst de vraag hoe turbulentiegemodelleerd moet worden binnen de linearisatie. Er zijn twee benaderingen:

1. Bevroren eddy-viscositeit (Frozen): De turbulentieparameters (zoals eddy-viscositeit $\nu_t$) worden als constant beschouwd tijdens de linearisatie; fluctuaties in turbulentie worden genegeerd.
2. Gestoorde eddy-viscositeit (Perturbed): De turbulentieparameters worden expliciet gelineariseerd, waardoor fluctuaties in de turbulentie (via de $k-\varepsilon$ vergelijkingen) meegenomen worden.

De huidige literatuur gebruikt vaak de "bevroren" benadering zonder sterke validatie, hoewel het onduidelijk is of dit leidt tot fysisch correcte resultaten voor gevoeligheidsanalyses, vooral voor vormoptimalisatie.

Methodologie

De auteurs ontwikkelen en valideren een fysisch kader op basis van de lineair gelineariseerde Reynolds-averaged Navier-Stokes (RANS) vergelijkingen met een standaard $k-\varepsilon$ turbulentiemodel.

• Configuratie: Het onderzoek focust op een Francis-turbine (gebaseerd op de Francis-99 benchmark) met een focus op de zuigbuis (draft tube) bij deelslast. De geometrie wordt gemodelleerd als een 2D axiale-symmetrische domein met een centraal lichaam (centerbody).
• Basisstroom: De basisstroom wordt berekend met RANS-simulaties. Om de bifurcatiepunten van de vortex rope te matchen met 3D URANS-simulaties en experimenten, worden de inlet-condities (turbulentie-intensiteit) aangepast (getuned) om de effecten van niet-opgeloste wake-structuren van de runner te simuleren.
• Linearisatie: De vergelijkingen worden gelineariseerd rondom de turbulente basisstroom. Er worden twee modellen vergeleken:
  • Het perturbed model: Linearisatie van alle variabelen, inclusief $k$, $\varepsilon$ en $\nu_t$.
  • Het frozen model: Linearisatie waarbij $\delta \nu_t = 0$, $\delta k = 0$ en $\delta \varepsilon = 0$.
• Gevoeligheidsanalyse: Er wordt een vormgevoeligheidsanalyse (shape sensitivity) uitgevoerd om te bepalen hoe kleine vervormingen van de wand (specifiek de dikte van het centerbody) de groeifactor (growth rate) van de instabiliteit beïnvloeden. De gevoeligheid wordt ontbonden in bijdragen van de basisstroom (base-flow contribution) en feedback-mechanismen (feedback contribution).
• Validatie: De resultaten worden vergeleken met experimentele data (drukmetingen en stochastische modelleermethodes voor groeifactoren) en 3D URANS-simulaties.

Belangrijkste Bijdragen

1. Vergelijking van Turbulentie-modellen: Het artikel biedt een systematische vergelijking tussen de "frozen" en "perturbed" eddy-viscositeit modellen binnen een RANS-LSA kader voor een complexe turbulente stroming.
2. Ontleding van Gevoeligheid: De auteurs tonen aan dat de totale vormgevoeligheid voornamelijk wordt gedreven door de bijdrage van de basisstroom (via de Hessian-operator), en niet door de directe feedback van de lineaire operator.
3. Fysische Mechanismen: Het werk identificeert specifieke fysische mechanismen die door het perturbed model worden gevangen (zoals diffusie en productie van turbulentiefluctuaties) maar volledig worden gemist door het frozen model.
4. Experimentele Validatie: Voor het eerst wordt de vormgevoeligheid van een turbulente vortex rope in een Francis-turbine experimenteel gevalideerd, waarbij wordt aangetoond dat alleen het perturbed model de correcte trends voorspelt.

Resultaten

• Eigenwaarden en Eigenmodes: De invloed van het perturbed eddy-viscositeit model op de eigenwaarden (frequentie en groeifactor) en de vorm van de eigenmodes is marginaal. Zowel het frozen als het perturbed model voorspellen vrijwel dezelfde stabiliteitsgrenzen en oscillatiefrequenties. Dit suggereert dat de lineaire operator zelf weinig gevoelig is voor de toevoeging van de gelineariseerde $k-\varepsilon$ vergelijkingen.
• Vormgevoeligheid (Shape Sensitivity): Er is een drastisch verschil in de voorspelde vormgevoeligheid tussen de twee modellen:
  • Het frozen model voorspelt dat het verdunnen van het centerbody op de plaats van maximale dikte de stroming stabiliseert.
  • Het perturbed model voorspelt juist dat het verdikken van het centerbody op die locatie de stroming stabiliseert.
  • De experimentele data bevestigt de voorspelling van het perturbed model: verdikken stabiliseert de stroming.
• Oorzaak van het verschil: Het verschil wordt veroorzaakt door de bijdrage van de basisstroom. In het perturbed model leidt een verdikking van het centerbody tot een toename van de schuifspanning en daarmee de eddy-viscositeit ($\nu_t$). Deze toename in $\nu_t$ dempt de coherente fluctuaties (stabiliserend effect), wat het destabiliserende effect van een grotere wake-achterstand (wake deficit) compenseert. Het frozen model mist dit dempingsmechanisme volledig en ziet alleen het destabiliserende effect van de wake-verandering.
• Dominantie van Basisstroom: De analyse toont aan dat de gevoeligheid van de eigenwaarde ten opzichte van vormveranderingen gedomineerd wordt door hoe de basisstroom verandert (via de Hessian), en niet door de directe verandering in de lineaire operator.

Significantie en Conclusie

De studie concludeert dat voor gevoeligheidsgebaseerde controle en vormoptimalisatie van turbulente globale instabiliteiten, het consistent lineariseren van turbulentiemodellen (het gebruik van een perturbed eddy-viscositeit model) essentieel is.

Hoewel het frozen model voldoende kan zijn om de stabiliteit (eigenwaarden) van een stroming te voorspellen, faalt het volledig bij het voorspellen van de richting van optimale vormveranderingen. Het negeren van turbulentiefluctuaties leidt tot fysisch onjuiste adviezen voor ontwerp (bijv. het verdunnen in plaats van verdikken van componenten).

De bevindingen onderstrepen dat de RANS-LSA framework, wanneer correct uitgerust met een perturbed turbulentiemodel, een krachtig en betrouwbaar instrument is voor het ontwerpen van Francis-turbines met een breder stabiel operationeel bereik, zelfs in volledig turbulente omstandigheden. Dit opent de weg voor gradient-based shape optimization in complexe turbulente stromingen.