🌱 ഹരിതഗൃഹ കാലാവസ്ഥാ നിയന്ത്രണം ലൈവ് മോഡൽ

അടിസ്ഥാനം — സംഖ്യകൾക്ക് പിന്നിലെ മോഡൽ

സിമുലേറ്റർ ഉപയോഗിക്കുന്ന ഓരോ ബന്ധവും അതിന്റെ ഉറവിടം സഹിതം. അനുമാനങ്ങളായി അടയാളപ്പെടുത്തിയ സ്ഥിരാങ്കങ്ങൾ ഉദാഹരണത്തിനുള്ള കാലിബ്രേഷനുകളാണ്.

വായുവിന്റെ താപനില ഊർജ്ജ ബാലൻസ് പിന്തുടരുന്നു: സൂര്യൻ അകത്തേക്ക്, ആവരണത്തിലൂടെയും വെന്റിലേഷൻ വായുവിനൊപ്പവും ചൂട് പുറത്തേക്ക്, വെള്ളം ബാഷ്പീകരിക്കാൻ ഉപയോഗിക്കുന്ന ചൂട്, കൂടാതെ ചൂടാക്കൽ.

C·dT/dt = τ·I·(1 − 0.6·0.75·screen) − λ·(E_crop + E_fog) + Q_heat + (U·Ac/Af + ρc_p·Q_v)·(T_out − T) + G_soil·(T_soil − T)[1]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

ആവരണത്തിലൂടെയുള്ള താപനഷ്ടം U × വിസ്തീർണ്ണം × താപനില വ്യത്യാസമാണ്; അടച്ച സ്ക്രീൻ U കുറയ്ക്കുന്നു, തെളിഞ്ഞ രാത്രി ആകാശം അത് കൂട്ടുന്നു.

U = 6.2 − (6.2 − 4.0)·screen W/m²K (single glass → with thermal blanket); ×(1 + 0.25·clear sky) at night[1]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

വെന്റിലേഷൻ എന്നത് ചോർച്ചയും തുറന്ന വെന്റുകളിലൂടെയുള്ള കാറ്റ് നയിക്കുന്ന പ്രവാഹവും ചേർന്നതാണ്.

Q_v = (N·H/3600 + vent·C_w·max(1, u_wind))·(1 − 0.4·screen), N = 0.75 h⁻¹[1][4]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

സാച്ചുറേറ്റഡ് നീരാവി മർദ്ദം താപനിലയനുസരിച്ച് കുത്തനെ കൂടുന്നു; വായു സാച്ചുറേഷനിൽ നിന്ന് എത്ര അകലെയാണെന്നാണ് VPD.

e_s(T) = 0.6108·exp(17.27·T / (T + 237.3)) kPa; VPD = e_s(T) − e_a; RH = e_a / e_s(T)[2]

ജലബാഷ്പ ബാലൻസ്: ട്രാൻസ്പിരേഷനും ഫോഗും ഈർപ്പം ചേർക്കുന്നു, വെന്റിലേഷനും തണുത്ത ഗ്ലാസിലെ ഘനീഭവിക്കലും അത് നീക്കം ചെയ്യുന്നു.

H·dχ/dt = E_crop + E_fog − Q_v·(χ − χ_out) − k_c·(Ac/Af)·(1 − 0.7·screen)·max(0, e_a − e_s(T_cover)), χ = e·M_w/(R·T)[6][2]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

ട്രാൻസ്പിരേഷന് ഒരു വികിരണ ഭാഗവും വരണ്ട വായു (VPD) നയിക്കുന്ന ഒരു ഭാഗവുമുണ്ട്.

E_crop = A·(1 − e^(−k·LAI))·I_in/λ + B·LAI·VPD (A 0.45, k 0.64, LAI 3, B 0.02 kg/m²·h·kPa by day)[3]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

ഫോഗ് ബാഷ്പീകരിക്കുന്നത്ര മാത്രമേ തണുപ്പിക്കൂ: ലിറ്ററിന് ഏകദേശം 2.45 MJ, ഈർപ്പമുള്ള വായുവിനെക്കാൾ വരണ്ട വായുവിൽ കൂടുതൽ.

E_fog = fog·0.4 kg/m²·h × min(1, VPD/0.8); each evaporated kg removes λ = 2.45 MJ of sensible heat; the rest wets leaves[1][2]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

CO₂ ബാലൻസ്: ഡോസിംഗ് അകത്തേക്ക്, വെന്റിലേഷനും വിളയുടെ ആഗിരണവും പുറത്തേക്ക്.

H·dC/dt = dose − Q_v·(C − C_out) − P_net; dosing ≈ 5.5 g/m²·h holds ≈ 1000 ppm at 1 air change per hour[1]

പ്രകാശസംശ്ലേഷണം പ്രകാശം, CO₂, താപനില, VPD എന്നിവയോട് പ്രതികരിക്കുന്നു; CO₂ പ്രതികരണം ഹരിതഗൃഹ ഫല വിളകൾക്കുള്ള പ്രസിദ്ധീകരിച്ച പ്രായോഗിക നിയമം പിന്തുടരുന്നു.

P = P_ref·(1 − e^(−I_in/300))·exp(150/400 − 150/C)·f(T: 1 at 17–27 °C)·f(VPD) (from the CO₂ rule dP/P = 1.5 %·(1000/C)² per 100 ppm)[5][7]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

വിളയുടെ സമ്മർദ്ദ കൗണ്ടറുകൾ: വളരെ ചൂടുള്ള, വളരെ തണുത്ത, അല്ലെങ്കിൽ ഗ്രേ മോൾഡിന് ഇടയാക്കുന്നത്ര ഈർപ്പമുള്ള മണിക്കൂറുകൾ.

heat hours: T > 35 °C · cold hours: T < 10 °C · disease-risk hours: RH ≥ 90 % or wet leaves · GDD = Σ max(0, T − 10)·Δt[1][7][8][10]

പുറത്തെ കാലാവസ്ഥ: ദിവസേനയുള്ള താപനില വക്രം, സ്ഥിരമായ മഞ്ഞുബിന്ദു, മേഘങ്ങൾ കുറച്ച സൂര്യപ്രകാശം, ക്രമരഹിതമായ മേഘ-കാറ്റ് ഏറ്റക്കുറച്ചിലുകൾ.

T_out = T̄ + A·cos(2π(h − 15)/24); e_out = e_s(T_dew); I = I_clear·(1 − 0.75·N^3.4); cloud N and wind fluctuate (normal, seeded)[9][2]അനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

മോഡൽ ഉപയോഗിക്കുന്ന മറ്റ് സ്ഥിരാങ്കങ്ങൾ.

height 4.5 m · cover/floor 1.2 · τ 0.8 · heat capacity 30 kJ/m²K · heating ≤ 180 W/m² · screen blocks 60 % of light (75 % of that as heat) and 40 % of air exchange, and hides 70 % of the cold glass (condensation) · vents 0.016 m³/s·m² per m/s wind · soil 2 W/m²K · fog ≤ 0.4 L/m²·h · CO₂ ≤ 8 g/m²·h · P_ref 4 g CO₂/m²·h · condensation 1.9·10⁻⁸ kg/m²·s·Paഅനുമാനം: ഹരിതഗൃഹത്തിന്റെ വലിപ്പം, താപ ശേഷി, വെന്റ്, സ്ക്രീൻ, ഫോഗ്, വിള ഗുണാങ്കങ്ങൾ എന്നിവ ആധുനിക ഗ്ലാസ്ഹൗസിനുള്ള ഉദാഹരണ മൂല്യങ്ങളാണ്; ഭൗതികശാസ്ത്രം നന്നായി കലർന്ന ഒരൊറ്റ വായു വ്യാപ്തമാണ്. യഥാർത്ഥ സൈറ്റുകൾ ഓരോ സ്ഥിരാങ്കവും കാലിബ്രേറ്റ് ചെയ്യുന്നു.

ക്രമരാഹിത്യം: സീഡ് ഉള്ള mulberry32 ജനറേറ്റർ; ഉപയോഗിക്കുന്ന വിതരണങ്ങൾ — യൂണിഫോം, എക്സ്പോണൻഷ്യൽ (ഇൻവേഴ്സ് CDF), നോർമൽ (Box–Muller), പോയിസൺ (Knuth). സീഡ് കാണിക്കുകയും പങ്കിടാൻ കഴിയുകയും ചെയ്യും.

ഉറവിടങ്ങൾ

1. ANSI/ASAE EP406.4 JAN2003 (R2008) — Heating, Ventilating and Cooling Greenhouses (eq. 1–4, Table 1 U-values, Table 2 infiltration: new glass or fiberglass 0.50–1.0 h⁻¹, §5.8.5 CO₂ injection, §7.1.3 temperature limits); since withdrawn, available from ASABE as a historical document — American Society of Agricultural and Biological Engineers, 2003
2. R. G. Allen, L. S. Pereira, D. Raes, M. Smith — Crop evapotranspiration, FAO Irrigation and Drainage Paper 56 (eq. 11 saturation vapour pressure, eq. 14 dew point, λ = 2.45 MJ/kg) — FAO, Rome, 1998
3. M. Baille, A. Baille, J. C. Laury — A simplified model for predicting evapotranspiration rate of nine ornamental species vs. climate factors and leaf area — Scientia Horticulturae 59:217–232, 1994
4. T. Boulard, A. Baille — Modelling of air exchange rate in a greenhouse equipped with continuous roof vents — Journal of Agricultural Engineering Research 61(1):37–48, 1995
5. E. M. Nederhoff — Effects of CO₂ concentration on photosynthesis, transpiration and production of greenhouse fruit vegetable crops (PhD thesis; CO₂ rule) — Wageningen Agricultural University, 1994
6. B. H. E. Vanthoor, C. Stanghellini, E. J. van Henten, P. H. B. de Visser — A methodology for model-based greenhouse design: Part 1, a greenhouse climate model for a broad range of designs and climates — Biosystems Engineering 110:363–377, 2011
7. R. R. Shamshiri et al. — Review of optimum temperature, humidity, and vapour pressure deficit for microclimate evaluation and control in greenhouse cultivation of tomato: a review — International Agrophysics 32:287–302, 2018
8. Identifying and Managing Three Major Disease Threats in Greenhouse Tomatoes: Leaf Mold, Gray Mold, and Bacterial Wilt (SPES-820) — Virginia Cooperative Extension
9. F. Kasten, G. Czeplak — Solar and terrestrial radiation dependent on the amount and type of cloud — Solar Energy 24:177–189, 1980
10. G. S. McMaster, W. W. Wilhelm — Growing degree-days: one equation, two interpretations — Agricultural and Forest Meteorology 87:291–300, 1997