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

Форма тока задавалась функцией Хейдлера [4]. где п = 10, 7’= 10"5, т = 2Т0"7, // = 1. Функция сдвинута на 0,5т, т.к. вначале растет слишком медленно. Расчеты в данной статье для упрощения проводились с диэлектрической проницаемостью без учета ее частотной зависимости (которая может быть учтена при помощи соответствующих методов [2,5]). Удельное сопротивление грунта равно 1500 Омм, относительная диэлектрическая проницаемость — 2 0 (что близко к диэлектрической проницаемости на частоте в несколько МГц по данным из [1]). Время рассчитанного процесса — 100 мкс. Заданный ток и рассчитанное напряжение показаны на рисунке 2. Рис.2. Заданный ток и рассчитанное напряжение для первых 3 мкс После расчета напряжения и тока, рассчитывается их отношение в частотной области, откуда выражаются значения удельного сопротивления и диэлектрической проницаемости [ 1 ]. Результаты расчета для сетки с размером ячейки 0,05 м показаны на рисунке 3. После этого были проведены расчеты для сеток с размерами 0,1 и 0,2 м. Результаты расчета представлены на рисунках 4 и 5. Из результатов видно, что рассчитанное удельное сопротивление для сетки 0 , 2 м заметно отличается от удельного сопротивления для меньшего шага. Т.е., по-видимому, расстояние в две ячейки между электродами мало. Проведены также расчеты для граничных условий UPML (uniaxial perfectly matched layer, [2]). Из результатов на рисунке 6 видно, что применение данных граничных условий приводит к недопустимо большой ошибке на низких частотах для диэлектрической проницаемости. Вероятно, увеличение объема уменьшило бы ошибку, но это привело бы к заметному увеличению времени расчета. Поэтому граничные условия CPML более предпочтительны. 58