Уравнение пуассона и лапласа для электростатического поля. Уравнение пуассона и математическая постановка задач электростатики

Я хотел бы в познавательных целях рассказать об уравнениях, которые применялись при выводе уравнения Дебая-Хюккеля. Это уравнение Пуассона и распределение Больцмана.

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

Мы выяснили, что плазма квазинейтральна в равновесном состоянии и что под действием электрического поля от движущихся зарядов, заряженные частицы смещаются на дебаевскую длину и поле в пределах этой длины затухает. В электростатике взаимодействие заряженных частиц описывается кулоновским уравнением:

Где – величины взаимодействующих точечных зарядов, – квадрат расстояния между зарядами. Коэффициент k является константой. Если мы используем систему в электростатических единицах СГС, обозначаемых СГСЭq, то k = 1. Если используется система СИ, то , где – диэлектрическая проницаемость среды, в которой расположены заряды, – электрическая постоянная, равная 8,86 ∙ .

В физике непосредственно силой не пользуются, а вводят понятие электростатического поля распределённых зарядов и измеряют поле величиной напряженности электрического поля . Для этого в каждую точку поля мысленно помещают единичный пробный заряд и измеряют силу, с которой поле зарядов действует на пробный заряд:


Отсюда, если подставить в это уравнение силу Кулона, то получим:
Но и этим физики не ограничиваются, для того чтобы описать полноценно электрическое поле. Рассмотрим единичный заряд, помещённый в электростатическое поле. Поле выполняет работу по перемещению этого заряда на элементарное расстояние ds из точки P1 в точку P2:
Величину называют разностью потенциалов или напряжением. Напряжение измеряется в Вольтах. Знак минус говорит нам о том, что само поле выполняет работу для переноса единицы положительного заряда. Силы, перемещающие заряды являются консервативными, так как работа по замкнутому пути равна всегда нулю, независимо от того, по какому пути перемещается заряд.

Отсюда следует глубокий смысл разности потенциалов. Если зафиксировать точку Р1 и перемещать заряд в переменную точку Р2, то работа зависит только от положения второй точки Р2. Таким образом мы можем ввести понятие потенциала. Потенциал – это силовая функция, показывающая какую необходимо выполнить работу полю, чтобы переместить заряд из бесконечности в данную точку P2, где условно принимают потенциал в бесконечности равным нулю.

Чтобы понять уравнение Пуассона, необходимо разбираться в «особой» векторной математике. Я вкратце расскажу про такие понятия как градиент поля и дивергенции (подразумевается, что читатель знаком с математическим анализом)
Пусть f(x,y,z) является некоторой непрерывной дифференцируемой функцией координат. Зная её частные производные в каждой точке пространства можно построить вектор, компоненты которого x, y, z равны соответствующим частным производным:


где – единичные векторы соответствующих осей x, y, z. Значок читается «набла» и является дифференциальным оператором
Этот оператор ввёл в математику Гамильтон. С набла можно выполнять обычные математические операции, такие как обычное произведение, скалярное произведение, векторное произведение и так далее.

Теперь вернёмся к электростатическому полю E. С одной стороны изменение потенциала при переходе из одной точки в другую имеет следующий вид:


С другой стороны, согласно формуле (*)
Применяя только что введённое понятие градиент, эта формула преобразуется в:
Теперь разберёмся с таким понятием, как дивергенция поля. Рассмотрим конечный замкнутый объем V произвольной формы (см. рис. ниже). Обозначим площадь этой поверхности S. Полный поток вектора F, выходящего из этого объема по определению равно
, где da является бесконечно малым вектором, величина которого равна площади малого элемента поверхности S, а направление совпадает с наружной нормалью к этому элементу.
Возьмём этот поток вектора F поделим на объём и найдём предел при стремящейся к нулю, т.е. будем стягивать объём в бесконечно малую точку.


Мы подошли к понятию дивергенции. Обозначается дивергенция символом div и является отношением потока вектора F к объёму V, при V стремящейся к нулю.

Прежде чем показать, как получается уравнение Пуассона, важно знать закон Гаусса и теорему Гаусса. Представим себе сферу, внутри которой находится заряд q. Заряд создаёт вокруг себя электрическое поле напряжённости E. Возьмём поток вектора E


где S площадь нашей сферы равная . Следовательно
Это и есть закон Гаусса, утверждающий, что поток электрического поля E через любую замкнутую поверхность равен произведению на полный заряд, охватываемый поверхностью:
где – плотность объёмного заряда, т.е. величина электрического заряда в единице объёма, и – элементарный объём, выделенный внутри нашего замкнутого объёма.

Теорема Гаусса (полное название теорема Гаусса-Остроградского) чисто математическая теорема о дивергенции. Перепишем полный поток вектора F следующим образом:


В пределе, когда N → ∞, →0 величина в скобках становится дивергенцией и сумма переходит в объёмный интеграл:
Это и есть теорема Гаусса, и является поистине самой важной формулой полевой теории. Применим эту теорему к электростатическому полю. С одной стороны, согласно закону Гаусса
А с другой стороны, согласно теореме Гаусса (только не путайте теорему с законом Гаусса):
Комбинируя два последних уравнения, получим:
Вспомним формулу (**) и подставим сюда вместо E потенциал поля
Дивергенция градиента это новый оператор, который в математике называют оператор Лапласа, или сокращённо лапласиан. Лапласиан обозначается значком набла следующим образом и равен
Перепишем предыдущую формулу в форме лапласиана:
Наконец мы получили уравнение Пуассона. В первой статье это уравнение было немного в другой форме, с учётом диэлектрической проницаемости среды. Вспомните силу Кулона в системе СИ, там константа . Соответственно в законе Гаусса будет не , а коэффициент . Таким образом получаем уравнение Пуассона в форме представленной в предыдущей статье
Таким образом по сути уравнение Пуассона – это закон Кулона (а точнее закон Гаусса) переписанный в другой форме, в обозначениях векторного дифференциального анализа.

В мы разберём важное распределение из математической статистики - распределение Больцмана.

Теги:

  • физика
  • электростатики
Добавить метки

К исследованию уравнений Лапласа и Пуассона приводит рассмотрение задач о стационарном процессе: это задачи гидродинамики, диффузии, распределения температуры, электростатики и др.

Эти уравнения относятся к уравнениям эллиптического типа.

Те задачи, которые приводят к уравнениям, содержащим время, называются нестационарными или динамическими задачами математической физики; задачи, приводящие к уравнениям, не содержащим время, называются стационарными или статическими.

Как было показано, уравнения математической физики имеют бесчисленное множество решений, зависящие от двух произвольных функций (речь идет об уравнениях второго порядка для функции двух переменных). Для того, чтобы из множества решений выделить определенное, характеризующее процесс, необходимо на искомую функцию наложить дополнительные условия, которые диктуются физическими соображениями. Таковыми условиями для уравнений в частных производных являются, чаще всего, начальные и граничные условия. Граничные условия - это условия заданные на границе рассматриваемой среды; начальные условия - это условия, относящиеся к какому-нибудь моменту времени, с которого начинается изучение данного физического явления. Дополнительные условия, так же как и само дифференциальное уравнение, выводятся на основе физических соображений, связанных с самим процессом. Вместе с тем дополнительные условия должны быть такими, чтобы обеспечить выделение единственного решения из всего множества решений. Число граничных и начальных условий определяются типом уравнения, а их вид - заданным исходным состоянием на границе объекта и внешней среды. Для рассматриваемых нами уравнений число начальных условий равно порядку старшей производной по времени, входящей в уравнение, а число граничных условий - порядку старшей производной по координате.

Совокупность дифференциального уравнения и дополнительных условий представляют собой математическую формулировку физической задачи, и называется задачей математической физики.

Итак, задача математической физики состоит в отыскании решений уравнений в частных производных, удовлетворяющих некоторым дополнительным условиям, скажем, граничным и начальным.

Задача математической физики считается поставленной корректно, если решение задачи, удовлетворяющее всем ее условиям, существует, единственно и устойчиво.

Колебания струны. Граничные и начальные условия. Постановка краевых задач

Пусть струна находится под действием сильного начального натяжения. Если вывести струну из положения равновесия и подвергнуть действию какой-либо силы, то струна начнет колебаться. Процесс колебания можно описать одной функцией, характеризующей вертикальное перемещение струны (отклонение от положения равновесия (рис. 2.2)). При каждом фиксированном значении график функции на плоскости дает форму струны в момент времени.

Функция удовлетворяет уравнению

описывает свободные колебания струны без воздействия внешних усилий.

Уравнение (2.69) является простейшим уравнением гиперболического типа и в то же время одним из важнейших уравнений матфизики.

Одного уравнения движения (2.69) или (2.70) при математическом описании физического процесса недостаточно. При рассмотрении задачи о колебании струны дополнительные условия могут быть двух видов: начальные и граничные (краевые).

Так как процесс колебаний струны зависит от ее начальной формы и распределения скоростей, то следует задать начальные условия:

Будем говорить о трех типах граничных условий:

где известные функции,

и известные постоянные.

Приведенные условия называют соответственно граничными условиями первого, второго, третьего рода. Условия I имеют место в том случае, когда концы объекта (струна, стержень и т.д.) перемещаются по заданному закону; условия II - в случае, когда к концам приложены заданные силы; условия III - в случае упругого закрепления концов.

Если функции, заданные в правой части равенства, равны нулю, то граничные условия называются однородными. Так, граничные условия (2.72) - однородные. Комбинируя различные перечисленные типы граничных условий, получим шесть типов простейших краевых задач.

В том случае, когда режим на концах не будет оказывать существенного влияния на ту часть струны, которая достаточно удалена от них, струну считают бесконечной. В силу этого вместо полной краевой задачи ставят предельную задачу - з а д а ч у К о ш и: найти решение уравнения (2.69) для при, удовлетворяющее начальным условиям

Если изучается процесс вблизи одной границы и влияние граничного режима на второй границе не имеет существенного значения на протяжении интересующего нас промежутка времени, то приходим к постановке задачи на полуограниченной прямой. В этом случае задаются начальные условия и одно из граничных условий I - III при.

Примеры решения задач

ПРИМЕР 2.42. Однородная струна длины совершает малые поперечные колебания. Поставить задача об определении отклонений точек струны от прямолинейного положения покоя, если в момент струна имела форму () и скорость каждой ее точки задается функцией. Рассмотреть случаи:

  • а) концы струны закреплены;
  • б) концы струны свободны;

в) к концам струны и, начиная с момента, приложены поперечные силы и соответственно;

г) концы струны закреплены упруго, т.е. каждый из концов испытывает сопротивление, пропорциональное отклонению конца.

Решение. Как известно, отклонения точек струны от положения равновесия удовлетворяют в отсутствии действующей внешней силы уравнению свободных колебаний (2.70)

Здесь, натяжение, линейная плотность, т.к. струна однородная.

Начальные условия имеют вид:

Займемся выводом граничных условий.

Случай а). Так как концы струны закреплены, то их отклонения в точках и должны быть равными нулю при любом, т.е.

Итак, физическая задача о колебаниях закрепленной на концах струны свелась к следующей математической задаче: найти функцию, определенную при и, являющуюся решением уравнения

и удовлетворяющую граничным условиям

и начальным условиям

Уравнение (10.2) устанавливает связь между потенциалом электростатического поля и напряженностью этого поля. Из этого уравнения можно получить соотношение между потенциалом и плотностью заряда. Для этого нужно образовать дивергенцию обеих частей этого уравнения и воспользоваться затем формулой (6.5):

Согласно правилам векторного анализа [см. уравнение (40]

так что уравнение (11.1) может быть записано так:

Это дифференциальное уравнение носит название уравнения Пуассона. В тех участках поля, где нет электрических зарядов

Уравнение это обращается в следующее:

Этот частный вид уравнения Пуассона носит название уравнения Лапласа.

Уравнение Пуассона дает возможность определить потенциал поля объемных зарядов, если известно расположение этих зарядов. Решение (интеграл) этого дифференциального уравнения (при определенных граничных условиях) должно, очевидно, совпадать с выведенной нами ранее формулой (8.8):

В дальнейшем мы докажем это непосредственным вычислением. Пока же отметим, что для решения некоторых задач удобнее исходить не из интеграла (8.8), а непосредственно из дифференциального уравнения (11.3).

Пример. Определить плотность термоионного тока между двумя бесконечными плоскими электродами в вакууме. Пример этот на применение уравнения Пуассона взят не из электростатики, а из учения о токе и имеет большое значение для теории катодных (усилительных) ламп.

Известно, что накаленные металлы испускают со своей поверхности в окружающее пространство поток свободных электронов. Если к двум металлическим электродам приложить определенную разность потенциалов и раскалить отрицательный электрод (катод), то непрерывно испускаемые накаленным катодом электроны будут притягиваться к поверхности положительного электрода (анода). Поток электронов, движущихся от катода к аноду, эквивалентен электрическому току. Ток этот называется термоионным.

Выберем оси декартовых координат так, чтобы начало их находилось на катоде, а ось х была перпендикулярна плоскости электродов и направлена к аноду. Примем потенциал катода равным нулю, а потенциал анода равным Из соображений симметрии явствует, что эквипотенциальные поверхности параллельны электродам, поэтому и уравнение Пуассона в пространстве между электродами принимает вид

Если обозначить через число электронов, приходящихся на единицу объема в пространстве между электродами на расстоянии х от катода, а через абсолютную величину заряда электрона, то плотность заряда на

этом расстоянии будет:

Предположим для простоты, что испускаемые катодом электроны при выходе из его поверхности не обладают никакой начальной скоростью. На пути от катода к аноду силы электрического поля будут совершать над электронами заряда работу - которая будет, очевидно, переходить в кинетическую энергию движения электронов. Обозначая через скорость электрона на расстоянии х от катода, а через потенциал на том же расстоянии, получим

где 771 - масса электрона. Наконец, плотность электрического тока, т. е. заряд, протекающий за единицу времени через перпендикулярную току (т. е. перпендикулярную оси площадку в равна, очевидно:

ибо есть число электронов, проходящих за единицу времени через эту площадку. В отличие от плотность тока есть величина постоянная, не зависящая от х, ибо по достижении стационарного состояния через любую параллельную электродам плоскость проходит, очевидно, одинаковое число электронов.

Исключим из уравнения (11.5) все неизвестные функции х, кроме Прежде всего

Но из (11.6) следует, что

стало быть,

Вводя обозначение А - получим

Как легко убедиться подстановкой, из решений этого дифференциального уравнения, которое, согласно условию задачи, обращается на катоде в нуль и, кроме того, удовлетворяет условию

Если обозначить расстояние от анода до катода через I, то при потенциал должен обращаться в Стало быть,

Таким образом, плотность термоионного тока не подчиняется закону Ома, а растет пропорционально степени 3/2 приложенного к электродам напряжения и обратно пропорционально квадрату расстояния между ними. Это отличие законов термоионного тока от законов тока в металлах обусловливается двоякого рода причинами. Во-первых, электроны в металлах соударяются с положительными ионами, образующими твердый скелет металла, и испытывают благодаря этому сопротивление своему движению, отсутствующее при движении в вакууме 1). Во-вторых, при термоионном токе в пространстве между электродами находятся лишь свободные электроны, заряд которых не компенсируется зарядом положительных ионов, как это имеет место в металлах, вследствие чего поле этого так называемого «пространственного заряда» искажает поле электродов.

Отметим, что формула (11.9) перестает быть справедливой при больших плотностях тока 2). При повышении потенциала анода наступает момент, когда все выделяемые катодом электроны немедленно же увлекаются к аноду. Дальнейшее повышение потенциала анода не может, очевидно, повести к увеличению плотности тока, которая, таким образом, достигает постоянного значения (ток насыщения).

Задача 10. Пусть означает расстояние данной точки пространства от некоторой произвольно выбранной начальной точки Показать, что скаляр

удовлетворяет уравнению Лапласа

Точка не рассматривается.

Задача 11. Бесконечная плоская пластина толщиной 2а равномерно заряжена электричеством с объемной плотностью Ось х перпендикулярна пластине, начало координат расположено в срединной плоскости, равноотстоящей от обеих поверхностей пластины. Показать, что потенциал поля внутри и вне пластины равен соответственно:

а вектор направлен вдоль оси х от срединной плоскости и численно равен:

Сравнить этот случай с предельным случаем бесконечной заряженной плоскости (§ 4).

Задача 12. Найти потенциал поля шара, равномерно заряженного по своему объему [формула (8.12)], исходя из уравнения Пуассона в сферических координатах.

Уравнения Пуассона и Лапласа являются основными уравнениями электростатики. Они вытекают из теоремы Гаусса в дифференциальной форме. Действительно, известно, что Е = - grad j . В то же время согласно теореме Гаусса

Подставим в (11.22) E из (11.7). Получим

.

Вынесем минус за знак дивергенции

.

Вместо того чтобы писать gradj, запишем его эквивалент Ñj. Вместо div напишем Ñ. Тогда

Уравнение (11.27) называется уравнением Пуассона. Частный вид уравнения Пуассона, когда ρ свб =0, называется уравнением Лапласа. Уравнение Лапласа запишется так:

Оператор называют оператором Лапласа или лапла­сианом и иногда обозначают еще символом D. Поэтому можно встретить иногда и такую форму записи уравнения Пуассона:

Раскроем в декартовой системе координат. С этой целью произведение двух множителей Ñ и запишем в развернутом виде

Произведем почленное умножение и получим

.

Таким образом, уравнение Пуассона в декартовой системе координат запишется следующим образом:

. (11.29)

Уравнение Лапласа в декартовой системе координат

. (11.30)

Приведем без вывода выражения Ñ 2 j в цилиндрической системе координат

, (11.31)

в сферической системе координат (11.32)

Уравнение Пуассона дает связь между частными производными второго порядка от j в любой точке поля и объемной плотностью свободных зарядов в этой точке поля. В то же время потенциал j в какой-либо точке поля зависит, разумеется, от всех зарядов, создающих поле, а не только от величины свободного заряда, находящегося в данной точке.

Уравнение Лапласа (1780 г.) первоначально было применено для описания потенциальных полей небесной механики и впоследствии было использовано для описания электрических полей. Уравнение Пуассона применяется к исследованию потенциальных полей (электрических и магнитных) с 1820 г.

Рассмотрим вопрос о том, как в общем виде может быть записано решение уравнения Пуассона. Пусть в объеме V есть объемные (r), поверхностные (s) и линейные (t) заряды. Эти заряды представим в виде совокупностей точечных зарядов rdV, sds, tdl; dV - элемент объема, ds -элемент заряженной поверхности, dl - элемент длины заряженной оси. Составляющая потенциала dj в некоторой точке пространства, удаленной от rdV на расстояние R , в соответствии с формулой (11.20) равна

Составляющие потенциала от поверхностного и линейного зарядов, рассматривая их как точечные, определим аналогичным образом:

Полное значение j определится как сумма (интеграл) составляющих потенциала от всех зарядов в поле:

. (11.33)

В формуле (11.33) r,s и t есть функции радиуса R . Практически формулой (11.33) пользуются редко, так как распределение s по поверхности, t по длине и r по объему сложным образом зависит от конфигурации электродов и, как правило, перед проведением расчета неизвестно. Другими словами, неизвестно, как r, s и t зависят от радиуса R .


Граничные условия

Под граничными условиями понимают условия, которым подчиняется поле на границах раздела сред с различными электрическими свойствами. При изучении раздела «переходные процессы» исключительно большое значение имел вопрос о начальных условиях и о законах коммутации. Начальные условия и законы коммутации позволяли определить постоянные интегрирования при решении задач классическим методом. В классическом методе они использовались в явном виде, в операторном методе - в скрытом. Без использования их нельзя решить ни одной задачи на переходные процессы.

Можно провести параллель между ролью граничных условий в электрическом (и в любом другом) поле и ролью начальных условий и законов коммутации при переходных процессах. При интегрировании уравнения Лапласа (или Пуассона) в решение войдут постоянные интегрирования. Их и определяют, исходя из граничных условий. Прежде чем перейти к подробному обсуждению граничных условий, рассмотрим вопрос о поле внутри проводящего тела в условиях электростатики.

ОПРЕДЕЛЕНИЕ

Описывает адиабатный процесс, протекающий в . Адиабатным называют такой процесс, при котором отсутствует теплообмен между рассматриваемой системой и окружающей средой: .

Уравнение Пуассона имеет вид:

Здесь – объем, занимаемый газом, – его , а величина называется показателем адиабаты.

Показатель адиабаты в уравнении Пуассона

В практических расчётах удобно помнить, что для идеального газа показатель адиабаты равен , для двухатомного – , а для трёхатомного – .

Как же быть с реальными газами, когда важную роль начинают играть силы взаимодействия между молекулами? В этом случае показатель адиабаты для каждого исследуемого газа можно получить экспериментально. Один из таких методов был предложен в 1819 году Клеманом и Дезормом. Мы наполняем баллон холодным газом, пока давление в нём не достигнет . Затем открываем кран, газ начинает адиабатически расширяться, а давление в баллоне падает до атмосферного . После того, как газ изохорно прогреется до температуры окружающей среды, давление в баллоне повысится до . Тогда показатель адиабаты можно рассчитать за формулой:

Показатель адиабаты всегда больше 1, поэтому при адиабатическом сжатии газа – как идеального, так и реального – до меньшего объема температура газа всегда возрастает, а при расширении газ охлаждается. Это свойство адиабатического процесса, называемое пневматическим огнивом, применяется в дизельных двигателях, где горючая смесь сжимается в цилиндре и воспламеняется от высокой температуры. Вспомним первый закон термодинамики: , где — , а А – выполняемая над ней работа. Поскольку то работа, осуществляемая газом, идёт только на изменение его внутренней энергии – а значит, температуры. Из уравнения Пуассона можно получить формулу для расчёта работы газа в адиабатном процессе:

Здесь n – количество газа в молях, R – универсальная газовая постоянная, Т – абсолютная температура газа.

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

Стоит помнить, что уравнение Пуассона точно описывает только равновесный адиабатный процесс, состоящий из непрерывно сменяющих друг друга состояний равновесия. Если же мы в реальности откроем кран в баллоне, чтобы газ адиабатически расширился, возникнет нестационарный переходной процесс с завихрениями газа, которые затухнут из-за макроскопического трения.

Примеры решения задач

ПРИМЕР 1

Задание Одноатомный идеальный газ адиабатически сжали так, что его объем увеличился в 2 раза. Как изменится давление газа?
Решение Показатель адиабаты для одноатомного газа равен . Однако его можно рассчитать и по формуле:

где R – универсальная газовая постоянная, а і – степень свободы молекулы газа. Для одноатомного газа степень свободы равен 3: это значит, что центр молекулы может совершать поступательные движения по трём координатным осям.

Поэтому показатель адиабаты:

Представим состояния газа в начале и конце адиабатного процесса через уравнение Пуассона:

Ответ Давление уменьшится в 3,175 раза.

ПРИМЕР 2

Задание 100 молей двухатомного идеального газа адиабатически сжали при температуре 300 К. При этом давление газа увеличилось в 3 раза. Как изменилась работа газа?
Решение Степень свободы двухатомной молекулы , так как молекула может двигаться поступательно по трём координатным осям, и вращаться вокруг двух осей.


Последние материалы раздела:

Развитие критического мышления: технологии и методики
Развитие критического мышления: технологии и методики

Критическое мышление – это система суждений, способствующая анализу информации, ее собственной интерпретации, а также обоснованности...

Онлайн обучение профессии Программист 1С
Онлайн обучение профессии Программист 1С

В современном мире цифровых технологий профессия программиста остается одной из самых востребованных и перспективных. Особенно высок спрос на...

Пробный ЕГЭ по русскому языку
Пробный ЕГЭ по русскому языку

Здравствуйте! Уточните, пожалуйста, как верно оформлять подобные предложения с оборотом «Как пишет...» (двоеточие/запятая, кавычки/без,...