Page 71 - dysertaciyahembara
P. 71
2.3.4. Вплив густини сітки розбиття та часового проміжку на
точність числового розв’язку. Розв’зувалась спочатку задача
теплопровідності, а потім дифузії водню у стальному зразку (сталь 22к) у
формі паралелепіпеда із розмірами 3х2х1 метра. Зразок був нагрітий до
початкової однорідно розподіленої температури 100 ⁰с, а далі
охолоджувався до 20 ⁰с. На рис. 2.5а зображено зміну температури у
центральній точці зразка з координатами (0,0,0).
а) b)
Рисунок 2.5. Часова залежність температури (а) та відносної
концентрації водню (б) у центральній точці зразка (суцільна лінія –
розв’язок згідно з (2.21); точки – розв’язок, отриманий МСЕ).
Аналогічно до температури досліджено розподіл концентрації водню
у цьому ж зразку під час дегазації водню, попередньо наводненого до
рівномірно розподіленої концентрації С0. На рис. 2.5б подано зміну
концентрації водню з часом у центральній точці тіла. Отримані числові
результати порівняно з аналітичними (2.21).
Оцінено вплив густоти сітки і величини часового кроку на точність
отриманих числових результатів. На рис. 2.6а показано залежність відносної
похибки від кількості скінченних елементів, а на рис. 2.6б – від величини
кроку по часу.
71