Это не официальный сайт wikipedia.org 01.01.2023

Уравнение диффузии — Википедия

Уравнение диффузии

(перенаправлено с «Теплопроводности уравнение»)

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

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

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

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

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


Общий видПравить

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

φ ( r , t ) t = [ D ( φ , r )   φ ( r , t ) ] ,  

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

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

φ ( r , t ) t = i = 1 3 j = 1 3 x i [ D i j ( φ , r ) φ ( r , t ) x j ] .  

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

ϕ ( r , t ) t = D 2 ϕ ( r , t ) ,  

также называемому уравнением теплопроводности.

История происхожденияПравить

Дифференциальное уравнение в частных производных было первоначально выведено Адольфом Фиком в 1855 году.[1]

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

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

Одномерный случайПравить

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

t c ( x , t ) = x D x c ( x , t ) + f ( x , t ) .  

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

t c ( x , t ) = D 2 x 2 c ( x , t ) + f ( x , t ) ,  

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

Трёхмерный случайПравить

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

t c ( r , t ) = ( , D c ( r , t ) ) + f ( r , t ) ,  

где = ( x , y , z )   — оператор набла, а ( , )   — скалярное произведение. Оно также может быть записано как

t c = d i v ( D g r a d c ) + f ,  

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

t c ( r , t ) = D Δ c ( r , t ) + f ( r , t ) ,  

где Δ = 2 = 2 x 2 + 2 y 2 + 2 z 2   — оператор Лапласа.

n-мерный случайПравить

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

= ( 1 , 2 , , n ) ,  
Δ = 2 = 1 2 + 2 2 + + n 2 .  

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

МотивацияПравить

A.Править

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

Φ = ϰ c x   (одномерный случай),
j = ϰ c   (для любой размерности),

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

c t + Φ x = 0   (одномерный случай),
c t + d i v j = 0   (для любой размерности),

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

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

B.Править

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

РешениеПравить

В одномерном случае фундаментальное решение однородного уравнения с постоянным — не зависящим от x   и t   — D   (при начальном условии, выражаемом дельта-функцией c f ( x , 0 ) = δ ( x )   и граничном условии c f ( , t ) = 0  ) есть

c f ( x , t ) = 1 4 π D t exp ( x 2 4 D t ) .  

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

x 2 = + x 2 c f ( x , t ) d x = 2 D t .  


В случае произвольного начального распределения c ( x , 0 )   общее решение уравнения диффузии представляется в интегральном виде как свёртка:

c ( x , t ) = + c ( x , 0 ) c f ( x x , t ) d x = + c ( x , 0 ) 1 4 π D t exp ( ( x x ) 2 4 D t ) d x .  

Физические замечанияПравить

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

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

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

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

( , D c ( r ) ) = f ( r ) .  
Δ c ( r ) = f ( r ) D ,  
Δ c ( r ) = 0.  

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

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

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

Найти решение уравнения теплопроводности в области x +   и t t 0  , удовлетворяющее условию u ( x , t 0 ) = φ ( x ) ( < x < + )  , где φ ( x )   — заданная функция.

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

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

Найти решение уравнения теплопроводности в области x +   и t t 0  , удовлетворяющее условиям

{ u ( x , t 0 ) = φ ( x ) , ( 0 < x < ) u ( 0 , t ) = μ ( t ) , ( t t 0 )  

где φ ( x )   и μ ( t )   — заданные функции.

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

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

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

{ u ( 0 , t ) = μ 1 ( t ) , u ( l , t ) = μ 2 ( t ) ,  

где μ 1 ( t )   и μ 2 ( t )   — заданные функции.

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

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

u t = a 2 u x x + f ( x , t ) , 0 < x < l , 0 < t T   — уравнение теплопроводности.

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

u ( x , 0 ) = φ ( x ) , 0 x l   — начальное условие в момент времени t = 0  , температура в точке x   задается функцией φ ( x )  .
u ( 0 , t ) = μ 1 ( t ) , u ( l , t ) = μ 2 ( t ) , } 0 t T   — краевые условия. Функции μ 1 ( t )   и μ 2 ( t )   задают значение температуры в граничных точках 0 и l   в любой момент времени t  .

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

α 1 u x ( 0 , t ) + β 1 u ( 0 , t ) = μ 1 ( t ) , α 2 u x ( l , t ) + β 2 u ( l , t ) = μ 2 ( t ) .  

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

Принцип максимумаПравить

Пусть функция u ( x , t )   в пространстве D × [ 0 , T ] , D R n  , удовлетворяет однородному уравнению теплопроводности u t a 2 Δ u = 0  , причем D   — ограниченная область. Принцип максимума утверждает, что функция u ( x , t )   может принимать экстремальные значения либо в начальный момент времени, либо на границе области D  .

ПримечанияПравить

  1. Fick A., Ueber Diffusion, Pogg. Ann. Phys. Chem.— 1855.— 170 (4. Reihe 94).— pp. 59-86.