2025-11-19T00:46:13.540157

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

Lavacot, Liu, Morgan et al.
While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier--Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semi-analytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh--Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with O(100) simulations rather than O(1000) simulations. We then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.
academic

एड़ी विसरणशीलता क्षणों के परिमाणीकरण में सुधारित सांख्यिकीय अभिसरण के लिए तकनीकें

मूल जानकारी

  • पेपर ID: 2503.06418
  • शीर्षक: Techniques for improved statistical convergence in quantification of eddy diffusivity moments
  • लेखक: Dana L. O.-L. Lavacot (Stanford University), Jessie Liu (Stanford University), Brandon E. Morgan (Lawrence Livermore National Laboratory), Ali Mani (Stanford University)
  • वर्गीकरण: physics.flu-dyn
  • प्रकाशन तिथि: 15 अक्टूबर, 2025 (arXiv प्रीप्रिंट)
  • पेपर लिंक: https://arxiv.org/abs/2503.06418v2

सारांश

यह पेपर मैक्रोस्कोपिक फोर्सिंग विधि (MFM) में एड़ी विसरणशीलता गुणांक के क्षणों की सांख्यिकीय अभिसरण समस्याओं को संबोधित करने के लिए सुधारित तकनीकें प्रस्तुत करता है। अनुसंधान से पता चलता है कि सभी MFM सिमुलेशन के लिए समान प्रत्यक्ष संख्यात्मक सिमुलेशन (DNS) का उपयोग एड़ी विसरणशीलता क्षणों की सांख्यिकीय अभिसरण को काफी हद तक सुधार सकता है। लेखकों ने विघटन विधि को जोड़ते हुए एक कार्यान्वयन योजना विकसित की है, जो MFM फोर्सिंग को अर्ध-विश्लेषणात्मक रूप से संभाल सकती है और सीमा शर्तों के सुसंगत उपचार को सुनिश्चित करती है। द्विविमीय Rayleigh-Taylor अस्थिरता केस अध्ययन में, यह विधि अभिसरण के लिए आवश्यक सिमुलेशन की संख्या को O(1000) से O(100) तक कम कर देती है।

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

समस्या की पृष्ठभूमि

अशांत प्रवाह मॉडलिंग में, Reynolds-औसत Navier-Stokes (RANS) समीकरणों को अनुबंधित पदों जैसे अदिश प्रवाह या Reynolds प्रतिबल को संभालने के लिए समापन मॉडल की आवश्यकता होती है। पारंपरिक Boussinesq सन्निकटन स्थानीय और समदैशिक एड़ी विसरणशीलता गुणांक मानता है, लेकिन जटिल प्रवाह के लिए अक्सर विफल होता है।

मुख्य समस्याएं

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

अनुसंधान प्रेरणा

मौजूदा मैक्रोस्कोपिक फोर्सिंग विधि (MFM) और Green फ़ंक्शन विधि अ-स्थानीय और अनिसोट्रॉपिक एड़ी विसरणशीलता की गणना कर सकती हैं, लेकिन उच्च-क्रम क्षणों की गणना करते समय धीमी सांख्यिकीय अभिसरण की समस्या होती है। यह पेपर सांख्यिकीय अभिसरण को तेज करने और कम्प्यूटेशनल लागत को कम करने के लिए सुधारित तकनीकें विकसित करने का लक्ष्य रखता है।

मुख्य योगदान

  1. सांख्यिकीय त्रुटि प्रवर्धन समस्या की पहचान और समाधान: स्वतंत्र दाता सिमुलेशन के उपयोग से उच्च-क्रम क्षणों में सांख्यिकीय त्रुटि प्रवर्धन होता है, एकल दाता का उपयोग करने का समाधान प्रस्तावित किया गया है
  2. विघटन MFM विधि का विकास: सामान्य अ-स्थिर और अराजक प्रवाह, साथ ही गति परिवहन के लिए विघटन विधि का विस्तार किया गया है
  3. कम्प्यूटेशनल दक्षता में उल्लेखनीय सुधार: 2D Rayleigh-Taylor अस्थिरता में लगभग 10 गुना कम्प्यूटेशनल त्वरण प्राप्त किया गया है
  4. संपूर्ण सैद्धांतिक और कार्यान्वयन ढांचा प्रदान करना: अदिश और गति परिवहन के लिए संपूर्ण गणितीय ढांचा और कार्यान्वयन विवरण प्रदान किए गए हैं

विधि विवरण

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

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

एड़ी विसरणशीलता गुणांक क्षणों का गणितीय विवरण

एड़ी विसरणशीलता गुणांक को Kramers-Moyal विस्तार के माध्यम से क्षणों के रूप में दर्शाया जा सकता है:

−⟨ui′c′⟩(x,t)=[Dij00(x,t)+Dijk10(x,t)∂∂xk+⋯+Dij01(x,t)∂∂t+⋯ ]∂⟨c⟩∂xj-\langle u'_i c' \rangle(x,t) = \left[D^{00}_{ij}(x,t) + D^{10}_{ijk}(x,t)\frac{\partial}{\partial x_k} + \cdots + D^{01}_{ij}(x,t)\frac{\partial}{\partial t} + \cdots\right]\frac{\partial\langle c \rangle}{\partial x_j}

जहां विभिन्न क्रम के क्षणों को इस प्रकार परिभाषित किया गया है:

  • Dij00(x,t)=∫∫Dij(x,x′,t,t′)dx′dt′D^{00}_{ij}(x,t) = \int\int D_{ij}(x,x',t,t')dx'dt' (शून्य-क्रम स्पेसटाइम क्षण)
  • Dijk10(x,t)=∫∫(xk′−xk)Dij(x,x′,t,t′)dx′dt′D^{10}_{ijk}(x,t) = \int\int (x'_k-x_k)D_{ij}(x,x',t,t')dx'dt' (प्रथम-क्रम स्थान क्षण)
  • Dij01(x,t)=∫∫(t′−t)Dij(x,x′,t,t′)dx′dt′D^{01}_{ij}(x,t) = \int\int (t'-t)D_{ij}(x,x',t,t')dx'dt' (प्रथम-क्रम समय क्षण)

मैक्रोस्कोपिक फोर्सिंग विधि (MFM)

MFM अदिश परिवहन समीकरण में मैक्रोस्कोपिक फोर्सिंग जोड़कर एड़ी विसरणशीलता गुणांक के क्षणों की गणना करता है:

∂c∂t+∂∂xi(uic)=DM∂2c∂xi∂xi+s\frac{\partial c}{\partial t} + \frac{\partial}{\partial x_i}(u_i c) = D_M \frac{\partial^2 c}{\partial x_i \partial x_i} + s

विभिन्न औसत अदिश प्रवणताओं (जैसे बहुपद रूप) को निर्दिष्ट करके, संबंधित एड़ी विसरणशीलता गुणांक के क्षणों को सीधे प्राप्त किया जा सकता है।

विघटन MFM विधि

विघटन विधि अदिश उतार-चढ़ाव c′c' को इस प्रकार विस्तारित करती है:

c′(x,t)=[cj00(x,t)+cjk10(x,t)∂∂xk+⋯+cj01(x,t)∂∂t+⋯ ]∂⟨c⟩∂xjc'(x,t) = \left[c^{00}_j(x,t) + c^{10}_{jk}(x,t)\frac{\partial}{\partial x_k} + \cdots + c^{01}_j(x,t)\frac{\partial}{\partial t} + \cdots\right]\frac{\partial\langle c \rangle}{\partial x_j}

यह MFM फोर्सिंग को अर्ध-विश्लेषणात्मक रूप से संभालने की अनुमति देता है, आवधिक सीमा शर्तों की संगतता समस्या को हल करता है।

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

  1. एकल दाता रणनीति: सभी receiver समीकरणों को प्रवाह क्षेत्र प्रदान करने के लिए समान DNS सिमुलेशन का उपयोग करके, स्वतंत्र दाता सिमुलेशन के बीच अंतर से उत्पन्न त्रुटि प्रवर्धन से बचा जाता है
  2. त्रुटि प्रसार विश्लेषण: सैद्धांतिक विश्लेषण से पता चलता है कि स्वतंत्र दाता का उपयोग करते समय, उच्च-क्रम क्षणों की सांख्यिकीय त्रुटि स्थानिक निर्देशांक की शक्तियों के अनुसार प्रवर्धित होती है: D1=F1−x1F0+x1(ε01−ε00)+ε11D^1 = F^1 - x_1 F^0 + x_1(\varepsilon^{01} - \varepsilon^{00}) + \varepsilon^{11}
  3. विघटन विधि का विस्तार: मूल रूप से स्थिर लामिनार प्रवाह के लिए उपयोग की जाने वाली विघटन विधि को सामान्य अ-स्थिर अराजक प्रवाह और गति परिवहन तक विस्तारित किया गया है

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

केस अध्ययन: द्विविमीय Rayleigh-Taylor अस्थिरता

  • कम्प्यूटेशनल डोमेन: 2049×2049 ग्रिड वाला वर्गाकार डोमेन
  • सीमा शर्तें: x₁ दिशा में आवधिक, x₂ दिशा में नो-स्लिप और नो-पेनिट्रेशन
  • भौतिक पैरामीटर:
    • Atwood संख्या: A = 0.05
    • अधिकतम Mach संख्या: Ma_max = 0.05
    • Grashof संख्या: Gr = 1
    • Schmidt संख्या: Sc = 1
    • Reynolds संख्या: Re_T = 54, Re_L = 8,000

मूल्यांकन मेट्रिक्स

  • सांख्यिकीय अभिसरण: विभिन्न कार्यान्वयन संख्याओं के तहत मूल माध्य वर्ग त्रुटि (RMSE) द्वारा मूल्यांकन
  • कम्प्यूटेशनल दक्षता: समान अभिसरण सटीकता प्राप्त करने के लिए आवश्यक सिमुलेशन संख्या की तुलना
  • मॉडलिंग प्रदर्शन: मोमेंट मैचिंग इनवर्जन (MMI) विधि द्वारा निर्मित समापन ऑपरेटर प्रदर्शन का मूल्यांकन

तुलनात्मक विधियां

  • मानक MFM (स्वतंत्र दाता): प्रत्येक क्षण के लिए स्वतंत्र दाता-receiver जोड़ी का उपयोग
  • विघटन MFM (एकल दाता): सभी receiver समीकरणों के लिए एक ही दाता सिमुलेशन का उपयोग

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

मुख्य परिणाम

सांख्यिकीय अभिसरण तुलना

एकल दाता का उपयोग करने वाली विघटन MFM स्वतंत्र दाता की मानक MFM की तुलना में सांख्यिकीय अभिसरण में उल्लेखनीय सुधार दिखाती है:

  1. शून्य-क्रम क्षण D⁰⁰: दोनों विधियां समान प्रदर्शन करती हैं, क्योंकि गणना विधि गणितीय रूप से समतुल्य है
  2. प्रथम-क्रम समय क्षण D⁰¹: विघटन MFM 100 कार्यान्वयनों में अच्छी अभिसरण प्राप्त करती है, जबकि मानक MFM को 1000 कार्यान्वयनों की आवश्यकता होती है
  3. प्रथम-क्रम स्थान क्षण D¹⁰: विघटन MFM की 100 कार्यान्वयनों में त्रुटि मानक MFM से स्पष्ट रूप से कम है

कम्प्यूटेशनल दक्षता में सुधार

  • त्वरण अनुपात: लगभग 10 गुना कम्प्यूटेशनल त्वरण (O(1000) से O(100) सिमुलेशन तक)
  • त्रुटि स्केलिंग: सैद्धांतिक विश्लेषण और संख्यात्मक परिणाम दोनों समय और स्थानिक निर्देशांक के अनुसार त्रुटि प्रवर्धन की पुष्टि करते हैं

विलोपन प्रयोग: अशांत चैनल प्रवाह सत्यापन

Re_τ = 180 के अशांत चैनल प्रवाह में विधि की प्रभावशीलता का सत्यापन किया गया:

  • मानक MFM (एकल दाता) और विघटन MFM परिणाम लगभग पूरी तरह से सुसंगत हैं
  • यह साबित करता है कि सुधार विघटन विधि के बजाय दाता के एकीकरण से आता है

मॉडलिंग प्रदर्शन मूल्यांकन

सुधारी गई अभिसरण एड़ी विसरणशीलता गुणांक के क्षणों का उपयोग करके निर्मित MMI समापन ऑपरेटर:

  • विघटन MFM: 200 कार्यान्वयनों के तहत निर्मित मॉडल DNS परिणामों के साथ अच्छी तरह से मेल खाता है
  • मानक MFM: समान कार्यान्वयन संख्या के तहत निर्मित मॉडल महत्वपूर्ण विचलन दिखाता है, जो गलत मॉडलिंग निष्कर्षों की ओर ले जा सकता है

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

एड़ी विसरणशीलता गुणांक गणना विधियां

  1. Green फ़ंक्शन विधि: Hamba आदि द्वारा विकसित सटीक लेकिन कम्प्यूटेशनल रूप से महंगी विधि
  2. तेज़ MFM: Bryngelson आदि द्वारा विकसित छिपी हुई विरलता का उपयोग करने वाली त्वरण विधि
  3. सहायक MFM: Liu आदि द्वारा प्रस्तावित सहायक विधि पर आधारित लक्ष्य गणना तकनीक

अ-स्थानीय मॉडलिंग सिद्धांत

  • Kraichnan का सटीक अ-स्थानीय अभिव्यक्ति
  • विसरण मॉडलिंग में क्षण-प्रवणता विधि का अनुप्रयोग
  • अशांत मॉडलिंग में Kramers-Moyal विस्तार का अनुप्रयोग

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

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

  1. सांख्यिकीय त्रुटि स्रोत की पहचान: स्वतंत्र दाता सिमुलेशन के बीच सूक्ष्म अंतर उच्च-क्रम क्षणों की गणना में O(1) त्रुटि में प्रवर्धित होते हैं
  2. प्रभावी समाधान: एकल दाता का उपयोग त्रुटि प्रवर्धन को समाप्त कर सकता है, सांख्यिकीय अभिसरण को काफी हद तक सुधार सकता है
  3. विधि की सार्वभौमिकता: विघटन MFM अदिश और गति परिवहन दोनों पर लागू होता है, आवधिक सीमा शर्तों की समस्या को संभाल सकता है

सीमाएं

  1. प्रवाह प्रकार प्रतिबंध: मुख्य रूप से 2D Rayleigh-Taylor अस्थिरता में सत्यापित, अधिक प्रवाह प्रकारों के परीक्षण की आवश्यकता है
  2. संख्यात्मक कार्यान्वयन निर्भरता: सुधार की डिग्री विशिष्ट संख्यात्मक सॉल्वर और समानांतर कार्यान्वयन पर निर्भर हो सकती है
  3. सैद्धांतिक विश्लेषण सरलीकरण: त्रुटि विश्लेषण सरलीकृत धारणाओं पर आधारित है, वास्तविक स्थिति अधिक जटिल हो सकती है

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

  1. अनुप्रयोग का विस्तार: अधिक जटिल प्रवाह में विधि की प्रभावशीलता का सत्यापन
  2. सैद्धांतिक सुधार: त्रुटि प्रसार का अधिक संपूर्ण सिद्धांत विकसित करना
  3. एल्गोरिदम अनुकूलन: कम्प्यूटेशनल दक्षता को और बढ़ाने के लिए अन्य त्वरण तकनीकों को जोड़ना

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

शक्तियां

  1. समस्या की पहचान सटीक: MFM में सांख्यिकीय त्रुटि प्रवर्धन समस्या की सटीक पहचान
  2. समाधान प्रभावी: प्रस्तावित एकल दाता रणनीति सरल और प्रभावी है
  3. सैद्धांतिक विश्लेषण गहन: त्रुटि प्रसार का गणितीय विश्लेषण प्रदान किया गया है
  4. प्रायोगिक सत्यापन पर्याप्त: विधि की प्रभावशीलता को कई मामलों में सत्यापित किया गया है
  5. व्यावहारिक मूल्य उच्च: कम्प्यूटेशनल लागत को काफी हद तक कम करता है, विधि की व्यावहारिकता में सुधार करता है

कमियां

  1. सीमित केस: मुख्य रूप से 2D RT अस्थिरता में सत्यापित, 3D और अन्य प्रकार के प्रवाह की कमी है
  2. सैद्धांतिक सीमाएं: त्रुटि विश्लेषण एक-आयामी सरलीकृत मामले पर आधारित है, वास्तविक बहु-आयामी स्थिति अधिक जटिल है
  3. निर्भरता समस्या: सुधार प्रभाव विशिष्ट संख्यात्मक कार्यान्वयन और समानांतर रणनीति पर निर्भर हो सकता है

प्रभाव

  1. शैक्षणिक योगदान: MFM विधि के सुधार के लिए महत्वपूर्ण तकनीकी पथ प्रदान करता है
  2. व्यावहारिक मूल्य: एड़ी विसरणशीलता गुणांक क्षणों की गणना की दक्षता में उल्लेखनीय सुधार
  3. विधि प्रचार: विघटन MFM विधि में अच्छी प्रचार संभावनाएं हैं

लागू परिदृश्य

  1. अशांत मॉडलिंग: सटीक एड़ी विसरणशीलता की आवश्यकता वाली RANS मॉडलिंग
  2. मिश्रण समस्याएं: Rayleigh-Taylor, Richtmyer-Meshkov आदि अस्थिरता अनुसंधान
  3. अ-स्थानीय प्रभाव अनुसंधान: अ-स्थानीय और अनिसोट्रॉपिक प्रभावों को परिमाणित करने की आवश्यकता वाली प्रवाह समस्याएं

संदर्भ

यह पेपर 37 महत्वपूर्ण संदर्भों का हवाला देता है, जो अशांत सिद्धांत, संख्यात्मक विधियों और MFM संबंधित अनुसंधान को कवर करते हैं, जो अनुसंधान के लिए एक ठोस सैद्धांतिक आधार प्रदान करते हैं।

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