Ин таҳқиқот реактори даври кушодаи нигоҳдории энергияи термохимиявиро дар асоси гидрати бромиди стронсий бо модели ададии сеченака меомӯзад, ки метавонад барои нигоҳдории дарозмуддати гармии офтобии ҳарорати паст ё гармии партов истифода шавад. Дар реактор зарраҳои сахти SrBr₂·H₂O зери таъсири ҷозиба ба поён ҳаракат мекунанд, дар ҳоле ки ҳавои намнок ба самти уфуқӣ аз дохили бистар мегузарад. Реаксияи экзотермии гидрататсия байни буғи об ва намак энергияи химиявии захирашударо дубора ба гармӣ табдил медиҳад ва ҳарорати ҳавои аз реактор хориҷшавандаро баланд мекунад.

Муҳаққиқон дар муҳити COMSOL Multiphysics модели сеченакае таҳия карданд, ки ҷараёни муҳити ковок, интиқоли буғи об, кинетикаи реаксия, интиқоли гармӣ дар фазаҳои сахт ва газ ва бистари ҳаракаткунандаи намакро муттаҳид мекунад. Ҳаракати поёнии зарраҳои сахт бо усули шабакаи ҳаракаткунанда бо номи Arbitrary Lagrangian–Eulerian (ALE) намояндагӣ карда шуд. Ҳамин тавр фронти реаксия, майдони ҳарорат, тақсимоти буғи об, талафоти фишор ва табдили маҳаллии намак дар дохили реактор вобаста ба вақт ва ҷой пайгирӣ карда шуданд.

Натиҷаи асосии модел ин аст, ки пурқувваттарин тағйирёбандаи амалиётӣ барои муайян кардани кори реактор миқдори намии ҳавои воридшаванда мебошад. Вақте намии воридшаванда аз 4 g/kg то 12 g/kg зиёд шуд, ҳарорати бистар тақрибан аз 26°C то 39°C ва гармии муфиди ҳаҷмие, ки ба ҳаво дода мешавад, тақрибан аз 2 kW/m³ то 5,5 kW/m³ зиёд шуд. Намии баландтар фишори қисмии буғи об ва қувваи термодинамикии пешбарандаи гидрататсияро зиёд карда, минтақаи реаксияи гармтар, васеътар ва равшантар ба вуҷуд овард.

Дар дебити ҳаво натиҷаи мутавозинтарин дар шароити 110 m³/соат ба даст омад. Дар 80 m³/соат ҳаво бо намак муддати дарозтар тамос дошт, аммо иқтидори умумии интиқоли гармӣ паст монд. Дар 140 m³/соат бошад ҳаво тезтар ҳаракат кард, аммо аз сабаби вақти кӯтоҳи тамос ва сардкунии қавии конвективӣ реаксия ва ҳарорати бистар маҳдуд шуданд. Суръати ҳаракати бистари сахт дар байни 0,5–2 cm/соат ҳарорат ва баромади фаврии гармиро танҳо ба андозаи маҳдуд тағйир дод; дар иваз, ҳаракати тезтари сахт вақти мондани намакро дар реактор кӯтоҳ карда, хатари гидрататсияи нопурраро зиёд намуд.

Модел бо маълумоти ҳарорат, намӣ ва баромади гармии таҳқиқоти қаблан нашршудаи реактори бистари ҳаракаткунанда муқоиса шуд. Тамоюлҳо ва бузургиҳои умумӣ бо таҷрибаҳо мувофиқ буданд. Бо вуҷуди ин, ҳисобкунии мустақими ҷадвали validation дар PDF нишон медиҳад, ки фарқҳои нисбӣ дар баъзе шароитҳо аз диапазони 1–4 дарсаде, ки дар abstract гузориш шудааст, гузашта, барои гармии муфид то 8 дарсад ва барои ҳарорати бистар тақрибан то 7 дарсад мерасанд.

Таҳқиқот реактори нави физикӣ ё маълумоти нави таҷрибавӣ пешниҳод намекунад. Натиҷаҳо ба моделсозии сеченакаи setup-и таҷрибавии дар адабиёт буда ва симулятсияҳои ададии параметрӣ такя мекунанд. Устувории дарозмуддати даврии мавод, agglomeration, wear, deliquescence, шикастани зарраҳо ва назорати воқеии plant омӯхта нашудаанд. Таҳқиқот preprint-и аз peer review нагузашта мебошад.

Чаро нигоҳдории энергияи термохимиявӣ муҳим аст?

Энергияи офтобӣ, гармии партови саноатӣ ва дигар манбаъҳои гармии барқароршаванда на ҳамеша дар ҳамон лаҳзае, ки тавлид мешаванд, истифода мешаванд. Хусусан дар гармидиҳии биноҳо, имкони истифодаи гармии дар тобистон ё рӯзона гирифташуда пас аз ҳафтаҳо ё моҳҳо усулҳои нигоҳдории дорои зичии баланди энергия ва талафоти ками интизориро талаб мекунад.

Технологияҳои нигоҳдории энергияи гармӣ умуман ба се гурӯҳ ҷудо мешаванд:

• Нигоҳдории гармии ҳисшаванда: Ба баланд кардани ҳарорати маводе мисли об, санг ё бетон такя мекунад.
• Нигоҳдории гармии ниҳонӣ: Энергияи гудозиш ё сахтшавии маводи phase-change-ро истифода мебарад.
• Нигоҳдории энергияи термохимиявӣ: Гармиро дар маҳсулоти як реаксияи химиявии баргардон нигоҳ медорад.

Дар нигоҳдории термохимиявӣ маҳсулоти реаксияи зарядшуда метавонанд дар ҳарорати муҳит аз ҳам ҷудо нигоҳ дошта шаванд, бинобар ин дар давраи дарози интизорӣ нигоҳ доштани зарфи доиман гарм лозим нест. Энергия ба ҷойи фарқи ҳарорат дар таркиби химиявӣ нигоҳ дошта мешавад. Аз ин рӯ, талафоти гармӣ дар нигоҳдории дарозмуддат ё мавсимӣ метавонад хеле паст бошад.

Гидратҳои намак чӣ гуна гармиро нигоҳ медоранд?

Ҷуфти реаксияе, ки дар таҳқиқот истифода шудааст, байни strontium bromide monohydrate ва strontium bromide hexahydrate мебошад. Дар PDF реаксия ба шакли зерин дода шудааст:

\[ SrBr_2\cdot H_2O(s) + 5H_2O(g) \rightleftharpoons SrBr_2\cdot 6H_2O(s) + 5\Delta H \]

Дар марҳилаи зарядкунӣ ба намаки гидратшуда ҳавои гарм ва хушк дода шуда, қисми оби пайвастшуда хориҷ карда мешавад. Барои ин dehydratation-и эндотермӣ гармӣ аз берун дода мешавад ва энергия дар сохтори намаки dehydrateшуда ҳамчун потенсиали химиявӣ нигоҳ дошта мешавад.

Дар марҳилаи discharge бошад ҳавои хунук ва намнок бо намаки dehydrateшуда тамос мегирад. Намак буғи оби ҳаворо боз пайваст карда hydrate мешавад. Азбаски hydration экзотермӣ аст, гармии озодшуда ҳаворо гарм мекунад. Ҳавои гарми хориҷӣ метавонад тавассути heat exchanger-и ҷудогона ба бино ё district-heating system энергия диҳад.

Ба амал омадани реаксия танҳо ба мавҷуд будани об дар ҳаво вобаста нест. Фишори воқеии буғи об дар ҳаво бояд аз фишори буғи equilibrium-и намак дар ҳарорати мавҷуда баландтар бошад. Аз ин рӯ намӣ, ҳарорат ва equilibrium pressure ба ҳам вобастаанд.

Фарқи системаҳои термохимиявии кушода ва пӯшида

Дар системаи термохимиявии пӯшида намак ва буғи тозаи об одатан зери vacuum реаксия мекунанд. Ин сохтор метавонад бештар назоратшаванда бошад; аммо evaporator, condenser, liquid-water tank, vacuum equipment ва sealing талаб мекунад.

Дар системаи кушодаи омӯхташуда ҳамчун gas reactant ҳавои намнок дар фишори атмосферӣ истифода мешавад. Бо ин роҳ ниёз ба гардиши ҷудогонаи буғи об кам мешавад. Дар иваз, система ба намии муҳити беруна хеле ҳассос мегардад. Намии аз ҳад зиёд метавонад боиси часпидани намак, agglomeration ё deliquescence гардад; намии паст бошад суръати нокифояи реаксияро ба вуҷуд меорад.

Реактори бистари ҳаракаткунанда чӣ гуна кор мекунад?

Аз қисми болоии реактор зарраҳои хушк ё камгидратшудаи SrBr₂ дода мешаванд. Бистари сахт зери таъсири ҷозиба амудӣ ба поён ҳаракат мекунад. Rotary valve дар поён суръати хориҷшавии сахт ва аз ин рӯ суръати сахтро назорат мекунад.

Ҳавои намнок аз як тарафи реактор ворид шуда, ба самти cross-flow аз бистари сахт гузашта, аз тарафи муқобил мебарояд. Дар ин cross-flow arrangement:

• Фазаи сахт аз боло ба поён,
• Фазаи газ аз чап ба рост,
• Буғи об аз ҳаво ба сахт,
• Гармии реаксия аз сахт ва минтақаи реаксия ба ҳаво

интиқол дода мешавад.

Дар схемаи реактори PDF, solid feed hopper дар боло, moving reactive bed дар марказ, solid discharge mechanism дар поён ва humid-air connections дар ду тарафи reactor body ҷой доранд. Бистари танги амудии марказӣ бо ду ҳаҷми паҳлуии васеъшавандаи тақсимоти ҳаво иҳота шудааст.

Андозаҳои асосии геометрияи реактор

Ду девори stainless steel дар болои reactive bed барои пешгирӣ аз он истифода шудаанд, ки ҳаво аз void region-и боло мустақиман ба outlet гузарад ва барои маҷбур кардани ҷараён ба дохили salt bed. Ин деворҳо ҳамчунин hydraulic resistance-и иловагӣ эҷод карданд.

Холигии илмие, ки таҳқиқот ҳадаф гирифтааст

Таҳқиқоти қаблии moving-bed асосан моделҳои дученака, lumped-parameter ё ба реаксияҳои ҳарорати баланд ва irreversible нигаронидашударо истифода бурдаанд. Ин равишҳо дар гирифтани ҳамзамони таъсирҳои зерин бо тафсил маҳдуд буданд:

• Ҳаракати пайвастаи solid bed,
• Тақсимоти сеченакаи ҳаво ва side leakage,
• Кинетикаи реаксия вобаста ба partial pressure-и буғи об,
• Пешравии reaction front дар фазо ва вақт,
• Таъсири local porosity ва permeability,
• Муносибати байни solid residence time ва conversion,
• Trade-off байни heat output ва total energy-storage efficiency.

Ҳадафи таҳқиқот якҷо кардани ин равандҳо дар як модели сеченакаи continuous-operation ва муқоисаи натиҷаҳои модел бо маълумоти қаблан нашршудаи experimental reactor мебошад.

Шароити кори истифодашуда дар модел

Фарзияҳои асосии модел

• Heat ва mass transport асосан аз ҷониби airflow кашонда мешаванд.
• Азбаски solid bed хеле суст ҳаракат мекунад, таъсири бозгаштии он ба air dynamics сарфи назар шудааст.
• Gas flow steady ва дар inlet cross-section uniform қабул шудааст.
• Solid particles бо суръати доимӣ ба поён ҳаракат мекунанд.
• Фарз шудааст, ки намак дар давоми reaction chemical reactivity ва mechanical structure-и худро нигоҳ медорад.
• Bed porosity, ба истиснои ҳолатҳое, ки bypass effect алоҳида таъриф шудааст, доимӣ қабул шудааст.
• Humid air ҳамчун ideal-gas mixture модел шудааст.
• Long-term cyclic degradation, agglomeration, particle breakage ва wear ба модел дохил нашудаанд.

Нигоҳдории буғи об ва табдили намак

Тавозуни буғи об дар porous bed аз accumulation, consumption by reaction ва transport by gas flow иборат аст. Шакли мухтасари notation-и PDF чунин аст:

\[ \varepsilon n_h\frac{\partial y_v}{\partial t} = \phi_v(1-y_v) - n_h\mathbf{u}_h\cdot\nabla y_v \]

• ε: Bed porosity,
• n_h: Терме, ки molar density-и humid air-ро намояндагӣ мекунад,
• y_v: Mole ё mass fraction-и water vapor,
• φ_v: Vapor source term, ки аз reaction ҳосил ё consume мешавад,
• u_h: Velocity vector-и humid air.

Тавозуни local hydration conversion-и намак:

\[ (1-\varepsilon)\frac{\partial X}{\partial t} = -\frac{\phi_v}{\nu} - \mathbf{u}_s\cdot\nabla X \]

• X: Local hydration conversion,
• ν: Stoichiometric coefficient,
• u_s: Downward velocity-и solid bed.

Дар equation-и дуюм term-и аввал conversion бо chemical reaction ва term-и дуюм transport-и converted ё unconverted solid-ро дар тӯли bed намояндагӣ мекунад.

Ҷараёни муҳити ковок: муодилаи Brinkman

Airflow дар reactor бо Brinkman equation модел шудааст, ки ба Darcy law viscous shear effects-ро илова мекунад:

\[ \rho(\mathbf{u}\cdot\nabla)\mathbf{u} = \nabla\cdot\left[-p\mathbf{I} +\mu\left(\nabla\mathbf{u}+(\nabla\mathbf{u})^T\right)\right] -\mu\alpha\mathbf{u} +\rho\mathbf{F} \]

• ρ: Зичии gas mixture, kg/m³,
• u: Velocity vector, m/s,
• p: Pressure, Pa,
• μ: Dynamic viscosity, Pa·s,
• α: Терми вобаста ба Darcy resistance-и porous medium, m⁻²,
• F: Body force, ки ба unit mass таъсир мекунад.

Equation барои ҳисоб кардани high-speed jets дар side-air inlet, circulation дар expansion regions, pressure drop дар salt bed ва viscous resistance дар outlet channel истифода шудааст.

Интиқоли гармӣ ва гармии реаксия

Energy equation дар gas phase convection, conduction ва volumetric heat sources-ро якҷо дар бар мегирад:

\[ \rho C_p\mathbf{u}\cdot\nabla T+\nabla\cdot\mathbf{q} = Q+Q_p+Q_{vd} \]

\[ \mathbf{q}=-k\nabla T \]

• C_p: Specific heat capacity, J/(kg·K),
• T: Temperature, K ё фарқи °C,
• q: Conductive heat flux, W/m²,
• Q: Volumetric heat source-и chemical reaction, W/m³,
• Q_p: Pressure work,
• Q_vd: Viscous dissipation term.

Вақте ҳаракат ва convection дар solid phase сарфи назар мешаванд, energy equation:

\[ \rho C_p\frac{\partial T}{\partial t}+\nabla\cdot\mathbf{q} = Q+Q_{ted} \]

ба ин шакл дода шудааст. Thermal expansion term Q_ted дар таҳқиқот сарфи назар шудааст.

Exothermic heat source-и hydration:

\[ Q=R_{react}\Delta H_{hydration} \]

ба ин шакл татбиқ шудааст. Бо зиёд шудани reaction rate volumetric heat-и дар ҳамон region озодшуда низ зиёд мешавад.

Кинетикаи реаксия вобаста ба намӣ ва ҳарорат

Барои local reaction rate equation-и Arrhenius form бо vapor-pressure driving force истифода шудааст:

\[ X_a=A_r\exp\left(-\frac{E_a}{RT}\right)(1-X) \left(1-\frac{p_v}{p_{eq}}\right) \]

• X_a: Local reaction rate, 1/s,
• A_r: Pre-exponential coefficient; дар PDF 1,63×10⁴ дода шудааст,
• E_a: Activation energy, 56,6 kJ/mol,
• R: Universal gas constant,
• T: Absolute temperature, K,
• X: Local conversion,
• p_v: Partial pressure-и water vapor,
• p_eq: Vapor pressure барои salt-hydrate equilibrium.

Ҷузъиёти татбиқи sign дар equation дар hydration direction дар PDF пурра шарҳ дода нашудааст. Матн мегӯяд, ки heat source танҳо дар region-ҳое фаъол аст, ки water-vapor pressure аз equilibrium threshold баландтар буда, air-salt contact вуҷуд дорад.

Фишори буғи мувозинатӣ

Equilibrium pressure бо Clausius–Clapeyron-type relation ҳисоб шудааст:

\[ P_{eq}=P_{ref}\exp\left(-\frac{\Delta H^0}{RT} +\frac{\Delta S^0}{R}\right) \]

• P_ref: Reference pressure,
• ΔH⁰: Standard reaction enthalpy; дар model table 67.400 J/mol,
• ΔS⁰: Standard reaction entropy; дар model table 175 J/(mol·K).

Ин relation муайян мекунад, ки бо баланд шудани temperature equilibrium condition-и намак бо water vapor чӣ гуна тағйир меёбад. Фарқ байни actual vapor pressure дар reactor ва P_eq таъсири асосии thermodynamic мебошад, ки hydration-ро пеш мебарад.

Хусусиятҳои самараноки гармии бистари ковок

Volumetric average thermal properties-и gas ва solid phase ба шакли зерин муттаҳид шудаанд:

\[ (\rho C_p)_{eff} = (1-\varepsilon)(\rho C_p)_s +\varepsilon(\rho C_p)_v \]

\[ k_{eff} = (1-\varepsilon)k_s+\varepsilon k_v \]

Хусусиятҳои намак вобаста ба conversion байни monohydrate ва hexahydrate хаттӣ тағйир дода шудаанд:

\[ (\rho C_p)_s = (1-X)(\rho C_p)_{s1} +X(\rho C_p)_{s0} \]

\[ k_s=(1-X)k_{s1}+Xk_{s0} \]

Дар PDF гуфта шудааст, ки subscripts s1 ва s0 мутаносибан hexahydrate ва monohydrate-ро нишон медиҳанд. Аммо дар equations маълум нест, ки X дар кадом endpoint кадом hydrate-ро ифода мекунад; ин definition барои independent reproduction бояд равшан карда шавад.

Гармии муфиди ба ҳаво додашуда

Гармии муфиде, ки outlet air мебарад, аз sensible enthalpy ва фарқи enthalpy-и вобаста ба water vapor ба шакли зерин ҳисоб шудааст:

\[ q_{heat} = \dot{m}_{a,in} \left[ (C_{pa}+w_{out}C_{pv})T_{out} - (C_{pa}+w_{in}C_{pv})T_{in} \right] \]

• \(\dot{m}_{a,in}\): Mass flow rate-и dry inlet air, kg/s,
• C_pa: Specific heat-и dry air,
• C_pv: Specific heat-и water vapor,
• w_in, w_out: Inlet ва outlet absolute humidities,
• T_in, T_out: Inlet ва outlet temperatures.

Дар result graphs ин қимат нисбат ба reactor volume бо kW/m³ гузориш шудааст.

Шартҳои сарҳадӣ

Консентратсияи water vapor ва partial pressure-и inlet:

\[ c_{in}=\frac{P_v}{RT_{in}} \]

\[ P_v=P_s\cdot RH \]

бо ин муносибатҳо сохта шудаанд. Дар ин ҷо P_s saturated water-vapor pressure дар inlet temperature ва RH relative humidity мебошад.

Усули ALE moving mesh

Ҳаракати поёнии moving bed ба ҷойи нигоҳ доштани geometry ҳамчун fixed solid-flow domain бо усули Arbitrary Lagrangian–Eulerian модел шудааст, ки ба mesh nodes имкон медиҳад displacement-и solid-ро пайгирӣ кунанд.

Дар татбиқи ALE:

• Дар solid outlet downward normal mesh velocity таъриф шуд,
• Дар side walls slip conditions бе tangential constraint истифода шуданд,
• Барои пешгирӣ аз deformation-и номатлуби mesh symmetry ё roller constraints татбиқ шуданд,
• Дар mesh deformation усули Yeoh ва stiffening бо қимати 10 истифода шуд.

Ин moving mesh ҳамзамон бо flow, water-vapor transport, reaction ва energy equations ҳал карда шуд. Ҳамин тавр ворид шудани намаки нав ба reactor ва ба поён интиқол ёфтани намаки қисман hydrateшуда модел гардид.

Мустақилияти шабака

Муҳаққиқон аз сабаби тавозуни байни computational cost ва accuracy шабакаи 644.783-element-ро истифода карданд. Бо вуҷуди ин, дар ҷадвал равшан таъриф нашудааст, ки deviations нисбат ба кадом mesh pair ё кадом final reference ҳисоб шудаанд.

Матн мегӯяд, ки дар three-dimensional computational domain “free triangular elements” истифода шудаанд. Барои three-dimensional volume mesh истилоҳи tetrahedral elements интизор меравад; mesh visual дар PDF ҳам ба tetrahedral volume mesh монанд аст. Ин terminology барои reproducibility бояд равшан карда шавад.

Муқоисаи модел бо маълумоти таҷрибавӣ

Модел аз ҷиҳати temperature, inlet-outlet humidity, useful heat ва pressure drop бо moving-bed experiments-и қаблан нашршуда муқоиса шуд. Афзоиши inlet humidity аз 0 то 9 g/kg дар тақрибан 0,22 соати аввал дар experiment ба модел бо time-dependent interpolation function, ки аз experimental graph рақамисозӣ шудааст, ворид карда шуд.

Дар temperature comparison graph ду measurement line дар дохили bed ва outlet temperature тақрибан дар тӯли ҳашт соат нишон дода шудаанд. Simulation temperature curves experimental trends-ро пайгирӣ мекунанд. Outlet temperature дар simulation тақрибан 32,7°C ва дар experiment тақрибан 31°C аст. Internal bed temperatures тақрибан дар диапазони 35–37°C қарор доранд.

Дар water-vapor graph inlet humidity ба тақрибан 0,009 kg/kg ва outlet humidity ба тақрибан 0,005 kg/kg мерасад. Пасттар будани outlet humidity аз inlet нишон медиҳад, ки water vapor аз ҷониби намак consume мешавад.

Арзёбии муфассали ҷадвали validation

Relative differences, ки аз рақамҳои ҷадвал ҳисоб шудаанд, барои useful heat тақрибан дар диапазони 1,3–8,0 дарсад ва барои bed temperature тақрибан дар диапазони 0,46–7,07 дарсад мебошанд. Аз ин рӯ, гарчанде модел general trends ва magnitudes-ро мегирад, ибораи “± 1–4 дарсад” дар abstract ҳамаи validation conditions-ро фаро намегирад.

Тақсимоти фишор ва сохтори ҷараёни паҳлуӣ

Streamlines ва pressure distribution нишон медиҳанд, ки airflow дар narrow inlet connection channel суръат мегирад ва ҳангоми расидан ба wide side volume large circulation regions ташкил медиҳад. Дар inlet side-и reactive bed фишор тақрибан то 1.600 Pa мерасад ва дар тӯли porous medium пайваста кам мешавад.

Total pressure drop байни model ва experiment тақрибан дар диапазони 1.500–1.700 Pa мувофиқ аст. Ин натиҷа нишон медиҳад, ки илова ба bed resistance, air-guiding walls ва expansion-contraction geometry низ барои hydraulic performance муҳиманд.

Streamlines нишон медиҳанд, ки тамоми side volumes ба таври яксон истифода намешаванд; дар баъзе regions circulation ва short-circuit flow метавонад ба вуҷуд ояд. Аз ин рӯ bypass factor дар таҳқиқот ҳамчун performance variable-и ҷудогона баррасӣ шудааст.

Таъсири суръати сахт

Вақте solid velocity аз 0,5 cm/соат то 2 cm/соат зиёд шуд, bed temperature тақрибан аз 35°C то 35,5°C баланд шуд. Таъмини тезтари намаки хунук ва камгидратшуда метавонад local reaction-ро зиёд кунад, аммо ҳамзамон residence time-и ҳар як solid mass-ро дар bed кам мекунад.

Outlet temperature ва useful heat curves-и ба ҳаво интиқолёфта барои чор суръат ба ҳам хеле наздиканд. Useful heat пас аз афзоиши тези ибтидоӣ тақрибан дар 3 kW/m³ устувор мешавад.

Дар solid velocities-и паст, азбаски намак дар reactor муддати дарозтар мемонад, hydration-и пешрафтатар ба ҳар unit solid mass имконпазир аст. Дар velocities-и баландтар гарчанде solid throughput бештар шуда метавонад, conversion-и ҳар particle пасттар мемонад. Ин байни power output ва recovery-и total stored chemical energy trade-off эҷод мекунад.

Бо вуҷуди ин, дар PDF graph-и мустақими total conversion ё energy efficiency бар зидди solid velocity дода нашудааст. Аз ин рӯ то чӣ андоза velocity-и баланд conversion efficiency-ро паст мекунад, quantitatively нишон дода намешавад.

Рушди вақти фронти реаксия

Дар оғоз тамоми bed тақрибан 22°C аст. Пас аз оғоз шудани humid-air inlet, reaction аввал дар region-и наздик ба inlet пайдо мешавад:

• t = 0 соат: Reactor хунук ва дар ҳолати пеш аз reaction аст.
• t = 0,2 соат: Дар inlet side yellow-green temperature region ба вуҷуд меояд.
• t = 0,3 соат: Reaction region дар тӯли bed васеъ мешавад ва hot core равшантар мегардад.
• t = 1 соат: Қисми муҳими reactive bed тақрибан ба диапазони 34–38°C мерасад.

Contour visuals нишон медиҳанд, ки temperature аз air inlet ба outlet direction баланд шуда, дар тӯли central vertical reactive bed пеш меравад. Region-е, ки reaction оғоз мешавад, он ҷойест, ки фарқ байни water-vapor partial pressure ва equilibrium pressure ба қадри кофӣ калон аст.

Дар матн ин initial state иштибоҳан ҳамчун “пеш аз оғоз шудани dehydration” таъриф шудааст. Азбаски марҳилаи омӯхташуда discharge ва hydration аст, дар ин ҷо ибораи “пеш аз оғоз шудани hydration” мутобиқтар мебошад.

Таъсири дебити ҳаво

Air flow rate барои 80, 110 ва 140 m³/соат муқоиса шудааст. Натиҷаҳо муносибати якнавохти “flow бештар беҳтар” нишон намедиҳанд.

80 m³/соат

Дар low flow rate ҳаво бо намак муддати дарозтар тамос мегирад. Гарчанде барои hydration contact time кофӣ аст, total air mass flow rate паст мебошад ва интиқоли reaction heat аз reactor маҳдуд мемонад. Useful heat тақрибан 3 kW/m³ аст.

110 m³/соат

Medium flow rate беҳтарин тавозунро байни contact time-и кофӣ ва heat-carrying capacity-и кофӣ ташкил дод. Bed ва outlet temperature баландтарин буданд ва useful heat тақрибан 4–4,5 kW/m³, яъне баландтарин байни се condition буд.

140 m³/соат

High flow rate ҳавои бештар ба reactor меорад, аммо residence time-и gas дар bed-ро кам мекунад. Strong convective cooling боиси bed-и сардтар ва uneven temperature distribution мегардад. Outlet temperature аз condition-и 110 m³/соат пасттар аст.

Velocity contours нишон медиҳанд, ки air-inlet jet бо flow rate қавитар мешавад; дар 140 m³/соат local velocities дар inlet connection тақрибан то 2,5 m/s мерасанд. High local velocity ҳамеша маънои conversion-и баландтарро надорад; зеро барои reaction gas-solid contact time низ лозим аст.

Таъсири намии воридшаванда

Намӣ пурқувваттарин result variable-и таҳқиқот аст. Вақте air flow rate 110 m³/соат ва solid velocity 0,5 cm/соат нигоҳ дошта шуд:

Бо зиёд шудани humidity water-vapor partial pressure баланд мешавад, hydration driving force нисбат ба equilibrium pressure меафзояд ва дар ҳамон time interval намаки бештар reaction мекунад. Дар натиҷа:

• Reaction front тезтар мешавад,
• High-temperature region васеъ мешавад,
• Water-vapor consumption зиёд мешавад,
• Outlet air гармтар мешавад,
• Heat transferred to air меафзояд.

Майдони консентратсияи буғи об

Дар water-vapor contours, дар inlet humidity-и 4 g/kg gas concentration тақрибан 0,2–0,3 mol/m³ аст. Дар ҳолати 9 g/kg тақрибан 0,5–0,7 mol/m³ ва дар ҳолати 12 g/kg тақрибан то 0,8 mol/m³ мерасад.

High vapor concentration дар inlet side ҳангоми ворид шудан ба reactive bed аз ҷониби salt consume мешавад. Кам шудани concentration ба тарафи outlet spatial trace-и hydration reaction-ро нишон медиҳад.

Табдили маҳаллии hydration

Conversion contours қимати X-ро дар диапазони 0–0,09 нишон медиҳанд. Баландтарин local conversion тақрибан 9 дарсад аст. Аз ин рӯ дар simulation time full hydration ба даст наомадааст.

• 4 g/kg: Reaction суст, танг ва пароканда аст.
• 9 g/kg: Reaction front равшантар шуда, дар тӯли bed дароз мешавад.
• 12 g/kg: Дар locations-и бештар local conversion-и баланд ба вуҷуд меояд, аммо X боз ҳам ба 1 наздик намешавад.

PDF мегӯяд, ки experimental total conversion дар 9 g/kg humidity тақрибан дар диапазони 1,7–8,7 дарсад аст ва numerical results дар ҳамон order of magnitude мемонанд. Ин conversion-и хеле паст нишон медиҳад, ки гарчанде stable heat power ба даст ояд, қисми калони chemical storage capacity-и solid дар як pass истифода намешавад.

Режими назоратшаванда бо vapor pressure чӣ маъно дорад?

Bed temperature ва useful heat нисбат ба solid velocity ба inlet humidity хеле қавитар вокуниш нишон доданд. Муҳаққиқон ин ҳолатро “vapor-pressure-controlled regime” меноманд.

Дар ин regime омили асосии reaction rate, пеш аз residence time-и salt дар bed, фарқ байни actual water-vapor pressure-и air ва equilibrium vapor pressure мебошад. Ҳатто агар solid муддати дарозтар дар bed монад, агар water-vapor driving force паст бошад, reaction суст мемонад.

Бо вуҷуди ин, барои full conversion танҳо high humidity кофӣ нест. Барои расидани reaction region ба тамоми solid volume вақти кофӣ, appropriate gas distribution ва low bypass лозиманд.

Таъсири bypass factor

Bypass factor қисми inlet air-ро ифода мекунад, ки метавонад бе effective contact бо reactive salt аз reactor гузарад. Дар air flow rate-и 110 m³/соат, inlet temperature-и 22°C, humidity-и 9 g/kg ва solid velocity-и 0,5 cm/соат:

Бо зиёд шудани bypass ҳавои камтар бо намак reaction мекунад, water-vapor consumption ва reaction heat кам мешаванд. Ин натиҷа нишон медиҳад, ки на танҳо chemical properties-и salt, балки inlet distributor, wall geometry ва bed homogeneity низ барои reactor performance муҳиманд.

Таъсири thermal conductivity-и намак

Вақте thermal conductivity-и намак байни 0,7–8 W/(m·K) тағйир дода шуд, bed ва outlet temperature curves қариб болои ҳам афтоданд. Муҳаққиқон инро бо он шарҳ медиҳанд, ки дар thin-bed structure-и reactor heat conduction нисбат ба convective air transport ва reaction heat дуюмдараҷа мемонад.

Ин натиҷа барои geometry ва operating conditions-и омӯхташуда эътибор дорад. Дар thicker beds, lower air velocities ё designs-и дигар бо local hot spots таъсири salt conductivity метавонад калонтар бошад.

Таъсири permeability

Вақте permeability байни 10⁻⁸–10⁻¹⁰ m² тағйир дода шуд, таъсир ба outlet temperature маҳдуд буд. Дар permeability-и пасттар airflow мушкилтар мешавад, аммо зиёд шудани gas contact бо reactive solid ва кам шудани internal bypass метавонад useful heat-ро зиёд кунад.

Дар PDF дида мешавад, ки дар condition-и 10⁻¹² m² useful heat ба таври равшан меафзояд. Аммо permeability-и ин қадар паст метавонад fan power ва pressure drop-и хеле баландтар талаб кунад. Азбаски дар parametric section pressure drop ё fan energy ба objective function дохил нашудаанд, net system-level benefit-и ин натиҷа нишон дода нашудааст.

Таъсири kinetic coefficient

Вақте kinetic factor тақрибан дар диапазони 0,8×10⁴–2×10⁴ s⁻¹ зиёд карда шуд, outlet temperature ва useful heat баланд шуданд. Reaction-и тезтар:

• Ба consume шудани water vapor дар масофаи кӯтоҳтар,
• Васеъ шудани active reaction region,
• Озод шудани heat тезтар,
• Гарм шудани outlet air бештар

меоварад.

Аммо модел нашудааст, ки kinetic coefficient-ро бо кадом material modification ё particle structure дар амал тағйир додан мумкин аст. Натиҷа potential-и composites ё structured salts бо improved kinetics-ро нишон медиҳад.

Муҳимтарин натиҷаҳои муҳандисии таҳқиқот

• Reactor performance нисбат ба solid velocity ба inlet-air humidity ҳассостар аст.
• Баландтарин heat output-и омӯхташуда дар inlet humidity-и 12 g/kg ба даст омад.
• Humidity-и 9 g/kg, ки ҳамчун practical winter condition пешниҳод шудааст, trade-off байни heat output ва achievable humidity мебошад.
• Барои air flow rate, 110 m³/соат condition-е буд, ки contact time ва heat-carrying capacity-ро беҳтарин мувозинат кард.
• High solid velocity метавонад continuous power output-ро нигоҳ дорад, аммо single-pass conversion-ро маҳдуд кунад.
• Кам кардани bypass design strategy-и муҳимест, ки метавонад reactor performance-ро бе тағйир додани chemical kinetics беҳтар кунад.
• Reaction front дар air-inlet side ба вуҷуд меояд ва бо зиёд шудани humidity дар тӯли bed васеъ мешавад.
• Дар simulation time full hydration ба амал наомад; highest local conversion тақрибан дар 9 дарсад монд.

Маънои натиҷаҳо барои гармидиҳии бино

Ба тақрибан 33–35°C расидани outlet air метавонад барои applications мисли low-temperature heating systems ё баланд кардани source temperature-и heat pump муфид бошад. Аммо ин таҳқиқот бевосита building heating load, radiator system, fan power, heat-exchanger performance ё annual climate-profile simulation анҷом надодааст.

Қимати kW/m³-и гузоришшуда useful heat per reactor volume мебошад. Total kW-и дар бинои воқеӣ таъминшаванда ба reactive volume, operating duration, salt amount, fan power, heat-exchanger efficiency ва cycling strategy вобаста хоҳад буд.

Ҷанбаҳои қавии таҳқиқот

• Open thermochemical moving-bed reactor-ро сеченака модел мекунад.
• Solid motion-ро ба ҷойи simple constant source term бо ALE moving mesh намояндагӣ мекунад.
• Porous flow, reaction, humidity transport ва heat transfer-ро ҳамзамон ҳал мекунад.
• Reaction front, water vapor, temperature ва conversion-ро spatially нишон медиҳад.
• Model results-ро барои temperature, humidity, useful heat ва pressure drop бо experimental literature муқоиса мекунад.
• Humidity, air flow rate, solid velocity, bypass, permeability, conductivity ва kinetic factor-ро алоҳида меомӯзад.
• Фарқи байни continuous heat power ва solid conversion-ро намоён мекунад.
• Circulation ва short-circuit flows дар reactor geometry-ро бо three-dimensional streamlines нишон медиҳад.

Маҳдудиятҳои асосии таҳқиқот

• Таҳқиқот new experimental data ё physical prototype пешниҳод намекунад.
• Модел ба як basic reactor geometry аз literature асос ёфтааст.
• Long-term cyclic stability-и salt модел нашудааст.
• Agglomeration, deliquescence, particle breakage, wear ва mechanical compaction ба назар гирифта нашудаанд.
• Porosity асосан доимӣ қабул шудааст; volume change with reaction ба таври муфассал ҳал нашудааст.
• Feedback-и solid motion ба airflow сарфи назар шудааст.
• Particle size distribution ва bed settling гузориш нашудаанд.
• Fan power ва pumping energy ба optimization дохил нашудаанд.
• Чӣ қадар аз outlet heat дар external heat exchanger истифода мешавад, ҳисоб нашудааст.
• Climate ва seasonal humidity variations ҳамчун time series модел нашудаанд.
• Deliquescence ё flowability loss дар high humidity омӯхта нашудааст.
• Conversion loss вобаста ба solid velocity бо direct conversion graph quantitatively нишон дода нашудааст.
• Simulation files, mesh, code ё COMSOL model мубодила нашудаанд.
• Solver settings, time step, convergence criteria ва computational hardware бо тафсил дода нашудаанд.
• Таҳқиқот аз peer review нагузаштааст.

Таҳқиқот чиро нишон медиҳад?

• Нишон медиҳад, ки moving-bed reactor дар асоси SrBr₂, ки бо humid air кор мекунад, метавонад дар three-dimensional ва continuous operation numerically модел шавад.
• Нишон медиҳад, ки inlet humidity ба temperature ва useful heat таъсири dominant дорад.
• Нишон медиҳад, ки medium air flow rate метавонад useful heat-и баландтар аз low ва high flow rates диҳад.
• Нишон медиҳад, ки reaction front дар inlet side ба вуҷуд меояд ва бо humidity васеъ мешавад.
• Нишон медиҳад, ки кам кардани bypass метавонад air-salt contact ва heat output-ро зиёд кунад.
• Нишон медиҳад, ки stable heat output ҳатто бе full salt conversion ба даст омада метавонад.
• Нишон медиҳад, ки ALE method метавонад moving solid bed-ро бо reaction ва transport processes муттаҳид кунад.

Таҳқиқот чиро исбот намекунад?

• Исбот намекунад, ки reactor дар бинои воқеӣ талаботи seasonal storage-ро қонеъ мекунад.
• Нишон намедиҳад, ки condition-и 12 g/kg humidity дар ҳамаи climates бехатар ё economic аст.
• Нишон намедиҳад, ки 9 g/kg барои ҳамаи winter climates universal optimum аст.
• Исбот намекунад, ки lower permeability ҳангоми ба ҳисоб гирифтани fan energy efficientтар аст.
• Нишон намедиҳад, ки strontium bromide дар давоми ҳазорҳо cycles structure-и худро нигоҳ медорад.
• Исбот намекунад, ки дар reactor agglomeration, deliquescence ё clogging ба вуҷуд намеояд.
• Нишон намедиҳад, ки model метавонад мустақиман ба different salt hydrates ё different geometries интиқол дода шавад.
• Нишон намедиҳад, ки useful-heat values net system efficiency ё storage-cycle efficiency мебошанд.
• Peer-reviewed scientific consensus ё commercial product validation пешниҳод намекунад.

Таҳқиқоти зарурии оянда

Барои validation-и модел бо independent experiments, дар ҳамон reactor бояд three-dimensional temperature sensors, local humidity probes, pressure measurements ва hydration analysis-и outlet solid якҷо истифода шаванд. Location-и reaction front метавонад, масалан, бо neutron imaging, X-ray methods ё solid samples аз heights гуногун тасдиқ карда шавад.

Long-term cycling studies бояд salt breakdown, volume change, porosity loss, agglomeration ва deliquescence behavior-ро муайян кунанд. Ин effects метавонанд ба модел ҳамчун cycle-dependent change-и permeability ва reaction kinetics илова шаванд.

Дар reactor optimization на танҳо useful heat, балки fan electricity consumption, pressure drop, total conversion, salt utilization, heat-exchanger performance ва capital cost низ бояд objective бошанд.

Барои кам кардани bypass air distributors, perforated plates, structured reactive materials, composite salts ё multi-zone bed designs омӯхта шаванд. Бо active learning ё surrogate models three-dimensional design space метавонад бо computational cost-и пасттар скан шавад.

Усул ва натиҷаҳои таҳқиқот

Хулосаи усули техникӣ

Натиҷаҳои асосии параметрӣ

Муқоисаи шароити намӣ

Trade-off-и дебити ҳаво

Фарқҳои нисбие, ки аз validation data ҳисоб шудаанд

Ин percentages аз experimental ва simulation values-и PDF ҳисоб шудаанд. Гарчанде модел general trends-ро мегирад, accuracy дар ҳамаи conditions дар диапазони 1–4 дарсад нест.

Қарорҳое, ки барои reactor design гирифта мешаванд

1. Operating humidity primary control variable барои reactor temperature ва power control мебошад.
2. Air flow rate набояд танҳо maximize шавад; бояд нисбат ба gas-solid contact time optimize шавад.
3. Ҳангоми танзими solid velocity мувофиқи power demand, single-pass conversion низ бояд ҷудогона monitor шавад.
4. Air-distribution geometry бояд барои кам кардани bypass ва pressure drop беҳтар карда шавад.
5. Thermal gain-и аз low permeability гирифташуда бояд якҷо бо fan electricity consumption арзёбӣ шавад.
6. Барои full energy-storage efficiency на танҳо outlet heat, балки hydration degree-и outlet solid низ бояд чен шавад.

Ёддошт оид ба манбаъ ва усул

Номи аслии таҳқиқот: 3D Multiphysics Modeling of a Moving-Bed Thermochemical Reactor: Performance Analysis and Optimization under Humid Air Conditions

Муаллифон ва тартиб: Sanaz Akbarzadeh; Abhishek K. Singh; Haoshui Yu; Ahmad Arabkoohsar.

Маълумоти саҳми баробар: Дар PDF equal contribution ё co-first authorship declaration нест.

Муаллифи масъул: Ahmad Arabkoohsar.

Почтаи электронии муаллифи масъул: ahmar@dtu.dk

Муассисаҳо:

1. Department of Civil and Mechanical Engineering, Technical University of Denmark, Kgs. Lyngby, Denmark.
2. Department of Thermal and Fluid Engineering, University of Twente, the Netherlands.
3. Department of Energy, Aalborg University, DK-9220 Aalborg East, Denmark.

Мувофиқати муаллиф–муассиса: Sanaz Akbarzadeh ва Ahmad Arabkoohsar Technical University of Denmark; Abhishek K. Singh University of Twente; Haoshui Yu Aalborg University.

Маълумоти ORCID: Дар PDF ORCID numbers-и муаллифон дода нашудаанд.

DOI: 10.2139/ssrn.6945146

Навъи манбаъ: Preprint-и three-dimensional numerical modeling ва parametric reactor-performance analysis.

Ҳолати peer review: Таҳқиқот аз peer review нагузаштааст.

Платформаи нашр: SSRN.

Санаи пешниҳоди SSRN: 15 июни 2026.

Шумораи саҳифаҳо: 46. Main text ва references дар 24 саҳифаи аввал, enlarged figures дар саҳифаҳои 25–46 ҷой доранд.

Маҷалла: Номи peer-reviewed journal ё final accepted publication information вуҷуд надорад.

Ношири аслии маҷалла: Барои ин version тасдиқ нашудааст. SSRN як open-access early-research ва preprint platform аст, ки Elsevier идора мекунад; таҳқиқот набояд ҳамчун peer-reviewed Elsevier journal article пешниҳод шавад.

Пайвандҳои расмӣ:Саҳифаи сабти SSRN ва Пайванди DOI.

Маблағгузорӣ: Таҳқиқот аз ҷониби European Union дастгирӣ шудааст. Дар PDF гуфта шудааст, ки он қисми Marie Skłodowska-Curie Action / Horizon Europe project бо project number 101147098 ва унвони “A Game-Changing Solar-Thermochemical Heat Storage and Intensification Technology for the Process Heat Sector” мебошад.

Ихтилофи манфиатҳо: Дар PDF declaration-и ҷудогонаи conflict of interest вуҷуд надорад.

Саҳми муаллифон: Дар PDF CRediT ё declaration-и муфассали монанди он вуҷуд надорад.

Дастрасии додаҳо ва код: COMSOL model file, geometric model, mesh file, solver settings ё analysis code мубодила нашудаанд. Experimental data-и дар validation истифодашуда аз previously published sources ва inlet-humidity data-и аз graph рақамисозишуда гирифта шудаанд.

Эъломияи истифодаи AI: Муаллифон гуфтаанд, ки барои беҳтар кардани grammar ва spelling-и матн Microsoft Copilot истифода шудааст; баъдан content-ро review ва edit карда, масъулияти пурраро бар дӯш гирифтаанд.

Ин мақолаи туркӣ бо баррасии пурраи матн, formula-ҳо, tables, reactor schematics, mesh visual, temperature-humidity-conversion contours, pressure distribution, streamlines ва parametric graphs-и PDF-и 46-саҳифагии боршуда омода шудааст. Ягон иддаои experimental success, economic feasibility, building-energy saving, long-term material life ё commercial applicability, ки дар PDF нест, илова нашудааст. External sources танҳо барои bibliographic verification-и title, authors, DOI, corresponding author, SSRN record date ва platform status истифода шудаанд.

Номувофиқиҳои дохилӣ ва масъалаҳои reproducibility дар PDF

• Диапазони air flow rate: Гарчанде дар methods section гуфта шудааст, ки inlet flow rate 25–100 m³/соат аст, дар results conditions-и 110 ва 140 m³/соат истифода шудаанд. Actual parametric range-и омӯхташуда равшан 80–140 m³/соат аст.
• Initial hydration degree: Гарчанде salt ҳамчун “monohydrate” таъриф шудааст, ифодаи SrBr₂·XH₂O, X≈0 истифода шудааст. Барои monohydrate X=1 интизор меравад; X≈0 ҳолати anhydrous ё near-anhydrous-ро ифода мекунад.
• Hydrate subscripts: Дар effective-property equations s1 ва s0 мутаносибан hexahydrate ва monohydrate шарҳ дода шудаанд, аммо маълум нест endpoint values-и X ба кадом phase мувофиқат мекунанд.
• Validation source: Матн дар як ҷой experimental validation-ро ба Reference [15] вобаста мекунад. Reference [15] numerical study аст; experimental reactor data-и истифодашуда мувофиқи figures ва tables аз Reference [28] омадааст.
• Validation error range: Abstract deviation-и ± 1–4 дарсад гузориш мекунад. Relative differences, ки аз Table 4 values ҳисоб шудаанд, барои useful heat тақрибан то 8 дарсад ва барои bed temperature то 7,07 дарсад мерасанд.
• Specific-heat unit: Дар Table 3 specific heat-и salt ҳамчун 970 J/(g·K) навишта шудааст. Ин қимат аз ҷиҳати физикӣ ғайриоддӣ аст ва эҳтимолан бояд 970 J/(kg·K) бошад.
• Table numbering: Ҳам reactor-dimension table ва ҳам mesh-independence table “Table 2” номида шудаанд.
• Mesh-element terminology: Барои three-dimensional volume model ибораи “free triangular element” истифода шудааст. Visual ва COMSOL terminology tetrahedral volume mesh-ро нишон медиҳанд.
• Reference-и mesh deviations: Шарҳ дода нашудааст, ки percentages дар mesh-independence table нисбат ба кадом mesh solution ҳисоб шудаанд.
• Reaction direction: Дар COMSOL implementation чӣ гуна signs-и hydration equation ва reaction-rate expression маҳдуд шудаанд, кофӣ нишон дода нашудааст.
• Equation-и bypass factor: F_b таъриф шуда, parametrically тағйир дода мешавад; аммо дар reaction-rate ё transport equations-и нишон додашуда чӣ гуна татбиқ шудааст, равшан навишта нашудааст.
• Dehydration/hydration expression: Caption-и Figure 7 дар discharge stage оғозро иштибоҳан “before dehydration starts” меномад. Раванди омӯхташуда hydration аст.
• Volumetric-heat units: Дар right axes-и permeability ва kinetic-factor graphs қиматҳои тақрибан 1.000–4.000 ҳамчун “kW/m³” label шудаанд. Ҳамаи натиҷаҳои дигар дар диапазони 2–5,5 kW/m³ мебошанд. Эҳтимол ин axes бояд W/m³ бошанд ё scale factor дошта бошанд.
• Solid velocity ва conversion: Results section мегӯяд, ки high solid velocity conversion efficiency-ро кам мекунад; аммо direct conversion data ё energy-efficiency graph бар зидди solid velocity пешниҳод нашудааст.
• Practical humidity optimum: 9 g/kg ҳамчун practical winter condition пешниҳод шудааст. Азбаски humidity profiles-и climates гуногун, deliquescence limit ё climatic availability analysis иҷро нашудааст, ин қимат набояд universal optimum тафсир шавад.
• Missing computational details: Time step, nonlinear solver, convergence thresholds, stabilization methods, hardware ва total computation time нишон дода нашудаанд.

Ин номувофиқиҳо general results-еро, ки таъсири қавии inlet humidity ва air distribution ба reactor performance-ро нишон медиҳанд, аз байн намебаранд. Аммо барои independent reproduction-и model, арзёбии дурусти error bounds ва интиқоли design recommendations ба real system, ин нуқтаҳо бояд ислоҳ шаванд.