Bu araşdırma turbulent kanal axınında divara verilən istilik selinin eksponensial şəkildə sürətlə artdığı halda suyun qaynama başlamazdan dərhal əvvəl necə davrandığını araşdırır. Tədqiqatın əsas diqqət nöqtəsi nüvəli qaynama başlanğıcı, yəni Onset of Nucleate Boiling (ONB) öncəsindəki birfazalı maye rejimidir. Bu mərhələdə axıcı hələ də mayedir; lakin divar yaxınlığındakı temperatur sahəsi, turbulentlik quruluşu və istilikötürmə mexanizmləri sürətlə dəyişir.

Elie Roumet və həmkarları bu prosesi yüksək ayırdetməli Large-Eddy Simulation (LES) üsulu ilə araşdırıblar. Simulyasiyalarda suyun temperaturdan asılı özlülüyü və istilikkeçirmə qabiliyyəti nəzərə alınıb. Bu çox vacibdir; çünki yüksək istilik seli altında divar yaxınlığındakı suyun temperaturu sürətlə yüksəlir, özlülük nəzərəçarpacaq dərəcədə azalır və istilikkeçirmə temperaturdan asılı olaraq qeyri-xətti dəyişir. Beləliklə, temperatur sahəsi yalnız passiv şəkildə daşınmır; eyni zamanda sürət sahəsinə, turbulentliyə və divar yaxınlığındakı axın quruluşuna əks təsir göstərir.

Araşdırma xüsusilə sürətli və yavaş isitmə rampalarını müqayisə edir. Xarakterik isitmə müddəti \(\tau = 5 \, \mathrm{ms}\) olan sürətli keçidlə \(\tau = 100 \, \mathrm{ms}\) olan quasi-steady, yəni demək olar ki, sabit keçid arasında mühüm fərqlər müəyyən edilib. Sürətli isitmə zamanı Stanton ədədi təxminən %30 artır, termal giriş uzunluğu qısalır və divar yaxınlığında demək olar ki, spanwise-homogen superqızdırılmış maye təbəqəsi yaranır. Bunun əksinə, yavaş isitmə zamanı superqızdırılmış maye turbulent kanal axınının xarakterik aşağı sürətli streak strukturları ilə əlaqəli olaraq axın istiqamətində zolaqlar şəklində qalır.

Araşdırmanın əsas mesajı budur: qaynama başlamazdan əvvəlki birfazalı proses yalnız temperaturun yüksəlməsindən ibarət deyil. Divar yaxınlığındakı turbulent strukturlar yenidən təşkil olunur, axın istiqamətində sürət artır, spanwise vortisite güclənir, buffer zonasında Couette-bənzər kəsilmə təbəqəsi yaranır və temperatur keçidi divara daha çox sıxışır. Buna görə sürətli istilik seli keçidlərində qaynama başlanğıcını anlamaq üçün yalnız orta divar temperaturuna baxmaq kifayət deyil; divar yaxınlığındakı üçölçülü turbulentlik, temperatur sahəsi, özlülük dəyişməsi və istilik seli mexanizmləri birlikdə qiymətləndirilməlidir.

Araşdırmanın əsas problemi nədir?

Bu araşdırma çox yüksək istilik selinə məruz qalan turbulent su axınlarında qaynama başlamazdan dərhal əvvəl nələrin baş verdiyini anlamağa çalışır. Mövzu xüsusilə nüvə tədqiqat reaktorları, qabaqcıl elektron soyutma sistemləri, raket mühərriki soyutma kanalları və gələcək füzyon reaktorlarının plazmaya baxan komponentləri üçün vacibdir.

Belə sistemlərdə divara düşən istilik seli qısa müddətdə çox sürətlə arta bilər. Məsələn, nüvə tədqiqat reaktorunda reaktivlik mənşəli qəza ssenarisi yanacaq ətrafındakı istilik yükünü eksponensial şəkildə artıra bilər. Elektron sistemlərdə qəfil güc artımı kiçik səthdə qaynar nöqtənin yaranmasına səbəb ola bilər. Raket mühərriklərində alışma və ya dartı dəyişməsi zamanı soyuducu axıcı çox sürətli istilik yüklərinə məruz qala bilər. Bu hallarda təhlükəsiz işi müəyyən edən kritik sual budur: Axıcı istiliyi divardan kifayət qədər sürətlə uzaqlaşdıra bilərmi, yoxsa qaynama başlayaraq sistem başqa istilikötürmə rejiminə keçəcək?

Qaynama həmişə pis deyil. Nüvəli qaynama bir çox halda çox güclü istilikötürmə mexanizmidir. Lakin qaynama nəzarətsiz davam edib film qaynamasına çevrilərsə, divarla maye arasında buxar təbəqəsi yarana bilər. Buxar maye ilə müqayisədə istiliyi xeyli zəif ötürdüyündən divar temperaturu sürətlə arta və səthin zədələnmə riski böyüyə bilər. Buna görə qaynamasının nə vaxt, harada və hansı birfazalı ilkin şərtlərdən sonra başladığını anlamaq kritik əhəmiyyət daşıyır.

Bu tədqiqat məhz qaynama başladıqdan sonrakı ikifazalı rejimi deyil, ondan dərhal əvvəlki birfazalı maye rejimini araşdırır. Yəni divar yaxınlığında maye yerli olaraq doyma temperaturunu aşsa belə, hələ qabarcıq yaranması başlamayıb. Bu qısa, lakin kritik dövr sistemin qaynamaya hansı yolla keçdiyini müəyyən edir.

Niyə “eksponensial divar istilik seli keçidi” vacibdir?

Araşdırmada divara verilən istilik seli zamanla eksponensial şəkildə artır. Bu, istilik yükünün tədricən xətti artmasından fərqlidir. Eksponensial artımda başlanğıcda dəyişmə nisbətən məhdud ola bilər; lakin zaman keçdikcə istilik seli çox sürətlə böyüyür. Bu hal nüvə təhlükəsizliyi təhlillərində və qəfil güc artımı ssenarilərində xüsusi əhəmiyyət daşıyır.

Araşdırmada istifadə olunan divar istilik seli \(q''_w(t)\) axıcıya keçən xalis divar istilik seli kimi nəzərdən keçirilir. Burada mühüm fərq budur: eksperimental qurğuda qızdırıcıya verilən elektrik gücü ilə axıcıya keçən xalis istilik seli eyni şey deyil. Qızdırıcı materialın, xüsusilə sapfir altlığın öz termal cavab müddəti var. Buna görə ən sürətli rampada, yəni \(\tau = 5 \, \mathrm{ms}\) halında, axıcıya keçən xalis istilik seli ideal eksponensial profilə tam uyğun gəlmir. Daha yavaş \(\tau = 100 \, \mathrm{ms}\) halında isə altlıq dəyişimi daha rahat izlədiyinə görə xalis divar istilik seli eksponensial profilə daha yaxın davranır.

Bu fərq araşdırmanın metodoloji dəqiqliyi baxımından vacibdir. Tədqiqatçılar simulyasiyalarda birbaşa elektrik gücünü deyil, eksperimental şəkildə emal edilmiş və axıcıya keçdiyi müəyyənləşdirilmiş zamandan asılı xalis divar istilik selini sərhəd şərti kimi istifadə ediblər. Beləliklə, simulyasiya real eksperimental istilik keçidini daha düzgün təmsil edir.

Araşdırma hansı fiziki sistemi araşdırır?

Araşdırma nazik düzbucaqlı kanallarda axan deionlaşdırılmış suyu öyrənir. Bu quruluş nüvə tədqiqat reaktorlarındakı dar soyutma kanallarını təmsil edən kiçildilmiş miqyaslı eksperimental qurğuya əsaslanır. İki sınaq hissəsindən istifadə olunub. Birincisi daha böyük kanal olan TS1-dir və ölçüləri \(3 \, \mathrm{cm} \times 1 \, \mathrm{cm}\)-dir. İkincisi daha dar kanal olan TS2-dir və \(18 \, \mathrm{mm} \times 1.95 \, \mathrm{mm}\) kanal aralığına malikdir.

Hər iki sınaq hissəsində qızdırıcı səthin uzunluğu 20 mm, eni isə 4 mm-dir. Qızdırıcı sapfir altlıq üzərində buxar üsulu ilə çökdürülmüş nano-hamar titan təbəqəsindən ibarətdir. Sapfirin infraqırmızı şüalanmanı qismən keçirməsi eksperimentlərdə divar temperaturunun infraqırmızı termoqrafiya ilə ölçülməsinə imkan verir. Sonra bu ölçmələrdə sapfir daxilində istilikkeçirmə təsirləri düzəldilir və nəticələr divar temperaturu ilə xalis divar istilik seli xəritələrinə çevrilir.

Axıcı təzyiqi 10 bar, giriş temperaturu isə \(353 \, \mathrm{K}\), yəni \(80^\circ \mathrm{C}\) qəbul edilib. Doyma temperaturuna nəzərən giriş subcooling-i \(\Delta T_{\mathrm{sub}} = 100 \, \mathrm{K}\) səviyyəsindədir. Bu, axıcının kanala daxil olarkən doyma temperaturundan 100 K daha soyuq olması deməkdir. Lakin divar yaxınlığında tətbiq edilən sürətli istilik seli mayenin yerli olaraq doyma temperaturunu aşmasına və superqızdırılmış maye zonalarının yaranmasına səbəb ola bilər.

Large-Eddy Simulation burada nə edir?

Araşdırmada istifadə olunan Large-Eddy Simulation turbulent axınların böyük miqyaslı burulğan strukturlarını birbaşa həll edən, çox kiçik miqyasları isə modelləşdirən yüksək dəqiqlikli CFD yanaşmasıdır. Birbaşa Ədədi Simulyasiya, yəni DNS bütün turbulentlik miqyaslarını həll etməyə çalışır; lakin bu problemdə divar yaxınlığında özlülük temperaturla sürətlə dəyişdiyindən DNS-in hesablama xərci çox yüksək olardı. Buna görə LES fiziki detallılıqla hesablama xərci arasında uyğun tarazlıq yaradır.

Simulyasiyalar TRUST/TrioCFD proqram təminatı ilə aparılıb və yüksək məhsuldarlıqlı hesablama resurslarından istifadə olunub. Ən aşağı Reynolds ədədi hallarında 10-40 milyon element, ən böyük hallarda isə \(2 \times 650\) milyon elementə qədər şəbəkə ölçülərindən istifadə edilib. Bu miqyasdakı həll araşdırmanın sadə korrelyasiya hesabı deyil, üçölçülü turbulentlik və istilikötürmə strukturunu ətraflı öyrənən yüksək ayırdetməli simulyasiya işi olduğunu göstərir.

PDF-nin 7-ci səhifəsindəki Figure 2.1 simulyasiya quruluşunu göstərir. Axın istiqaməti \(x\) oxu boyunca irəliləyir. Spanwise, yəni \(z\) istiqamətində periodik sərhəd şərtləri tətbiq edilir. Sınaq hissəsindən əvvəl recirculation box yerləşir; bu bölmə sınaq hissəsinə tam inkişaf etmiş turbulent axının daxil olmasını təmin edir. Sonra sınaq hissəsində aşağı divara zamandan asılı istilik seli verilir.

İstifadə olunan əsas tənliklər hansılardır?

Araşdırma sıxılmayan Navier-Stokes tənlikləri ilə temperatur tənliyini birlikdə həll edir. Temperatur passiv kəmiyyət kimi deyil, aktiv skalyar kimi nəzərdən keçirilir. Çünki temperatur artdıqca suyun özlülüyü \(\mu(T)\) və istilikkeçirmə əmsalı \(k(T)\) dəyişir; bu dəyişiklik də axın sahəsinə əks təsir göstərir.

Araşdırmada verilən əsas tənliklər belə yazıla bilər:

\[ \nabla \cdot \mathbf{u} = 0 \]

Bu, kəsilməzlik tənliyidir. Sıxılmayan axında sürət sahəsinin divergensiyası sıfırdır. Sadə dillə desək, axıcı hesablama sahəsində sıxılıb-genişlənən və sıxlığı dəyişən qaz kimi davranmır; kütlənin qorunması bu şərtlə təmin olunur.

Momentum tənliyinin quruluşu belədir:

\[ \frac{\partial \mathbf{u}}{\partial t} + \nabla \cdot (\mathbf{u}\mathbf{u}) = -\frac{1}{\rho}\nabla P + \nabla \cdot \left(\nu(T)\nabla \mathbf{u}\right) - \nabla \cdot \tau^{SGS} \]

Burada \(\mathbf{u}\) sürət vektorunu, \(t\) zamanı, \(P\) təzyiqi, \(\rho\) sıxlığı, \(\nu(T)\) temperaturdan asılı kinematik özlülüyü, \(\tau^{SGS}\) isə LES alt-şəbəkə miqyası gərginliyini ifadə edir. Tənlik axıcının sürət sahəsinin təzyiq, özlü yayılma və modelləşdirilən kiçik turbulentlik miqyasları ilə necə dəyişdiyini göstərir.

Temperatur tənliyi isə belədir:

\[ \frac{\partial T}{\partial t} + \mathbf{u}\cdot\nabla T = \frac{1}{\rho c_p}\nabla \cdot \left(k(T)\nabla T\right) + \nabla \cdot \left(\alpha_t \nabla T\right) \]

Bu tənlik temperaturun zamana görə dəyişməsini, axınla daşınmasını, molekulyar istilikkeçirmə ilə yayılmasını və alt-şəbəkə turbulent termal diffuziyanın töhfəsini göstərir. Burada \(T\) temperaturu, \(c_p\) xüsusi istilik tutumunu, \(k(T)\) temperaturdan asılı istilikkeçirməni, \(\alpha_t\) isə turbulent termal diffuziyanı ifadə edir.

Niyə temperaturdan asılı özlülük və istilikkeçirmə vacibdir?

Su yüksək istilik seli altında divar yaxınlığında çox sürətlə qızdıqda axıcının xüsusiyyətləri sabit qalmır. PDF-nin 9-cu səhifəsindəki Figure 2.2 10 bar təzyiqdə suyun temperaturdan asılı özlülüyünü, istilikkeçirməsini və Prandtl ədədini göstərir.

Bu qrafikə görə özlülük temperatur artdıqca təxminən eksponensial şəkildə azalır və araşdırılan temperatur intervalında təxminən üç dəfə aşağı düşür. Bu çox güclü təsirdir. Özlülüyün azalması divar yaxınlığındakı axının daha asan sürətlənməsinə, kəsici gərginliyin dəyişməsinə və turbulentlik strukturlarının yenidən təşkil olunmasına səbəb olur.

İstilikkeçirmə isə monoton davranmır. Təxminən \(416 \, \mathrm{K}\) ətrafında maksimuma çatır və ümumilikdə təxminən %10 dəyişir. Bu, daha kiçik təsir kimi görünə bilər; lakin divar yaxınlığındakı istilikötürmə birbaşa keçiriciliklə başladığından, istilikkeçirmədəki yerli azalma və ya artım divar istilikötürməsinə təsir göstərə bilər.

Prandtl ədədi isə belə müəyyən edilir:

\[ Pr = \frac{\mu c_p}{k} \]

Burada \(Pr\) momentum diffuziyası ilə istilik diffuziyasının nisbi əhəmiyyətini göstərir. \(\mu\) dinamik özlülük, \(c_p\) xüsusi istilik tutumu, \(k\) isə istilikkeçirmə əmsalıdır. Araşdırmada Prandtl ədədinin divar yaxınlığında 1-dən aşağı düşə bildiyi göstərilib. Bu hal ən kiçik termal və hidrodinamik miqyasların xarakterini dəyişir və divar yaxınlığındakı termal sərhəd təbəqəsinin həllini daha həssas edir.

Simulyasiyalar eksperimentlərlə necə doğrulanıb?

Araşdırmada simulyasiyalar eksperimental divar temperaturu ölçmələri ilə müqayisə edilərək doğrulanıb. Bunun üçün ölçüsüz divar temperaturu və ölçüsüz zaman istifadə olunub.

Ölçüsüz divar temperaturu belə müəyyən edilir:

\[ T_w^* = \frac{T_w - T_{\mathrm{sat}}}{\Delta T_{\mathrm{sub}}} \]

Burada \(T_w^*\) ölçüsüz divar temperaturudur. \(T_w\) divar temperaturunu, \(T_{\mathrm{sat}}\) doyma temperaturunu, \(\Delta T_{\mathrm{sub}}\) isə giriş subcooling-ni ifadə edir. Əgər \(T_w^* = 0\)-dırsa, divar temperaturu doyma temperaturuna bərabərdir. Müsbət dəyərlər divarın doyma temperaturundan yuxarı qalxdığını göstərir.

Ölçüsüz zaman isə belə müəyyən edilib:

\[ t^* = \frac{t - t_{\mathrm{ONB}}}{\tau} \]

Burada \(t_{\mathrm{ONB}}\) nüvəli qaynama başlanğıcı anıdır. \(t^* = 0\) qaynama başlanğıcı anını təmsil edir. Mənfi \(t^*\) dəyərləri qaynama öncəsini, müsbət dəyərlər isə eksperimentdə qaynama başladıqdan sonrakı ikifazalı bölgəni ifadə edir.

PDF-nin 10-cu səhifəsindəki Figure 3.1 eksperimental və ədədi \(T_w^*\) əyrilərini müqayisə edir. \(\tau = 5 \, \mathrm{ms}\), \(20 \, \mathrm{ms}\) və \(100 \, \mathrm{ms}\) hallarında simulyasiyaların eksperimentlərlə kifayət qədər yaxşı uyğunlaşdığı görünür. Birfazalı bölgədə xəta ümumilikdə %5-dən aşağı qalır. \(t^* > 0\) bölgəsində xəta artır; bu gözləniləndir, çünki simulyasiyalar birfazalıdır, eksperimentlərdə isə artıq qaynama başlayıb.

Bu doğrulama araşdırmanın etibarlılığı baxımından vacibdir. Simulyasiyalar qaynama öncəsi divar temperaturunun təkamülünü düzgün tuta bildiyinə görə, tədqiqatçılar daha sonra eksperimentlərdə birbaşa ölçülə bilməyən üçölçülü sürət, temperatur, turbulentlik quruluşu və istilik seli sahələrini təhlil edə biliblər.

Sürətli isitmə divar temperaturuna necə təsir edir?

Araşdırma sürətli isitmə rampalarının qaynama başlanğıcında daha yüksək divar superqızması yaratdığını göstərir. Quasi-steady, yəni daha yavaş rampalarda divar superqızması təxminən \(0.2\Delta T_{\mathrm{sub}}\), yəni təxminən 20 K olduğu halda, ən sürətli rampalarda bu dəyər \(0.65\Delta T_{\mathrm{sub}}\), yəni təxminən 65 K səviyyəsinə çatır.

Bu nəticə çox vacibdir. Eyni axıcı və kanal daxilində istilik yükünün nə qədər sürətlə artması qaynama başlamazdan əvvəl divarın nə dərəcədə superqıza biləcəyini dəyişir. Sürətli keçidlərdə maye və turbulentlik strukturu istilik yükünə uyğunlaşmaq üçün kifayət qədər zamana malik olmur. Buna görə divar yaxınlığında daha yüksək temperaturlar və daha güclü termal qeyri-tarazlıq yaranır.

Divar yaxınlığındakı turbulentlik quruluşu necə dəyişir?

PDF-nin 12-ci səhifəsindəki Figure 3.2 \(t^* = 0\) anında temperatur sahəsini və temperaturun zamana görə dəyişməsini göstərir. Şəkildə divar yaxınlığında yüksək temperaturlu filamentlərin axın istiqamətində uzandığı görünür. Bu quruluş turbulent kanal axınındakı aşağı sürətli streak nümunələri ilə əlaqəlidir. Bu streak-lərin spanwise aralığı təxminən 60-120 özlü vahiddir; bu da klassik divar turbulentliyi ədəbiyyatı ilə uyğun gəlir.

Eyni şəkil temperatur sahəsinin vahid təbəqə olmadığını göstərir. Soyuq, subcooled zonalar bulk axından daha böyük ləkələr şəklində gəldiyi halda, superqızdırılmış zonalar divara bağlı turbulentlik strukturları daxilində yaranır. Temperaturun zamana görə törəməsi \(\partial T/\partial t\) isti və soyuq zonaların harada sürətlə dəyişdiyini göstərir. Bu sahə turbulentliyin isti və soyuq mayeni qatlayıb qarışdırdığı yerlərdə ən intensivdir.

PDF-nin 13-cü səhifəsindəki Figure 3.3 \(\lambda_2\) meyarı ilə müəyyənləşdirilən turbulent burulğan strukturlarını göstərir. \(t^* = -2\) anında sürətli və yavaş rampalar demək olar ki, eynidir. Lakin \(t^* = 0\) anına çatdıqda iki hal ayrılır. Sürətli rampa \(\tau = 5 \, \mathrm{ms}\) halında divar yaxınlığındakı bütün zona doyma temperaturunu aşan qırmızı superqızdırılmış təbəqəyə çevrilir. Quasi-steady \(\tau = 100 \, \mathrm{ms}\) halında isə superqızdırılmış zonalar axın istiqamətində zolaqlar şəklində qalır.

Bu nəticə qaynama başlanğıcının məkan nümunəsini anlamaq baxımından çox vacibdir. Sürətli isitmə zamanı qaynama daha geniş və homogen səthdə başlaya bilər. Yavaş isitmə zamanı isə qaynama axın istiqamətində zolaqlı və lokal sahələrdən başlaya bilər. Yəni isitmə zaman miqyası yalnız “nə vaxt qaynayacaq?” sualına deyil, “harada və necə qaynayacaq?” sualına da təsir edir.

Divar yaxınlığındakı axın niyə sürətlənir?

Araşdırmanın dinamik təhlilində isitmə irəlilədikcə divar yaxınlığındakı axının sürətləndiyi göstərilib. PDF-nin 15-ci səhifəsindəki Figure 3.5 ölçüsüz axın istiqaməti sürəti \(U^+\) və spanwise vortisite \(\omega_z^+\) profillərini verir.

Divar turbulentliyində klassik özlü alt təbəqədə sürət profili təxminən belədir:

\[ U^+ = y^+ \]

Burada \(U^+ = U_x/u_\tau\) ölçüsüz axın istiqaməti sürətini, \(y^+\) isə divardan məsafənin özlü vahidlərlə ifadəsini göstərir. Araşdırmada \(t^* = 0\) anında bu klassik xətti profilin meylinin təxminən %30 artdığı bildirilir. Bu, divar yaxınlığındakı axının sürətləndiyini göstərir.

Bu sürətlənmənin əsas səbəbi temperatur artdıqca özlülüyün azalmasıdır. Özlülük azaldıqca divar yaxınlığındakı momentum tarazlığı dəyişir. Spanwise vortisite demək olar ki, iki dəfə artır. Bundan əlavə, \(y^+ = 60\) ilə \(140\) arasında vortisite platosu və buna uyğun olaraq axın istiqamətində demək olar ki, xətti sürət zonası yaranır. Müəlliflər bunu Couette-like shear layer, yəni Couette-bənzər kəsilmə təbəqəsi kimi şərh edirlər.

Bu nəticə çox vacibdir, çünki qaynama öncəsi proses yalnız temperatur sahəsinin inkişafı deyil. İstilik axının özünü yenidən təşkil edir; yenidən təşkil olunan axın isə istiliyin daşınma formasını dəyişir. Bu ikitərəfli qarşılıqlı təsir araşdırmanın əsas elmi töhfələrindən biridir.

Özlü gərginlik tenzoru nə göstərir?

Araşdırmada özlü gərginlik tenzorunun komponentləri də öyrənilib. İstifadə olunan ümumi ifadə belədir:

\[ \tau_{ij} = \left(\mu(T) + \mu_t\right) \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right) \]

Burada \(\tau_{ij}\) özlü gərginlik tenzorunun komponentini, \(\mu(T)\) temperaturdan asılı dinamik özlülüyü, \(\mu_t\) turbulent eddy özlülüyünü, \(u_i\) sürət komponentini, \(x_i\) isə koordinat komponentini göstərir. Bu düstur temperaturdan asılı özlülüyün momentum ötürülməsini birbaşa dəyişdirdiyini aydın göstərir.

PDF-nin 17-ci səhifəsindəki Figure 3.6 \(\tau_{xy}\) kəsici gərginlik komponentini göstərir. İsitmə irəlilədikcə divar yaxınlığındakı ölçüsüz kəsici gərginlik azalır. Bunun səbəbi temperatur artdıqca yerli özlülüyün azalmasıdır. Sürətli keçiddə \(y^+ \approx 60-140\) arasında aydın kəsici gərginlik platosu yaranır. Bu plato Figure 3.5-də görünən Couette-bənzər sürət profili ilə uyğun gəlir.

PDF-nin 18-ci səhifəsindəki Figure 3.7 normal gərginlik komponentləri \(\tau_{yy}\) və \(\tau_{zz}\) dəyərlərini göstərir. İsitmə irəlilədikcə \(\tau_{yy}\) daha mənfi olur; bu, divara normal istiqamətdə sıxılma deməkdir. \(\tau_{zz}\) isə daha müsbət olur; bu spanwise istiqamətdə yayılma deməkdir. Müəlliflərin şərhi belədir: divar yaxınlığında isti ləkələr low-speed streak strukturları daxilində inkişaf edir, spanwise istiqamətdə yayılır və divara normal istiqamətdə sıxılır. Beləliklə, istilik divar yaxınlığında daha çox sıxışıb qalır.

Temperatur profilləri nə göstərir?

Araşdırmada ölçüsüz maye temperaturu belə müəyyən edilir:

\[ T^* = \frac{T - T_{\mathrm{sat}}}{\Delta T_{\mathrm{sub}}} \]

Burada \(T^*\) ölçüsüz temperaturdur. \(T = T_{\mathrm{sat}}\) olduqda \(T^* = 0\) olur. \(T^* > 0\) mayenin yerli olaraq doyma temperaturunu aşdığını, yəni superqızdırılmış olduğunu göstərir. \(T^* < 0\) isə mayenin hələ subcooled olduğunu bildirir.

PDF-nin 19-cu səhifəsindəki Figure 3.8 sürətli və quasi-steady keçidlərdə divara normal temperatur profillərini müqayisə edir. Hər iki halda temperatur artımı əsasən divar yaxınlığında məhdud qalır. Təxminən \(y^+ \gtrsim 100\) bölgəsində bulk axın demək olar ki, izotermik qalır. Bu, keçici isitmə müddətinin qısa olduğunu və istiliyin kanal nüvəsinə tam nüfuz edə bilmədiyini göstərir.

Lakin sürətli və yavaş keçid arasında mühüm fərq var. \(t^* = 0\) anında sürətli keçiddə divar yaxınlığında \(T^* \approx 0.75\) dəyərinə çatıldığı halda, quasi-steady keçiddə \(T^* \approx 0.30\) ətrafındadır. Bu, sürətli keçiddə divar yaxınlığındakı mayenin xeyli daha çox superqızdırıldığını göstərir.

Eyni nəticə qaynamasının harada başlaya biləcəyi baxımından da vacibdir. Sürətli keçiddə bütün qızdırıcı səth doyma temperaturunun üstünə çıxa bilər. Quasi-steady keçiddə isə yalnız orta və çıxış zonası doyma səviyyəsinə çatdığı halda, giriş zonası subcooled qala bilər.

Stanton ədədi nəyi ölçür?

Araşdırmanın mühüm istilikötürmə göstəricilərindən biri Stanton ədədidir. Stanton ədədi divardan axıcıya keçən istilikötürmənin axıcının daşıma qabiliyyətinə nisbətini göstərən ölçüsüz kəmiyyətdir. Araşdırmada belə müəyyən edilir:

\[ St(x^*,t^*) = \frac{q''_w(t^*)} {\rho c_p U\left(T_w(x^*,t^*) - T_{\mathrm{in}}\right)} \]

Burada \(St\) Stanton ədədidir. \(q''_w\) divardan axıcıya keçən xalis istilik selini, \(\rho\) sıxlığı, \(c_p\) xüsusi istilik tutumunu, \(U\) orta axın sürətini, \(T_w\) divar temperaturunu, \(T_{\mathrm{in}}\) isə giriş temperaturunu göstərir.

PDF-nin 20-ci səhifəsindəki Figure 3.9 Stanton ədədinin axın istiqamətində necə dəyişdiyini göstərir. Bütün hallarda Stanton ədədi girişdə yüksəkdir və axın istiqamətində azalaraq plato dəyərinə yaxınlaşır. Bu plato termal cəhətdən tam inkişaf etmiş bölgəni göstərir.

Sürətli keçiddə \(\tau = 5 \, \mathrm{ms}\) Stanton ədədi \(t^* = -2\) ilə \(t^* = -1\) arasında təxminən %30 artır və sonra bu yüksək səviyyəyə yaxın qalır. Quasi-steady halda isə Stanton ədədi yalnız təxminən %8 dəyişir. Bu nəticə sürətli isitmə prosesinin istilikötürməni müvəqqəti olaraq gücləndirdiyini göstərir.

Lakin bu artım sadəcə “axın hər şeyi daha yaxşı soyudur” demək deyil. Araşdırmanın digər təhlilləri sürətli keçiddə divara normal konvektiv daşınmanın tam inkişaf etmədiyini, istiliyin daha çox divar yaxınlığında toplandığını və diffuziya/keçici enerji saxlanması təsirlərinin daha mühüm hala gəldiyini göstərir.

Termal giriş uzunluğu niyə dəyişir?

Termal giriş uzunluğu qızdırılan kanal boyunca temperatur sərhəd təbəqəsinin inkişaf edib müəyyən tam inkişaf etmiş quruluşa çatması üçün lazım olan məsafəni ifadə edir. Figure 3.9-a görə sürətli keçiddə Stanton ədədi platosu daha tez yaranır. Sürətli \(\tau = 5 \, \mathrm{ms}\) keçidində termal cəhətdən tam inkişaf etmiş zona təxminən \(x^* \approx 2\pi\) ətrafında başladığı halda, quasi-steady \(\tau = 100 \, \mathrm{ms}\) halında bu, yalnız çıxışa yaxın \(x^* \approx 5\pi\) ətrafında baş verir.

Bu fərq sürətli isitmə zamanı istilikötürmənin daha qısa məsafədə divar yaxınlığında cəmləşdiyini göstərir. Quasi-steady halda isə turbulent strukturlar istiliyi axın istiqamətində və divara normal istiqamətdə daha uzun məsafə boyu daşıya bilər; buna görə termal inkişaf daha uzun çəkir.

Temperatur profilləri niyə sigmoidal qanuna uyğunlaşır?

Araşdırmanın diqqətçəkən nəticələrindən biri termal cəhətdən tam inkişaf etmiş bölgədə divara normal temperatur profillərinin iki parametrli sigmoidal qanuna yaxşı uyğun gəlməsidir. Əvvəlcə ölçüsüz temperatur belə müəyyən edilir:

\[ \Theta = \frac{T - T_w}{T_{\mathrm{in}} - T_w} \]

Burada \(\Theta\) divar temperaturu ilə giriş temperaturu arasında normallaşdırılmış temperaturdur. Divara yaxın bölgədə \(T\) dəyəri \(T_w\)-yə yaxınlaşdığından \(\Theta\) kiçik olur. Bulk bölgəyə doğru \(T\), \(T_{\mathrm{in}}\)-ə yaxınlaşdıqca \(\Theta\) 1-ə yaxınlaşır.

Araşdırmada istifadə olunan sigmoidal ifadə belədir:

\[ \Theta(y^+) = \frac{(y^+)^k}{(y_0^+)^k + (y^+)^k} = \frac{1}{1+\exp\left[-k\left(\ln y^+ - \ln y_0^+\right)\right]} \]

Bu düsturda \(y_0^+\) termal keçidin orta nöqtəsidir; yəni \(\Theta = 1/2\) olduğu divar məsafəsini ifadə edir. \(k\) isə keçidin kəskinliyini müəyyən edir. Böyük \(k\) daha dik və divara daha çox sıxışmış temperatur keçidi deməkdir. Kiçik \(k\) daha geniş yayılmış və daha qalın termal keçidi göstərir.

PDF-nin 21-ci səhifəsindəki Figure 3.10 və 22-ci səhifədəki Table 3.1 bu sigmoidal uyğunlaşdırmanın \(R^2 > 0.99\) ilə çox yaxşı işlədiyini göstərir. Sürətli keçiddə \(k\) daha böyük, \(y_0^+\) isə daha kiçikdir. Bu, sürətli isitmə temperatur keçidini daha kəskin və divara daha yaxın etdiyini göstərir. Quasi-steady keçiddə isə temperatur perturbasiyası divardan daha uzağa yayıla bilir.

Yerli Péclet ədədi və \(\xi\) göstəricisi nə göstərir?

Araşdırma yerli olaraq daşınma və diffuziyanın nisbi əhəmiyyətini anlamaq üçün Péclet ədədinə əsaslanan göstəricidən istifadə edir. Yerli Péclet ədədi belə müəyyən edilir:

\[ Pe = \frac{\left|\rho c_p \mathbf{u}\cdot\nabla T\right|} {\left|\nabla\cdot(k\nabla T)\right|} \]

Burada surət advektiv, yəni axınla daşınan istilikötürmə terminini, məxrəc isə keçiricilik/diffuziya terminini təmsil edir. \(Pe\) böyükdürsə adveksiya üstünlük təşkil edir; kiçikdirsə diffuziya dominantdır.

Lakin \(Pe\) çox böyük dəyərlərə çata bildiyindən araşdırma bounded, yəni məhdudlaşdırılmış bir göstərici müəyyən edir:

\[ \xi = \frac{\left|\rho c_p \mathbf{u}\cdot\nabla T\right|} {\left|\rho c_p \mathbf{u}\cdot\nabla T\right| + \left|\nabla\cdot(k\nabla T)\right|} = \frac{Pe}{1+Pe} \]

Bu göstəricidə \(\xi = 0\) tam diffuziya dominant, \(\xi = 1\) tam adveksiya dominant, \(\xi = 1/2\) isə hər iki mexanizmin bərabər töhfə verdiyi haldır.

PDF-nin 23-cü səhifəsindəki Figure 3.11 superqızdırılmış sərhəd təbəqəsinin yuxarıdan görünüşünü \(\xi\) ilə rəngləndirərək göstərir. \(t^* = -1\) anında superqızdırılmış təbəqə əsasən diffuziya dominant streak-lər şəklində görünür. \(t^* = 0\) anında isə termal giriş uzunluğundan sonra maye demək olar ki, tamamilə doyma temperaturunun üstünə çıxıb və spanwise istiqamətdə daha homogen superqızdırılmış təbəqə yaranıb.

Bu şəklin əhəmiyyəti budur: qaynama öncəsi superqızdırılmış təbəqədə diffuziya dominant, qarışıq və adveksiya dominant bölgələr eyni anda mövcud olur. Yəni birfazalı istilikötürmə tək mexanizmin hökm sürdüyü sadə proses deyil.

İstilik tənliyinin terminləri necə ayrılıb?

Araşdırma istilik tənliyinin terminlərini enerji saxlanması, diffuziya və adveksiya töhfələrinə ayıraraq təhlil edir. PDF-nin 24-cü səhifəsindəki Figure 3.12 bu terminlərin divara normal profillərini göstərir. Terminlər \(q''_w/e\) ilə normallaşdırılıb; burada \(e\) kanalın yarım qalınlığıdır.

Təhlilə görə termal sərhəd təbəqəsini təxminən üç bölgəyə ayırmaq olar:

• Termal alt təbəqə: \(y^+ \lesssim 2\). Burada molekulyar diffuziya dominantdır.
• Buffer təbəqəsi: \(2 \lesssim y^+ \lesssim 25\). Burada diffuziya və adveksiya oxşar böyüklükdədir.
• Xarici təbəqə: \(y^+ \gtrsim 25\). Burada adveksiya üstünlük təşkil edir; lakin divar mənşəli termal perturbasiya bulk bölgəyə hələ zəif çatır.

Sürətli keçiddə operatorların minimumları divara daha yaxın yerləşir. Quasi-steady halda isə istilikötürmə quruluşu divardan daha uzağa yayıla bilir. Bu, sigmoidal temperatur uyğunlaşdırmasındakı nəticələrlə üst-üstə düşür.

İstiqamətli istilik selinin ayrılması nə göstərir?

Araşdırma istilik tənliyini sel formasında da yazır:

\[ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot \left(\mathbf{j}_q + \mathbf{j}_{adv}\right) \]

Diffuziv istilik seli və advektiv istilik seli belə müəyyən edilir:

\[ \mathbf{j}_q = -k\nabla T \]

\[ \mathbf{j}_{adv} = -\rho c_p \mathbf{u}(T - T_{\mathrm{in}}) \]

Burada \(\mathbf{j}_q\) istilikkeçirmə/diffuziya selini, \(\mathbf{j}_{adv}\) isə axınla daşınan istilik selini ifadə edir. Bu ayrılma istiliyin hansı istiqamətdə və hansı mexanizmlə daşındığını anlamaq üçün istifadə olunub.

PDF-nin 25-ci səhifəsindəki Figure 3.13 \(x\), \(y\) və \(z\) istiqamətlərində advektiv və diffuziv istilik sellərini müqayisə edir. İstilik seli böyüklükləri dörd tərtib fərqlənə bilər. Ümumi iyerarxiya belədir:

\[ j^{adv}_x > j^{adv}_z > j^{adv}_y \approx j^q_y > j^q_z > j^q_x \]

Bu sıralama ən böyük töhfənin axın istiqamətindəki advektiv daşınma olduğunu göstərir. Lakin bu töhfə istiliyi birbaşa divardan uzaqlaşdırmaqdan daha çox, istiliyi kanal boyunca aşağı axına daşıyır. Spanwise adveksiya ikinci ən böyük töhfədir və quasi-streamwise burulğanların aşağı sürətli streak-lərdən istiliyi yanlara yaymasında mühüm rol oynayır.

Divara normal istiqamətdə advektiv və diffuziv sellər oxşar böyüklükdədir. Bu xüsusilə vacibdir, çünki divardan bulk axına birbaşa istilik uzaqlaşdırılması divara normal istiqamətdəki ötürmə ilə bağlıdır. Quasi-steady halda divara normal advektiv sel daha güclüdür. Sürətli halda isə termal sahə divara daha çox sıxışmış qalır və diffuziya/keçici enerji saxlanması daha mühüm olur.

Sürətli və quasi-steady keçidlər qaynamaya eyni yolla gedirmi?

Araşdırmanın ən güclü nəticələrindən biri sürətli və quasi-steady isitmə proseslərinin ONB-yə eyni birfazalı yol üzərindən getmədiyini göstərməsidir. Quasi-steady keçiddə turbulent strukturlar istiliyi axın istiqamətində zolaqlar şəklində daşımağa davam edir. Superqızdırılmış maye streak strukturları ilə əlaqəli uzununa zolaqlar daxilində qalır. Bu halda qaynama başlanğıcı da daha lokal və zolaqlı ola bilər.

Sürətli keçiddə isə istilik yükü o qədər tez artır ki, turbulentliyin divara normal daşınması tam inkişaf edə bilmir. Stanton ədədi artır, termal giriş uzunluğu qısalır və divar yaxınlığında daha homogen superqızdırılmış təbəqə yaranır. Bu təbəqə metastabil ola bilər; yəni kiçik perturbasiya geniş sahədə nüvəli qaynamanı başlada bilər.

Bu fərq praktiki mühəndislik baxımından çox vacibdir. Eyni maksimum istilik selinə baxmaq kifayət deyil. İstilik selinin hansı zaman miqyasında artması qaynama başlanğıcı nümunəsini, divar superqızmasını və istilikötürmə mexanizmini dəyişə bilər.

Araşdırmanın gündəlik həyat baxımından mənası nədir?

Bu araşdırma çox texniki kanal axını və nüvə termohidravlikası problemi kimi görünsə də, gündəlik həyatda qarşılığı olan sadə fiziki həqiqəti izah edir: səthə çox sürətlə istilik versəniz, maye və axın sistemi bu istiliyi həmişə eyni sürətlə uzaqlaşdıra bilməz. İsitmə yavaşdırsa, axın strukturları istiliyi daşımağa vaxt tapır. İsitmə çox sürətlidirsə, istilik divar yaxınlığında toplanır, lokal superqızma yaranır və qaynama fərqli şəkildə başlayır.

Qazanı yavaş-yavaş qızdırmaqla onu qəfil çox güclü istilik mənbəyinə məruz qoymaq arasında necə fərq varsa, mikrokanal soyutması, reaktor soyutması və yüksək güclü elektronikanın soyudulmasında da oxşar fərq var. Lakin burada proses millisaniyələr ərzində və mikrometr miqyasında baş verdiyi üçün onu birbaşa müşahidə etmək çox çətindir. LES simulyasiyaları bu görünməyən keçidi üçölçülü şəkildə araşdırmağa imkan verir.

Keçmiş, bu gün və gələcək baxımından araşdırmanın əhəmiyyəti

Keçmiş baxımından keçici qaynama və yüksək istilik seli problemləri uzun müddətdir eksperimental şəkildə araşdırılır. Lakin bir çox eksperiment divar temperaturu və ya ümumi istilikötürmə əmsalı kimi orta kəmiyyətlər təqdim edir. Bu dəyərlər vacib olsa da, divar yaxınlığındakı turbulent temperatur strukturlarının necə inkişaf etdiyini birbaşa göstərmir.

Bu gün baxımından yüksək istilik selinə malik sistemlər getdikcə daha kritik olur. Elektron qurğular kiçildikcə güc sıxlığı artır. Nüvə təhlükəsizliyi təhlillərində sürətli güc keçidlərinin düzgün modelləşdirilməsi tələb olunur. Füzyon və raket texnologiyalarında səthlər çox intensiv və dəyişkən istilik yüklərinə məruz qalır. Buna görə qaynama öncəsi birfazalı sərhəd təbəqəsinin inkişafını düzgün anlamaq təhlükəsiz layihələndirmə üçün əsas bilikdir.

Gələcək baxımından bu araşdırma qabaqcıl CFD modelləri və yüksək məhsuldarlıqlı hesablama vasitəsilə keçici istilikötürmədə təhlükəsizlik sərhədlərinin daha düzgün proqnozlaşdırıla biləcəyini göstərir. Əgər sürətli istilik seli keçidlərində qaynama başlanğıcı daha dəqiq proqnozlaşdırılsa, soyutma sistemlərini daha təhlükəsiz və səmərəli layihələndirmək olar. Araşdırma həmçinin klassik steady-state istilikötürmə korrelyasiyalarının sürətli keçici rejimlərdə kifayət etməyə biləcəyini xatırladır.

Araşdırmanın güclü tərəfləri hansılardır?

Araşdırmanın ən güclü tərəflərindən biri eksperimental divar temperaturu ölçmələri ilə doğrulanmış yüksək ayırdetməli LES simulyasiyalarından istifadə etməsidir. Bunun sayəsində yalnız orta divar temperaturu deyil, üçölçülü sürət, temperatur, vortisite, özlü gərginlik və istilik seli strukturları araşdırıla bilib.

İkinci güclü tərəf temperaturdan asılı özlülüyü və istilikkeçirməni nəzərə almasıdır. Yüksək istilik seli altında suyun xüsusiyyətləri sabit deyil. Bu dəyişən xüsusiyyətlər nəzərə alınmasa, divar yaxınlığındakı sürətlənmə, turbulentlik quruluşu və istilikötürmə natamam təmsil oluna bilər.

Üçüncü güclü tərəf sürətli və quasi-steady isitmə rampalarını açıq şəkildə müqayisə etməsidir. Beləliklə, isitmə zaman miqyasının yalnız temperatur səviyyəsini deyil, turbulentlik strukturunu və qaynama başlanğıcı nümunəsini də dəyişdirdiyi göstərilib.

Dördüncü güclü tərəf temperatur profillərini iki parametrli sigmoidal qanunla təmsil etməsidir. Belə kompakt təqdimatlar gələcəkdə daha aşağı xərcli mühəndislik modelləri və ya divar funksiyası təkmilləşdirmələri üçü