Steady, Quasi-Periodic, and Chaotic Regimes of Natural Convection Around a Heated Cylinder in a Square Cavity

Technische Samenvatting: Stabiele, Quasi-periodieke en Chaotische Regimes van Natuurlijke Convectie rond een Verwarmde Cilinder in een Vierkante Holte

Probleemstelling
Deze studie onderzoekt de onstabiele natuurlijke convectie van lucht ($Pr = 0,71$) binnen een koude vierkante behuizing met een centraal verwarmde cirkelvormige cilinder met een straal $R = 0,2L$. Hoewel de differentieel verwarmde vierkante holte een standaard benchmark is, introduceert de toevoeging van een interne verwarmde cilinder complexe reorganisatie van de stroming. Eerdere literatuur heeft zich grotendeels gericht op stationaire oplossingen of beperkte onstabiele bereiken, waarbij vaak kunstmatige symmetriebeperkingen werden opgelegd (bijv. het afdwingen van een nul streamfunctie op de cilinder), wat asymmetrische stromingsmodi onderdrukt. Bijgevolg bleven het kritische Rayleigh-getal ($Ra_c$) voor het begin van onstabiliteit, de specifieke route naar chaos en de volledige sequentie van stromingsregimes voor deze configuratie onopgelost. Dit werk beoogt deze hiaten te overbruggen door de overgang van stationaire laminaire stroming naar chaos te simuleren zonder symmetrie af te dwingen, waardoor asymmetrische modi natuurlijk kunnen groeien.

Methodologie
De auteurs gebruikten een drukgebaseerde eindelementenmethode (FEM) op een lichaam-conforme mesh. Belangrijke methodologische kenmerken zijn:

• Beherende Vergelijkingen: De onstabiele 2D Navier-Stokes en energievergelijkingen onder de Boussinesq-benadering werden opgelost in primitieve variabelen. Er werden geen symmetriecondities opgelegd op de verticale middellijn.
• Mesh-generatie: Een lichaam-conforme mesh werd gegenereerd vanuit een Cartesiaans achtergrondrooster via een ray-cast projectie gevolgd door harmonische smoothing. Deze aanpak lost de gebogen cilindergrens op zonder coördinatentransformatie of immersed boundary interpolatie.
• Numeriek Schema: Een expliciet fractional-step (projectie) algoritme met een tweede-orde Adams-Bashforth behandeling voor convectie-diffusie termen werd gebruikt. De druk-Poisson vergelijking werd opgelost met een warm-gestarte Conjugate Gradient methode met Element-By-Element (EbE) matrix-vrije producten om ongestructureerde meshes efficiënt te verwerken.
• Validatie en Grid-onafhankelijkheid: De solver werd gevalideerd tegen de De Vahl Davis benchmark ($Ra = 10^6$). Een Grid Convergence Index (GCI) studie bevestigde grid-onafhankelijkheid met een fijne-grid onzekerheid van 0,84% op het Nusselt-getal van de cilinder.
• Simulatieomvang: Acht Rayleigh-getallen werden gesimuleerd, variërend van $10^4$ tot $5 \times 10^7$. Tijdshistorieën van het cilinder-gemiddelde Nusselt-getal ($Nu_{cyl}$) werden opgenomen, en energie-spectrale dichtheden (ESD) en Poincaré-kaarten werden gegenereerd met behulp van data van elf meetpunten verdeeld over de holte.

Belangrijkste Resultaten
De studie identificeert drie verschillende stromingsregimes en karakteriseert de overgang tussen hen:

1. Stationair Regime ($10^4 \le Ra \le 5 \times 10^6$):

   • De stroming blijft stationair en symmetrisch rond de verticale middellijn, gekenmerkt door twee spiegelbeeldige recirculatiecellen en een enkele thermische pluim die van de cilinder naar de bovenwand stijgt.
   • Bij $Ra = 5 \times 10^6$ verschijnt het eerste fysieke signaal van symmetriebreking: de magnitude van de negatieve streamfunctie-kern ($|\psi_{min}|$) overtreft de positieve kern ($\psi_{max}$) met ongeveer 3,6%, hoewel de stroming visueel symmetrisch blijft.
   • Warmteoverdracht wordt gedomineerd door de bovenwand (die ~55% van de totale wandwarmte absorbeert), terwijl de onderwand minder dan 2% accountt door thermische stratificatie in de onderste holte.

2. Begin van Onstabiliteit ($5 \times 10^6 < Ra_c < 7,5 \times 10^6$):

   • Het kritische Rayleigh-getal voor het begin van onstabiliteit wordt ingeklemd tussen $5 \times 10^6$ en $7,5 \times 10^6$. Dit wordt bepaald door een scherpe sprong (meer dan drie grootteordes) in de standaarddeviatie van $Nu_{cyl}$.
   • Deze kritische waarde is hoger dan wat door eerdere stationaire studies werd geïmpliceerd die symmetrie afdwongen, wat bevestigt dat het onderdrukken van asymmetrische modi de transitie kunstmatig vertraagt.

3. Onstabiel Regime ($Ra \ge 7,5 \times 10^6$):

   • Limit Cycle ($Ra = 7,5 \times 10^6$): De stroming vertoont een single-frequency periodieke oscillatie (limit cycle). De wanorde vindt oorsprong in de pluim boven de cilinder.
   • Quasi-periodiek ($Ra = 10^7$): Een tweede incommensurabele frequentie ontstaat, waardoor de stroming op een twee-torus terechtkomt. De Poincaré-secties tonen een gesloten curve, en de onderste gestratificeerde pool begint de secundaire frequentie te registreren.
   • Chaos ($Ra = 5 \times 10^7$): De stroming gaat over naar een breedbandige strange attractor. De Poincaré-secties onthullen een verspreide wolk van punten, en de energie-spectrale dichtheid wordt continu. De wanorde plant zich voort van de pluim naar de zijwanden en bereikt uiteindelijk de onderste gestratificeerde pool, die actief wordt over een breed scala aan tijdschalen.
   • Symmetrie in Gemiddelden: Ondanks dat instantane snapshots sterk asymmetrisch zijn, herstellen de tijdgemiddelde velden de bilaterale symmetrie. Dit is een statistisch effect waarbij het systeem de spiegelbeeldige toestanden even vaak bezoekt als de andere.

Warmteoverdracht Correlaties

• Een enkele machtswet, $Nu_{cyl} = 0,764 Ra^{0,226}$, past op de tijdgemiddelde cilinder Nusselt-getallen over het gehele bereik ($10^4 \le Ra \le 5 \times 10^7$) met een afwijking van minder dan 6%.
• De exponent voor de stationaire tak is $n \approx 0,224$, en voor de onstabiele tak is $n \approx 0,243$. Beide waarden liggen tussen de klassieke laminaire grenslaag voorspelling ($n=0,25$) en de Castaing soft-turbulence schaling ($n \approx 0,286$), wat suggereert dat de stroming niet volledig turbulent is maar beïnvloed wordt door fluctuerende velden.
• Het Nusselt-getal van de onderwand ($Nu_{bot}$) neemt af naarmate $Ra$ toeneemt, een contra-intuïtieve trend die wordt toegeschreven aan de versterkende thermische stratificatie en het vangen van koude jets in de onderste hoek-eddy's.

Betekenis en Claims
Het artikel claimt de eerste uitgebreide numerieke karakterisering te bieden van de volledige transitievolgorde van stationaire naar chaotische natuurlijke convectie in een centraal verwarmde cilinderconfiguratie zonder kunstmatige symmetriebeperkingen. Belangrijke bijdragen zijn:

• Kritisch Rayleigh-getal: Het vaststellen van een nauwkeuriger bereik voor $Ra_c$ door asymmetrische modi toe te laten, wat eerdere overschattingen door symmetriebeperkingen corrigeert.
• Route naar Chaos: Bevestiging dat de transitie de Ruelle-Takens-Newhouse route volgt (Limit Cycle $\to$ Two-Torus $\to$ Strange Attractor) in plaats van een periodedubbelingscascade.
• Ruimtelijke Voortplanting van Wanorde: Het aantonen dat chaos initieert in de thermische pluim en naar beneden propageert, waarbij de onderste gestratificeerde pool relatief rustig blijft tot de hoogste Rayleigh-getallen.
• Referentiedata: Het leveren van betrouwbare benchmarkdata voor het kritische Rayleigh-getal, regime-transities en globale warmteoverdrachtcorrelaties voor configuraties waar interne geometrie en onstabiele dynamica interageren.

De auteurs merken op dat de studie beperkt is tot 2D, een enkele cilinderstraal ($R=0,2L$), en $Pr=0,71$, wat suggereert dat uitbreidingen naar 3D, uit het midden geplaatste posities en variërende Prandtl-getallen toekomstige richtingen zijn.