Справочник функций

Ваш аккаунт

Войти через: 
Забыли пароль?
Регистрация
Информацию о новых материалах можно получать и без регистрации:

Почтовая рассылка

Подписчиков: -1
Последний выпуск: 19.06.2015

Решение параболического уравнения

74K
13 ноября 2011 года
saabru
4 / / 01.10.2011
Задача:Разработать программу для решения параболического уравнения
Uxx+Ut=0, используя явную схему, для следующих условий: x (0, 10), t (0, 100);
граничное значение: U(0, t)=100; U(10, t)=200;
начальное условие


Проверьте, кто может, правильность моего решения, а то меня терзают сомнения.:confused:

И ещё,может кто подскажет, как лучше визуализировать решение под borland c++ 3.1?
278
13 ноября 2011 года
Alexander92
1.1K / / 04.08.2008
Не смотрел ваше решение, у меня сразу вопрос к условию: как вы можете даже пытаться решать уравнение, у которого начальные и граничные условия рассогласованы? Смотрите, у вас

* из граничных условий: U(0, t) = 100 => U(0, 0) = 100;
* из начального условия: U(x, 0) = sin(Pi*x/5) => U(0, 0) = 0.
74K
14 ноября 2011 года
saabru
4 / / 01.10.2011
Вот да. Вот в этом вся фишка задачи. Если у нас синусоида, то у нас граничные условия должны быть нулевыми. Вот что получается у меня,если они нулевые Но преподаватель утверждает, что всё задано правильно.
74K
14 ноября 2011 года
saabru
4 / / 01.10.2011
Если выкинуть начальное условие (предположить, что начальное условие постоянно по x) , то задача решается. Вот решение на 81 слое 100 109.6389769 119.3124282 129.0548284 138.887484 148.8317019 158.8874823 169.0548256 179.3124253 189.6389752 200
А как решать с такими условиями, я не имею понятия. У кого какие есть ещё варианты решения задачи?
278
14 ноября 2011 года
Alexander92
1.1K / / 04.08.2008
Ну, для нулевых граничных условий задача имеет решение в замкнутом виде. :)
 
Код:
U(x, t) = sin(pi*x/5) * exp(pi^2*t/25)

Можете сравнить это с тем, что вы получили.

[QUOTE=saabru]
преподаватель утверждает, что всё задано правильно.
[/QUOTE]
А вы ему показывали то, что я написал в предыдущем посте? Чем он это объясняет?

Да, и еще вопрос такой... Честно говоря, у вас немного странное уравнение. По картинке, я смотрю, вы ищете какую-то температуру, т.е. предполагаете, насколько я понимаю, что ваша система должна прийти в термодинамическое равновесие. А уравнение описывает, напротив, экспоненциальный уход от состояния равновесия. Проверьте знак перед вторым членом в самом уравнении, уравнение Uxx - Ut = 0 для температуры гораздо логичнее.

[QUOTE=saabru]
А как решать с такими условиями, я не имею понятия.
[/QUOTE]
А никак. Бред - он и есть бред.
74K
14 ноября 2011 года
saabru
4 / / 01.10.2011
Да это уравнения теплопроводности. Минус точно должен быть,т.к. uxx = ut. Да ничем он не объяснил, ссылался на отсутствие времени. Сказал только, что решение меня сходится и dt dx я выбрал правильно, т.к. сигма=dt/dx*dx <0.5 :facepalm:
278
14 ноября 2011 года
Alexander92
1.1K / / 04.08.2008
Цитата: saabru
сигма=dt/dx*dx <0.5


Это к чему?

Знаете кого-то, кто может ответить? Поделитесь с ним ссылкой.

Ваш ответ

Реклама на сайте | Обмен ссылками | Ссылки | Экспорт (RSS) | Контакты
Добавить статью | Добавить исходник | Добавить хостинг-провайдера | Добавить сайт в каталог