Перейти к содержанию
    

Уравнение Пуассона

Всем привет!

 

В книге [1] рассматривается численное решение динамики одномерной плазмы.

 

Одним из этапов моделирования является вычисление потенциала поля Fi(x) по заданной плотности распределения заряда Rho(x), то есть, решение уравнения Пуасона: d2Fi/dx2 = -Rho(x)/esilon0.

 

Уравнение Пуассона решается довольно муторно, через Фурье преобразование разностной схемы, потом прогонки и обратного фурье преобразования.

 

А зачем??!! Одномерный случай дает нам ОДУ, которое вполне решаемо, например Рунге-Куттом. Граничные условия сами напрашиваются, потенциал на удалении равен нулю, его производная – поле, тоже равно нулю подальше от заряда.

 

Однако, не могли авторы [1] ошибиться в таком критичном для программы месте. Где я ошибаюсь?

 

1. Бэдсел Ч., Ленгдон.А. Физика плазмы и численное моделирование.М.: Энергоатомиздат, 1989.

 

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Дошло, Рунге-Куттом не получается удовлетворить граничным условиям с ‘дальней’ стороны.

Граничные условия имеют вид Fi(0)=0, Fi(L)=0. А у Рунге-Кутта Fi(L) – какое получится.

 

 

Поделиться сообщением


Ссылка на сообщение
Поделиться на другие сайты

Присоединяйтесь к обсуждению

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

Гость
К сожалению, ваш контент содержит запрещённые слова. Пожалуйста, отредактируйте контент, чтобы удалить выделенные ниже слова.
Ответить в этой теме...

×   Вставлено с форматированием.   Вставить как обычный текст

  Разрешено использовать не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отображать как обычную ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставлять изображения напрямую. Загружайте или вставляйте изображения по ссылке.

×
×
  • Создать...