Зарегистрироваться

Сообщество «Физика»: Обсуждение

Столкновение бризеров, Какие начальные условия нужны для этого?

01 Сентябрь 2010 | Автор: Коновалов Константин Александрович

Численно смоделировал бризероподобное решение, по начальным условиям: z(x,t)=1/ch(x/a) z'(x,t)=iw/ch(x/a). z' - производная по t. Обе при t=0. Применил преобразование Лоренца, и получил движущийся бризероподобное решение. Вопрос: Как смоделировать - столкновение двух бризеров? Поставил вместо вышеуказанного нач. условия (для z'), это: z'(x,t)=1/ch(x+V*tao)+1/ch(x-V*tao) (V скорость, tao - шаг по времени, для разностной схемы 2-го порядка точности), но получаю - один бризер, который стоит на месте, и "дышит" как и прежде. (вычисления производил в Фортране PS. 4.0, визуализация в Матлабе. 7.5) Спасибо зараннее.

Только авторизованные пользователи могут участвовать в обсуждении.

А в чём смысл так-то менять? Вам надо тоже сложить его из двух бегущих волн, так же, как вы и z' сложили.
Получилось! Большое спасибо. Получил: 1. Лобовое столкновение; 2. Столкновение движущегося на неподвижный; 3. "Обгоняющее". кажется для одномерного бризера это все варианты, если не ошибаюсь. Совсем скоро хочу перейти на столкновения двумерных бризеров. Я моделировал двумерный но не перемещающийся бризер (Не смог "передвинуть). Хотя то решение не очень "чистое", в смысле - наблюдается затухание, незаметно, но затухает.
Нет. Стоп. Проблема в другом. Вы взяли два солитона, и наложили их пространственно в одном месте. Но это ошибка, потому что когда солитоны взаимодействуют, они складываются нелинейно. Сохраняют свою формул они только вдалеке от взаимодействия (на асимптотической бесконечности). Поэтому вам надо взять два начальных условия вид и разнести их на большое расстояние, то есть, например, взять для одного большое положительное X, а для другого противоположное X. И преобразовать каждое по Лоренцу, чтобы они двигались друг к другу. И вот их уже сложить.
А в чём смысл так-то менять? Вам надо тоже сложить его из двух бегущих волн, так же, как вы и z' сложили.
z - я тоже менял (просто умножил на 2) следующим образом: z(x,t)=2/ch(x/a) (Сорри, что не упоминул. Также, в выражении z' пропустил "/a", т.е. правильный вариант: z'(x,t)=1/ch((x+V*tao)/a)+1/ch((x-V*tao)/a)) В инете искал, и на английском, не нашел решения. Что можно еще сделать?
А почему вы z' поменяли, а z нет?