знаходять значення температури у внутрішніх вузлах сітки.
Температуру в межах програми слід виражати в кельвінах. Так як коефіцієнти теплопровідності футерування обчислюються в точках, що лежать на гранях контрольних обсягів між вузлами сітки, то в розрахункові формули підставляють середньоарифметичні значення температури в сусідніх вузлових точках:
В
де а, b - числові коефіцієнти.
Згідно прийнятої тут нумерації точок на гранях контрольних обсягів, формули для розрахунку коефіцієнтів дискретних рівнянь будуть представлені наступним чином:
В
Коефіцієнт теплопровідності шару гарнісажу, утвореного на внутрішньої поверхні футеровки застиглим клінкерною розплавом, приймають рівним 1 Вт/(м • К).
У циклі зворотного прогонки обчислюють у вузлах сітки відносні різниці значень температури в поточній і попередньої ітераціях:
В
і вибирають з них максимальну різницю. У кінці загального ітераційного циклу роблять оцінку збіжності ітерацій, порівнюючи абсолютну величину із заданим малим числом. Якщо | | більше ніж, наприклад, 0,00001, то ітерації повторюються, якщо ж менше, то ітерації завершуються. Щоб повторити знову вибір максимальної відносної різниці температур, на початку кожної ітерації величину встановлюють рівною нулю.
В кінці програми передбачають виведення на екран і на друк вихідних даних і результатів розрахунку. При цьому температуру представляють в градусах Цельсія. p> Друга програма для розрахунку стаціонарного температурного поля в футеровке повинна враховувати теплообмін футерування з дифузійним факелом і повітрям у печі, згідно граничним умовам третього роду. З цією метою в попередню програму вносяться відповідні зміни.
Коефіцієнти прогонки Рj, Sj на внутрішній поверхні футеровки або гарнісажу визначають тепер, попередньо обчисливши коефіцієнти дискретного рівняння. При цьому використовуються значення ступеня чорноти і поглинальної здатності газоподібного середовища, наведеної ступеня чорноти і коефіцієнта конвективної тепловіддачі, знайдені при розрахунку параметрів дифузійного факела. Температуру на внутрішній поверхні футеровки або гарнісажу визначають зворотного прогоном.
Третю програму складають для розрахунку нестаціонарного теплообміну, маючи на увазі, що при обертанні печі температура внутрішньої поверхні футеровки змінюється. При контакті з технологічним матеріалом вона дорівнює температурі цього матеріалу, що відповідає граничним умовам першого роду, а при нагріванні дифузійним факелом залежить від умов радіаційної та конвективної тепловіддачі, згідно граничним умовам третього роду.
Для обліку нестаціонарності в програмі організують цикл по інтервалах часу, який є зовнішнім по відношенню до ітераційного циклу. Цикл по інтервалах часу виконується протягом декількох обертів печі, так щоб розрахунковий час прогріву футерування виявилося достатнім для стабілізації змін температурного поля в н...
