Труды КНЦ вып.8 (ЭНЕРГЕТИКА вып.4 1/2012(8)

E Г +0’5 = C E Г ° ’5 + C, ( я Г - Я Г + Я Г 'li, j+0,5,k +0,5 a 'li, j+0,5,k+0,5 b \ z \i,j+ 1,k+0,5 ZU,j,k +0,5 У|, - Я } I ,j+0,5,k y li,j+0,5,k+1/ n JJ in i —я z li—1,j+ 1,k+0,5 z li,j+1,k+0,5 n n 0,5,j+0,5,k +1 + C b V y \ i , j+0,5,k+1 _ Я Ai — 1,j+0,5,k+1 + “ 'li —0,5,/,k+1 " ' Я Г = Я Г + D - ^ 1 ......... + E |n+0,5 - E |n+0,5 'li —0,5,j+1,k+1 ■'li—0,5,j+1,k+1 ^ \ y li—0,5,j+1,k+1,5 y li—0,5,j+1,k+0,5 Я Г = я Г + D ( e Г+0,5 — E Г+0,5 + E Г‘+0,5 — E . y |i,j+0,5,k+1 y |i,j+0,5,k+1 b V z li+0,5,j+0,5,k+1 z li—0,5,j+0,5,k+1 ' li,j+0,5,k+0,5 ' li,j+0,5,5,k+1,5 я Г+1 = н Г z li,j+1,k+0,5 z li,j+1,k+0,5 E l " ' 05 = C E [ " + C ( я Г —Я Г + Я Г Я | y li—0,5,j+1,k+0,5 a y li—0,5,j+1,k+0,5 b ' 'li —0,5,j,k+1 'li —0,5,j,k E | n+0,5 = c e |n—0,5 + C ( я \" — Я \n + Я Г — Я \ I л > \ z li—0,5,j+0,5,k+1 a z li— j \ У li, y li ,j ' ,j 'li —0,5,j+1,k+1/ . ( 2 ) |n+1 tj in 'li , , , n J ( z li—0,5,j+0,5,k+1 z li—0,5,j+1,5,k+1/ n+0,5 in+0,5 \ n +0,5 n +0,5 n +0,5 n +0,5 ь E —E + E —E I b \'li,j+1,5,k+0,5 '\i,j+0,5,k+0,5 y li—0,5,j+1,k+0,5 y li+0,5,j+1,k+0,5/ В (2) верхний индекс указывает на момент времени в который происходит вычисление, а нижние - на положение вычисляемой величины в пространстве. Рис.1. Ячейка Yee Рис.2. Расчетная область Константы уравнений Максвелла Са,Сь и D описывают параметры среды в соответствующей точке пространства: с =[1 |[ 1+ст' dt 2 s 2 s (3) с - S If . (4) (5) Выбор длины ребра зависит от характерных размеров моделируемых объектов. В данном случае шаг расчета принят равным 0.25 м. Шаг расчета по времени At вычислялся из условия устойчивости: 47