🔭 astrophysics

Nonlinear evolution of unstable solar inertial modes: The case of viscous modes on a differentially rotating sphere

यह शोध पत्र एक विभेदक घूर्णन वाले गोले पर सूर्य के सबसे प्रमुख उच्च-अक्षांश जड़त्व मोड ( $m=1$ ) के गैररेखीय विकास की जांच करता है, जो प्रत्यक्ष संख्यात्मक सिमुलेशन के माध्यम से यह प्रदर्शित करता है कि यह एक सुपरक्रिटिकल हॉपफ बाइफर्केशनेशन (supercritical Hopf bifurcation) से गुजरता है जहाँ रेनॉल्ड्स स्ट्रेस (Reynolds stresses) विभेदक घूर्णन को सुचारू करते हैं ताकि अस्थिरता को सौर अवलोकनों के तुलनीय आयामों पर संतृप्त किया जा सके।

मूल लेखक: Muneeb Mushtaq, Damien Fournier, Rama Ayoub, Peter J. Schmid, Laurent Gizon

प्रकाशित 2026-03-10

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

CC BY 4.0

मूल लेखक: Muneeb Mushtaq, Damien Fournier, Rama Ayoub, Peter J. Schmid, Laurent Gizon

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

एक बड़ी तस्वीर: सूर्य का "डगमगाता" नृत्य

कल्पना कीजिए कि सूर्य एक स्थिर आग का गोला नहीं है, बल्कि एक विशाल, घूमता हुआ, तरल गुब्बारा है। क्योंकि यह गैस और तरल से बना है, इसलिए यह एक ठोस चट्टान की तरह नहीं घूमता (जहाँ ऊपर और नीचे एक ही गति से चलते हैं)। इसके बजाय, सूर्य का भूमध्य रेखा (equator) ध्रुवों (poles) की तुलना में तेज़ी से घूमता है। इसे विभेदक घूर्णन (differential rotation) कहा जाता है।

इस तरल नृत्य में, अदृश्य लहरें होती हैं जिन्हें इनर्शियल मोड्स (inertial modes) कहा जाता है। इन्हें आप उस लहरों की तरह समझें जो तब दिखाई देती हैं जब आप जेली (Jell-Jello) के एक कटोरे को हिलाते हैं। एक विशिष्ट लहर, जिसे $m=1$ मोड कहा जाता है, सबसे तेज़ और सबसे ऊर्जावान है जो हम सूर्य पर देखते हैं। यह एक उच्च-अक्षांश वाली लहर है जो ध्रुवों के पास आगे-पीछे चलती रहती है।

समस्या: यह लहर बड़ी और बड़ी क्यों नहीं होती जाती?

सरल भौतिकी (रैखिक सिद्धांत/linear theory) के अनुसार, यह विशिष्ट लहर अस्थिर (unstable) होनी चाहिए। क्योंकि सूर्य भूमध्य रेखा पर तेज़ी से और ध्रुवों पर धीरे घूमता है, यह एक 'शीयर फोर्स' (shear force) की तरह है जो इस लहर में ऊर्जा भरता रहता है।

यदि आप एक झूले को बिल्कुल सही समय पर धक्का देते हैं, तो वह ऊँचा और ऊँचा जाता है। सूर्य के मामले में, विभेदक घूर्णन लगातार इस लहर को "धक्का" दे रहा है। सैद्धांतिक रूप से, लहर को तब तक बढ़ना चाहिए जब तक कि वह सूर्य को फाड़ न दे या एक विशाल तूफान न पैदा कर दे। लेकिन वास्तव में, लहर एक स्थिर, प्रबंधनीय आकार (लगभग 10 मीटर प्रति सेकंड) में रहती है।

प्रश्न: क्या चीज़ लहर को हमेशा के लिए बढ़ने से रोकती है? वह "ब्रेक" क्या है जो इसे नियंत्रण में रखता है?

प्रयोग: एक डिजिटल सैंडबॉक्स

लेखकों ने इसका अनुकरण (simulate) करने के लिए एक कंप्यूटर मॉडल बनाया। उन्होंने सूर्य को एक 2D घूमते हुए गोले में सरल बना दिया (फिलहाल जटिल 3D ताप और चुंबकीय क्षेत्रों को अनदेखा करते हुए) और देखा कि क्या होता है जब वे इस "अस्थिरता" को चालू करते हैं।

उन्होंने पाया कि लहर वास्तव में बढ़ती तो है, लेकिन फिर यह एक सीमा तक पहुँचती है और स्थिर हो जाती है। इसे संतृप्ति (saturation) कहा जाता है।

उपमा: पहाड़ी पर कार

इस लहर को एक पहाड़ी पर चढ़ने वाली कार की तरह समझें।

इंजन (अस्थिरता): विभेदक घूर्णन इंजन है, जो कार को पहाड़ी पर ऊपर धकेल रहा है।
पहाड़ी (संतृप्ति): जैसे-जैसे कार तेज़ होती है, वह ज़मीन पर अपने पैर घसीटने लगती है (घर्षण/friction)।
परिणाम: अंततः, इंजन का धक्का घर्षण के खिंचाव के बराबर हो जाता है। कार का त्वरण रुक जाता है और वह एक स्थिर गति से चलती रहती है।

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

गणित: "लैंडौ" समीकरण (The "Landau" Equation)

यह शोध पत्र इस व्यवहार का वर्णन करने के लिए एक प्रसिद्ध गणितीय उपकरण, स्टुअर्ट-लैंडौ समीकरण (Stuart-Landau equation) का उपयोग करता है।

सरल उपमा: कल्पना कीजिए कि एक सूत्र कहता है:
- विकास दर (Growth Rate) = (आप कितना बढ़ना चाहते हैं) - (आप पहले से कितना बढ़ चुके हैं)।
यदि आप छोटे हैं, तो आप तेज़ी से बढ़ते हैं।
यदि आप बहुत बड़े हो जाते हैं, तो "ब्रेकिंग" शब्द (braking term) ज़ोर से काम करता है, और आप बढ़ना बंद कर देते हैं।

लेखकों ने सिद्ध किया कि इस प्रकार का "रुकना और चलना" वाला व्यवहार एक सुपरक्रिटिकल बिफर्केशन (Supercritical Bifurcation) है।

सुपरक्रिटिकल (अच्छे प्रकार का): यदि आप सिस्टम को थोड़ा सा हिलाते हैं, तो यह धीरे से एक नए, स्थिर आकार तक बढ़ता है। यह एक स्प्रिंग की तरह है जो एक नई लंबाई तक खिंचता है और वहीं रहता है।
सबक्रिटिकल (बुरे प्रकार का): यदि आप इसे हिलाते हैं, तो यह अराजकता (chaos) में बदल सकता है या पूरी तरह से अलग स्थिति में कूद सकता है।

सूर्य की लहरें "सुपरक्रिटिकल" हैं। वे सुव्यवस्थित हैं। वे बढ़ती हैं, एक सीमा तक पहुँचती हैं, और सुचारू रूप से स्थिर हो जाती हैं।

"गूँज": हार्मोनिक्स (Harmonics)

जब मुख्य लहर ( $m=1$ ) मज़बूत होती है, तो वह केवल अकेली नहीं रहती। यह "गूँज" या हार्मोनिक्स ( $m=2, m=3$ , आदि) बनाना शुरू कर देती है।

उपमा: गिटार के तार के बारे में सोचें। जब आप मुख्य स्वर (fundamental note) को बजाते हैं, तो आप उच्च-पिच वाले ओवरटोन्स (overtones) भी सुनते हैं।
लेखकों ने पाया कि ये ओवरटोन्स मुख्य लहर की तुलना में बहुत छोटे हैं। दूसरी गूँज लगभग 8 गुना छोटी है, और तीसरी 25 गुना छोटी है।
महत्वपूर्ण रूप से, लेखकों ने दिखाया कि ये "गूँज" यादृच्छिक शोर (random noise) नहीं हैं; ये पूरी तरह से व्यवस्थित पैटर्न हैं जिन्हें मुख्य लहर अस्तित्व में लाती है।

"वास्तविक दुनिया" से संबंध

टीम ने अपना सिमुलेशन एक श्यानता (viscosity/चिपचिपाहट) के साथ चलाया जो हमारे सूर्य की सतह पर मौजूद होने की संभावना है (जो सुपरग्रैन्यूल्स द्वारा उत्पन्न होती है, जो उठते हुए गैस के विशाल बुलबुलों की तरह हैं)।

परिणाम: लहर 28 मीटर प्रति सेकंड की गति पर स्थिर हुई।
यह उल्लेखनीय रूप से वही है जो खगोलविद वास्तव में सूर्य पर देखते हैं (जो लगभग 10-20 मीटर/सेकंड के आसपास है)। यह सुझाव देता है कि "लहरों द्वारा स्पिन को सुचारू करने" का सरल भौतिकी ही एक प्रमुख कारण है कि सूर्य की लहरें अनियंत्रित क्यों नहीं होती हैं।

एक चेतावनी (लेकिन...)

लेखक सावधान करते हैं कि: "इसका अत्यधिक अर्थ न निकालें।"
उनका मॉडल 2D (समतल) है। वास्तविक सूर्य 3D है (गहराई वाला एक गोला)। 3D में, भौतिकी अधिक जटिल है क्योंकि ताप (heat) और एंट्रॉपी (entropy) एक बड़ी भूमिका निभाते हैं। 2D मॉडल एक सरलीकृत "खिलौना मॉडल" (toy model) है जो व्यवहार के सार को पकड़ता है, लेकिन वास्तविक सूर्य एक अधिक जटिल नृत्य है।

एक वाक्य में सारांश

सूर्य की सबसे बड़ी लहरें तब तक बढ़ती हैं जब तक कि वे उस गति के अंतर को "सुचारू" नहीं कर देतीं जो उन्हें ईंधन देते हैं, जिससे एक आदर्श, स्थिर संतुलन बनता है जिसे एक सरल गणितीय नियम द्वारा अनुमानित किया जा सकता है, ठीक वैसे ही जैसे एक कार एक समतल सड़क पर अपनी शीर्ष गति खोज लेती है।

1. समस्या का विवरण (Problem Statement)

यह शोध पत्र सूर्य के उच्च-अक्षांश जड़त्व मोड (high-latitude inertial mode) के अनुदैर्ध्य वेव-संख्या (longitudinal wavenumber) $m=1$ के गैररेखीय विकास (nonlinear evolution) की जांच करता है।

संदर्भ: अवलोकन दर्शाते हैं कि यह मोड सूर्य पर सबसे बड़े आयाम (RMS वेग $>10$ m/s) वाला है। रैखिक सिद्धांत (linear theory) भविष्यवाणी करता है कि यह अक्षांशीय विभेदक घूर्णन (latitudinal differential rotation - तीव्र भूमध्य रेखा, धीमे ध्रुव) के कारण होने वाली कतरनी अस्थिरता (shear instability) के कारण अस्थिर है।
अंतराल (Gap): जबकि रैखिक स्थिरता विश्लेषण इस अस्थिरता की पहचान करता है, यह देखे गए परिमित आयामों (finite amplitudes) की व्याख्या करने में असमर्थ है। सौर संवहन के पूर्णतः गैररेखीय 3D सिमुलेशन गणनात्मक रूप से अत्यधिक महंगे हैं और विशिष्ट मोड अंतःक्रियाओं के संबंध में व्याख्या करना कठिन है।
उद्देश्य: एक सरलीकृत 2D टॉरॉइडल मॉडल में इस अस्थिरता के संतृप्ति तंत्र (saturation mechanism) को संख्यात्मक और सैद्धांतिक रूप से लक्षणबद्ध करना, संतुलन आयाम, उच्च हार्मोनिक्स की भूमिका और दुर्बल रूप से गैररेखीय (weakly nonlinear - WNL) सिद्धांत की वैधता निर्धारित करना।

2. कार्यप्रणाली (Methodology)

लेखक प्रत्यक्ष संख्यात्मक सिमुलेशन (Direct Numerical Simulations - DNS) और दो विचलित (perturbative) दुर्बल रूप से गैररेखीय (WNL) सिद्धांतों के संयोजन का उपयोग करते हैं।

A. भौतिक मॉडल और शासी समीकरण (Physical Model & Governing Equations)

डोमेन: एक 2D विभेदक घूर्णन गोला (सतही परत)।
प्रवाह: विशुद्ध रूप से टॉरॉइडल गतियाँ (असंपीड्य श्यान प्रवाह/incompressible viscous flow)।
समीकरण: रेडियल भंवर समीकरण (radial vorticity equation) को संवेग समीकरण से व्युत्पन्न किया गया है, जिसमें शामिल हैं:
- कोरिओलिस बल (Coriolis force)।
- HMI हेलियोसेस्मिक डेटा से अनुमानित अक्षांशीय विभेदक घूर्णन ( $\Omega_0(\theta)$ )।
- श्यानता (viscosity) (Ekman number $E$ के माध्यम से मॉडल की गई)।
- $\Lambda$ -प्रभाव (गैर-विसारी कोणीय संवेग परिवहन) ताकि प्रारंभिक विभेदक घूर्णन प्रोफाइल को बनाए रखा जा सके।
नियंत्रण पैरामीटर: Ekman number $E = \nu / (r^2 \Omega_{ref})$ । अध्ययन $10^{-3} \lesssim E < E_c \approx 1.5 \times 10^{-3}$ की सीमा पर केंद्रित है।

B. संख्यात्मक सिमुलेशन (DNS)

सॉल्वर: समय डोमेन में डायरेक्ट न्यूमेरिकल सिमुलेशन का उपयोग करते हुए 3rd-order Adams–Bashforth टेम्पोरल स्कीम और स्थानिक विविक्तकरण (spatial discretization) के लिए गोलाकार हार्मोनिक अपघटन (degree $L_{max}=200$ तक) का उपयोग किया गया है।
आरंभीकरण (Initialization): सौर जैसा विभेदक घूर्णन प्रोफाइल + छोटा यादृच्छिक शोर (random noise)।
प्रक्रिया: सिमुलेशन तब तक चलते हैं जब तक कि प्रणाली एक स्थिर संतुलन अवस्था (संतृप्ति) तक नहीं पहुँच जाती।

C. दुर्बल रूप से गैररेखीय (WNL) सिद्धांत

दो शास्त्रीय विचलन विधियों का उपयोग स्टुअर्ट-लैंडौ आयाम समीकरण (Stuart-Landau amplitude equation) को व्युत्पन्न करने के लिए किया जाता है:
$\frac{d|A|}{dt} = \sigma_I |A| + \beta_I |A|^3$

बहु-पैमाना विस्तार (Multiscale Expansion): एक छोटे पैरामीटर $\epsilon$ का उपयोग करता है जो क्रिटिकल Ekman number ( $E_c - E$ ) से दूरी को दर्शाता है। समय को तेज़ (दोलन) और धीमी (आयाम विकास) स्केल्स में विभाजित किया गया है।
आयाम विस्तार (Amplitude Expansion): स्वयं मोड आयाम $A$ को एक छोटे पैरामीटर के रूप में उपयोग करता है।

दोनों विधियों का लक्ष्य रैखिक विकास दर ( $\sigma_I$ ) और गैररेखीय संतृप्ति गुणांक ( $\beta_I$ ) की गणना करना है।

3. मुख्य योगदान (Key Contributions)

सुपरक्रिटिकल हॉफ द्विभाजन (Supercritical Hopf Bifurcation) का सत्यापन: अध्ययन पुष्टि करता है कि $m=1$ मोड एक सुपरक्रिटिकल हॉफ द्विभाजन से गुजरता है। जैसे ही $E$ , $E_c$ से नीचे गिरता है, प्रणाली एक स्थिर अवस्था से एक परिमित-आयाम सीमा चक्र (limit cycle) में सुचारू रूप से परिवर्तित होती है।
संतृप्ति स्केलिंग का व्युत्पन्न: लेखक यह व्युत्पन्न करते हैं कि संतृप्त आयाम $\sigma_I^{1/2}$ (या $(E_c - E)^{1/2}$ ) के रूप में स्केल करता है, जो लैंडौ समीकरण के अनुरूप है जहाँ $\beta_I < 0$ है।
हार्मोनिक्स का लक्षण वर्णन: शोध पत्र प्रदर्शित करता है कि मौलिक $m=1$ $m = 1$ मोड की गैररेखीय अंतःक्रिया उच्च हार्मोनिक्स ( $m=2, 3, \dots$ $m = 2, 3, \dots$ ) उत्पन्न करती है। इन्हें एसिम्प्टोटिक कूपमैन मोड (asymptotic Koopman modes) के रूप में पहचाना गया है।
- $m=2$ हार्मोन उत्तर-दक्षिण (NS) प्रति-सममित (antisymmetric) है।
- $m=3$ हार्मोन उत्तर-दक्षिण (NS) सममित (symmetric) है (मौलिक मोड के समान)।
- इनके आयाम $|A|^m$ के रूप में स्केल होते हैं।
रेनॉल्ड्स स्ट्रेस तंत्र: संतृप्ति तंत्र को रेनॉल्ड्स स्ट्रेस (Reynolds stress) के रूप में पहचाना गया है, जो विभेदक घूर्णन पर फीडबैक देकर उत्पन्न होता है, जिससे कतरनी प्रोफाइल (shear profile) को सुचारू बनाया जाता है जब तक कि रैखिक विकास दर शून्य न हो जाए।

4. मुख्य परिणाम (Key Results)

A. संख्यात्मक निष्कर्ष (DNS)

द्विभाजन प्रकार (Bifurcation Type): $E < E_c$ के लिए, प्रणाली प्रारंभिक विक्षोभ की शक्ति की परवाह किए बिना एक अद्वितीय संतुलन अवस्था की ओर अभिसरित होती है, जो एक सुपरक्रिटिकल द्विभाजन की पुष्टि करती है।
आयाम: $E = 4 \times 10^{-4}$ (सौर सुपरग्रैनुलर श्यानता के अनुरूप) पर, संतृप्त RMS वेग 28 m/s तक पहुँच जाता है, जो देखे गए सौर मानों ( $\sim 10$ m/s औसत, 20 m/s तक) के तुलनीय है।
हार्मोनिक्स: संतृप्ति शासन में, $m=1$ $m = 1$ मोड हावी रहता है। हालाँकि, $m=2$ $m = 2$ और $m=3$ $m = 3$ में महत्वपूर्ण ऊर्जा स्थानांतरित होती है।
- $E = 4 \times 10^{-4}$ पर, $m=2$ का आयाम $m=1$ से $\sim 8$ गुना छोटा है, और $m=3$ का आयाम $\sim 25$ गुना छोटा है।
- जैसे-जैसे $E$ और कम होता है (कम श्यानता), कई मोड अस्थिर हो जाते हैं, और सरल एकल-मोड WNL सिद्धांत विफल हो जाता है।

B. WNL सिद्धांत के साथ तुलना

सहमति: WNL सिद्धांत (दोनों बहु-पैमाना और आयाम विस्तार) क्रिटिकल पॉइंट ( $E \approx E_c$ ) के पास संतृप्ति आयाम और स्केलिंग नियमों की सटीक भविष्यवाणी करते हैं।
गुणांक: DNS से प्राप्त लैंडौ गुणांक ( $\sigma_I, \beta_I$ ) सैद्धांतिक भविष्यवाणियों से मेल खाते हैं।
स्थानिक संरचना: WNL सिद्धांत मौलिक मोड और पहले दो हार्मोनिक्स के स्थानिक आइगेनफंक्शंस (eigenfunctions) को सफलतापूर्वक पुनरुत्पादित करता है। $E$ के $E_c$ से दूर जाने पर WNL और DNS के बीच विसंगति बढ़ती है।
विभेदक घूर्णन: सिद्धांत सही ढंग से रेनॉल्ड्स स्ट्रेस द्वारा प्रेरित बैकग्राउंड रोटेशन प्रोफाइल ( $\Delta \Omega$ ) के संशोधन की भविष्यवाणी करता है, जो DNS के साथ अच्छी सहमति दिखाता है।

5. महत्व और निहितार्थ (Significance and Implications)

भौतिक अंतर्दृष्टि: यह अध्ययन पूर्ण 3D संवहन सिमुलेशन की आवश्यकता के बिना सौर जड़त्व मोड के संतृप्ति के लिए एक स्पष्ट हाइड्रोडायनामिक स्पष्टीकरण प्रदान करता है। यह कतरनी अस्थिरता और रेनॉल्ड्स स्ट्रेस फीडबैक को प्राथमिक चालक के रूप में अलग करता है।
अवलोकन संबंधी प्रासंगिकता:
- अनुमानित स्केलिंग ( $|A| \propto \sigma_I^{1/2}$ ) देखे गए मोड आयामों की व्याख्या करने के लिए एक मात्रात्मक ढांचा प्रदान करती है।
- हार्मोनिक्स का अस्तित्व यह सुझाव देता है कि सौर पावर स्पेक्ट्रा में मौलिक $m=1$ आवृत्ति के $2\times$ और $3\times$ आवृत्तियों पर देखे गए शिखर (peaks), जो पहले गिज़ोन एट अल. (Gizon et al.) द्वारा रिपोर्ट किए गए थे, गैररेखीय हार्मोनिक्स हो सकते हैं, न कि अलग रैखिक मोड।
सीमाएँ और भविष्य का कार्य:
- 2D मॉडल में बैरोक्लिनिक प्रभाव (baroclinic effects) और एंट्रॉपी परिवहन की कमी है, जो 3D मॉडलों में महत्वपूर्ण हैं (जहाँ अस्थिरता बैरोक्लिनिक होती है)।
- यह मॉडल अभी तक मोड आयामों के सौर चक्र मॉड्यूलेशन की व्याख्या नहीं करता है।
- भविष्य के कार्यों में कई अस्थिर मोड वाले मामलों के लिए इन उपकरणों को 3D सिमुलेशन या डायनेमिक मोड डिकंपोजिशन (DMD) पर लागू किया जाना चाहिए।

निष्कर्ष (Conclusion)

यह शोध पत्र स्थापित करता है कि सौर $m=1$ जड़त्व मोड का गैररेखीय संतृप्ति एक सुपरक्रिटिकल हॉफ द्विभाजन द्वारा नियंत्रित होता है जिसे दुर्बल रूप से गैररेखीय सिद्धांत द्वारा अच्छी तरह से वर्णित किया जा सकता है। संतृप्ति आयाम, रैखिक विकास और विभेदक घूर्णन के सुचारू बनाने (smoothing) के बीच रेनॉल्ड्स-स्ट्रेस-प्रेरित संतुलन द्वारा निर्धारित होता है। यह 2D ढांचा रैखिक सिद्धांत और देखे गए सौर आयामों के बीच सफलतापूर्वक सेतु बनाता है, जो मोड हार्मोनिक्स और उनके स्केलिंग व्यवहार के लिए एक भविष्य कहनेवाला मॉडल प्रदान करता है।