गैर-आदर्श संपीड्य द्रव गतिविज्ञान के लिए जालक बोल्ट्जमैन मॉडल

मूल जानकारी

• पेपर ID: 2510.14712
• शीर्षक: गैर-आदर्श संपीड्य द्रव गतिविज्ञान के लिए जालक बोल्ट्जमैन मॉडल
• लेखक: S. A. Hosseini, M. Feinberg, I. V. Karlin (ETH Zurich)
• वर्गीकरण: physics.flu-dyn
• प्रकाशित पत्रिका: Journal of Fluid Mechanics
• पेपर लिंक: https://arxiv.org/abs/2510.14712

सारांश

यह पेपर संपीड्य गैर-आदर्श द्रव प्रवाह के अनुकरण के लिए एक जालक बोल्ट्जमैन विधि प्रस्तावित करता है। यह विधि प्रथम निकटतम जालक का उपयोग करती है, अर्ध-संतुलन आकर्षक के माध्यम से सुसंगत सुधार पदों का एक समूह प्रस्तुत करती है, जो सकारात्मक निश्चितता और गैलिलियन-अपरिवर्तनीय Navier-Stokes विसरण दर सुनिश्चित करती है। यह निर्माण विस्तारित टेम्पलेट या अस्थायी नियमितीकरण की आवश्यकता से बचता है, जबकि प्रवाह स्थितियों की एक विस्तृत श्रृंखला में संख्यात्मक स्थिरता और ऊष्मागतिक सुसंगतता बनाए रखता है। परिणामी मॉडल यूलर और Navier-Stokes पदानुक्रम स्तरों पर द्रव गतिविज्ञान को सटीकता से पुनः प्रस्तुत करता है। कठोर सत्यापन के रूप में, लेखकों ने पहली बार जालक बोल्ट्जमैन ढांचे के भीतर मैक संख्या 1.47 तक की बूंद-आघात परस्पर क्रिया के मात्रात्मक रूप से सटीक अनुकरण का प्रदर्शन किया है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

गैर-आदर्श संपीड्य द्रव गतिविज्ञान द्रव यांत्रिकी की एक उभरती और तेजी से विकसित होने वाली शाखा है, मुख्य रूप से निकट-महत्वपूर्ण, अतिक्रांतिक और अतिसंकट अवस्थाओं में काम करने वाली विधियों और तकनीकों के उदय के कारण। इस प्रकार के द्रव में शामिल हैं:

1. अतिसंकट द्रव: महत्वपूर्ण बिंदु से ऊपर दबाव और तापमान स्थितियों में द्रव
2. सघन वाष्प: संपीड़न कारक जो आदर्श गैस व्यवहार से विचलित होते हैं
3. BZT द्रव: नकारात्मक मौलिक गैस गतिविज्ञान व्युत्पन्न वाले उच्च बहुलक जटिलता द्रव

महत्व

1. औद्योगिक अनुप्रयोग की आवश्यकता: जैविक रैंकिन चक्र, अतिसंकट CO₂ टर्बाइन आदि ऊर्जा प्रौद्योगिकियां
2. प्रायोगिक डेटा की कमी: प्रायोगिक सेटअप जटिल हैं, डेटा प्राप्त करना कठिन है
3. संख्यात्मक उपकरण की आवश्यकता: जटिल भौतिक घटनाओं को समझने के लिए सुसंगत और कुशल संख्यात्मक उपकरणों की आवश्यकता है

मौजूदा विधियों की सीमाएं

1. पारंपरिक CFD विधियां: ऊष्मागतिक रूप से अस्थिर स्पिनोडल क्षेत्र के भीतर, रुद्धोष्म ध्वनि गति वर्ग ऋणात्मक हो जाता है, जिससे विकास समीकरण अतिपरवलयिक से परवलयिक में बदल जाता है
2. मौजूदा LBM मॉडल: मुख्य रूप से उबलते अनुप्रयोगों के लिए, संपीड्य गैर-आदर्श प्रवाह के लिए व्यवस्थित मॉडल की कमी
3. संख्यात्मक स्थिरता समस्याएं: उच्च मैक संख्या और मजबूत गैर-आदर्श प्रभावों के तहत स्थिरता चुनौतियां

मुख्य योगदान

1. नया द्वि-वितरण फलन LBM मॉडल प्रस्तावित किया: गैर-आदर्श संपीड्य प्रवाह अनुकरण के लिए शास्त्रीय प्रथम निकटतम जालक (D3Q27) का उपयोग करके
2. अर्ध-संतुलन आकर्षक तंत्र प्रस्तुत किया: fᵢ और gᵢ सुधार पदों के माध्यम से सकारात्मक निश्चितता और गैलिलियन-अपरिवर्तनीय विसरण दर सुनिश्चित करने के लिए
3. स्वतंत्र आयतन श्यानता नियंत्रण प्राप्त किया: BGK संरचना में आयतन श्यानता संभवतः नकारात्मक होने की समस्या को हल किया
4. पहली बार उच्च मैक संख्या बूंद-आघात परस्पर क्रिया अनुकरण प्राप्त किया: LBM ढांचे के भीतर Ma=1.47 का मात्रात्मक रूप से सटीक अनुकरण
5. ऊष्मागतिक सुसंगतता बनाए रखी: प्रवाह स्थितियों की विस्तृत श्रृंखला में संख्यात्मक स्थिरता

विधि विवरण

कार्य परिभाषा

Navier-Stokes-Korteweg (NSK) समीकरण समूह को संतुष्ट करने वाले संपीड्य गैर-आदर्श द्रव का अनुकरण:

• द्रव्यमान संरक्षण: ∂ₜρ + ∇·(ρu) = 0
• गति संरक्षण: ∂ₜ(ρu) + ∇·(ρu⊗u) + ∇P + ∇·T^NS + ∇·T^K = 0
• ऊर्जा संरक्षण: ∂ₜ(ρE) + ∇·((ρE+P)u) + u·(∇·T^K) + ∇·(u·T^NS) + ∇·q = 0

जहां T^K Korteweg सतह तनाव टेंसर है, T^NS श्यान प्रतिबल टेंसर है।

मॉडल आर्किटेक्चर

1. द्वि-वितरण फलन डिजाइन

दो वितरण फलन समूह {fᵢ, gᵢ}, i=1,...,27 का उपयोग:

• fᵢ: द्रव्यमान और गति संरक्षण को संभालता है
• gᵢ: ऊर्जा संरक्षण को संभालता है

2. जालक बोल्ट्जमैन समीकरण

fᵢ(x+cᵢδt, t+δt) = fᵢ + 2β(fᵢᵉᵠ - fᵢ) + (1-β)(fᵢ* - fᵢᵉᵠ)
gᵢ(x+cᵢδt, t+δt) = gᵢ + 2β(gᵢᵉᵠ - gᵢ) + (1-β)(gᵢ* - gᵢᵉᵠ)

जहां β = Pδt/(2μ + Pδt) शिथिलन पैरामीटर है।

3. संतुलन वितरण फलन

उत्पाद रूप अपनाया:

fᵢᵉᵠ = ρ ∏_α Ψᵢα(uα, P/ρ + uα²)

जहां Ψᵢα आधार फलन है:

Ψᵢα(ξα, ζαα) = 1 - cᵢα² + ½[(3cᵢα² - 2)ζαα + cᵢαξα]

4. अर्ध-संतुलन सुधार पद

f वितरण फलन सुधार:

• बल पद: Fα = -κρ∂α∇²ρ (Korteweg बल)
• प्रतिबल सुधार: Φαα विकर्ण प्रतिबल टेंसर सुधार और स्वतंत्र आयतन श्यानता नियंत्रण सम्मिलित करता है

g वितरण फलन सुधार:

• स्थानांतरित वेग: u* = u + δt F/ρ
• स्थानांतरित तापमान: T* = T - δt²F·F/(2ρ²cᵥ)
• ऊष्मा प्रवाह सुधार: गैर-आदर्श गैस के एन्थैल्पी प्रवणता प्रभाव पर विचार

तकनीकी नवाचार बिंदु

1. अर्ध-संतुलन आकर्षक तंत्र

fᵢ और gᵢ प्रस्तुत करके, प्रणाली केवल स्थानीय संतुलन की ओर शिथिल नहीं होती, बल्कि संशोधित अर्ध-संतुलन अवस्था की ओर शिथिल होती है, जो सुनिश्चित करता है:

• सही मैक्रोस्कोपिक गतिविज्ञान
• संख्यात्मक स्थिरता
• गैलिलियन अपरिवर्तनीयता

2. स्वतंत्र आयतन श्यानता नियंत्रण

पारंपरिक BGK मॉडल में आयतन श्यानता η' = ((D+2)/D - ρcs²/P)μ नकारात्मक हो सकता है। सुधार पद Φ' के माध्यम से स्वतंत्र नियंत्रण प्राप्त किया:

Φ' = P((D+2)/D - ρcs²/P - η/μ)(∇·u)

3. ऊष्मागतिक सुसंगतता

सावधानीपूर्वक डिजाइन किए गए स्थानांतरित तापमान T* और ऊर्जा सुधार के माध्यम से, गैर-रैखिक अवस्था समीकरण के तहत ऊष्मागतिक सुसंगतता सुनिश्चित की।

प्रायोगिक सेटअप

सत्यापन केस

1. मूल सत्यापन:
   • ध्वनि गति माप
   • कतरनी श्यानता सत्यापन
   • तापीय Couette प्रवाह
   • सामान्य मोड विसरण
2. बहु-चरण प्रवाह सत्यापन:
   • द्रव-वाष्प सह-अस्तित्व घनत्व
   • इंटरफेस सुसंगतता और अभिसरण
3. संपीड्य प्रवाह सत्यापन:
   • एक-आयामी गैर-आदर्श आघात ट्यूब
   • आघात-द्रव स्तंभ परस्पर क्रिया

पदार्थ पैरामीटर

मुख्य रूप से नाइट्रोजन (N₂) का उपयोग करके सत्यापन:

• महत्वपूर्ण दबाव: Pc = 3.4×10⁶ Pa
• महत्वपूर्ण घनत्व: ρc = 241.96 kg/m³
• महत्वपूर्ण तापमान: Tc = 126.2 K
• R/cv = 0.4

अवस्था समीकरण

van der Waals अवस्था समीकरण अपनाया:

P(ρ,T) = ρRT/(1-bρ) - aρ²

प्रायोगिक परिणाम

मुख्य परिणाम

1. मूल भौतिक राशि सत्यापन

• ध्वनि गति: संतृप्त द्रव और वाष्प शाखाओं पर, Tr∈0.7,1 श्रेणी में विश्लेषणात्मक समाधान के साथ पूर्ण मेल
• कतरनी श्यानता: विभिन्न मैक संख्याओं पर गैलिलियन अपरिवर्तनीयता बनाए रखता है, त्रुटि <1%
• तापीय Couette प्रवाह: Pr∈{0.6,1.2,4.9} और Ma∈{0.8,1.2,1.6} पर विश्लेषणात्मक समाधान के साथ सुसंगत

2. बहु-चरण प्रवाह गुण

• सह-अस्तित्व घनत्व: Tr∈0.3,0.99 श्रेणी में Maxwell समान क्षेत्र निर्माण के साथ पूर्ण सुसंगतता
• इंटरफेस अभिसरण: विभिन्न जालक संकल्पों पर अच्छी अभिसरण विशेषताएं प्रदर्शित करता है

3. गैर-आदर्श आघात ट्यूब

तीन विभिन्न गैर-आदर्श आघात ट्यूब कॉन्फ़िगरेशन का सफल अनुकरण:

• केस I: Γ<0 क्षेत्र में विरल आघात
• केस II: शास्त्रीय संपीड़न तरंग व्यवहार
• केस III: Γ=0 रेखा को गतिशील रूप से पार करने वाली मिश्रित विरल तरंग

सभी केस संदर्भ डेटा के साथ उत्कृष्ट मेल दिखाते हैं।

4. आघात-बूंद परस्पर क्रिया

• पहली बार LBM ढांचे के भीतर Ma=1.47 का उच्च मैक संख्या अनुकरण प्राप्त किया
• घटना आघात, संचारित तरंग, परावर्तित आघात, मैक रॉड आदि जटिल तरंग संरचनाओं को सटीकता से पकड़ा
• बूंद विरूपण विकास प्रायोगिक डेटा के साथ अत्यधिक सुसंगत

विलोपन प्रयोग

विभिन्न जालक संकल्पों के साथ इंटरफेस परीक्षण के माध्यम से मॉडल की अभिसरण सत्यापित:

• δx से 5μm से 0.1μm तक चार संकल्प स्तर
• अच्छी जालक अभिसरण विशेषताएं प्रदर्शित करता है

मुख्य निष्कर्ष

1. संख्यात्मक स्थिरता: उच्च मैक संख्या और मजबूत गैर-आदर्श प्रभावों के तहत स्थिर रहता है
2. कम्प्यूटेशनल दक्षता: मानक प्रथम निकटतम जालक का उपयोग, नियंत्रणीय कम्प्यूटेशनल लागत
3. भौतिक निष्ठा: गैर-आदर्श द्रव की जटिल भौतिक घटनाओं को सटीकता से पुनः प्रस्तुत करता है

संबंधित कार्य

पारंपरिक LBM बहु-चरण प्रवाह मॉडल

1. रंग प्रवणता मॉडल (Gunstensen et al. 1991)
2. छद्म-संभावित मॉडल (Shan & Chen 1993)
3. मुक्त ऊर्जा मॉडल (Swift et al. 1996)

संपीड्य LBM विकास

• He & Doolen (2002): संपीड्य गैर-आदर्श प्रवाह का पहला प्रयास
• मौजूदा मॉडल मुख्य रूप से उबलते अनुप्रयोगों के लिए
• संपीड्य गैर-आदर्श प्रवाह मॉडल की व्यवस्थित कमी

इस पेपर के सापेक्ष लाभ

1. व्यवस्थित विधि: संपूर्ण द्वि-वितरण फलन ढांचा
2. सैद्धांतिक आधार: कठोर बहु-पैमाने विश्लेषण
3. सत्यापन पूर्णता: मूल भौतिक राशियों से जटिल प्रवाह घटनाओं तक व्यापक सत्यापन

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

1. संपीड्य गैर-आदर्श द्रव के लिए उपयुक्त LBM मॉडल का सफल विकास
2. अर्ध-संतुलन आकर्षक तंत्र के माध्यम से संख्यात्मक स्थिरता और भौतिक सुसंगतता प्राप्त
3. पहली बार LBM ढांचे के भीतर उच्च मैक संख्या बूंद-आघात परस्पर क्रिया अनुकरण
4. मॉडल प्रवाह स्थितियों की विस्तृत श्रृंखला में ऊष्मागतिक सुसंगतता बनाए रखता है

सीमाएं

1. एकल शिथिलन समय मॉडल: वर्तमान में सरल BGK टकराव मॉडल का उपयोग, अशांत अनुकरण के लिए पर्याप्त मजबूत नहीं हो सकता
2. अवस्था समीकरण प्रतिबंध: मुख्य रूप से van der Waals समीकरण सत्यापित, अन्य जटिल अवस्था समीकरणों को आगे सत्यापन की आवश्यकता
3. त्रि-आयामी सत्यापन सीमित: अधिकांश सत्यापन एक-आयामी और द्वि-आयामी कॉन्फ़िगरेशन पर केंद्रित

भविष्य की दिशाएं

1. अधिक मजबूत टकराव मॉडल तक विस्तार: अशांत कॉन्फ़िगरेशन अनुकरण के लिए
2. अधिक जटिल अवस्था समीकरण: जैसे Peng-Robinson आदि
3. बड़े पैमाने पर त्रि-आयामी अनुप्रयोग: वास्तविक इंजीनियरिंग समस्याओं का अनुप्रयोग

गहन मूल्यांकन

शक्तियां

1. सैद्धांतिक कठोरता: संपूर्ण बहु-पैमाने विश्लेषण और गणितीय व्युत्पत्ति
2. विधि नवाचार: अर्ध-संतुलन आकर्षक तंत्र महत्वपूर्ण नवाचार है
3. सत्यापन व्यापकता: मूल भौतिकी से जटिल अनुप्रयोगों तक व्यवस्थित सत्यापन
4. व्यावहारिक मूल्य: महत्वपूर्ण इंजीनियरिंग अनुप्रयोग समस्याओं को हल करता है
5. संख्यात्मक दक्षता: मानक जालक का उपयोग कम्प्यूटेशनल जटिलता से बचता है

कमियां

1. टकराव मॉडल सरल: BGK मॉडल जटिल प्रवाह के लिए सीमित प्रयोज्यता
2. पैरामीटर निर्भरता: कुछ सुधार पैरामीटर के चयन में व्यवस्थित मार्गदर्शन की कमी
3. त्रि-आयामी सत्यापन अपर्याप्त: अधिक त्रि-आयामी जटिल ज्यामिति सत्यापन की आवश्यकता

प्रभाव

1. शैक्षणिक योगदान: LBM को संपीड्य गैर-आदर्श प्रवाह क्षेत्र में नई दिशा खोलता है
2. इंजीनियरिंग मूल्य: अतिसंकट द्रव प्रौद्योगिकी के लिए महत्वपूर्ण संख्यात्मक उपकरण
3. पद्धति महत्व: अर्ध-संतुलन आकर्षक अवधारणा अन्य संख्यात्मक विधियों को प्रभावित कर सकती है

प्रयोज्य परिदृश्य

1. अतिसंकट द्रव प्रौद्योगिकी: CO₂ टर्बाइन, जैविक रैंकिन चक्र
2. बहु-चरण प्रवाह: बूंद विभाजन, चरण परिवर्तन प्रक्रियाएं
3. आघात गतिविज्ञान: उच्च गति प्रवाह में जटिल तरंग परस्पर क्रिया
4. प्रत्यक्ष संख्यात्मक अनुकरण: उच्च सटीकता की आवश्यकता वाली मूल भौतिकी अनुसंधान

संदर्भ

1. Guardone, A. et al. (2024). गैर-आदर्श संपीड्य द्रव गतिविज्ञान सघन वाष्प और अतिसंकट द्रव का। Annual Review of Fluid Mechanics, 56, 241-269.
2. Hosseini, S.A. & Karlin, I.V. (2023). गैर-आदर्श द्रव के लिए जालक बोल्ट्जमैन: मौलिक सिद्धांत और व्यवहार। Physics Reports, 1030, 1-137.
3. Shan, X. & Chen, H. (1993). कई चरणों और घटकों के साथ प्रवाह का अनुकरण करने के लिए जालक बोल्ट्जमैन मॉडल। Physical Review E, 47(3), 1815.