2025-11-14T05:04:11.127420

A combined thermal-resistance-capacity and finite-element model for very fast and accurate short- and medium-term simulations of single U-tube borehole heat exchangers

Zanchini, Zanchini, Naldi
An accurate design of a ground-coupled heat pump system requires the knowledge of the outlet fluid temperature from the borehole heat exchangers (BHEs), both in the short and long term. This paper fucuses on the short and medium term. In this time range, either 3D finite-element simulations or Thermal Resistance Capacity Models (TRCMs) can be applied. The former can yield very accurate results but require long computation times. The latter are much faster but cannot be fully precise, because they require simplifying assumptions. In this paper, we present a new method for the short-term and medium-term simulation of single U-tube BHEs, which combines the speed of TRCMs and the accuracy of finite-elements simulations. The method uses a TRCM to estimate the thermal response of the BHE, then corrects the results by interpolation with a dataset, which has been produced by running 54 finite-element simulations in various configurations. The model is implemented in a C++ program, available at the open-source online data repository of the University of Bologna. The program provides, within two seconds, the time evolution of the inlet, outlet and mean fluid temperature, of the mean BHE surface temperature, of the 3D and the effective borehole thermal resistance. It can be easily connected to long-term simulation tools to obtain the full-time-scale thermal response of a bore field.
academic

एकल यू-ट्यूब बोरहोल हीट एक्सचेंजर के अत्यंत तीव्र और सटीक अल्पकालीन और मध्यकालीन सिमुलेशन के लिए संयुक्त तापीय-प्रतिरोध-क्षमता और परिमित-तत्व मॉडल

मूल जानकारी

  • पेपर आईडी: 2510.14421
  • शीर्षक: एकल यू-ट्यूब बोरहोल हीट एक्सचेंजर के अत्यंत तीव्र और सटीक अल्पकालीन और मध्यकालीन सिमुलेशन के लिए संयुक्त तापीय-प्रतिरोध-क्षमता और परिमित-तत्व मॉडल
  • लेखक: एंजो ज़ांचिनी, फ्रांसेस्को ज़ांचिनी, क्लाउडिया नाल्डी (बोलोग्ना विश्वविद्यालय)
  • वर्गीकरण: physics.app-ph
  • प्रकाशन समय: 2024
  • पेपर लिंक: https://arxiv.org/abs/2510.14421

सारांश

भूतापीय ऊष्मा पंप प्रणालियों के सटीक डिजाइन के लिए बोरहोल हीट एक्सचेंजर (बीएचई) के अल्पकालीन और दीर्घकालीन आउटलेट द्रव तापमान को समझना आवश्यक है। यह पेपर अल्पकालीन और मध्यकालीन विश्लेषण पर केंद्रित है। इस समय सीमा में, 3डी परिमित-तत्व सिमुलेशन या तापीय प्रतिरोध-क्षमता मॉडल (टीआरसीएम) का उपयोग किया जा सकता है। पहला उच्च सटीकता लेकिन लंबी गणना समय प्रदान करता है, जबकि दूसरा तीव्र गणना लेकिन सरलीकृत मान्यताओं के कारण सीमित सटीकता प्रदान करता है। यह पेपर एकल यू-ट्यूब बोरहोल हीट एक्सचेंजर के अल्पकालीन और मध्यकालीन सिमुलेशन के लिए टीआरसीएम की गति और परिमित-तत्व सिमुलेशन की सटीकता को जोड़ने वाली एक नई विधि प्रस्तावित करता है। यह विधि टीआरसीएम का उपयोग करके बीएचई की तापीय प्रतिक्रिया का अनुमान लगाती है, फिर 54 विभिन्न कॉन्फ़िगरेशन परिमित-तत्व सिमुलेशन द्वारा उत्पन्न डेटासेट के माध्यम से इंटरपोलेशन सुधार करती है। मॉडल को सी++ प्रोग्राम के रूप में लागू किया गया है और बोलोग्ना विश्वविद्यालय के ओपन-सोर्स डेटाबेस में उपलब्ध है, जो दो सेकंड में इनलेट और आउटलेट तापमान, औसत द्रव तापमान, औसत बीएचई सतह तापमान, 3डी और प्रभावी बोरहोल तापीय प्रतिरोध का समय विकास प्रदान कर सकता है।

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

समस्या विवरण

भूतापीय ऊष्मा पंप प्रणालियां उच्च दक्षता वाली भवन हीटिंग और शीतलन तकनीक हैं, जिनके डिजाइन के लिए बोरहोल हीट एक्सचेंजर की तापीय प्रतिक्रिया विशेषताओं का सटीक पूर्वानुमान आवश्यक है। मुख्य चुनौतियाँ हैं:

  1. गणना सटीकता और दक्षता का विरोधाभास: 3डी परिमित-तत्व सिमुलेशन उच्च सटीकता लेकिन घंटों का समय लेता है; टीआरसीएम तीव्र लेकिन सीमित सटीकता प्रदान करता है
  2. अल्पकालीन और मध्यकालीन सिमुलेशन की आवश्यकता: 10⁻²·⁵ घंटे से 10³ घंटे की समय सीमा में मौजूदा विधियों में सटीकता या दक्षता की समस्याएं हैं
  3. व्यावहारिक इंजीनियरिंग अनुप्रयोग की आवश्यकता: प्रणाली डिजाइन और अनुकूलन के लिए तीव्र और सटीक सिमुलेशन उपकरण की आवश्यकता है

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

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

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

टीआरसीएम की गति और परिमित-तत्व सटीकता दोनों को संयोजित करने वाली एक संकर सिमुलेशन विधि विकसित करना, जो इंजीनियरिंग व्यावहारिक आवश्यकताओं को पूरा करे।

मुख्य योगदान

  1. एचयूएसटी मान्यता के प्रभाव की खोज और परिमाणीकरण: उच्च तापीय चालकता परत के साथ और बिना परिमित-तत्व सिमुलेशन की तुलना करके, यह साबित किया कि एचयूएसटी मान्यता T_fm, T_in, T_out, R_b3D और R_beff को कम आंकती है
  2. सुधारित टीआरसीएम मॉडल का विकास: द्रव, पाइप, ग्राउट और भूस्तर नोड्स के साथ स्तरीकृत तापीय नेटवर्क का उपयोग करके, सिमुलेशन सटीकता में सुधार
  3. सुधार गुणांक डेटाबेस की स्थापना: 54 विभिन्न कॉन्फ़िगरेशन परिमित-तत्व सिमुलेशन के माध्यम से, समय-निर्भर सुधार गुणांक डेटासेट का निर्माण
  4. कुशल संकर सिमुलेशन विधि का कार्यान्वयन: टीआरसीएम आधार गणना और इंटरपोलेशन सुधार को संयोजित करके, 2 सेकंड में परिमित-तत्व सिमुलेशन सटीकता प्राप्त करना
  5. ओपन-सोर्स सिमुलेशन उपकरण प्रदान करना: सी++ प्रोग्राम बीएचई के तीव्र और सटीक अल्पकालीन और मध्यकालीन सिमुलेशन के लिए उपयोग किया जा सकता है

विधि विवरण

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

इनपुट: बीएचई ज्यामितीय पैरामीटर (लंबाई L, त्रिज्या r_b, पाइप दूरी s आदि), तापीय गुण पैरामीटर, प्रवाह दर, तापीय भार आउटपुट: समय विकास के साथ इनलेट और आउटलेट तापमान (T_in, T_out), औसत द्रव तापमान (T_fm), बीएचई सतह तापमान (T_b), 3डी तापीय प्रतिरोध (R_b3D), प्रभावी तापीय प्रतिरोध (R_beff) समय सीमा: 10⁻²·⁵ घंटे से 10³ घंटे तक

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

1. परिमित-तत्व आधार मॉडल

  • कॉमसोल मल्टीफिजिक्स के ऊष्मा स्थानांतरण ठोस इंटरफेस का उपयोग करके ठोस डोमेन ऊष्मा चालन का अनुकरण
  • द्रव ऊर्जा संतुलन का अनुकरण करने के लिए पाइप प्रवाह मॉड्यूल का उपयोग
  • जाल स्वतंत्रता सत्यापन सटीकता सुनिश्चित करता है
  • ऊर्ध्वाधर समन्वय स्केलिंग (कारक 20) जाल गुणवत्ता में सुधार करता है

2. सुधारित टीआरसीएम मॉडल

स्तरीकृत संरचना: बीएचई और भूस्तर को m क्षैतिज परतों में विभाजित किया जाता है, प्रत्येक परत की ऊंचाई l = L/m है

तापीय नेटवर्क: प्रत्येक परत में निम्नलिखित नोड्स शामिल हैं

  • द्रव नोड्स: f1 (नीचे की ओर), f2 (ऊपर की ओर)
  • पाइप नोड्स: p1, p2
  • ग्राउट नोड्स: gt1, gt2
  • बोरहोल नोड: b
  • भूस्तर नोड्स: g1, g2, ..., gn

तापीय प्रतिरोध गणना:

R₁^Δ = 2R_b                    (द्रव से बोरहोल दीवार तक)
R₁₂^Δ = 4R_b R_a/(4R_b - R_a)   (द्रव के बीच तापीय प्रतिरोध)
R_conv = 1/(2πr_pi h)           (संवहन तापीय प्रतिरोध)

ऊर्जा संतुलन समीकरण: प्रत्येक नोड के लिए क्षणिक ऊर्जा संतुलन समीकरण स्थापित किए जाते हैं, जो रैखिक समीकरण समूह बनाते हैं।

3. सुधार गुणांक विधि

सुधार गुणांक परिभाषा:

C_coeff = R_b3D,COMSOL / R_b3D,TRCM

सुधार प्रक्रिया:

  1. टीआरसीएम आधार तापीय प्रतिक्रिया की गणना करता है
  2. पैरामीटर इंटरपोलेशन के आधार पर सुधार गुणांक प्राप्त करता है
  3. R_b3D को सुधारता है और तापमान क्षेत्र की पुनः गणना करता है

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

  1. एचयूएसटी प्रभाव का मात्रात्मक विश्लेषण: उच्च तापीय चालकता परत को शामिल करके समान सीमा तापमान को लागू करके, एचयूएसटी मान्यता के सिमुलेशन सटीकता पर प्रभाव का मात्रात्मक विश्लेषण
  2. स्तरीकृत तापीय नेटवर्क डिजाइन: पारंपरिक टीआरसीएम की तुलना में, पाइप और ग्राउट नोड्स जोड़े गए हैं, बीएचई के अंदर ऊष्मा स्थानांतरण का अधिक सटीक विवरण
  3. पश्च-सत्यापन सुधार रणनीति: टीआरसीएम की T_b को सटीक रूप से भविष्यवाणी करने की विशेषता का उपयोग करके, R_b3D को सुधारकर अन्य पैरामीटर को अप्रत्यक्ष रूप से सुधारना
  4. बहु-पैरामीटर इंटरपोलेशन तकनीक: बोरहोल त्रिज्या, पाइप दूरी, ग्राउट तापीय चालकता, पाइप बाहरी व्यास चार मुख्य पैरामीटर के आधार पर परवलयिक इंटरपोलेशन

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

सिमुलेशन कॉन्फ़िगरेशन

ज्यामितीय पैरामीटर सीमा:

  • बोरहोल त्रिज्या r_b: 64-76 मिमी
  • पाइप दूरी s: 27-47 मिमी
  • पाइप बाहरी व्यास r_pe: 16, 20 मिमी
  • ग्राउट तापीय चालकता k_gt: 1.0-2.2 डब्ल्यू/(मी·के)

निश्चित पैरामीटर:

  • बीएचई लंबाई: 100 मीटर
  • दफन गहराई: 1.8 मीटर
  • भूस्तर तापीय चालकता: 1.8 डब्ल्यू/(मी·के)
  • तापीय भार: 50 डब्ल्यू/मीटर

डेटासेट निर्माण

  • 54 विभिन्न कॉन्फ़िगरेशन बीएचई पैरामीटर संयोजन
  • प्रत्येक समूह के लिए कॉमसोल और टीआरसीएम तुलनात्मक सिमुलेशन
  • समय बिंदु: 56 लॉगरिदमिक समान दूरी वाले बिंदु (10⁻²·⁵ से 10³ घंटे तक)
  • इंटरपोलेशन के लिए सुधार गुणांक डेटाबेस उत्पन्न करना

सत्यापन मामले

तीन विशिष्ट बीएचई कॉन्फ़िगरेशन का चयन सत्यापन के लिए:

  • बीएचई 1: मानक कॉन्फ़िगरेशन, सुधार गुणांक गणना के लिए उपयोग किया जाता है
  • बीएचई 2: पैरामीटर जिन्हें इंटरपोलेशन की आवश्यकता है
  • बीएचई 3: छोटी त्रिज्या कॉन्फ़िगरेशन जिसे एक्सट्रापोलेशन की आवश्यकता है

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

मुख्य परिणाम

एचयूएसटी मान्यता प्रभाव विश्लेषण

उच्च तापीय चालकता परत के साथ और बिना सिमुलेशन परिणामों की तुलना करके:

  • T_fm, T_in, T_out को कम आंका जाता है, T_b लगभग अप्रभावित है
  • R_b3D और R_beff में महत्वपूर्ण कम आंकलन
  • कम आंकलन की डिग्री समय के साथ बढ़ती है, मध्यकालीन (>4 घंटे) में विशेष रूप से स्पष्ट

सुधार प्रभाव सत्यापन

बीएचई 1 परिणाम:

  • T_fm रूट मीन स्क्वायर त्रुटि: 0.037°C
  • T_in, T_out रूट मीन स्क्वायर त्रुटि: <0.038°C
  • R_beff अंतिम समय त्रुटि: 0.066%

बीएचई 2 परिणाम (इंटरपोलेशन की आवश्यकता):

  • T_fm रूट मीन स्क्वायर त्रुटि: 0.051°C
  • R_b3D अंतिम समय त्रुटि: 0.13%
  • R_beff अंतिम समय त्रुटि: 0.23%

बीएचई 3 परिणाम (एक्सट्रापोलेशन की आवश्यकता):

  • T_fm रूट मीन स्क्वायर त्रुटि: 0.052°C
  • R_b3D अंतिम समय त्रुटि: 0.20%
  • R_beff अंतिम समय त्रुटि: 0.40%

गणना दक्षता

  • सुधारित टीआरसीएम: 2 सेकंड में पूर्ण
  • कॉमसोल सिमुलेशन: घंटों
  • सटीकता सुधार: परिमित-तत्व सिमुलेशन सटीकता प्राप्त करता है

तुलनात्मक विश्लेषण

पास्क्वियर-मार्कोट कृत्रिम तंत्रिका नेटवर्क के साथ तुलना:

  • दोनों विधियों में मध्यकालीन कम आंकलन की समस्या है
  • यह पेपर सुधार गुणांक के माध्यम से इस समस्या को प्रभावी ढंग से हल करता है

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

विश्लेषणात्मक समाधान विधि

  • लामार्चे और ब्यूचैम्प: अल्पकालीन विश्लेषणात्मक समाधान, लेकिन द्रव के बीच ऊष्मा विनिमय पर विचार नहीं कर सकते
  • बीयर: द्रव ऊर्जा संतुलन पर विचार करने वाला विश्लेषणात्मक मॉडल, लाप्लास व्युत्क्रम परिवर्तन की आवश्यकता है

संख्यात्मक सिमुलेशन विधि

  • ली और झेंग: 3डी ठोस डोमेन + 1डी द्रव ऊर्जा संतुलन
  • मार्कोट और पास्क्वियर: 3डी ठोस और द्रव डोमेन पूर्ण युग्मन सिमुलेशन

टीआरसीएम विधि

  • बाउर आदि: मानक त्रिकोणीय तापीय प्रतिरोध नेटवर्क, एचयूएसटी मान्यता का उपयोग करते हैं
  • पास्क्वियर और मार्कोट: सुधारित टीआरसीएम और अर्ध-3डी मॉडल का विकास
  • रुइज़-कैल्वो आदि: बोरहोल-भूस्तर मॉडल, समान तापमान मान्यता से बचते हैं लेकिन सीमित सटीकता

मशीन लर्निंग विधि

  • पास्क्वियर आदि: 15000 टीआरसीएम सिमुलेशन पर प्रशिक्षित कृत्रिम तंत्रिका नेटवर्क

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

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

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

सीमाएं

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

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

  1. अन्य प्रकार के बीएचई (दोहरे यू-ट्यूब, समाक्षीय) तक विस्तार
  2. भूजल पारगम्यता जैसी जटिल ऊष्मा स्थानांतरण तंत्र पर विचार करना
  3. दीर्घकालीन जी-फ़ंक्शन विधि के साथ युग्मन, पूर्ण समय-स्केल सिमुलेशन प्राप्त करना
  4. स्व-अनुकूली सुधार रणनीति विकसित करना, पूर्व-निर्मित डेटाबेस पर निर्भरता को कम करना

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

शक्तियां

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

कमियां

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

प्रभाव

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

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

  1. भूतापीय ऊष्मा पंप प्रणाली डिजाइन: बीएचई के अल्पकालीन और मध्यकालीन तापीय प्रदर्शन का तीव्र और सटीक मूल्यांकन
  2. तापीय प्रतिक्रिया परीक्षण विश्लेषण: क्षेत्र परीक्षण डेटा व्याख्या और पैरामीटर पहचान में सहायता
  3. प्रणाली अनुकूलन अनुसंधान: बीएचई पैरामीटर अनुकूलन के लिए उच्च दक्षता सिमुलेशन उपकरण प्रदान करना
  4. शिक्षण और अनुसंधान: बीएचई ऊष्मा स्थानांतरण विश्लेषण के लिए मानक उपकरण के रूप में

संदर्भ

पेपर 45 संबंधित संदर्भों का हवाला देता है, जो बीएचई सिमुलेशन की मुख्य विधियों और नवीनतम प्रगति को कवर करते हैं, इस अनुसंधान के लिए एक मजबूत सैद्धांतिक आधार प्रदान करते हैं। मुख्य संदर्भों में शामिल हैं:

  • क्लेसन और जेवेड (2011): जी-फ़ंक्शन विश्लेषण विधि
  • बाउर आदि (2011): टीआरसीएम मानक विधि
  • पास्क्वियर और मार्कोट (2014, 2020): सुधारित टीआरसीएम और मशीन लर्निंग विधि
  • ज़ांचिनी (2024): बीएचई सिमुलेशन विधि तुलनात्मक विश्लेषण

समग्र मूल्यांकन: यह एक उच्च गुणवत्ता वाला अनुप्रयुक्त भौतिकी पेपर है, जो भूतापीय ऊष्मा पंप बीएचई सिमुलेशन की व्यावहारिक इंजीनियरिंग आवश्यकताओं के लिए, एक नई संकर मॉडलिंग विधि प्रस्तावित करता है। पेपर में गहन सैद्धांतिक विश्लेषण, पर्याप्त प्रायोगिक सत्यापन, और महत्वपूर्ण शैक्षणिक मूल्य और इंजीनियरिंग अनुप्रयोग संभावनाएं हैं। ओपन-सोर्स प्रोग्राम का प्रावधान इसके प्रभाव और व्यावहारिक मूल्य को और बढ़ाता है।