Решение осесимметричных задач моделирования нестационарных электромагнитных полей

Министерство образования и науки РФ

Новосибирский Государственный Технический Университет

Кафедра ПМт

Лабораторная работа №3

по дисциплине: технологии математического моделирования физических процессов

Решение осесимметричных задач моделирования нестационарных электромагнитных полей

Факультет: ПМИ

Группа: ПМ-23

Студенты:   Чусовкова А.В.

Ефимова А.А.

Преподаватели:    Рояк М.Э.

Чернышев А.В.

Новосибирск 2006

Цель работы

изучение проблем аппроксимации параболической задачи.

Постановка задачи

Решить краевую задачу, описываемую уравнением:

,                                         

где:

 - электрический потенциал,

 - удельная электрическая проводимость материала.

Расчетная область для данной задачи:

Краевые условия:

Вид сетки для данной задачи:

Физические параметры задачи:

Выполнение работы

1. Решение задачи вычисления потенциала от одного источника:  

На верхней границе расчетной области были определены значения потенциала в некоторых точках, со сгущением вблизи источника. По координатам данных точек и значениям потенциала в них был построен интерполяционный линейный сплайн.

2. Проверка на оптимальность построенной сетки

Значения потенциала в некоторых точках на верхней границе области:

4. Построение электрического поля на поверхности Земли.

С помощью интерполяционного сплайна было построено распределение напряжения между регулярно расположенными на поверхности Земли парами точек M и N. Считаем, что измерение разности потенциалов осуществляется между точками M и N, отстоящими на 5 метров от равномерной сетки точек О. Расстояние между точками A и B взято равным 1000м, между точками M и N – 10 метров.

Разность потенциалов между точками M и N будет равна:

, где  и , где:

- расстояние между A и М,

- расстояние между B и М,

- расстояние между A и N,

- расстояние между B и N.

Построенное поле изображено на рисунке:

5.Определение глубинности.

Необходимо определить глубинность – влияние слоев разной проводимости, находящихся на разной глубине, на поле на поверхности Земли. Будем менятиь электрическую проводимость слоя, если при этом величина  окажется более 5%, считаем слой видимым в данной точке.

Ниже приведены графики  при проверке видимости всех слоев. Выступающая часть первого графика для слоя – видимые точки. На втором графике для каждого слоя изображен вид сверху.

Рис. 1 Слой 1

Рис. 2 Слой 2

Рис. 3 Слой 3

Рис. 4 Слой 4

Рис. 5 Слой 5

На следующих трех графиках изображен вид сверху для графика видимости слоев 6 - 8:

Заметно, что хорошо видны слои 1 - 4. Далее видимость значительно ухудшается. Распространение идет от источников в стороны.

Выводы

В рамках данной математической модели возможно найти полезные ископаемые по  поведению электрического поля от горизонтальной линии в горизонтально-слоистой среде до глубины приблизительно 4000 метров.