Стационарное уравнение диффузии в ограниченной среде. Теплопроводность; уравнение диффузии

Явлением диффузии называется процесс установления внутри фаз равновесного распределения концентраций.

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

Уравнение диффузии в одномерном случае

Уравнение диффузии в одномерном случае () в двухкомпонентной системе - это первый закон Фика:

где dm – масса первого компонента, которая переносится за время dt через элементарную площадку dS в направлении нормали x к рассматриваемой площадке в сторону убывания плотности первого компонента, – градиент плотности, D – коэффициент диффузии.

Если в однокомпонентной системе выделить группу молекул, выравнивание концентрации выделенных частиц по объёму сосуда называется самодиффузией. Самодиффузия тоже описывается уравнением диффузии (первым законом Фика), в котором коэффициент D- называется коэффициентом самодиффузии.

Уравнение диффузии в трехмерном случае

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

где D- коэффициент диффузии, t- время. Если D не зависит от концентрации, то уравнение диффузии будет иметь вид:

Уравнение (3) еще называют вторым законом Фика, где - дифференциальный оператор Лапласа.

В том случае, если перенос вещества вызван лишь градиентом его концентрации уравнение диффузии можно записать и в следующем виде:

где c(x, t) - концентрация вещества в точке среды в момент времени t, D – коэффициент диффузии, q - коэффициент поглощения, a F - интенсивность источников вещества. Величины D, q и F обычно являются функциями координат и времени, а также могут зависеть от концентрации с(x, t). B последнем случае, уравнение диффузии (4) становится нелинейным. В анизотропной среде коэффициент диффузии D является тензорным полем. В случае, когда величины D и q постоянны уравнение (4) является уравнением параболического типа. Для такого типа уравнений в математической физике разработаны методы решения. Допущение о постоянстве коэффициента диффузии справедливо в большинстве случаев реализуемых на практике. Уравнения диффузии не содержат ни каких сведений о механизмах этого процесса. Основная цель решения уравнения - найти распределение примеси c(x,t) после диффузии в течение определенного времени при различных условиях осуществления процесса.

Решение уравнения диффузии

Для выделения единственного решения для уравнения (4) необходимо задать начальные и граничные условия. Обычно, рассматривают следующие граничные условия:

1) на границе поверхности S поддерживается заданное распределение вещества

2)на границе поверхности S поддерживается заданная плотность потока вещества, входящего в V через границу S:

где n – внутренняя нормаль к поверхности S

3) S- полупроницаема, и диффузия во внешнюю среду с заданной концентрацией через поверхность S происходит по линейному закону:

В простейшем случае, когда диффузия происходит только вдоль одной прямой и c=c(x,t)уравнение (3) запишется в виде:

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

Тогда уравнение (5) имеет решение вида:

Текущая координата интегрирования.

Выражение (6) называется фундаментальным решением уравнения диффузии в случае (5).

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

ПРИМЕР 1

Задание Найти массу газа ( с молярной плотностью прошедшего вследствие диффузии через площадку за время , если градиент плотности в направлении, перпендикулярном площадке, равен . Температура газа T, средняя длинна свободного пробега молекулы .
Решение Запишем первый закон Фика в терминах условий задачи:

Знак минус означает, направление вектора плотности. Возьмем модуль от правой части выражения (1.1):

Зная, что , где - средняя длина свободного пробега молекулы, - средняя скорость молекулы газа и она равна: .

Соответственно преобразуем (1.2), найдем искомую массу газа:

Ответ Искомая масса газа может быть найдена по формуле:

· Ньютоновская жидкость · Неньютоновская жидкость · Поверхностное натяжение

См. также: Портал:Физика

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

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

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

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

  • Ближайшим формальным, а во многом и содержательным, аналогом уравнения диффузии является уравнение Шрёдингера , отличающееся от уравнения диффузии множителем мнимая единица перед производной по времени. Многие теоремы о решении уравнения Шрёдингера и даже некоторые виды формальной записи его решений прямо аналогичны соответствующим теоремам об уравнении диффузии и его решениях, однако качественно их решения различаются очень сильно.

Общий вид

Уравнение обычно записывается так:

\frac{\partial\varphi(\mathbf{r},t)}{\partial t} = \nabla \cdot \big[ D(\varphi,\mathbf{r}) \ \nabla\varphi(\mathbf{r},t) \big],

где φ(r , t ) - плотность диффундирующего вещества в точке r и во время t и D (φ, r ) - обобщённый диффузионный коэффициент для плотности φ в точке r ; ∇ - оператор набла . Если коэффициент диффузии зависит от плотности - уравнение нелинейно, в противном случае - линейно.

Если D - симметричный положительно определённый оператор , уравнение описывает анизотропную диффузию:

\frac{\partial\varphi(\mathbf{r},t)}{\partial t} = \sum_{i=1}^3\sum_{j=1}^3 \frac{\partial}{\partial x_i}\left

Если D постоянное, то уравнение сводится к линейному дифференциальному уравнению:

\frac{\partial\phi(\mathbf{r},t)}{\partial t} = D\nabla^2\phi(\mathbf{r},t),

История происхождения

Нестационарное уравнение

Нестационарное уравнение диффузии классифицируется как параболическое дифференциальное уравнение . Оно описывает распространение растворяемого вещества вследствие диффузии или перераспределение температуры тела в результате теплопроводности .

Одномерный случай

В случае одномерного диффузионного процесса с коэффициентом диффузии (теплопроводности) D уравнение имеет вид:

\frac{\partial}{\partial t}c(x,\;t)=\frac{\partial}{\partial x}D\frac{\partial}{\partial x}{c(x,\;t)}+f(x,\;t).

При постоянном D приобретает вид:

\frac{\partial}{\partial t}c(x,\;t)=D\frac{\partial^2}{\partial x^2}{c(x,\;t)}+f(x,\;t),

где c(x,\;t) - концентрация диффундирующего вещества, a f(x,\;t) - функция, описывающая источники вещества (тепла).

Трёхмерный случай

В трёхмерном случае уравнение приобретает вид:

\frac{\partial}{\partial t} c(\vec{r},\;t)=(\nabla,\;D\nabla c(\vec{r},\;t))+f(\vec{r},\;t),

где \nabla=(\partial_x,\;\partial_y,\;\partial_z) - оператор набла , а (\;,\;) - скалярное произведение. Оно также может быть записано как

\partial_t c=\mathbf{div}\,(D\,\mathbf{grad}\,c)+f,

а при постоянном D приобретает вид:

\frac{\partial}{\partial t} c(\vec{r},\;t)=D\Delta c(\vec{r},\;t)+f(\vec{r},\;t),

где \Delta=\nabla^2=\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2} - оператор Лапласа .

n -мерный случай

n-мерный случай - прямое обобщение приведенного выше, только под оператором набла, градиентом и дивергенцией, а также под оператором Лапласа надо понимать n-мерные версии соответствующих операторов:

\nabla=(\partial_1,\;\partial_2,\;\ldots,\;\partial_n), \Delta=\nabla^2=\partial_1^2+\partial_2^2+\ldots+\partial_n^2.

Это касается и двумерного случая n=2.

Мотивация

A.

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

\Phi=-\varkappa\frac{\partial c}{\partial x} (одномерный случай), \mathbf j=-\varkappa\nabla c (для любой размерности),

в сочетании с уравнением непрерывности, выражающим сохранение вещества (или энергии):

\frac{\partial c}{\partial t}+\frac{\partial\Phi}{\partial x}=0 (одномерный случай), \frac{\partial c}{\partial t}+\mathrm{div}\,\mathbf j=0 (для любой размерности),

с учетом в случае уравнения теплопроводности ещё теплоёмкости (температура = плотность энергия / удельная теплоемкость).

  • Здесь источник вещества (энергии) в правой части опущен, но он, конечно же, может быть легко туда помещён, если в задаче есть приток (отток) вещества (энергии).

B.

Кроме того, оно естественно возникает как непрерывный предел аналогичного разностного уравнения, возникающего в свою очередь при рассмотрении задачи о случайном блуждании на дискретной решётке (одномерной или n-мерной). (Это простейшая модель; в более сложных моделях случайных блужданий уравнение диффузии также возникает в непрерывном пределе). Простейшей интерпретацией функции c в этом случае служит количество (или концентрация) частиц в данной точке (или вблизи неё), причём каждая частица движется независимо от остальных без памяти (инерции) своего прошлого (в несколько более сложном случае - с ограниченной по времени памятью).

Решение

c(x,\;t)=\int\limits_{-\infty}^{+\infty}c(x",\;0)c_f(x-x",\;t)\,dx"=\int\limits_{-\infty}^{+\infty}c(x",\;0)\frac{1}{\sqrt{4\pi Dt}}\exp\left(-\frac{(x-x")^2}{4Dt}\right)\,dx".

Физические замечания

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

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

Стационарное уравнение

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

-(\nabla,\;D\nabla c(\vec{r}))=f(\vec{r}).

  • При D, не зависящем от \vec{r}, стационарное уравнение диффузии становится уравнением Пуассона (неоднородное), или уравнением Лапласа (однородное, то есть при f=0):
\Delta c(\vec{r})=-\frac{f(\vec{r})}{D}, \Delta c(\vec{r})=0.

Постановка краевых задач

  • Задача с начальными условиями (задача Коши) о распределении температуры на бесконечной прямой

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

и t\geqslant t_0, удовлетворяющее условию u(x,\;t_0)=\varphi(x)\quad(-\infty, где \varphi(x) - заданная функция.

  • Первая краевая задача для полубесконечного стержня

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

Найти решение уравнения теплопроводности в области -\infty\leqslant x\leqslant +\infty и t\geqslant t_0, удовлетворяющее условиям

\left\{\begin{array}{l}

u(x,\;t_0)=\varphi(x),\quad(0 где \varphi(x) и \mu(t) - заданные функции.

  • Краевая задача без начальных условий

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

Найти решение уравнения теплопроводности в области 0\leqslant x\leqslant l и -\infty, удовлетворяющее условиям

\left\{\begin{array}{l}

u(0,\;t)=\mu _1(t), \\ u(l,\;t)=\mu _2(t), \end{array}\right. где \mu_1(t) и \mu_2(t) - заданные функции.

  • Краевые задачи для ограниченного стержня

Рассмотрим следующую краевую задачу:

u_t=a^2 u_{xx}+f(x,\;t),\quad 0 - уравнение теплопроводности.

Если f(x,\;t)=0, то такое уравнение называют однородным , в противном случае - неоднородным .

u(x,\;0)=\varphi(x),\quad 0\leqslant x\leqslant l - начальное условие в момент времени t=0, температура в точке x задается функцией \varphi(x). \left.\begin{array}{l}

u(0,\;t)=\mu_1(t), \\ u(l,\;t)=\mu_2(t), \end{array}\right\}\quad 0\leqslant t\leqslant T - краевые условия. Функции \mu_1(t) и \mu_2(t) задают значение температуры в граничных точках 0 и l в любой момент времени t.

В зависимости от рода краевых условий, задачи для уравнения теплопроводности можно разбить на три типа. Рассмотрим общий случай (\alpha_i^2+\beta_i^2\ne 0,\;(i=1,\;2)).

\begin{array}{l}

\alpha_1 u_x(0,\;t)+\beta_1 u(0,\;t)=\mu_1(t), \\ \alpha_2 u_x(l,\;t)+\beta_2 u(l,\;t)=\mu_2(t). \end{array}

Если \alpha_i=0,\;(i=1,\;2), то такое условие называют условием первого рода , если \beta_i=0,\;(i=1,\;2) - второго рода , а если \alpha_i и \beta_i отличны от нуля, то условием третьего рода . Отсюда получаем задачи для уравнения теплопроводности - первую, вторую и третью краевую.

Принцип максимума

Пусть функция u(x,\;t) в пространстве D\times,\;D\in\R^n, удовлетворяет однородному уравнению теплопроводности \frac{\partial u}{\partial t}-a^2\Delta u=0, причем D - ограниченная область. Принцип максимума утверждает, что функция u(x,\;t) может принимать экстремальные значения либо в начальный момент времени, либо на границе области D.

Напишите отзыв о статье "Уравнение диффузии"

Примечания

Рассмотрим баланс нейтронов в единице объема dV при заданных Ф(r ), Ss.

Баланс нейтронов

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

Рождение – утечка – поглощение.

Рождение нейтронов обусловлено источником: S(r ) -число нейтронов, рождающихся в единицу времени в единице объема вблизи r . Поглощение нейтронов определяется числом реакций в единицу времени в единице объема . Нужно найти выход реакции в элементе объёма

.

Найдем утечку нейтронов, зная вектор плотности J из закона Фика

Если известенвектор J в каждой точке поверхности элементарного объема dV, то утечка равна divJ - число нейтронов, пересекающих поверхность единичного объема в единицу времени. Причем

div /D= const/= – D DФ

где

Таким образом, имеем уравнение

В стационарном случае

Замечания:

При выводе данных уравнений пользовались законом Фика, который справедлив, если распределение потока по координатам является линейным на расстоянии в несколько . Значит, эти уравнения плохо работают вблизи границы источника. Коэффициент D здесь уже учитывает возможную несферичность рассеяния(см. ранее).

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

1) поток Ф нейтронов конечен и неотрицателен в области, где применимо уравнение диффузии;

2) на границе двух сред, отличающихся хотя бы одной характеристикой взаимодействия нейтронов с ядрами.

Взаимодействие нейтронов с ядрами

В точке а :

Нормаль к поверхности;

Ток нейтронов.

Так как сама граница не поглощает нейтроны, то сколько нейтронов уходит из среды А, столько и приходит в среду В, т.е. проекции на нормаль

т. е. поток на границе неразрывен.

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

Итак, имеем условия на границе

Условия на границе

Условия на границе

3) на границе среды с вакуумом (это условие необходимо при решении задач о конечном реакторе) нет потока внутрь среды из вакуума. Это условие можно выразить, если задать функцию Ф(r, E, W). На границе имеем:

функция Ф(r, E, W).

Видно, что это граничное условие нельзя записать, зная только зависимость Ф от r . Используем следующий прием: изобразим Ф(r ) в плоском реакторе. Очевидно, поток на границе меньше, чем в центре активной зоны, но не равен 0, т.е. . Уравнение наиболее просто решается при нулевых граничных условиях.

Поток на границе

Решение уравнения диффузии особенно просто, когда на какой-либо границе поток равен 0. Будем считать, что поток образуется в 0 не на физической, а на некоторой экстраполированной границе реактора (экстраполяция линейная).

Длина экстраполяции d – величина неопределенная, но вносящая малую поправку в уравнение диффузии. Оценка d была сделана как теоретически, так и экспериментально. Оказалось, что при d = 0,71λ tr наблюдается наилучшее совпадение теории с опытом.

Изучение процессов диффузии велось также в направлении создания на основе экспериментальных результатов более точных моделей, которые давали бы возможность предсказывать протекание процесса диффузии путем теоретического анализа. Конечная цель исследования процесса диффузии – возможность расчетным путем определять электрические характеристики полупроводниковых приборов на основе технологических параметров процесса. Диффузионные модели развивались с позиции двух основных приближений: 1) теории сплошных сред с использованием уравнения диффузии Фика и 2) атомистической теории, которая принимает во внимание взаимодействие между точечными дефектами (вакансиями и межузельными атомами), с одной стороны, и примесными атомами – с другой. Теория сплошных сред описывает явление диффузии исходя из диффузионного уравнения Фика с учетом соответствующих коэффициентов диффузии. Коэффициенты диффузии легирующих элементов могут быть определены путем экспериментальных измерений поверхностной концентрации, глубины р п -перехода или профиля концентрации и из решения уравнения диффузии Фика.

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

В 1855 Фик предложил теорию диффузии. В основу этой Теории положена аналогия между процессами переноса в жидких растворах и тепла за счет теплопроводности. Фик предложил следующие уравнение, получившее название I закона Фика :

Для одномерного случая:

, (3.2)

Здесь: j – поток атомов диффундирующего вещества через единичную площадку (например, через см 2) за единицу времени (с), N – количество таких атомов в единице объема, t (с)– время диффузии, а D –коэффициент пропорциональности, связывающий j и grad N, имеющий размерность см 2 /с, D (см 2 /с)называют коэффициентом диффузии . Знак (-) отражает тот факт, что поток атомов идет в направлении уменьшения их концентрации. Диффузия идет всегда, но направленный поток имеет место только в случае неоднородного по пространству распределения диффундирующих частиц и исчезает, когда система становится однородной.

Выполнив дифференцирование по координате, в одномерном варианте получим:

(3.3)

Легко устанавливается факт, что

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

В частных производных. Бывает нестационарным и стационарным.

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

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

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

  • Ближайшим формальным, а во многом и содержательным, аналогом уравнения диффузии является уравнение Шрёдингера , отличающееся от уравнения диффузии множителем мнимая единица перед производной по времени. Многие теоремы о решении уравнения Шрёдингера и даже некоторые виды формальной записи его решений прямо аналогичны соответствующим теоремам об уравнении диффузии и его решениях, однако качественно их решения различаются очень сильно.

Общий вид

Уравнение обычно записывается так:

История происхождения

Одномерный случай

В случае одномерного диффузионного процесса с коэффициентом диффузии (теплопроводности) D уравнение имеет вид:

\frac{\partial}{\partial t}c(x,\;t)=\frac{\partial}{\partial x}D\frac{\partial}{\partial x}{c(x,\;t)}+f(x,\;t).

При постоянном D приобретает вид:

\frac{\partial}{\partial t}c(x,\;t)=D\frac{\partial^2}{\partial x^2}{c(x,\;t)}+f(x,\;t),

где c(x,\;t) - концентрация диффундирующего вещества, a f(x,\;t) - функция, описывающая источники вещества (тепла).

Трёхмерный случай

В трёхмерном случае уравнение приобретает вид:

\frac{\partial}{\partial t} c(\vec{r},\;t)=(\nabla,\;D\nabla c(\vec{r},\;t))+f(\vec{r},\;t),

где \nabla=(\partial_x,\;\partial_y,\;\partial_z) - оператор набла , а (\;,\;) - скалярное произведение. Оно также может быть записано как

\partial_t c=\mathbf{div}\,(D\,\mathbf{grad}\,c)+f,

а при постоянном D приобретает вид:

\frac{\partial}{\partial t} c(\vec{r},\;t)=D\Delta c(\vec{r},\;t)+f(\vec{r},\;t),

где \Delta=\nabla^2=\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2} - оператор Лапласа .

n -мерный случай

n-мерный случай - прямое обобщение приведенного выше, только под оператором набла, градиентом и дивергенцией, а также под оператором Лапласа надо понимать n-мерные версии соответствующих операторов:

\nabla=(\partial_1,\;\partial_2,\;\ldots,\;\partial_n), \Delta=\nabla^2=\partial_1^2+\partial_2^2+\ldots+\partial_n^2.

Это касается и двумерного случая n=2.

Мотивация

A.

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

\Phi=-\varkappa\frac{\partial c}{\partial x} (одномерный случай), \mathbf j=-\varkappa\nabla c (для любой размерности),

в сочетании с уравнением непрерывности, выражающим сохранение вещества (или энергии):

\frac{\partial c}{\partial t}+\frac{\partial\Phi}{\partial x}=0 (одномерный случай), \frac{\partial c}{\partial t}+\mathrm{div}\,\mathbf j=0 (для любой размерности),

с учетом в случае уравнения теплопроводности ещё теплоёмкости (температура = плотность энергия / удельная теплоемкость).

  • Здесь источник вещества (энергии) в правой части опущен, но он, конечно же, может быть легко туда помещён, если в задаче есть приток (отток) вещества (энергии).

B.

Кроме того, оно естественно возникает как непрерывный предел аналогичного разностного уравнения, возникающего в свою очередь при рассмотрении задачи о случайном блуждании на дискретной решётке (одномерной или n-мерной). (Это простейшая модель; в более сложных моделях случайных блужданий уравнение диффузии также возникает в непрерывном пределе). Простейшей интерпретацией функции c в этом случае служит количество (или концентрация) частиц в данной точке (или вблизи неё), причём каждая частица движется независимо от остальных без памяти (инерции) своего прошлого (в несколько более сложном случае - с ограниченной по времени памятью).

Решение

c(x,\;t)=\int\limits_{-\infty}^{+\infty}c(x",\;0)c_f(x-x",\;t)\,dx"=\int\limits_{-\infty}^{+\infty}c(x",\;0)\frac{1}{\sqrt{4\pi Dt}}\exp\left(-\frac{(x-x")^2}{4Dt}\right)\,dx".

Физические замечания

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

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

Стационарное уравнение

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

-(\nabla,\;D\nabla c(\vec{r}))=f(\vec{r}).

  • При D, не зависящем от \vec{r}, стационарное уравнение диффузии становится уравнением Пуассона (неоднородное), или уравнением Лапласа (однородное, то есть при f=0):
\Delta c(\vec{r})=-\frac{f(\vec{r})}{D}, \Delta c(\vec{r})=0.

Постановка краевых задач

  • Задача с начальными условиями (задача Коши) о распределении температуры на бесконечной прямой

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

и t\geqslant t_0, удовлетворяющее условию u(x,\;t_0)=\varphi(x)\quad(-\infty, где \varphi(x) - заданная функция.

  • Первая краевая задача для полубесконечного стержня

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

Найти решение уравнения теплопроводности в области -\infty\leqslant x\leqslant +\infty и t\geqslant t_0, удовлетворяющее условиям

\left\{\begin{array}{l}

u(x,\;t_0)=\varphi(x),\quad(0 где \varphi(x) и \mu(t) - заданные функции.

  • Краевая задача без начальных условий

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

Найти решение уравнения теплопроводности в области 0\leqslant x\leqslant l и -\infty, удовлетворяющее условиям

\left\{\begin{array}{l}

u(0,\;t)=\mu _1(t), \\ u(l,\;t)=\mu _2(t), \end{array}\right. где \mu_1(t) и \mu_2(t) - заданные функции.

  • Краевые задачи для ограниченного стержня

Рассмотрим следующую краевую задачу:

u_t=a^2 u_{xx}+f(x,\;t),\quad 0 - уравнение теплопроводности.

Если f(x,\;t)=0, то такое уравнение называют однородным , в противном случае - неоднородным .

u(x,\;0)=\varphi(x),\quad 0\leqslant x\leqslant l - начальное условие в момент времени t=0, температура в точке x задается функцией \varphi(x). \left.\begin{array}{l}

u(0,\;t)=\mu_1(t), \\ u(l,\;t)=\mu_2(t), \end{array}\right\}\quad 0\leqslant t\leqslant T - краевые условия. Функции \mu_1(t) и \mu_2(t) задают значение температуры в граничных точках 0 и l в любой момент времени t.

В зависимости от рода краевых условий, задачи для уравнения теплопроводности можно разбить на три типа. Рассмотрим общий случай (\alpha_i^2+\beta_i^2\ne 0,\;(i=1,\;2)).

\begin{array}{l}

\alpha_1 u_x(0,\;t)+\beta_1 u(0,\;t)=\mu_1(t), \\ \alpha_2 u_x(l,\;t)+\beta_2 u(l,\;t)=\mu_2(t). \end{array}

Если \alpha_i=0,\;(i=1,\;2), то такое условие называют условием первого рода , если \beta_i=0,\;(i=1,\;2) - второго рода , а если \alpha_i и \beta_i отличны от нуля, то условием третьего рода . Отсюда получаем задачи для уравнения теплопроводности - первую, вторую и третью краевую.

Принцип максимума

Пусть функция u(x,\;t) в пространстве D\times,\;D\in\R^n, удовлетворяет однородному уравнению теплопроводности \frac{\partial u}{\partial t}-a^2\Delta u=0, причем D - ограниченная область. Принцип максимума утверждает, что функция u(x,\;t) может принимать экстремальные значения либо в начальный момент времени, либо на границе области D.

{{#ifeq: Image:Wiki_letter_w.svg|none||Шаблон:!class ="ambox-image"Шаблон:! }}