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




Скачать 69.83 Kb.
Название1. Постановка задачи
Дата публикации23.02.2013
Размер69.83 Kb.
ТипРешение
shkolnie.ru > Математика > Решение
1. Постановка задачи:
Необходимо найти функцию , в области , удовлетворяющую уравнению переноса
(1.1)
, начальному
(1.2)
и левому краевому условию.
(1.3)
2. Метод конечных разностей.
Метод конечных разностей, или метод сеток - такой приближённый метод решения задач, связанных с дифференциальными уравнениями, при котором область непрерывного изменения аргументов, в которой ищется решение уравнения, заменяется дискретным множеством точек (узлов), называемым сеткой; вместо функций непрерывного аргумента рассматриваются функции дискретного аргумента, определяемые в узлах сетки и называемые сеточными функциями; производные, входящие в уравнение, краевые и начальные условия, аппроксимируются разностными отношениями; интегралы аппроксимируются квадратурными формулами; при этом исходное уравнение заменяется системой (линейных, если исходная задача была линейной) алгебраических уравнений (системой сеточных уравнений, а применительно к дифференциальным уравнениям — разностной схемой).
- сеточная функция
- равномерная сетка


- конечные разности

Если полученная таким образом система сеточных уравнений разрешима, по крайней мере, на достаточно мелкой сетке, и её решение при неограниченном измельчании сетки сходится к решению исходного уравнения, то полученное на любой фиксированной сетке решение и принимается за приближённое решение исходного уравнения.
Метод сеток является одним из наиболее распространённых приближённых методов решения задач, связанных с дифференциальными уравнениями. Широкое применение Сеток метод объясняется его большой универсальностью и сравнительной простотой реализации на ЭВМ.
Решение линейной разностной задачи сходится к решению дифференциальной задачи, если разностная задача устойчива и аппроксимирует дифференциальную задачу на ее решении. При этом порядок аппроксимации совпадает с порядком сходимости.
2.1. Аппроксимация
Для того чтобы уравнения были как-то связаны с приближаемым дифференциальным оператором нужно, чтобы уравнения удовлетворяли многочленам. Если уравнения выполняются для всех многочленов степени не выше r с точностью , то говорят, что разностная схема имеет r-тый порядок аппроксимации.
2.2. Устойчивость
Условия аппроксимации не достаточно для того, чтобы результат разностной схемы приближался к точному ответу при h → 0. В случае схем, коэффициенты которых не зависят от решения дифференциального уравнения, нужно выполнение условия устойчивости. Такие схемы можно представить как некоторый линейный оператор, который преобразует значения функции в момент t в значения функции в момент t + τ. Условие устойчивости требует, чтобы собственные числа (вообще говоря, комплексные) этого оператора не превосходили по модулю 1+ch, где с — некоторая константа, при h→0. Если это условие не выполнено, то погрешности схемы быстро возрастают и результат тем хуже, чем меньше шаг. Если выполнены как условие аппроксимации, так и условие устойчивости, то результат разностной схемы сходится к решению дифференциального уравнения (теорема Филиппова - Рябенького).
^ 3. Явные и неявные разностные схемы по времени.
Явными схемами называются такие разностные схемы для эволюционных уравнений, когда данные на следующем слое по времени находятся непосредственно из данных на предыдущем слое без решения алгебраических систем уравнений. Если же на верхнем временном слое для определения значений сеточной функции необходимо решать систему алгебраических уравнений, то схема называется неявной.
Простейшая явная по времени схема для уравнения (1.1) имеет вид:
(3.1)
Простейшая неявная по времени схема имеет вид:
(3.2)
Ниже показаны шаблоны схемы для (3.1) и (3.2), т.е. конфигурации расчетных узлов в области интегрирования, используемых на каждом элементарном шаге вычислений

Явная схема по времени.


Неявная схема по времени.
Направление "уголка" (левый, правый) будет рассмотрено позднее, в связи с его влиянием на устойчивость решения.
Решение разностной задачи будет условно устойчивым, если σ < 1, где
- число Куранта (3.3)
При σ > 1 сходимости разностной задачи ожидать не следует.
^ 4. Схема с центральными разностями.
Схемы с центральными разностями имеют симметричный шаблон. Например, схема со вторым порядком точности по координате и времени имеет вид:
(4.1)
Второй порядок точности означает, что погрешность, обусловленная заменой дифференциалов конечными разностями, имеет величину второго порядка малости, т.е.
(4.2)

(4.3)
Вычитаем (4.3) из (4.2) и находим:
(4.4)
Аналогичное выражение будет и с производной по времени.
Для аппроксимации начальных данных необходимо задать не только u(m,0), но и u(m,1), иначе временной слой u(m,2) уже не вычислить. Происходит это из-за того, что исходное дифференциальное уравнение имеет первый порядок, а разностная задача - второй.
^ 5. Схема против потока.
Схема против потока, или схема Куранта - Изаксона - Риса (КИР) имеет вид:
(5.1)

(5.2)
Пусть для определённости a > 0.

Тогда схему можно записать как:
(5.3)
, где σ определена ранее (3.3)

Схема КИР условно устойчива при выполнении условия Куранта, т.е. σ < 1. Преимущество схем против потока состоит в более точном учете области зависимости решения.
^ 6. Практическая часть.
Найдём численное решение задачи в области :
(6.1)
, где (6.2)
, начальное условие:

(6.3)
, левое граничное условие:

(6.4)
В качестве разностной схемы выберем явную по времени схему "против потока" (5.1)
Уравнение будем рассматривать на сетке 101 х 101, т.е с шагом

(6.5)
Число Куранта выходит σ = 1/2 < 1, что позволяет рассчитывать на сходимость решения.
^ 7. Листинг программы.
(Программа написана в среде Turbo Pascal 7.0)
Начальные условия задаются функцией F_N(i,j), неоднородность - функцией F_F(i,j)

Выходные даные записываются в файл "data.txt"


Program kurant;
Function F_N(i : integer) : real;
var

x : real;

s : real;
begin

x := 0.01*i;

s := sin(6.283185*x);

F_N := 0.0;

if (i < 50) then

{F_N := 8.0*x*(0.5 - x);}

F_N := s*s*s;

end; {function F_N}
Function F_F(i,j : integer) : real;
var

x : real;

t : real;
begin

x := 0.01*i;

t := 0.01*j;

{F_F := 0.0;}

F_F := x*(1-x)*t;

end; {function F_F}
var

lett : char;

u : array [0..100,0..100] of real;

a : real;

m,n,i,k : integer;

t,h,s : real;

fout : text;
const

tab = #9;
begin


assign(fout,'data.txt');

rewrite(fout);
a := 0.5;

h := 0.01;

t := 0.01;

s := a*t/h;
{nachalnoe uslovie}

for i := 0 to 100 do

begin

u[i,0] := F_N(i);

end;
{granichnoe uslovie}

for i := 0 to 100 do

begin

u[0,i] := 0.0;

end;
{reshenie}

for n := 1 to 100 do

begin
for m := 1 to 100 do

begin

u[m,n] := (1 - s) * u[m,n-1] + s * u[m-1,n-1] + t*F_F(m,n);

end;

end;
for i := 0 to 100 do

begin

writeln(fout,i,tab,u[i,0]:2:4,tab,u[i,50]:2:4,tab,u[i,100]:2:4);

end;
close(fout);
WriteLn('');

writeln('Enter a letter to end the program');

Readln(lett);

end.

8. Результаты.
Для начала проверим численное решение однородного уравнения переноса, т.е. f = 0.
Очевидно, что решением дифференциального уравнения (6.1) будет сдвиг функции из начальных условий вправо вдоль оси x со скоростью 1/2
Кривая из чёрных элементов - начальное условие, из красных - решение через t = 0.5, из зелёных - при t = 1.0


Рис. 8.1
В действительности наблюдается некое уменьшение "амплитуды", что связано с погрешностью усечения, т.е. замены дифференциалов на конечные разности. Условная устойчивость разностной задачи позволяет рассчитывать, что данная ошибка будет уменьшаться при дальнейшем измельчении сетки.

Численное решение при заданной неоднородности (6.2) изображено ниже.


Рис. 8.2
Цветовые обозначения соответствую рис. 8.1.

Список использованной литературы:
1. А.И.Лобанов, И.Б.Петров "Численные методы решения уравнений в частных производных";
2. Н.С.Бахвалов "Численные методы";
3. В.С.Рябенький, А.Ф.Филиппов "Об устойчивости разностных уравнений";
4. В.В.Фаронов "Turbo Pascal 7.0 Начальный курс"




Похожие:

1. Постановка задачи iconПостановка задачи Описание входных и выходных данных. Математическая...
Создать сайты по предоставлению следующих сервисов. Сервисы должны быть доступны только для зарегистрированных пользователей. Информация...
1. Постановка задачи iconВозможная схема самоанализа урока
Постановка задач урока ( особенно ценны исследовательские задачи, задачи в увлекательной форме)
1. Постановка задачи iconПостановка задачи
Постановка задачи 3
1. Постановка задачи iconБилет Этапы развития теории программирования
Программирование велось в основном в машинном коде. Решались в главным образом научно-технические задачи. Была точная постановка...
1. Постановка задачи iconИнструкция пользователю 24
Постановка задачи 3
1. Постановка задачи iconЗадание № Одномерные массивы Постановка задачи
Анализ постановки задачи показывает, что в алгоритме задачи должны быть элементы циклической и разветвляющейся структур. Исходные...
1. Постановка задачи iconВсероссийский интернет конкурс педагогического творчества ( 2013/14 учебный год)
Задачи: повторить, что такое бсп; проверить, как 9-классники усвоили правила постановки запятой и точки с запятой в бсп; вспомнить...
1. Постановка задачи iconУрок по теме "Отечественная война 1812 года"
Первый этап – подготовительный: постановка цели, озвучивание темы урока, постановка проблемы
1. Постановка задачи iconНазвание занятий
...
1. Постановка задачи iconТема: Разработка реляционной базы данных
Постановка задачи: упорядочить информацию о клиентах фирмы, ассортименте продукции и сформированных заказах
Вы можете разместить ссылку на наш сайт:
Школьные материалы


При копировании материала укажите ссылку © 2014
shkolnie.ru
Главная страница