Главная | Обратная связь | Поможем написать вашу работу!
МегаЛекции

Метод переменных направлений

Рассмотрим двумерное уравнение теплопроводности:

 

dU = LU + f (x,t), x Î G02, t Î [0,t0]

 dt

U | г = m (x,t)                                                             (1)

U (x,0) = U0 (x)

     


LU = LU = (L 1 +L2) U, где L a U = d2U   , a =1,2

dx2

Область G0 a =G0 = {0<= x a <= l a, a =1,2} -прямоугольник со сторонами l1 и l2, Г - граница G0 = G0 + Г.

В G0 построили равномерную по xa сетку v h с шагами h1 = l1 /N1, h2 = l2 /N2. Пусть n h - граница сеточной области w h, содержащая все узлы на сторонах прямоугольника, кроме его вершин, v h = w h + n h.

Оператор L a заменим разностным оператором L a:

 

L a y = yx a x a, L = L 1 + L 2

В случае одномерного уравнения теплопроводности неявная схема на каждом слое приводит к разностной краевой задаче вида:

 

Aiyi-1 - Ciyi + Biyi+1 = -F, i=1,...,N-1

y0= m 1                                                                                                                            (2)

yn= m N

Ai > 0, Bi > 0, Ci > Ai + Bi

которая решается методом прогонки.

Рассмотрим теперь нашу двимерную задачу в прямоугольнике. Сетку v h можно представить как совокупность узлов, расположенных на строках i2=0,1,2,...,N2, или как совокупность узлов расположенных на столбцах i1=1,2,...,N1. Всего имеется N1+1 столбцов и N2+1 строк. Число узлов в каждой строке равно N1+1, а в каждом столбце N2+1 - узлов.

Если на каждой строке (или столбце) решать задачу вида (2) методом прогонки при фиксированом i2 (или i1), то для отыскания решения на всех строках (или столбцах), т.е. во всех узлах сетки, понадобится О(N1N2) арифметических действий. Основная идея большинства экономичных методов и состоит в сведении перехода со слоя на слой к последовательному решению одномерных задач вида (2) вдоль строк и вдоль столбцов.

Наряду  с основными значениями искомой  сеточной  функции   y(x,t),  т.е.  с   y = yn и y` = yn+1 вводится промежуточное значение y = yn+½, которое можно формально рассматривать как значение при t = tn+½ = t n+½. Переход от слоя n на слой n+1 совершается в два этапа с шагами 0.5 t.

 

yn+½ - yn = L 1 yn+½ + L 2 yn + j n                                                         (3)

T

yn+1 - yn+½  = L 1 yn+½ + L 2 yn+1 + j n                                (4)

T

 

Эти уравнения пишутся во всех внутренних узлах x = xi сетки v h и для всех t=th > 0.

Первая схема неявная по направлению х1 и явная по х2, вторая схема явная по х1 и неявная по х2. К уравнениям (3),(4) надо добавить начальные условия:

 


y(x,0) = U 0(x),   x Î v h                                                       (5)

и разностно краевые условия, например, в виде:

 

yn+1 = m n+1 при i1=0, i2=N2                            (6)

yn+½  = m при  i1=0, i2=N1                           (7)

где m = 1 (m n+1 + m n) -  t  L2(m n+1 - m n)                                                              (8)

 2                  4

Т.о., разностная краевая задача (3)-(8) соответствует задаче (1). Остановимся на методе решения этой задачи. Пререпишем (3) и (4) в виде:

 


2 y - L 1 y = F, F = 2 y + L 2 y + j

t                         t                                    (9)

             


2 y` - L 2 y` = F’, F = 2 y + L 1 y + j

t                         t

Введём обозначения:

 

xi = (i1h1, i2h2)

F = Fi1,i2

y = yi1,i2

при этом, если в уравнении один из индексов фиксирован, то его не пишем. Тогда (9) можно записать в виде (2), т.е.:

         
   


1   yi1-1 - 2 1   + 1 yi1 + 1  yi1+1 = - Fi1

h21            h21 t          h21

i1 = 1,...,N1-1                                                                (10)

y = m при i1 = 0,N1

 


1   y`i2-1 - 2 1   + 1 y`i2 + 1  y`i2+1 = - Fi2

h22             h22 t            h22

i2 = 1,...,N2-1                                                                (11)

y` = m ` при i2 = 0,N2

 

Пусть задано у=у n. Тогда вычисляем   ò F, затем методом прогонки вдоль строк i2=1,...,N2-1 решаем задачу (10) и определим y’ во всех узлах сетки w h, после чего вычисляем F и решаем задачу (11) вдоль столбцов i1=1,...,N1-1, определяя y`=yn+1. При переходе от слоя n+1 к слою n+2 процедура повторяется, т.е. происходит всё время чередование направлений.

 

Построение разностных схем

 

Для каждой области МДП - структуры построим консервативную разностную схему, учитывая при этом заданные условия.

Разобьём данную МДП - структуру на несколько областей следующим образом:

 

 

L                                M       N

 y

 

K 0

 

K 1

 

 

x

 


I: jk0,y = Un

t. j k+½ i-1,y + 1 +      t   +      t    . j k+½ ij -      t. j k+½ i+1y = Y ij

2h*ihi                         2h*ihi+1    2h*i2hi                2h*ihi+1

j k1,y = Un

где Y ij = j k ij +  t (L y j k ij + f k ij)

2

L y = 1 j kij+1 - j kij - j kij - j kij-1

   r*j     rj+1              rj

     

 


II: j ij =U3

   t   . j k+½ i-1,j + 1 +     t         +     t       . j k+½ ij -       t j k+½ i+1,j =

2h*ihi                              2h*ihi+1  2h*ihi               2h*ihi+1

= j k ij +  t L y j k ij

           2                   ,           0 < i < k0-1 L < j < M

 eok  . j k+½ i-1,j + -  enn  -  eok  . j k+½ ij +   en    . j k+½ i+1,j = Y * ij, i=k0

h*i-1                        h*hi h*hi-1                h*ihi

t   . j k+½ i-1,j + 1 +     t       +    t  . j k+½ ij -     t   . j k+½ i+1,j =

 2h*ihi                              2h*ihi  2h*ihi              2h*ihi+1

=   j k ij +  t L y j k ij - f kij                         ,k0+1 < i < k1

         2

j k 1,j = Un

 

...

 

III: j k0,j =Uc

      t   . j k+½ i-1,j + 1 +     t         +     t       . j k+½ ij -       t j k+½ i+1,j =

2h*ihi                            2h*ihi+1  2h*ihi               2h*ihi+1

= j k ij +  t L y (j k ij - f kij),                      M+1 < j < N

 2

j k1,j = Un

Разностные схемы (I) -(III) решаются методом прогонки в направлении оси OX.

 

y

 


K0

 

K1

       x

 

()

Разностные схемы (IV) -(VI) также решаются методом прогонки в направлении оси OY.

ЛИТЕРАТУРА

1. Годунов С.К.,Рыбинский В.С.: ”Разностные схемы”

2. Кобболд Р.: “Теория и приминение транзисторов”

3. Самарский А.М.: “Теория разностных схем”

4. Самарский А.М.,Николаев Е.С.: “Методы решения сеточных уравнений”

5. Самарский А.А.,Андреев В.Б.: “Разностные методы решения эллиптических уравнений”

6. Калиткин Н.Н.: ”Численные методы”

Поделиться:





Воспользуйтесь поиском по сайту:



©2015 - 2024 megalektsii.ru Все авторские права принадлежат авторам лекционных материалов. Обратная связь с нами...