Будь умным!


У вас вопросы?
У нас ответы:) SamZan.net

Нанотехнологий и новых материалов Сафронов А

Работа добавлена на сайт samzan.net:

Поможем написать учебную работу

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

Предоплата всего

от 25%

Подписываем

договор

Выберите тип работы:

Скидка 25% при заказе до 7.11.2024

МИНИСТЕРСТВО ОБРАЗОВАНЯ И НАУКИ

Тольяттинский государственный университет

Физико-технический институт

Кафедра «Нанотехнологий и новых материалов»

Сафронов А.И.

Методические указания

по выполнению лабораторной работы № 5

РЕШЕНИЕ КРАЕВОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ПАРАБОЛИЧЕСКОГО ТИПА РАЗНОСТНЫМ МЕТОДОМ

Тольятти 2010

УДК

ББК

М

Рецензент: А.А. Викарчук – профессор, заведующий кафедрой «Нанотехнологий и новых материалов», д. ф.-м.н.

Методические указания к лабораторной работе № 5 «Решение краевой задачи для уравнения параболического типа разностным методом»,  Сафронов А.И., 2010 – 10с.

Лабораторная работа предназначена для ознакомления студентов с разностным методом решения краевой задачи для уравнения параболического типа c граничными условиями 2-го рода. 

УДК

ББК

© Тольяттинский государственный университет, 2010

© Авторы пособия

Цель работы: усвоить численный метод решения  краевой задачи для уравнения параболического типа

Приборы и принадлежности: миллиметровая бумага.

ТЕОРЕТИЧЕСКИЙ МАТЕРИАЛ

Рассмотрим смешанную задачу для однородного уравнения теплопроводности. Задача состоит в отыскании функции и (х, t), удовлетворяющей в областиD ={ (х, t)|0 < х < а,0 <t < Т}уравнению                                        (1) начальному условию  u(x,t) = f(x)                         (2) и краевым условиям второго рода

   .                                          (3)

К задаче (1)—(3) приводит, а частности, задача о распространении тепла в однородном стержне длины а, на левом конце которого задаётся тепловой поток в стержнь (q), а правый конец теплоизолирован.

В данной лабораторной работе рассматриваются краевые условия второго рода. Замена переменных t kt приводит уравнение (1) к виду

поэтому в дальнейшем будем считать k =1.

Построим в области D = { (х, t) )|0 < х < а, 0 <t < Т} равномерную прямоугольную сетку а шагом h в направлении х и шагом τ — в направлении t (рис. 1).  Обозначим  узлы  cетки  через (x1,t1), …(xi,tj), …, а приближенные значения функции u(x, t) в узлах – uij. Тогда xi=ih, i=0,1,…,n, h=a/n, tj= jτ, j=0,1,…,m, τ=T/m.

Аппроксимируем уравнение (1) на четырехточечном шаблоне, который изображен на рис. 1 жирными линиями. В результате получаем неявную двухслойную разностную схему:

λui+1,j(1+2λ)ui, j + λui-1,j = – ui,j –1         (4)

которая аппроксимирует уравнение (1) с погрешнос тью О(τ + h2). Здесь λ = τ / h2.

Рис. 1     Рис. 2

Схема (4) аппроксимирует уравнение (1) только во внутренних узлах сетки, поэтому число уравнений в схеме (4) меньше числа неизвестных uij. Недостающие уравнения получают из краевых условий

.                                         ( 5 )

Схема (4)—(5) неявная, поэтому значения uij находят как решение системы линейных уравнений (4). Для решения системы (4) можно применять любой алгоритм решения систем линейных уравнений, однако система (4) обладает трехдиагональной матрицей и рациональнее всего решать ее методом прогонки. Таким образом, решив систему разностных уравнений, найдем значения функции и на временном слое j, если известно решение на временном слое j – 1.

Итак, алгоритм численного решения задачи имеет следующий вид. На нулевом временном слое (j = 0) решение известно из начального условия ui0 = f (xi). На каждом следующем слое искомая функция определяется как решение системы (4)—(5).

Метод прогонки

Уравнение (4) в общем виде записывается

ai ui-1 — 2bi ui + ci ui+1 =  fi ,   i= 1, 2, … , n-1             (6)

на левой границе области

c0u1 – 2b0u0 = f0                                                       (7)        

на правой границе области

 aNuN-1-2bNuN=fN.                                                    (8)

Имеем  u0,  u1  , u2  ,...

Будем искать решение в виде:

ui-1=Рiui+Qi                                                            ( * )

u0=c0/2b0·u1+(-f0/2b0)                                               (9)

 Р1 = c0 /2b0 , Q1 = - f0 / 2b0 ;

Для всех остальных Рi , Qi получим рекуррентные формулы, подставив (*) в (6):

ai (Рi  ui + Qi) — 2bi ui + ci ui+1 =  fi , 

(ai Рi  –2bi ) ui + ci ui+1 = fi – ai Qi

ui =

ui = Рi+1ui+1 + Qi+1 ,

                                                               

Если известны Р1, Q1, то можно найти Р2, Q2, и т.д.

uN-1 = PN uN +Q N ;                     

uN-1 = 2bN /aN ·uN + fN /aN ;

Имеем систему двух уравнений.                                       

Отсюда:      uN =  

Зная uN по формуле  (*) вычисляем все  ui от i = N-1 до 0.

Таким образом, находятся все ui на верхнем слое.

Замечательным свойством неявной схемы (4) является ее устойчивость при любых значениях параметра λ = τ / h2 > 0. Преимущества схемы (4) особенно ощутимы при сравнении с явной схемой, которая
получается при аппроксимации уравнения (1) на шаблоне, изображенном на рис. 2.

Явная схема оказывается устойчивой только при λ < 1/2 , т. е. при τ < h2 /2. Это означает, что вычисления по явной схеме придется вести с очень малым шагом по τ, что конечно может привести к большим затратам машинного времени. В неявной схеме вычисления на одном шаге требуют больше операций, чем в явной схеме, но зато величину шага τ можно выбирать, как угодно большой без риска нарушить устойчивость схемы. Все это позволяет значительно уменьшить машинное время, необходимое для решения задачи. Схема (4) обладает сходимостью. Это означает, что при h, τ → 0 решение разностной задачи (4) — (5) стремится к точному решению смешанной задачи (1)—(3).

Для численного решения задачи (1)—(3) предназначе на подпрограмма CALC.

Подпрограмма CALC решает систему (4) методом прогонки и вычисляет решение на каждом слое по значениям решения с предыдущего слоя.

Входные параметры: НХ — шаг сетки h по переменной х; НТ — шаг сетки τ его переменной t; N — количество шагов сетки по х; U — массив из N + 1 чисел, содержащий перед обращением к CALC начальное значение решения в узлах первого слоя сетки; Р, Q — рабочие массивы из N — 1 чисел каждый /это прогоночные коэффициенты/.

Выходные параметры: U — массив из N + 1 чисел (он же входной), содержащий значение решения в узлах очередного временного слоя сетки.

Перед обращением к программе CALC  необходимо:

1)описать массивы U, P, Q:

2)присвоить значения параметрам НХ, НТ, N;

3)присвоить элементам U начальное значение решения в узлах первого слоя сетки.

Задание. Используя программу CALC, решить краевую задачу с граничными условиями 2-го рода для однородного уравнения теплопроводности в области  D = { (х, t)|0 < х < а, 0 <t < Т}. Результаты печатать на каждом временном слое.

      Порядок   выполнения   лабораторной  работы

1. Составить  основную   программу,  содержащую  цикличес кое  обращение  к подпрограмме решения  параболического уравнения CALC с  печатью  результатов  на   каждом  шаге. Образец дан в приложении1.

2. Составить     подпрограмму CALC  решения  параболичес кого уравнения методом прогонки.

3. Составить   подпрограмму  -  функцию,  вычисляющую   значение   функции   F(x)  для  начального  значения  u( x, 0 ).

4. Провести   вычисления   на  компьютере.

5. Результаты  решения  нанести на мимиллиметровку в  виде  графика  изменения  функции u по времени.

Примечание:  для  более  плавного   поведения   функции  увеличить  начальное  количество  точек  разбиения  отрезка  X [ 0 , 1 ]  в два раза.

Варианты заданий

Решить краевую задачу для уравнения теплопровод ности с начальным условием  u(х, 0) = f (х)   и граничными условиями    . Во всех вариантах q по вариантам.

№ вари анта

f (x)

a

b

c

d

1

2

3

4

5

6

7

8

9

10

1,1

1.2

1,3

1,4

1,5

1,6

1,7

1,8

1,9

2,0

3,0

3,5

4,0

4,5

5,0

5,5

6,5

7,0

7,5

8,0

0,05

0,15

0,25

0,35

0,45

0,55

0,65

0,75

0,85

0,95

11

12

13

14

15

16

17

18

19

20

8,0

9,0

10,0

11,0

12,0

13,0

14,0

15,0

16,0

17,0

3,0

4,0

5,0

6,0

7,0

8,0

9,0

10,0

11,0

12,0

20,0

21,0

22,0

23,0

24,0

25,0

26,0

27,0

28,0

29,0

0,05

0,15

0,25

0,35

0,45

0,55

0,65

0,75

0,85

0,95

РЕКОМЕНДУЕМАЯ ЛИТЕРАТУРА

1. А.А.Самарский, А.В.Гулин. Численные  методы.  М.: Главная  ред. Физ.- мат. лит. - 1989 г., 430 с.

2. Бахвалов Н.С. Численные методы / Н. С.Бахвалов, Н. П. Жидков, Г.М. Кобельков. – изд.4-е; Гриф МО. – М.: Бином. Лаборатория знаний. 2006 г. – 636 с.

Приложение1

Листинг основного блока программы на Pascal

program PAR_5;

const ht=0.02;

     n=10;

     hx=1/n;

var j,i,i1,N1,N2:integer;

   x,y,t,T1,al,B1,A1:real;

   u:array [0..n] of real;

   P:array [0..n] of real;

   Q:array [0..n] of real;

   XP:array [0..n] of real;

begin

XP[0]:=0.0;

XP[10]:=1.0;

u[0]:=10;

u[10]:=20;

T1:=9*HT;

T:=0;

for i:=1 to n do

   begin

   x:=i*HX;XP[i]:=x;

   u[i]:=f(x);

   end;

writeln(T:8:5);

   for i:=0 to n do write(u[i]:6:3);

   writeln;

for j:=1 to 40 do

   begin

calc; T:=T+HT;

if (T>T1) then begin writeln(T:8:5);

   for i:=0 to n do write(u[i]:6:3);

   writeln;

   T1:=T1+10*HT;end;

   end;

end.

Методические указания

по выполнению лабораторной работы № 5

РЕШЕНИЕ КРАЕВОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ПАРАБОЛИЧЕСКОГО ТИПА РАЗНОСТНЫМ МЕТОДОМ

Сафронов Александр Иванович

Подписано в печать 10.09.11.Формат 60х84/16

Печать оперативная. Уч.п.л. 0,5. Уч.-изд.л.0,5.

Тираж 50 экз.

Тольяттинский государственный университет

Тольятти, Белорусская, 14




1. Реферат. Идея хеширования двойного блока состоит в том чтобы построить функцию сжатия на 2nбитах используя.
2. Я була іграшкою багатих чоловіків за іншимі ~ даним терміном вперше позначається колаж Р
3. Мне лучше знать. Пусть миссис Пирибингл твердит хоть до скончания века что она не может сказать кто начал пе
4.  Наблюдение является одним из наиболее часто используемых в психологии исследовательских методов
5. УлГПУ им ИН Ульянова СМК 021092013 ПОЛОЖЕНИЕ о конкурсе Science Slm стр
6. придуманы чтобы общая сумма материального обеспечения неработающего пенсионера не была ниже величины пр
7. Деятельность третейских судов
8. Налоги на прибыль организаций
9. Организация работы вагонного участка1
10. 1 Таблица 2 Рис
11. в качестве субъекта гражданин иностранного государства иностранная организация или иностранное государ
12. Yesterdy из альбома Help 1965 Посвящаю Ют и Наату Рэмпл сущ
13. Реферат- Институт возмещения вреда в XIX - начале XX века
14. Величины характеризующие вращательное движение
15. Архивная реформа
16. .1 История становления договора поставки в России
17. зима бродила по пустым улицам заходила в распахнутые двери и окна в явном раздражении разбрасывала по проме
18. ОреховоЗуевский медицинский колледж УЧЕБНАЯ СЕСТРИНСКАЯ ИСТОРИЯ БОЛЕЗН
19. тематический сбор и обработка информации об общественном мнении и его тенденциях о появляющихся проблемах и
20. Нильсен Протокольная часть Сектофоб