close

Вход

Забыли?

вход по аккаунту

?

О непрерывной зависимости решений от данных задачи для системы уравнений Максвелла в квазистационарном электрическом приближении.

код для вставкиСкачать
Математическое моделирование. Оптимальное управление
Вестник Нижегородского
им. Н.И.задачи
Лобачевского,
2011, уравнений
№ 5 (1), с. 169–173
О непрерывной
зависимостиуниверситета
решений от данных
для системы
Максвелла
169
УДК 517.9
О НЕПРЕРЫВНОЙ ЗАВИСИМОСТИ РЕШЕНИЙ ОТ ДАННЫХ ЗАДАЧИ
ДЛЯ СИСТЕМЫ УРАВНЕНИЙ МАКСВЕЛЛА В КВАЗИСТАЦИОНАРНОМ
ЭЛЕКТРИЧЕСКОМ ПРИБЛИЖЕНИИ
 2011 г.
А.А. Жидков
Нижегородский госуниверситет им. Н.И. Лобачевского
Artem.Zhidkov@gmail.com
Поступила в редакцию 04.04.2011
Рассматривается квазистационарная система уравнений Максвелла, находящая применение при
решении математических задач атмосферного электричества. Доказывается теорема об устойчивости
решения при малых изменениях исходных данных (коэффициентов задачи, внешних источников, начальных данных).
Ключевые слова: атмосферное электричество, корректность, устойчивость, уравнения в частных
производных, функциональное пространство.
Введение
Исследование электромагнитных процессов
в атмосфере Земли в настоящее время вызывает
высокий интерес. В первую очередь это связано
с развитием климатических и метеорологических моделей. Хорошо известно такое атмосферное явление, как глобальная электрическая
цепь [1–5] – поддерживаемая достаточно стабильная разность потенциалов между поверхностью земли и верхними слоями атмосферы.
Попытки математического и численного моделирования поведения электрического поля в
атмосфере Земли при известных возмущающих
токах описаны в работах [1, 2, 5–10].
В последнее время интерес представляют
обратные задачи об определении источников
электрических токов, локализованных в грозовых образованиях, по результатам измерений
электрического поля, создаваемого этими источниками.
Одним из необходимых условий разрешимости обратных задач является устойчивость решения прямой задачи при возмущении исходных данных.
В настоящей работе обсуждаются вопросы
непрерывной зависимости решения прямой задачи от внешних источников электрического
тока, проводимости и начальных условий. Данная работа является продолжением работ [9, 10].
Постановка задачи
Пусть t  (0, T ) , x  Ω  R 3 ; Ω – открытое
ограниченное подмножество в R 3 , диффео-
морфное шаровому слою, с границей  
   1  2 ( Γ1 , Γ2 – две компоненты связности границы Ω , каждая из которых диффеоморфна сфере в R 3 ). В соответствии с физическим смыслом задачи, Γ1 – поверхность земли,
Γ2 – верхняя граница атмосферы. Предполагается, что в каждой точке x  Ω определён

единичный вектор внешней нормали n (x ) .
При отсутствии высокоскоростных процессов, связанных с молниевыми разрядами, исследование электрических процессов в атмосфере Земли обычно проводится в рамках квазистационарного электрического приближения
для системы уравнений Максвелла:


4π
rot H ( x, t ) 
σ ( x ) E ( x, t ) 
c

(1)

4π ст
1 E ( x, t )

J ( x, t ) 
,
c
c t

rot E ( x, t )  0,
(2)


div H ( x, t )  0, div E ( x, t )  4πρ( x, t ). (3)

3
Здесь J ст  C [0, T ]; L2 (Ω )
– объёмная


плотность сторонних токов, σ  L (Ω ) – удельная проводимость атмосферы. Данные функции
считаются заданными при решении прямых задач. Также предполагается, что существуют
постоянные σ  , σ  такие, что
0  σ   σ( x)  σ  для любого x  Ω .
Замечание: Определения основных функциональных пространств, используемых в настоящей работе, могут быть найдены, например,
в [9–12].
170
А.А. Жидков
Учитывая, что проводимость земли существенно превышает проводимость приземных
слоёв атмосферы и что проводимость атмосферы возрастает с ростом высоты по экспоненциальному закону [3–5, 9], будем считать, что границы рассматриваемой области  являются
идеально проводящей средой, что соответствует
заданию на границе  нулевых условий для
тангенциальной компоненты напряжённости
электрического поля [12]

E τ ( x, t )  0 , x  Γ
(4)

(для произвольной вектор-функции u , опреде 
лённой на  , приняты обозначения u n  (u  n ) ,

 
 

u n  u n  n , u τ  u  u n , x  Γ , где n – единичный вектор внешней нормали).
Дополнительно предполагается, что функция

E ( x, t ) удовлетворяет начальному условию


E ( x ,0 )  E 0 ( x ) , x  Ω .
(5)
Равенство (2) означает, что векторное поле

E ( x, t ) является потенциальным и представимо
в виде

E ( x, t )   grad ( x, t ) .
Тогда уравнение (1) преобразуется к виду

4
rot H ( x, t )   ( x) grad ( x, t ) 
c
(6)
4  ст
1 

J ( x, t ) 
grad ( x, t ).
c
c t
Не нарушая общности рассуждений, можем
считать, что электрический потенциал удовлетворяет граничным условиям, аналогичным (4):
( x, t ) x  0 , ( x, t ) x  C (t ) .
(7)
1
2
Подчеркнём, что в данном случае потенциал
верхних слоёв атмосферы C (t ) считается неизвестной функцией от времени и (7) констатирует, что ( x, t ) x не зависит от x  Γ2 .
2
Также из начального условия (5) для электрического поля можем получить начальное
условие для электрического потенциала
φ( x,0)  φ 0 ( x ) .
(8)
В работе [10] показано, что задача (6)–(8)
допускает замкнутую постановку в терминах
скалярного электрического потенциала:
Найти функцию φ  C 1 ([0, T ]; V (Ω )) , удовлетворяющую интегральному тождеству
d
grad ( x, t )  grad ( x) dx 
dt 

 4 ( x) grad ( x, t )  grad  ( x ) dx 



 4 J ст ( x, t )  grad ( x) dx,



(9)
при всех ψ  V (Ω) , и начальному условию
φ t  0  φ 0  V (Ω ) ,
(10)
где V () – гильбертово пространство

V ()  u  H 1 ( ) : u x  0, u x  const
1
2

со скалярным произведением

(u  v) V ( Ω )  (grad u ( x)  grad v( x))dx .
Ω
Включение φ(, t )  V (Ω) обеспечивает выполнение граничных условий (7) при каждом
фиксированном t  (0, T ) .
В работе [10] доказана теорема о существовании и единственности решения задачи (9),
(10) в функциональном пространстве V (Ω) .
Непрерывная зависимость решения
задачи (9), (10) от параметров
Теорема. Пусть
φ
(1)
0

, σ (1) , J (ст1) ,

J (ст2 ) , ( φ (0i )  V (Ω) , σ (i )  L (Ω ) ,

  
(2)
0
,  (2) ,
J  C 0, T ;
ст
(i)
L2 () ) ,
3
i  1, 2 ) – два различных набора исходных данных в задаче (9), (10). Пусть
φ 1 , φ 2  C 1 ([0, T ]; V (Ω )) – решения, соответствующие каждому набору. Тогда справедлива
оценка
 2  1 C ([ 0,T ];V (  )) 


 C    (02 )   (01)
 T J (ст2 )  J (ст1)
 (11)
V ()
L (  )


 ,
  ( 2 )   (1)    (02 )
 T J (ст2 )

V
(

)
L (  )  

где C  0 зависит от коэффициента σ * , вида
области Ω и значения T .
Доказательство. Введём вспомогательные
обозначения
a (φ, ψ)  4 π σ( x)(grad φ( x, t )  grad ψ( x )) dx ,
1
3
2
3
2

Ω

f [ψ]  4π J ст ( x, t )  grad ψ( x) dx .

Ω


Из условий на функции σ(x ) и J ст ( x, t ) следует, что a(,) – билинейная, симметричная,
ограниченная, коэрцитивная форма, а f – линейный ограниченный функционал. Тогда, на
основании леммы Рисса об общем виде линейного непрерывного функционала [13], существуют элементы A[φ] и F (t ) , такие, что
a(φ, ψ)  ( A[φ]  ψ) V ( Ω ) , f [ψ]  ( F (t )  ψ)V ( Ω ) .
Тогда уравнение (9) эквивалентно следующему
 d 
 
 ( A[]  )V (  )  ( F (t )   )V (  ) .

 dt
V (  )
171
О непрерывной зависимости решений от данных задачи для системы уравнений Максвелла
Поскольку равенство справедливо для любой функции ψ  V (Ω) , следовательно, φ удовлетворяет следующему абстрактному операторному уравнению
d(t )
 A[(t )]  F (t ) .
(12)
dt
Оператор A : V (Ω)  V (Ω) определяется
формулой
A[φ]  v ,
где функция v  V (Ω) является решением интегрального тождества


(grad v  grad ψ) dx  σ (grad φ  grad ψ) dx
Ω
для любого ψ  V (Ω) .
Тогда норма оператора A по определению
запишется в виде
A[φ] V ( Ω )
(13)
A  sup
 σ L ( Ω) .

φ V (Ω )
Аналогично, оператор
F (t ) : L2 ( ) 
 V () определяется из соотношения

F (t )[u ]  v ,
где v  V (Ω) – решение интегрального тождества

(grad v  grad ψ)dx  (u  grad ψ)dx
3

Ω
Ω
для любого ψ  V (Ω) . Тогда


F (t )[u ] V ( Ω )  u
L2 ( Ω )3
.
Получаем, что оператор A : V (Ω)  V (Ω)
является самосопряжённым, ограниченным,
положительно определённым. Тогда решение
уравнения (12) с начальным условием (10) запишется в виде [14]:
t
φ(t )  e  At [φ 0 ]  e A( τ t ) [ F ( τ)]dτ ,

2
2
где e u 
(14)
u
n
 n!
1
0
C ([ 0 ,T ];V (  ))
V ()
t


e A ( t )  e A ( t )  F2 ( ) V (  ) d.
2
(19)
1
Пользуясь определением операторной экспоненты и свойствами операторов Ai , легко
доказать неравенство
e  A2t  e  A1t  e α  A2  A1 ,
где α  T  σ * .
Следовательно, на основании (13), неравенство (19) примет вид
φ2  ~
φ  e α φ (02 )  T F2  A2  A1 




 e α φ (02 )  T F2  σ 2  σ1 ,
где F2  F2
C ( 0,T ;V ( Ω ))
, остальные нормы берут-
ся в соответствующих функциональных пространствах.
Рассмотрим задачи (15), (17). Функция
~  φ удовлетворяет следующей задаче:
Φφ
1
dΦ
 A1 [Φ]  F (t ) , Φ t 0  Φ 0 ,
dt
где F (t )  F2 (t )  F1 (t ) , Φ 0  φ (02 )  φ (01) .
Решение данной задачи существует, единственно и записывается формулой (14):
t
0

2
2
~ представляются формулой (14). Тогда, исхоφ
дя из вида задач (16), (17), на основании линейности операторов Ai , получаем оценку
~
 
 e  A t  e  A t  (2)

0
Ω

Данная задача также будет иметь единственное решение.
Тогда справедлива следующая оценка
~  ~
φ 2  φ1  φ 2  φ
φ  φ1 .
(18)
Рассмотрим норму φ  ~
φ . Функции φ и
– операторная экспонента.
n0
Равенство (12) позволяет заключить, что решения φ 1 , φ 2 задачи (9), (10) удовлетворяют
абстрактным задачам
dφ 1
 A1 [φ 1 ]  F1 (t ) , φ 1 t  0  φ (01) ,
(15)
dt
dφ 2
 A2 [φ 2 ]  F2 (t ) , φ 2 t 0  φ (02 ) . (16)
dt
Здесь операторы Ai порождаются функциями σ ( i ) ( x) , а функции Fi (t ) – функциями

J (стi ) ( x, t ) .
Рассмотрим вспомогательную задачу
d~
φ
~  φ ( 2 ) . (17)
~ ]  F (t ) , φ
 A1 [φ
2
0
t 0
dt
Φ  e  A1t [Φ 0 ]  e A1 ( τ t ) [ F ( τ )]dτ .

0
Поскольку оператор A1 ограничен и положительно определён, справедлива оценка
*
*
Φ  eσ T  Φ 0  e σ T T  F .
Подставляя полученные результаты в неравенство (18), получаем


φ 2  φ1 C ([0,T ];V (Ω ))  eα φ(02)  φ (01)  T J (ст2)  J (ст1) 



 e α φ (02 )  T J (ст2 )  σ 2  σ 1 .


Аналогично оценим
d
d
d ~
(φ 2  φ 1 ) 
(φ 2  ~
φ) 
(φ  φ 1 ) .
dt
dt
dt
Из задач (16), (17) и соотношения (14) получаем
172
А.А. Жидков
d
~ ]  A [φ ] 
(φ 2  ~
φ)  A1 [φ
2
2
dt
t
 A1 e  A1t [φ (02 ) ] 
Ae
A1 ( τ t )
1
[ F2 ( τ )]dτ 
0
t
 A2 e  A2 t [φ (02 ) ] 
A e
2
A2 ( τ  t )
[ F2 ( τ)]dτ .
0
Аналогично предыдущим рассуждениям,
используя оценку
A2 e  A2t  A1e  A1t  Te α  A2  A1 ,
получаем
d
( 2  ~
)  Te   (02 )  T F2   2  1 .
dt
Таким же образом может быть получена
оценка
d ~
(  1 )   * e   (02 )   (01)  *Te  F2  F1 .
dt
Тогда
φ 2  φ 1 C 1 ( 0 ,T ;V ( Ω )) 


 max  φ 2  φ 1

C ( 0,T ;V ( Ω ))
*

,
d
(φ 2  φ 1 )
dt


.
C ( 0 ,T ;V ( Ω )) 


Обозначая C  e σ T max 1, σ * , T , получаем
оценку (11) из утверждения теоремы.
Доказанная теорема позволяет заключить,
что при достаточно малых изменениях параметров задачи (9), (10) решение также меняется
достаточно мало.
Формулировка задачи в операторной форме
(12) с ограниченным оператором A позволяет
построить решение задачи (9), (10) на ограниченном отрезке времени. В частности, такая
постановка позволяет решать обратную задачу
финального наблюдения.
Автор выражает благодарность доценту кафедры математической физики А.В. Калинину
за внимание к работе и высказанные замечания.
Работа выполнена при финансовой поддержке в
рамках Аналитической целевой ведомственной программы «Развитие научного потенциала высшей
школы» (2009–2011 годы ) Минобрнауки РФ (регистрационный номер 2.1.1/3927), Федеральной целевой
программы «Научные и научно-педагоги-ческие
кадры
инновационной
России»
на
2009–2013
годы (шифр проекта НК-13П-13), гранта РФФИ (код
проекта 09-01-97019-р_поволжье_а).
Список литературы
1. Hays P.B., Roble R.G. A Quasi-Static Model of
Global Atmospheric Electricity. 1. The Lower Atmosphere // J. of Geophysical Research. 1979. Vol. 84,
No A7. P. 3291–3305.
2. Roble R.G., Hays P.B. A Quasi-Static Model of
Global Atmospheric Electricity. 2. Electrical Coupling
between the Upper and Lower Atmophere // J. of Geophysical Research. 1979. Vol. 84, No A12. P. 7247–
7256.
3. Mareev E.A., Anisimov S.V. Global Electric Circuit as an Open Dissipative System // Proc. 12th Int.
Conf. on Atmospheric Electricity. 2003. P. 797–800.
4. Mareev E.A., Anisimov S.V. Lifetime of the
Energy in the Global Electric Circuit // Proc. 13th Int.
Conf. on Atmospheric Electricity. 2007. P. 17–21.
5. Давыденко С.С., Беспалов П.А. Проявление
широтного хода проводимости атмосферы в распределении электрических полей и токов глобальной
электрической цепи // Геомагнетизм и аэрономия.
2000. Т. 40, № 2. С. 71–77.
6. Морозов В.Н. Модель нестационарного электрического поля в нижней атмосфере // Геомагнетизм и аэрономия. 2005. Т. 45, № 2. С. 268–278.
7. Морозов В.Н. Распределение электрического
поля, создаваемого нестационарным током заряжения грозового облака в атмосфере с неоднородной
электрической проводимостью // Прикладная метеорология. 2006. Вып. 7 (555). С. 51–67.
8. Browning G.L., Tzur I., Roble R.G. A Global
Time-Dependent Model of Thunderstorm Electricity.
Part I: Mathematical Properties of the Physical and Numerical Models // J. of the Atmospheric Sciences. 1987.
Vol. 44, No 15. P. 2166–2177.
9. Жидков А.А., Калинин А.В. Корректность одной математической задачи атмосферного электричества // Вестник Нижегородского госуниверситета
им. Н.И. Лобачевского. 2009. № 4. С. 123–129.
10. Жидков А.А., Калинин А.В. Некоторые вопросы математического и численного моделирования
глобальной электрической цепи в атмосфере // Вестник Нижегородского госуниверситета им. Н.И. Лобачевского. 2009. № 6. С. 150–158.
11. Темам Р. Уравнения Навье–Стокса. Теория и
численный анализ. М.: Мир, 1981. 408 с.
12. Дюво Г., Лионс Ж.-Л. Неравенства в механике
и физике. М.: Наука, 1980. 384 с.
13. Иосида К. Функциональный анализ. М.: Мир,
1967. 624 с.
14. Крейн С.Г. Линейные дифференциальные
уравнения в банаховом пространстве. М.: Наука,
1967. 464 с.
О непрерывной зависимости решений от данных задачи для системы уравнений Максвелла
ON A CONTINUOUS DEPENDENCE OF PROBLEM SOLUTIONS
ON PROBLEM PARAMETERS FOR A SYSTEM OF MAXWELL EQUATIONS
IN THE QUASI-STATIONARY ELECTRIC APPROXIMATION
A.A. Zhidkov
A quasi-stationary system of Maxwell equations is considered, which is used in solving mathematical problems
of atmospheric electricity. A stability theorem for the solution at small variations of initial problem parameters (coefficients, external sources, initial data) is proved.
Keywords: atmospheric electricity, correctness, stability, partial differential equations, functional space.
173
1/--страниц
Пожаловаться на содержимое документа