Многоуровневые вычислительные схемы для решения трехмерных векторных уравнений Гельмгольца в неоднородных областях тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Нечаев, Олег Валентинович

• Нечаев, Олег Валентинович
• кандидат физико-математических наук
• 2005

Оглавление диссертации кандидат физико-математических наук Нечаев, Олег Валентинович

ВВЕДЕНИЕ

ГЛАВА 1. МЕТОДЫ МОДЕЛИРОВАНИЯ ЭЛЕКТРОМАГНИТНЫХ

ПОЛЕЙ

1.1. Математические модели.

1.2. Современные методы решения задач электромагнетизма. Векторный метод конечных элементов.

1.3. Многосеточные, многоуровневые методы и методы декомпозиции

ГЛАВА 2. ВЕКТОРНЫЕ ВАРИАЦИОННЫЕ ФОРМУЛИРОВКИ И ИХ

ДИСКРЕТНЫЕ АНАЛОГИ

2.1. Математическая модель.

2.2. Векторная вариационная постановка.

2.3. Дискретные подпространства.

2.4. Дискретные аналоги вариационных задач.

ГЛАВА 3. МНОГОУРОВНЕВЫЕ МЕТОДЫ РЕШЕНИЯ СЛАУ

3.1. Двухуровневый итерационный решатель.

3.2. Алгоритм решения СЛАУ, использующий ядро rot оператора . 59 ф 3.3. Многосеточный алгоритм

3.4. Мультипликативный алгоритм.

3.5. Особенности построения и реализации матриц перехода

ГЛАВА 4. СТРУКТУРА ПРОГРАММНОГО КОМПЛЕКСА. ЧИСЛЕННЫЕ РЕЗУЛЬТАТЫ

4.1. Структура программного комплекса.

4.2. Верификация программного комплекса.

4.3. Тестирование многоуровневых алгоритмов.

4.4. Моделирование работы высокочастотного каротажного зонда в неоднородной среде.

Введение диссертации (часть автореферата) на тему «Многоуровневые вычислительные схемы для решения трехмерных векторных уравнений Гельмгольца в неоднородных областях»

Моделирование электромагнитных полей играет важную роль при разработке различных приборов и интерпретации данных физических экспериментов. Часто эти две задачи дополняют друг друга, например, в электроразведке необходимо разработать прибор, обладающий заданными свойствами, и методы интерпретации данных, полученных с его помощью. Методы моделирования должны, с одной стороны, обеспечивать возможность учета особенностей конструкции прибора (различающиеся на порядки геометрические фрагменты прибора), а с другой — сложные по строению вмещающие среды, обладающие разрывными физическими свойствами. Часто возникают ситуации, когда из-за наличия геометрических и физических неоднородностей модель невозможно свести к одномерному, двумерному или осесимметрично-му случаю.

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

Одним из недостатков векторного метода конечных элементов являются плохие спектральные свойства матриц систем линейных алгебраических уравнений, полученных после дискретизации непрерывной задачи. Это связано с наличием большого ядра у rot оператора, что приводит к отсутствию свойств эллиптичности у дифференциального оператора и неопределенности матрицы. Поэтому применение многосеточных и многоуровневых методов, не учитывающих этой особенности, не приносит ожидаемого эффекта.

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

Цель работы. Разработка и реализация вычислительных схем на базе векторного метода конечных элементов, которые позволили бы выполнять многовариантные расчеты векторного электромагнитного поля (гармониче-екая зависимость от времени) в трехмерных областях с резко контрастными по физическим свойствам материалами. Для решения дискретного аналога векторного уравнения Гельмгольца, построенного при помощи элементов Неделека 1-го и П-го типов, разработать специальные процедуры предобу-словливания.

Методы исследования. Методы вычислительной математики. Сравнительный анализ результатов моделирования и имеющегося аналитического Ш решения. Расчеты на последовательности сгущающихся сеток с последующим анализом сходимости к аналитическим решениям.

Защищаемые положения и научная новизна:

• На базе векторного метода конечных элементов разработан и реализован алгоритм моделирования гармонических по времени электромаг

Ф нитных полей в трехмерных неоднородных по физическим свойствам областях.

• Разработана технология учета ядра rot оператора, при использовании векторных базисных функций Неделека низкого порядка 1-го и Н-го типов на параллелепипеидальной и тетраэдральной сетках.

• С использованием построенного алгоритма учета ядра rot оператора разработан и реализован многосеточный алгоритм решения системы линейных алгебраических уравнений, полученных при дискретизации векторного уравнения Гельмгольца на параллелепипеидальной сетке. Численно показана его эффективность.

• С использованием алгоритма учета ядра разработан мультипликативный алгоритм решения системы линейных алгебраических уравнений, полученных при дискретизации векторного уравнения Гельмгольца на тетраэдральной сетке, с использованием элементов Неделека 2-го типа. Численно показана его эффективность.

Значимость работы. Предложены, реализованы и исследованы многоуровневые алгоритмы решения систем линейных алгебраических уравнений, полученных после дискретизации векторного уравнения Гельмгольца при помощи векторного метода конечных элементов в неоднородных трехмерных областях. Для дискретизации непрерывной задачи использовались векторные конечные элементы 1-го и П-го типов первого порядка на параллелепипеидальной и тетраэдральных сетках. На ряде тестов показана эффективность предложенных алгоритмов по сравнению со стабилизированным методом бисопряженных градиентов. При помощи разработанного комплекса программ было выполнено моделирование процесса высокочастотного каротажного зондирования для различных конструкций зонда и различных типов сред.

Апробация работы. Основные результаты работы были представлены и докладывались на:

• Международная конференция молодых ученых по математическому моделированию и информационным технологиям (Новосибирск 2002),

• Региональная научная конференция студентов, аспирантов и молодых ученых «Наука. Техника. Инновации» (Новосибирск 2002),

• Международная конференция «Математические методы в геофизике» (Новосибирск 2003),

• Международной конференции по вычислительной математике (Новосибирск 2004).

• Объединенном семинаре института вычислительной математики и математической геофизики и кафедры вычислительной математики Новосибирского государственного университета. (Новосибирск, 2005).

Публикации. По теме диссертационной работы опубликовано 5 работ [Ю]-[12], [23], [26].

Структура работы. Диссертация состоит из введения, четырех глав, заключения и списка использованной литературы (134 наименования). Работа изложена на 118 страницах, включая 26 иллюстраций и 12 таблиц.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Нечаев, Олег Валентинович

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

ЗАКЛЮЧЕНИЕ

В работе предложены и реализованы методы решения задач моделирования гармонических по времени электромагнитных полей в неоднородных по физическим и геометрическим свойствам трехмерных областях. Для векторного уравнения Гельмгольца, описывающего поведение гармонического по времени электрического поля, сформулирована векторная вариационная постановка, учитывающая скачки поля на границах материалов с различными физическими свойствами. Показано, что электрическое поле, удовлетворяющее данной вариационной постановке, также удовлетворяет закону сохранения зарядов в слабой форме. Построены дискретные вариационные постановки, использующие векторные конечные элементы первого порядка 1-го на параллелепипеидальных сетках и первого порядка 1-го и Н-го типов на тетраэдральных сетках. Базисные функции первого порядка П-го типа построены в виде иерархического базиса.

Предложен двухуровневый алгоритм решения систем линейных алгебраических уравнений, полученных при дискретизации векторного уравнения Гельмгольца на параллелепипеидальной и тетраэдральной сетках, учитывающий ядро rot оператора. На базе данного двухуровневого алгоритма (в качестве сглаживателя) разработан многосеточный алгоритм решения СЛАУ на параллелепипеидальных сетках, и мультипликативный алгоритм решения СЛАУ, являющихся аппроксимациями векторного уравнения Гельмгольца на тетраэдральных сетках.

Проведенное тестирование разработанного программного комплекса на задаче, имеющей аналитическое решение, показало, что использование векторных элементов Н-го типа на тетраэдральной сетке является более экономичным по сравнению с элементами 1-го типа.

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

Реализованный программный комплекс был использован для моделирования процесса высокочастотного каротажного зондирования при помощи зондов ВИКИЗ и MPR в среде, состоящей из скважины, вмещающей среды и пропластков. Для зонда ВИКИЗ в случае однородной среды было проведено сравнение результатов моделирования с данными, полученными аналитически. Программный комплекс предоставляет возможность задания конструктивных особенностей зонда: корпус зонда различных конфигураций, количество и расположение генераторных и приемных катушек, что позволяет применять комплекс не только для моделирования процесса зондирования в сложных геометрических средах, но и для оптимизации конструкции зонда.

Список литературы диссертационного исследования кандидат физико-математических наук Нечаев, Олег Валентинович, 2005 год

1. Гилева J1.B Каскадный многосеточный алгоритм в методе конечных элементов для трехмерной задачи Дирехле //Сиб. журн. вычисл. математики / РАН. Сиб. отд-ние — Новосибирск, 1998 — Т.1; №3.— с.217— 226.

2. Годунов С.К., Рябенький B.C. Разностные схемы,— М.: Наука, 1977.— 439с.

3. Давыдычева С.Н., Друскин B.J1. Прямая задача расчета электромагнитного поля на оси наклонной скважины в анизотропной слоистой среде //Труды Международной Конференции «Горногеологической службе России 300 лет»,- С.-Петербург.- 2000.- с. 19-31.

4. Дворецкий П.И., Ярмахов И.Г. Электромагнитные и гидродинамические методы при освоении нефтегазовых месторождений.— М.: Недра, 1998.- 318 с.

5. Ильин В.П. Методы конечных разностей и конечных объемов для эллиптических уравнений.— Новосибирск: Изд-во Ин-та математики, -2000.- 345с.

6. Лаврентьев М.М., Романов В.Г., Васильев В.Г. Многомерные обратные задачи для дифференциальных уравнений — Новосибирск: Наука, 1969.- 67 с.

7. Лаевский Ю.М. Метод конечных элементов (основы теории, задачи).— Новосибирск: Новосиб. гос. ун-т, 1999.— 166 с.

8. Лаевский Ю.М. О некоторых итогах развития современной вычислительной математики //Вычислительные технологии.- 2002.— T.7;No2.— с.74—83.

9. Нечаев О.В. Шурина Э.П. Многосеточный алгоритм решения векторным методом конечных элементов трехмерного уравнения Гельмголь-ца// Математическое моделирование. 2005. - том. 17; No 6. — с. 92-102.

10. Нечаев О. В., Шурина Э.П., Федору к М.П. Использование метода конечных элементов для численного решения квазистационарных уравнений Максвелла// Вычислительные технологии. — 2004. — Том 9; No5.— с.73—81.

11. Потаов А.П., Кнеллер JI.E. Численное решение задачи становления магнитного диполя в скважинах много колонной конструкции //НТВ «Каротажник»,— Тверь, 1999.— Вып. 52.— с.76—83.

12. Самарский А.А., Гулин В.А. Численные методы.- М.: Наука, 1989.— 432с.

13. Технология исследования нефтегазовых скважин на основе ВИКИЗ. Методическое руководство/ Ред. Эпов М.И., Антонов Ю.Н.- Новосибирск: НИЦ ОИГГМ СО РАН, Издательство СО РАН, 2000.- 121 с.

14. Тихонов А.Н., Самарский А.А. Уравнения математической физики.— М.: Наука, 1977,- 735 с.

15. Федоренко Р.П. Релаксационный метод решения разностных эллиптических уравнений// ЖВМ и МФ.- 1961,- Т.1; №5.- с.922-927.

16. Федоренко Р.П. О скорости сходимости одного итерационного процесса// ЖВМ и МФ,- 1964,- Т.4; №5.- с.559-564.

17. Шайдуров В. В. Многосеточные методы конечных элементов.— М.: • Наука 1989, 288с.

18. Шурина Э.П., Гельбер М.А. О векторном методе конечных элементов для решения задач электромагнетизма //Сибирский журнал вычислительной математики /РАН. Сиб. Отделение — Новосибирск, 2004.— T.7,Nol с.79—95.

19. Шурина Э.П., Гельбер А.В., Гельбер М.А., Эпов М.И. Математическое моделирование нестационарного электромагнитного поля дефектоскопа обсадных колонн// Вычислительные технологии,— 2002,— Т.7; No6.-С.114-129.

20. Шурина Э.П., Гельбер М.А. Применение векторного метода конечных элементов для моделирования гармонических электромагнитных полей в неоднородных областях// Вестник НГУ.- Новосибирск.— 2004.

21. Шурина Э.П., Соловейчик Ю.Г., Рояк М.Э. Решение трехмерных нелинейных магнитостатических задач с использованием двух потенциалов: Препринт 1070. ВЦ СО РАН, 1996,- 26 с.

22. Adaptive multilevel methods for edge element discretization of Maxwell's equations / Beck R., Deuflhard P., Hiptmair R., Wohlmuth B. ; Preprint SC-97 66 Conrad-Zuse-Centrum, Berlin, 1997.- 45p.

23. Albanese R., Rubinacci G. Analysis of three-dimensional electromagnetic fields using edge elements //J. Comput. Phys. 1993. - Vol.108. -P. 236-245.

24. A mixed face-edge finite element formulation for the 3D magnetostatic problems / Alotto P., Delefino F., Molfino P., Nervi M., Perugeia I. //IEEE Trans. Magn.- 1998.- Vol.34.- p.2445-2448.

25. A nonlinear inversion method for 3D electromagnetic imaging using ajoint fields / Dorn 0., Bertete-Aguirre H., Berryman J.G., Papanicolaou G.C. //Inverse Problems.- 1999.- Vol.15.- P.1523-1558.

26. Assous F., Ciarlet P., Segre J. Numerical solution to the timedependent Maxwell equations in two-dimensional singular domains: the singular complement method //J. Comput. Phys.- 2000.- Vol.161.— P.218-250.

27. Arnold D. N. , Falk R.S. and Winther R. Multigrid in H(div) and H(curl)// Numerische Mathematik.- 2000.- Vol.85.- p. 197-218.

28. Beilennhoff K., Heinrich W. Treatment of field singularities in the finite difference approximation //Int. Micriwave Symp. Dig. 1993. - Vol.2. -p.979-982.

29. Bernardo Cockburn, Fengyan Li, Chi Wang Shu Locally divergence-free discontinuos Galerkin methods for Maxwell equations //Journal of computational physics.- 2004.- Vol.194.- p.588-610.

30. Bokil V.A., Bukas M.W. A 2D mixed finite element formulation of the unaxial perfectly matching layer //Ргерг. CP1677802, University of Houston.- 2002,- 18p.

31. Boner S., Fredette M., Clark В., Mills R., Williams R. Resistivity while drilling — images from the string //Slumberger olilfield review. 1996.-Vol. 8;Nol - p.4-19

32. Bossavit A. Computational electromagnetism.— Academic Press, San Diego, 1998.- 324p.

33. Bossavit A. Solving Maxwell equation in a closed cavity, and the question of spurious modes //IEEE Trans. Magn.- 1990.- Vol.26;No2 p.702-707.

34. Bourgeois В., Suignard K., Perrusson G. Electric and magnetic dipoles for geometric interpretation of three-component electromagnetic data in geophysics //Inverse problems.— 2000.— Vol.16.— P.1225—1261.

35. Braess D. Finite elements. Theory, fast solvers, and applications in solid mechanics.— Cambridge University Press.— 2001.— 370p.

36. Brenner S.C., Scott L.R. The mathematical theory of finite elements methods;— Springer-Verlag.— New York.— 2002.— 361p.

37. Brezzi F., Fortin M. Mixed and hybrid finite element methods.— Springer--Verlag, New York, 1991.- 351 p.

38. Bouillault F., Buffa A., Maday Y., Rapetti F. The mortar edge element method in three dimensions:application to magnetostatics //Publication No 1208,— Istitiuto di analisi numerica 2001.

39. Cenders Z.J. Vector finite elements for electromagnetic field calculations // IEEE Trans. Magn.- 1991,- Vol.27;No5.- P. 3958-3966.

40. Chen Z., Du Q., Zou J. Finite element methods with matching and nonmatching meshes for Maxwell equations with discontinuous coefficients //SIAM J. Numer. Anal.- 2000,- Vol.37;No5.- P.1542-1570.

41. Cioni J.-P., Fezoui L., Issautier D. High order upwind schemes for solving timedomain Maxwell equation //La Recherche Aerospatiale.- 1994.— No5 p.319-328.

42. Computation of optical modes inside axisimmetric open cavity resonators / Chinellato 0., Arbenz P., Streiff M., Geus R. ; Preprint: FGSC.- Zurich, 2003.- 20p.

43. Computation of 3D current driven skin effect using current vector potential / Biro O., Preis K., Reinhart W., Vrisk G., Richter F. //IEEE Trans. Magn.- 1993,- Vol.29.- P.1325-1332.

44. Cooray F., Costache G. An overview of the absorbing boundary conditions //J. Electromagn. Waves Appl 1991.- Vol.5;NolO - P.1041-1055.

45. Cpuliette D.L., Koch M. On the difficulties and remedies in enforcing the div=0 condition in the finite element analysis of thermal plumes with strongly temperaturedependent viscosity //Int. j. numer. methods fuids.-1994,- Vol.18.- p.189—214.

46. Deuflhard P. Cascadic conjugate gradient methods for eleptic partial differential equations. I. Algorithm and numerical results/ Tecnical Report SC 93-23- Berlin: Konrad-Zuse Zentrum.- 1993.

47. Deufhard P., Friesse Т., Schmidt F. A nonlinear multigrid eigenproblem solver for the complex Helmholtz equation // Preprint SC97 55,- Conrad-Zuse-Centrum, Berlin, 1997.- 26p.

48. Dietrich Braess, Wolfgang Dahmen A cascadic multigrid algorithm for the Stokes equations//Numer.Math.— 1999,- Vol.82.- p.179-191.

49. Douglas С.С., Douglas J. A unified convergence theory for abstract multigrid or multilevel algoritms, srial and parallel// SIAM J. Numer. Anal.- 1993,- Vol.30.- p.136-158.

50. DyczijEdlinger R., Peng G., Lee J.F. A fast vector potential method using tangentially continuous vector finite elements //IEEE Trans. Microwave Theory Tech.- 1998.- Vol.46.- P.863-868.

51. Eddy current computations using adaptive grids and edge elements / Liu Y.Q., Bondeson R., R. Bergstrom, Johnson C., Larson M.G., Samuelson K. //IEEE Trans. Magn.- 2002.- Vol.38;No2.- P.449-452.

52. Electromagnetics via the Taylor-Galerkin finite element method on unstructured grids / Ambrosiano J., Brandon S. Lohner R., Devore C.R. //J. Comput. Phys.- 1994,- Vol.110.- P.546—569.

53. Engquist В., Majda A. Absorbing boundary conditions for the numerical simulation of waves //Math. Comput 1977.-Vol.31.- P.629-651.

54. Fast simulation of 3D electromagnetic problems using potentials / Haber E., Ascher U.M., Aruliah D.A., Oldenburg D.W. //J. Comput. Phys.- 2000.-Vol.163.- p.150-171.

55. Garry Rodrigue and Daniel White A vector finite element time—domain method for solvig Maxwell's equations on unstructured hexahedral grids //SIAM J. Sci. Comput.- 2001.- Vol.23.- p.683-706.

56. Givoli D. Numerical methods for problems in infinite domains/Elsevier.-Amsterdam, 1992.- 293p.

57. Heise В., Kuhn M., Langer U. A mixed variational formulation for 3D linear and nonlinear magnetostatics in the space Ho(rot) Ho(div) //HEJ Manuscript No. ANM981030A.- 1998,- 16 p.

58. Hermeline F. Two coupled particlefinite volume methods using Delaunay-Voronoi meshes for the approximation of VlasovMaxwell equations //J. Comput. Phys.- 1993.- Vol.106.- p.1-18.

59. Hesthaven J.S., Warburton, Т. Nodal high—order methods on unstructured grids. I. time—domain solution of Maxwell's equations, Journal of computatioanl physisc.- 2002.- Vol.181.- p.186—221.

60. Hiptmair R. Finite elements in computational electromagnetism// Acta Numerica.- 2002,- p.237-339.

61. Hiptmair R. Multigrid methods for Maxwell's equations// SIAM J. Nymer. Anal.- 1998,- Vol.36;Nol p.204-225.

62. Hoppe R.H.W. Adaptive multigrid and domain decomposition methods in the computation of electromagnetic fields// Journal of Computational and Applied Mathematics 2004.— Vol.168.- p.245—254.

63. Hybrid time domain solvers for the Maxwell equations in 2D / Arbenius E., Andersson U., Edelvik F., Eriksson L., Ledfeld G. //Int. J. Numer. Meth. Engng.- 2002,- Vol.53.- P.2185-2199.

64. Igarashi H. On the property of the curl—curl matrix in finite element analysis with edge element //IEEE Transaction on magnetics.— 2001.— Vol.37;No5. p.3129-3132

65. Igarashi H., Honma T. On convergence of ICCG applied to finite element equation for quasystatic fields // IEEE Trans. Magn.- 2002 Vol. 38; No2-P. 565-568.

66. Ise K., Inoue К. , M. Koshiba Three—dimensional finite—element method with edge—elements for electromagnetic waveguide discontinuities// IEEE Trans.- MTT—39;No8.— p.1289-1295.

67. Jackson J. Classical electrodinamics New York, Wiley, 1962.- 839 p.

68. Jean Piquet and Xavier Vasseur Multigrig preconditioned Krylov subspace methods for three—dimensional numerical solutions of the incompressible Navier—Stokes equations// Numerical Algorithm.— 1998.— Vol.17.— p.l— 32.

69. Jiang B.N., Wu J., Povinelli L.A. The origin of spurious solutions in computational electromagnetics //J. Comput. Phys.- 1996.- Vol.125.-P. 104-123.

70. Johnson S.G., Joannopoulos J.D. Block-iterative frequencydomain methods for Maxwell's equations in a plainwave basis //Optics express.-2000,- Vol.8;No3 P. 173-190.

71. Kameari A. Calculation of transient 3D eddy current using edge-elements// IEEE Trans. Magn.- Vol.26;No2 p.466-469.

72. Kanayama H. Largescale magnetic field analysis //Annual Report ADV-991,- 1999.- 12p.

73. Kikuhi F. Mixed and penalty formulations for finite element analysis of an eigenvalue problem in electromagnetism //Comput. Meth. Applied. Mech. Engin.- 1987.- Vol.64. P.509-521.

74. Klaus Stuben A Review of Algebraic Multigrid// Institute for Algorithm and Scientific Computing, Gremany.— 1999.— GMD—Report 69.

75. Kunze M., Heinrich W. Efficient FD formulation for lossy waveguide analysis based on quasy-static field characteristics //IEEE Microwave and Guid. Wave Lett.- 1999,- Vol.9;Nol2 p.499-501.

76. Lacoste P. Solution of Maxwell equation in axisymmetric geometry by Fourier series decomposition and by use of H(rot) conforming finite element //Numer. Math.- 2000.- Vol.84. p.577-609.

77. Lakhtakia A., Weiglhofer W.S. Timeharmonic electromagnetic fields in source regions in a simple uniaxial bianisotropic medium //Int. J. Applied Electromagnetics in Materials.— 1994,— No5.— P. 101-108.

78. Ledger P.D. , Peraire J., Morgan K., Hassan O. and Weatherill N.P. Adaptive hp finite element comutations of the scattering width output of Maxwell's equations// Int. J. Numer. Meth. Fluids.— 2003 — Vol.43.— p. 953-978.

79. Lee J.F., Lee R., Cangellaris A. Time domain finite elements method //IEEE Ant. Prop.— 1997,— Vol.45.— P.430-441.

80. Lee R.L., Madsen N.K. A mixed finite element formulation for Maxwell's equations in the time domains //J. Comput. Phys.— 1990.— Vol.88.— p. 284-304.

81. Lee J.F., Sun D.K., Cendes Z.J. Full-wave analysis of dielectric waveguides using tangential vector finite element method //IEEE Trans. Microwave Theor. Tech.— 1991.— Vol.39;No8.— P.1262—1271.

82. Lynch D.R., Paulsen K.D. Origin of vector parasites in numerical Maxwell solutions //IEEE Trans. Microwave Theor. Tech.— 1991.— Vol.39;No3.— p.383—394.

83. Madsen N., Ziolkowski R. A 3 dimensional modified finite volume technique for Maxwell's equations //Electromagnetics.— 1990.— Vol.l0;Nol.— p.147-161.

84. Manges J.B. and Cendes Z.J. Tree—cotree decompositions for first-order complete tangential vector finite elements// International journal for numerical methods in engineerin.— 1997.— Vol.40.— p. 1667—1685.

85. Mahadevan K., Mittra R. Radar cross section computation of inhomogeneous scatteres using edge based finite element method in time and frequency domains //Radio Science.— 1993.— Vol.8;No6.— P.1181--1193.

86. Mark Ainsworth and Joe Coyle Hierarchic finite element bases on untstructured tetrahedral meshes// Int. J. Numer. Meth. Engng.— 2003,— Vol.58.- p.2103—2130.

87. Moulton David J. , Dendy Joel E. Jr. and James M. Hyman The Black Box Multigrid Numerical Homogenization Algorithm// Journal of computatioanal physics — 1998 Vol.142.— p. 80—108.

88. Miniowitz R., Webb J.P. Covariant—projection quadrilateral elements for the analysis of waveguides with sharp edges, IEEE Trans, on Microwave Theory and Technique, MTT-39 No3.- p.505.-1991.

89. Mifune Т., Iwashita T. and Shimasaki M. New algebraic multigrid preconditioning for iterative solvers in electromagnetic finite edge—element analyses// IEEE transaction on magnetics.— 2003.— Vol.39;No3.

90. Monk P. An analysis of Nedelec's method for spatial discretization of Maxwell's equations //J. Comput. and Applied Math.— 1993.— Vol.47.— P.101—121.

91. Munz C.D., Schneider R., Voss U. A finite-volume method for the Maxwell equations in the time domain //SIAM J. Sci. Comput.— 2000.— Vol.22;No2.— p.449—475.

92. Nedelec J. C. Mixed Finite Elements in R3 //Numerische Mathematik.— 1980.- Vol.35; No3.- p.315-341.

93. Nedelec J. C. A New Family of Mixed Finite Elements in R3 //Numerische Mathematik.- 1986.- Vol.50.- p.57-81.

94. Newman G.A., Alumbaugh D.L. Frequency-domain modelling of air-borne electromagnetic responses using staggered finite differences // Geophys. Prospecting.— 1995,— Vol.43.— p.1021—1042.

95. Newman G.A., Hoversten G.M. Solution strategies for two and three-dimensional electromagnetic inverse problems //Inverse Problems.— 2000.- Vol.16.— P.1357—1366.

96. On a finite—element method for solving the three-dimensional Maxwell equations / Assous F., Degond P., Heintze E., Raviart P.A., Segre J. //J. Comput. Phys.— 1993. — Vol.109.— p.222-237.

97. Perugia I. A mixed formulation for 3D magnetostatic problems: theoretical analysis and face-edge finite element approximation //Numer. Math.— 1999,— Vol.84.— P.305—326.

98. Perugia I., SchotzauD., Monk P. Stabilized interior penalty methods for the time-harmonic Maxwell equations //Сотр. Meth. Appl. Mech. Engrg.-- 2002,— Vol.191.— p.4675—4697.

99. Piperno S., Fezoui L. A centered discontinuous Galerkin finite volume scheme for the 3D heterogeneous Maxwell equations on unstructured meshes //Rapport de recherche, No4733.— 2003.

100. Polstyank Sergey V. о and Jin—Fa Lee Two Level Hierarchical FEM Method for Modeling Passive Microwave Devices// J. Comput. Phy.— 1998,- Vol.140.- p.400—420.

101. Quasi-static conductor loss calculations in transmission lines using a new conforming mapping technique / Tuncer E., Lee B.T., Islam M.S., Neikirk D.P. //IEEE Trans. Microwave Theor. Tech.— 1994.— Vol.42.— P. 1807-1815.

102. Randolph E. Bank Hierarchical bases an the finite element methods// Acta Numerica.- 1996.- p. 1-43.

103. Rudolf Beck and Ralf Hiptmair Multilevel solution of the time—harmonic Maxwell's equations based on edge elements// Int. J. Numer. Meth. Engng.- 1999.- Vol.45.- p.901-920.

104. Riley D.J., Turner C.D. Volmax A solidmodelbased, transient volumetric Maxwell solver using hybrid grids //IEEE Antennas Propagat.— 1997.— Vol.39;Nol.— p.20—33.

105. Rylander Т., Bondeson A. Application of Stable FEMFDTD Hybrid to Scattering Problems //IEEE Trans. Antennas Propagat.— 2002.— Vol.50;No2.— p.141—144.

106. Saad Y. Iterative Methods for Sparse Linear Systems.— PWS Publishing Company, 1996.

107. Scales J.A., Snieder A. The anatomy of inverse problems //Geophysics.— 2000,— Vol.65;No6.— P. 1708-1710.

108. Soln P. Scalar and vector—valued finite element of variable order— TICAM Report 02-36.

109. S. Reitzinger, J. Schoberl An algebraic multigrid method for finite element discretizations with edge elements// Numer. Linear Algebra Appl.— 2002.— Vol.9.- p.223—238.

110. Shaidurov V.V. Some estimates of the rate of convergence for the cacscadic conjugate—gradient method— Magdeburg.— 1994.— (Preprint/ Otto—von— Guericke—Universitat; 4).

111. Shurina E.P., Soloneneko O.P., Voitovich T.V. Technologies of finite volume-finite element method for the solution of convection-diffusion problems on unstructured grids //Вычислительные технологии.— 2002.— T.7;No3.— c.98—120.

112. Studen K. Algebraic multigrid (AMG): Experiences and comparisons//Appl. Math. Сотр.- 1983 — Vol.13.- p419-452.

113. Tiao Lu, Pingwen Zhang, Wei Cai Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions// Journal of computational physics.— 2004.— Vol. 200,— p.549—580.

114. The mortar edge element methods in three dimensions: application to magnetostatics / Bouillault F., Buffa A., Maday Y., Rapetti F. //SIAM J. Sci. Comput.— 2002,— Vol.24.— P.1303—1327.

115. Tony F.: Chan Domain decomposition algorithms// Acta numerica.— 1994,- p.61—143.

116. Wang Т., Hohman G.W. A finite difference time domain method for three-dimensional electromagnetic modelling //Geophysics.— 1993.— Vol.58.— p.797-809.

117. Webb J.P. Edge elements and what they can do for you // IEEE Trans, on Magn.- 1993.- Vol.29;No2.

118. Webb J.P. Hierarchical vector basis function of arbitrary order for triangular and tetrahedral elements// IEEE Trans. Ant. Prop.— 1999.— Vol.47.- p. 1244-1253.

119. Webb J.P., Miniowitz R. Analiysis of 3D microwave resonators using covariant—projection elements// IEEE Trans, on Microwave Theory and Technique.- 1991.- MTT-39; Noll.- p.1895-1899.

120. Wesseling P. An introduction to multigrid method. — Chichester—New York—Brisbane—Toronto—Singapore: John Wiley and S—s.— 1992.

121. White D. Numerical modeling of optical gradient traps using the vector finite element method //J. Comput. Phys.— 2000,— Vol.159.— P.13—37.

122. Winslow A.M. Numerical solution of the quasilinear Poisson equation in a nonuniform triangle mesh// J. Comput. Phys.— 1967.— Vol.2.— p.149--172.

123. Wojcik G., Mould J., West L.C. Time—domain finite element modeling of 3D integrated optical devices //Optical Society of America: technical digest series (integrated photonic research).— 1993.— Vol.10.— P.112-115.

124. Yioultsis T.V., Tsiboukis T.D. Vector finite element analysis of waveguide discontinuities involving anisotropic media //IEEE Trans. Magn.— 1995.— Vol.31 ;No3.— P.1550-1553.

125. Yuguo Li A finite—element algorithm for electromagnetic induction in two dimensional anisotropic conductivity structures// Geophys. J. Int.— 2002.— Vol.148.- p.389—401.