Труды КНЦ вып.32 (ГЕЛИОГЕОФИЗИКА вып. 6/2015(32))

где волновое число к = к0 (го) = royč^jJ^ j с е Размеры области моделирования и поглощающего слоя были І?1=4.3А,, Л2=4.8л, где X =2п1к - длина волны. Шаги пространственной сетки составляли А г= У N x = У 160 и А(3 = nl Np = n l 360 рад = 0.5 град., т. е. число шагов было 7ѴГ=768, 7Ѵр=360 и на поглощающий слой приходится 80 шагов по г. Шаг по времени т должен удовлетворять условию Куранта (сеточная скорость Vh больше скорости света в среде а]), которое в рассматриваемом случае удобно представить в виде: rmnVh = т іп ((і?0А р /т);(А г/т ))>а1= с /А/в 1ц1 . (4) Рассматривались два варианта значения 11, и шага по времени т = О/Л',,: і?оі=0.5)ц Ті= Ѳ/720,7Ѵѳі=720 и R02= 1 5 \ т2= Ѳ/360,ТѴѲ2=360. ИМ Рис.1. Зависимость для каждой компоненты поля рассчитанной по формуле (5) усредненной установившейся относительной ошибки (§р(г,р = 0)) от г в плоскости поперек диполя {|3=0}: слева а - при і?02= 1-5^, т2= Ѳ/360 и справа б - при і?м=0.5)ц іі= Ѳ/720 На рисунке 1 для присутствующих компонент поля Ег, Ер, Ва представлены оценки зависимости от расстояния г до источника усредненной установившейся относительной ошибки ^8^(г,Р = 0)^ численного решения /'/, = 0,//. j по сравнению с аналитическим решением І 'а Р = 0, //- j в плоскости поперек диполя {|3=0}, (где F = Е г,Е ^ ,В а , гр = рАг, тк =кт ), которая рассчитывалась по формуле: /f k.+k' ^ (5) f ку+К Л I f ку К (Ър {гр $ = 0))= £ |F,,(^,p = 0 , ^ ) - F a ( ^ ,p = 0 ,^ )| / £ |Fa ( ^ ,p = 0 ,^ )| I k=k, +1 к=кл+\ Усреднение (в установившемся режиме) проводилось по интервалу \Т\, / 2|. где 77і=700Ѳ и 7^=19000, т. е. £і=700 N e, А=1200 Ne. На рис. la представлены результаты для і?02= 1-5^, і 2= 0/360, а на рис. 16 при /^„і=0.5/.. і,= Ѳ/720. Видно, что размер R0 окрестности источника, в которой задается аналитическое решение, оказывает влияние на величину ошибки: для меньшего і?оі ошибка больше, чем для R(a, несмотря на меньший шаг по времени. 122