🔬 physics

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

एक वर्गाकार गुहा (कैविटी) में गर्म बेलन के चारों ओर प्राकृतिक संवहन का यह संख्यात्मक अध्ययन यह प्रकट करता है कि जैसे-जैसे रेले संख्या (रेले नंबर) $10^4$ से बढ़कर $5 \times 10^7$ हो जाती है, प्रवाह 'रुएल-टेकेन्स-न्यूहाउस' मार्ग के माध्यम से एक स्थिर अवस्था से अराजकता (केओस) में परिवर्तित हो जाता है, जबकि बेलन के नुसेल्ट संख्या के लिए एक विशिष्ट पावर-लॉ सहसंबंध स्थापित करता है और तापीय स्तरीकरण के कारण ऊष्मा स्थानांतरण में ऊपरी दीवार की प्रमुख भूमिका को रेखांकित करता है।

मूल लेखक: Mannu Yadav, Anil Lal Sadasivan

प्रकाशित 2026-07-07

📖 6 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Mannu Yadav, Anil Lal Sadasivan

मूल पेपर CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

एक वर्गाकार कमरे की कल्पना करें जिसकी दीवारें ठंडी हैं और ठीक बीच में एक गर्म, चमकता हुआ हीटर रखा है। यह अध्ययन में मौजूद "वर्गाकार गुहा" (square cavity) है। कमरे के अंदर की हवा केवल स्थिर नहीं है; क्योंकि हीटर गर्म है और दीवारें ठंडी हैं, इसलिए हवा हिलना चाहती है। गर्म हवा ऊपर उठती है, ठंडी हवा नीचे गिरती है, जिससे एक प्राकृतिक परिसंचरण (circulation) बनता है। यह प्राकृतिक संवहन (natural convection) है।

शोधकर्ता यह देखना चाहते थे कि जब वे गर्मी बढ़ाते हैं (जिसे रेले संख्या/Rayleigh number कहा जाता है), तो इस हवा के नृत्य के साथ क्या होता है। उन्होंने केवल अनुमान नहीं लगाया; उन्होंने हवा की गति को वास्तविक समय में देखने के लिए एक कंप्यूटर प्रोग्राम का उपयोग करके एक सुपर-सटीक डिजिटल मॉडल बनाया, बिना उसे किसी विशिष्ट तरीके से व्यवहार करने के लिए मजबूर किए।

यहाँ उनके द्वारा की गई खोज की कहानी है, जिसे सरल चरणों में विभाजित किया गया है:

1. शांत शुरुआत (स्थिर अवस्था - Steady State)

जब गर्मी कम या मध्यम होती है, तो हवा एक अच्छी तरह से प्रशिक्षित बैले (ballet) की तरह व्यवहार करती है।

पैटर्न: हीटर के दोनों ओर हवा के दो विशाल, दर्पण-प्रतिबिंब वाले लूप बनते हैं। हवा हीटर के किनारों से ऊपर जाती है, ठंडी छत से टकराती है, दीवारों के सहारे नीचे फिसलती है और नीचे की ओर वापस आती है।
प्लूम (Plume): गर्म हवा का एक एकल, स्थिर स्तंभ (एक "प्लूम") हीटर के ऊपरी हिस्से से सीधे छत की ओर उठता है।
परिणाम: सब कुछ अनुमानित और सममित (symmetrical) है। यदि आप कमरे के बीच में एक रेखा खींचते हैं, तो बायां हिस्सा बिल्कुल दाएं हिस्से जैसा दिखता है। यह एक निश्चित गर्मी के स्तर (लगभग 5 मिलियन यूनिट गर्मी) तक होता है।

2. परेशानी का पहला संकेत (टिपिंग पॉइंट)

जैसे ही उन्होंने गर्मी को थोड़ा और बढ़ाया, कुछ दिलचस्प हुआ। पूर्ण समरूपता (symmetry) टूटने लगी, हालांकि आप इसे नग्न आंखों से देख नहीं सकते थे।

बदलाव: बाईं और दाईं ओर के हवा के लूप अब पूरी तरह से समान नहीं रहे। एक तरफ का लूप दूसरी तरफ की तुलना में थोड़ा अधिक मजबूत होने लगा।
खोज: शोधकर्ताओं ने पाया कि वह "टिपिंग पॉइंट" जहाँ शांत, स्थिर प्रवाह अराजकता (chaos) में बदल जाता है, वह 5 मिलियन और 7.5 मिलियन यूनिट गर्मी के बीच कहीं है। इससे पहले हवा शांत थी; इसके बाद, इसमें थरथराहट शुरू हो गई।

3. नृत्य जटिल होता जाता है (क्वासी-पीरियडिक - Quasi-Periodic)

एक बार जब गर्मी उस टिपिंग पॉइंट को पार कर गई, तो हवा ने एक शांत बैले करना बंद कर दिया और एक जटिल, लयबद्ध नृत्य करना शुरू कर दिया।

लिमिट साइकिल (Limit Cycle): शुरुआत में, हवा एक एकल, पूर्ण लय (जैसे एक पेंडुलम) में आगे-पीछे झूलने लगी।
टू-टोरस (Two-Torus): जैसे-जैसे उन्होंने गर्मी बढ़ाई, नृत्य और जटिल हो गया। अब यह केवल आगे-पीछे झूलना नहीं था; यह झूलते हुए डगमगाना भी शुरू हो गया। एक घूमते हुए लट्टू की कल्पना करें जो आगे-पीछे भी डगमगा रहा है। हवा अब दो अलग-अलग लयों पर एक साथ चल रही थी जो एक-दूसरे से मेल नहीं खाती थीं। इसे "क्वासी-पीरियडिक" कहा जाता है।

4. अराजक तूफान (Chaos)

उच्चतम गर्मी के स्तरों (50 मिलियन यूनिट) पर, व्यवस्थित नृत्य पूरी तरह से बिखर गया।

स्ट्रेंज अट्रैक्टर (Strange Attractor): हवा अराजक हो गई। यह यादृच्छिक शोर (random noise) नहीं था, बल्कि अप्रत्याशित था। गर्म प्लूम केवल सीधा ऊपर नहीं उठता था; यह कभी बाईं ओर झुकता, कभी दाईं ओर, और फिर अलग होकर छोटे-छोटे गोलों में घूमने लगता।
अराजकता का प्रसार: दिलचस्प बात यह है कि अराजकता हर जगह एक साथ नहीं हुई। यह हीटर के ठीक ऊपर, प्लूम से शुरू हुई। कमरे के निचले हिस्से की हवा लंबे समय तक शांत और स्तरित (stratified) रही। केवल तभी जब गर्मी अत्यधिक उच्च हो गई, तब अराजक भंवर अंततः कमरे के निचले हिस्से तक पहुँचे।

5. "भूतिया" समरूपता (The "Ghost" Symmetry)

यहाँ एक पेचीदा हिस्सा है जिसे पेपर में उजागर किया गया है:

तत्काल बनाम औसत: यदि आप अराजक चरण के दौरान किसी भी क्षण हवा की एक तस्वीर लेते हैं, तो यह अव्यवस्थित और असममित (लॉपसाइड) दिखाई देती है। हालाँकि, यदि आप एक "टाइम-लैप्स" फोटो लेते हैं (एक लंबे समय के दौरान औसत निकालते हैं), तो वह अव्यवस्था समाप्त हो जाती है, और छवि फिर से पूरी तरह से सममित दिखाई देती है, ठीक शांत शुरुआत की तरह।
रूपक (Metaphor): एक भीड़ की कल्पना करें जो अराजक, जंबल में दौड़ रही है। यदि आप एक सेकंड के लिए उनकी फोटो लेते हैं, तो वह एक धुंधला सा दिखता है। यदि आप एक घंटे में फोटो लेते हैं, तो वह औसत निकलकर एक पूर्ण वृत्त की तरह दिखता है। औसत शांत दिखता है, लेकिन वास्तविकता अराजक है।

6. ऊष्मा स्थानांतरण (गर्मी कहाँ जाती है?)

शोधकर्ताओं ने इस बात को भी ट्रैक किया कि दीवारों ने कितनी गर्मी अवशोषित की।

छत विजेता है: चाहे हवा कितनी भी अराजक क्यों न हो जाए, ऊपरी दीवार ने हमेशा लगभग 55% गर्मी अवशोषित की। बगल की दीवारों ने बाकी हिस्सा (लगभग 22% प्रत्येक) साझा किया।
फर्श को अनदेखा किया गया: निचली दीवार ने लगभग कुछ भी (2% से कम) अवशोषित नहीं किया।
क्यों? जैसे-जैसे हवा गर्म होती गई, नीचे की ठंडी हवा एक परत में "फंसी" रह गई। हीटर से उठने वाली गर्म हवा और किनारों से नीचे गिरने वाली ठंडी हवा ने एक बाधा बना दी जिसने फर्श को अलग कर दिया। यह कमरे के निचले हिस्से के एक "डेड ज़ोन" बनने जैसा है जहाँ गर्मी पहुँच नहीं सकती, भले ही हीटर पहले से कहीं अधिक काम कर रहा हो।

सारांश

यह पेपर एक बॉक्स में हीटर के साथ हवा की कहानी बताता है।

कम गर्मी: शांत, सममित लूप।
मध्यम गर्मी: लूप डगमगाने लगते हैं और समरूपता टूटने लगती है।
उच्च गर्मी: हवा दो आवृत्तियों (frequencies) के साथ एक जटिल, लयबद्ध नृत्य में प्रवेश करती है।
अत्यधिक गर्मी: हवा पूरी तरह से अराजक हो जाती है, जो ऊपर से शुरू होकर धीरे-धीरे नीचे की ओर फैलती है।
सीख: भले ही हवा अराजक हो, औसत ऊष्मा वितरण शांत दिखता है, और कमरे का निचला हिस्सा आश्चर्यजनक रूप से ठंडा रहता है क्योंकि गर्मी ऊपरी परतों में फंसी रहती है।

शोधकर्ताओं ने यह सुनिश्चित करने के लिए एक विशेष कंप्यूटर विधि (फाइनाइट एलीमेंट मेथड) का उपयोग किया कि उनका डिजिटल मॉडल सटीक हो, और अपने नए प्रयोग चलाने से पहले इसे ज्ञात बेंचमार्क के साथ जांचा। उन्होंने पाया कि पिछले अध्ययनों ने इस अराजकता की शुरुआत को मिस कर दिया होगा क्योंकि उन्होंने हवा को सममित रहने के लिए मजबूर किया था, जिससे "डगमगाहट" स्वाभाविक रूप से नहीं हो सकी।

तकनीकी सारांश: एक वर्गाकार कैविटी में गर्म सिलेंडर के चारों ओर प्राकृतिक संवहन की स्थिर, अर्ध-आवधिक और अराजक अवस्थाएँ

समस्या विवरण
यह अध्ययन एक ठंडे वर्गाकार एनक्लोजर के भीतर एक केंद्रीय रूप से गर्म गोलाकार सिलेंडर (त्रिज्या $R = 0.2L$ ) में वायु ($Pr = 0.71 $) के अनस्थिर प्राकृतिक संवहन की जांच करता है। जबकि विभेदित रूप से गर्म वर्गाकार कैविटी एक मानक बेंचमार्क है, आंतरिक गर्म सिलेंडर का समावेश जटिल प्रवाह पुनर्गठन पेश करता है। पिछला साहित्य मुख्य रूप से स्थिर अवस्था समाधानों या सीमित अनस्थिर श्रेणियों पर केंद्रित रहा है, जो अक्सर कृत्रक समरूपता बाधाओं (जैसे, सिलेंडर पर स्ट्रीम फंक्शन को शून्य पर क्लैंप करना) को लागू करता है जो असममित प्रवाह मोड को दबा देते हैं। फलस्वरूप, अनस्थिरता के प्रारंभ के लिए महत्वपूर्ण रेले संख्या ($ Ra_c$), अराजकता का विशिष्ट मार्ग, और इस विन्यास के लिए प्रवाह व्यवस्थाओं का पूर्ण अनुक्रम अनसुलझा रहा। यह कार्य समरूपता को लागू किए बिना, असममित मोड को स्वाभाविक रूप से विकसित होने की अनुमति देकर, स्थिर लैमिनार प्रवाह से अराजकता (chaos) तक के संक्रमण का अनुकरण करने के माध्यम से इन अंतरालों को भरने का लक्ष्य रखता है।

कार्यप्रणाली
लेखकों ने एक बॉडी-कन्फॉर्मिंग मेश पर प्रेशर-बेस्ड फाइनाइट एलीमेंट मेथड (FEM) का उपयोग किया। प्रमुख कार्यप्रणाली संबंधी विशेषताएं शामिल हैं:

शासी समीकरण (Governing Equations): बूसिनेस्क सन्निकटन (Boussinesq approximation) के तहत प्रिमिटिव वेरिएबल्स में अनस्थिर 2D नेवियर-स्टोक्स और ऊर्जा समीकरणों को हल किया गया। ऊर्ध्वाधर मध्य रेखा पर कोई समरूपता शर्तें लागू नहीं की गईं।
मेश जनरेशन: एक रे-कास्ट प्रोजेक्शन के बाद हार्मोनिक स्मूथिंग का उपयोग करके एक कार्टेशियन बैकग्राउंड ग्रिड से बॉडी-कन्फॉर्मिंग मेश तैयार किया गया। यह दृष्टिकोण बिना किसी कोऑर्डिनेट ट्रांसफॉर्मेशन या इमर्स्ड बाउंड्री इंटरपोलेशन के घुमावदार सिलेंडर सीमा को हल करता है।
संख्यात्मक योजना (Numerical Scheme): कन्वेक्शन-डिफ्यूजन टर्म्स के लिए सेकंड-ऑर्डर एडम्स-बाशफोर्थ उपचार के साथ एक एक्सप्लिसिट फ्रैक्शनल-स्टेप (प्रोजेक्शन) एल्गोरिदम का उपयोग किया गया। प्रेशर पॉइसन समीकरण को एलीमेंट-बाय-एलीमेंट (EbE) मैट्रिक्स-फ्री उत्पादों के साथ एक वॉर्म-स्टार्टेड कंजुगेट ग्रेडिएंट विधि का उपयोग करके हल किया गया ताकि अनस्ट्रक्चर्ड मेश को कुशलतापूर्वक संभाला जा सके।
सत्यापन और ग्रिड स्वतंत्रता: सॉल्वर को डे वाल डेविस बेंचमार्क ( $Ra = 10^6$ ) के विरुद्ध सत्यापित किया गया। एक ग्रिड कन्वर्जेंस इंडेक्स (GCI) अध्ययन ने सिलेंडर नसेलट संख्या पर 0.84% की फाइन-ग्रिड अनिश्चितता के साथ ग्रिड स्वतंत्रता की पुष्टि की।
सिमुलेशन दायरा: आठ रेले संख्याओं का सिमुलेशन किया गया, जो $10^4$ से $5 \times 10^7$ तक की सीमा में थीं। सिलेंडर-औसत नसेलट संख्या ( $Nu_{cyl}$ ) के समय इतिहास को रिकॉर्ड किया गया, और कैविटी में वितरित ग्यारह प्रोब बिंदुओं का उपयोग करके ऊर्जा स्पेक्ट्रल डेंसिटी (ESD) और पोइनकेरे मैप्स (Poincaré maps) बनाए गए।

प्रमुख परिणाम
अध्ययन तीन विशिष्ट प्रवाह व्यवस्थाओं की पहचान करता है और उनके बीच के संक्रमण को स्पष्ट करता है:

स्थिर अवस्था ( $10^4 \le Ra \le 5 \times 10^6$ ):
- प्रवाह स्थिर रहता है और ऊर्ध्वाधर मध्य रेखा के चारों ओर सममित होता है, जिसमें दो मिरर-इमेज पुनरावृत्ति सेल और सिलेंडर से ऊपरी दीवार तक उठने वाला एक एकल थर्मल प्लूम (thermal plume) होता है।
- $Ra = 5 \times 10^6$ पर, समरूपता टूटने का पहला भौतिक संकेत दिखाई देता है: नकारात्मक स्ट्रीम फंक्शन कोर का परिमाण ( $|\psi_{min}|$ ) सकारात्मक कोर ( $\psi_{max}$ ) से लगभग 3.6% अधिक हो जाता है, हालांकि प्रवाह दृश्य रूप से सममित रहता है।
- ऊष्मा स्थानांतरण ऊपरी दीवार द्वारा संचालित होता है (जो कुल दीवार ऊष्मा का ~55% अवशोषित करती है), जबकि निचली दीवार निचले कैविटी में थर्मल स्ट्रैटिफिकेशन के कारण 2% से भी कम योगदान देती है।
अनस्थिरता का प्रारंभ ( $5 \times 10^6 < Ra_c < 7.5 \times 10^6$ ):
- अनस्थिरता के प्रारंभ के लिए महत्वपूर्ण रेले संख्या $5 \times 10^6$ और $7.5 \times 10^6$ के बीच स्थित है। यह $Nu_{cyl}$ के मानक विचलन में एक तीव्र उछाल (तीन आदेशों से अधिक) द्वारा निर्धारित किया गया है।
- यह महत्वपूर्ण मान उन पिछले स्थिर-केवल अध्ययनों द्वारा सूचित मान से अधिक है जिन्होंने समरूपता लागू की थी, जो पुष्टि करता है कि असममित मोड को दबाने से संक्रमण कृत्रिम रूप से विलंबित हो जाता है।
अनस्थिर अवस्था ( $Ra \ge 7.5 \times 10^6$ ):
- लिमिट साइकिल ( $Ra = 7.5 \times 10^6$ ): प्रवाह एक सिंगल-फ्रीक्वेंसी आवधिक दोलन (लिमिट साइकिल) प्रदर्शित करता है। विकार (disorder) सिलेंडर के ऊपर प्लूम में उत्पन्न होता है।
- अर्ध-आवधिक/क्वासी-पीरियॉडिक ( $Ra = 10^7$ ): एक दूसरी असंगत आवृत्ति उभरती है, जो प्रवाह को एक टू-टोरस (two-torus) पर रखती है। पोइनकेरे सेक्शन एक बंद वक्र दिखाते हैं, और निचला स्ट्रैटिफाइड पूल माध्यमिक आवृत्ति को दर्ज करने लगता है।
- अराजकता/केओस ( $Ra = 5 \times 10^7$ ): प्रवाह एक ब्रॉड-बैंड स्ट्रेंज अट्रैक्टर की ओर संक्रमण करता है। पोइनकेरे सेक्शन बिखरे हुए बिंदुओं का एक क्लाउड दिखाते हैं, और ऊर्जा स्पेक्ट्रल डेंसी निरंतर (continuous) हो जाती है। विकार प्लूम से नीचे की ओर साइड दीवारों तक फैलता है, और अंततः निचले स्ट्रैटिफाइड पूल तक पहुँच जाता है, जो विभिन्न समय-पैमानों (timescales) पर सक्रिय हो जाता है।
- औसत में समरूपता: हालांकि तात्कालिक स्नैपशॉट अत्यधिक असममित होते हैं, समय-औसत क्षेत्र द्विपक्षीय समरूपता को पुनः प्राप्त करते हैं। यह एक सांख्यिकीय प्रभाव है जहाँ सिस्टम समय के साथ मिरर-इमेज अवस्थाओं का समान रूप से दौरा करता है।

ऊष्मा स्थानांतरण सहसंबंध (Heat Transfer Correlations)

एक एकल पावर लॉ, $Nu_{cyl} = 0.764 Ra^{0.226}$ , पूरे रेंज ( $10^4 \le Ra \le 5 \times 10^7$ ) में समय-औसत सिलेंडर नसेलट नंबरों के लिए फिट बैठता है, जिसमें विचलन 6% से कम है।
स्थिर शाखा के लिए घातांक $n \approx 0.224$ है, और अनस्थिर शाखा के लिए $n \approx 0.243$ है। दोनों मान शास्त्रीय लैमिनार बाउंड्री लेयर भविष्यवाणी ( $n=0.25$ ) और कास्टांग सॉफ्ट-टर्बुलेंस स्केलिंग ( $n \approx 0.286$ ) के बीच स्थित हैं, जो सुझाव देते हैं कि प्रवाह पूरी तरह से टर्बुलेंट नहीं है बल्कि उतार-चढ़ाव वाले क्षेत्रों (fluctuating fields) से प्रभावित है।
निचली दीवार की नसेलट संख्या ( $Nu_{bot}$ ) $Ra$ बढ़ने के साथ घटती है, जो एक प्रति-सहज (counter-intuitive) प्रवृत्ति है जिसे मजबूत थर्मल स्ट्रैटिफिकेशन और निचले कॉर्नर एडीज में कोल्ड जेट्स के फंसने के कारण जिम्मेदार माना गया है।

महत्व और दावे
पेपर यह दावा करता है कि यह कृत्रक समरूपता बाधाओं के बिना, केंद्रीय रूप से गर्म सिलेंडर विन्यास में स्थिर से अराजक प्राकृतिक संवहन के पूर्ण संक्रमण अनुक्रम का पहला व्यापक संख्यात्मक लक्षण वर्णन प्रदान करता है। प्रमुख योगदानों में शामिल हैं:

महत्वपूर्ण रेले संख्या: असममित मोड को विकसित होने की अनुमति देकर $Ra_c$ के लिए एक अधिक सटीक रेंज स्थापित करना, जो समरूपता बाधाओं के कारण पिछले अति-अनुमानों को सुधारता है।
अराजकता का मार्ग: यह पुष्टि करना कि संक्रमण रुएले-टाकेन्स-न्यूहाउस मार्ग (लिमिट साइकिल $\to$ टू-टोरस $\to$ स्ट्रेंज अट्रैक्टर) का अनुसरण करता है, न कि पीरियड-डबलिंग कैस्केड का।
विकार का स्थानिक प्रसार: यह प्रदर्शित करना कि अराजकता थर्मल प्लूम में शुरू होती है और नीचे की ओर प्रसारित होती है, जहाँ निचला स्ट्रैटिफाइड पूल उच्चतम रेले संख्याओं तक अपेक्षाकृत शांत रहता है।
संदर्भ डेटा: उन विन्यासों के लिए महत्वपूर्ण रेले संख्या, शासन संक्रमण और वैश्विक ऊष्मा स्थानांतरण सहसंबंधों के लिए विश्वसनीय बेंचमार्क डेटा प्रदान करना जहाँ आंतरिक ज्यामिति और अनस्थिर गतिशीलता परस्पर क्रिया करती है।

लेखक नोट करते हैं कि यह अध्ययन 2D, एकल सिलेंडर त्रिज्या ( $R=0.2L$ ), और $Pr=0.71$ तक सीमित है, जो भविष्य के दिशाओं के रूप में 3D, ऑफ-सेंटर पोजीशन और विविध प्रैंड्टल नंबरों की ओर विस्तार का सुझाव देता है।