📦 কোল্ড-চেইন লজিস্টিক কেন্দ্র লাইভ মডেল

ভিত্তি — সংখ্যার পেছনের মডেল

সিমুলেটরের ব্যবহৃত প্রতিটি সম্পর্ক, তার উৎসসহ। অনুমান হিসেবে চিহ্নিত ধ্রুবকগুলো দৃষ্টান্তমূলক ক্যালিব্রেশন।

কোল্ড রুমের তাপ ভারসাম্য: দেয়াল ও খোলা দরজা তাপ আনে, রেফ্রিজারেশন তা বের করে।

C · dT/dt = UA·(T_out − T) + Q_door − Q_refrig[3][2]অনুমান: তাপীয় ভর, দেয়াল ও দরজার পরিবাহিতা, কম্প্রেসারের আকার ও সেবা-সময় দৃষ্টান্তমূলক মান।

দরজা দিয়ে অনুপ্রবেশ খোলা দরজার সংখ্যা ও তাপমাত্রার পার্থক্যের সঙ্গে বাড়ে; এয়ার কার্টেন তা কমায়।

Q_door = n_open · 0.6 · 1.2 kW/K · (T_out − T) · (1 − 0.7·curtain)[2]অনুমান: তাপীয় ভর, দেয়াল ও দরজার পরিবাহিতা, কম্প্রেসারের আকার ও সেবা-সময় দৃষ্টান্তমূলক মান।

কার্নো সীমার ভগ্নাংশ হিসেবে রেফ্রিজারেশন দক্ষতা।

COP = 0.45 · T_evap / (T_cond − T_evap) [K][4]

পণ্যের কেন্দ্রীয় তাপমাত্রা দুই ঘণ্টা পিছিয়ে বাতাসের সঙ্গে চলে।

T_product(t+Δt) = T_room + (T_product − T_room)·e^(−Δt/2h)[3]

বিচ্যুতি: ২–৮ °C সংরক্ষণ পরিসরের বাইরে ডিগ্রি-ঘণ্টা (বেশি গরম বা জমে যাওয়া)।

excursion = Σ max(0, T − 8, 2 − T) · Δt [°C·h][1]

মডেলের ব্যবহৃত অন্যান্য পরিচালন ধ্রুবক।

cold room 60 MJ/K, walls 0.6 kW/K · door open 60 % of unloading · air curtain −70 % infiltration, +3 kW fan per busy door · 2 × 40 kW compressors, −1 %/K above 32 °C · 4 kW per busy door (levellers, lights) · failed reefer τ = 3 h · scenario site settings: cold-room pre-alarm 7 °C, overtime budget 68 door-hours per 14 h day · booked trucks miss 5 % of slots · unloading Erlang-2অনুমান: তাপীয় ভর, দেয়াল ও দরজার পরিবাহিতা, কম্প্রেসারের আকার ও সেবা-সময় দৃষ্টান্তমূলক মান।

বুকিং ছাড়া ট্রাক এলোমেলোভাবে আসে, বুকিং করা ট্রাক তাদের নির্ধারিত স্লটে আসে (কয়েকটি স্লট মিস করে), এবং দরজাগুলো সমান্তরালে সেবা দেয়। অপেক্ষার সময় লিটলের সূত্র মানে: সারিকে গত এক ঘণ্টায় ঘণ্টাপ্রতি প্রকৃতপক্ষে খালাস হওয়া ট্রাকের সংখ্যা দিয়ে ভাগ, তাই আগমন হার বদলালে অপরিবর্তিত সারি লাফ দেয় না।

arrivals ~ Poisson(λΔt) + booked slots (5 % no-show), service ~ Erlang-2 (mean 45 min), c staffed doors; X = n_busy / 45 min (smoothed, τ = 1 h); W_q = L_q / X[6][5]

কর্মী-সময়: কর্মীযুক্ত প্রতিটি দরজার জন্য অর্থ দিতে হয়, সেখানে ট্রাক থাকুক বা না থাকুক।

crew time = Σ c_staffed · Δt [door-hours]

এলোমেলোতা: সিডযুক্ত mulberry32 জেনারেটর; ব্যবহৃত বণ্টন — সুষম, সূচকীয় (বিপরীত CDF), স্বাভাবিক (Box–Muller), পয়সঁ (Knuth)। সিড দেখানো হয় এবং শেয়ার করা যায়।

উৎস

1. PAHO — Cold chain: vaccines are stored and transported between +2 °C and +8 °C — Pan American Health Organization
2. ASHRAE Handbook — Refrigeration, ch. 24 Refrigerated-Facility Loads (door infiltration, Gosney–Olama model; air curtains) — ASHRAE, 2022
3. F. P. Incropera et al. — Fundamentals of Heat and Mass Transfer: lumped capacitance — Wiley, 2011
4. Y. A. Çengel, M. A. Boles — Thermodynamics: Carnot COP of refrigerators — McGraw-Hill, 2019
5. J. D. C. Little — A Proof for the Queuing Formula L = λW — Operations Research 9(3), 1961
6. D. Gross, C. M. Harris — Fundamentals of Queueing Theory (M/M/c) — Wiley, 2008