Главная
Новости
Строительство
Ремонт
Дизайн и интерьер

















Яндекс.Метрика

Уравнение непрерывности

Уравнения непрерывности — (сильная) локальная форма законов сохранения. Ниже приведены примеры уравнений непрерывности, которые выражают одинаковую идею непрерывного изменения некоторой величины.

Дифференциальная форма

Дифференциальная форма общего уравнения непрерывности такова:

где

∇ ⋅ {displaystyle abla cdot } — дивергенция, ρ {displaystyle ho } — количество величины q {displaystyle q} на единицу объёма (плотность величины q {displaystyle q} ), t {displaystyle t} — время, j {displaystyle j} — плотность потока величины q {displaystyle q} (см. ниже), σ {displaystyle sigma } — добавление q {displaystyle q} на единицу объёма в единицу времени. Члены, которые добавляют ( σ > 0 {displaystyle sigma >0} ) или удаляют ( σ < 0 {displaystyle sigma <0} ) q {displaystyle q} , называются «источниками» и «стоками» соответственно.

Это общее уравнение может быть использовано для вывода любого уравнения непрерывности, начиная с простого уравнения неразрывности и до уравнения Навье — Стокса.

Если q {displaystyle q} — сохраняющаяся величина, которая не может быть создана или уничтожена (например, энергия), тогда σ = 0 {displaystyle sigma =0} , и уравнение непрерывности принимает вид

∂ ρ ∂ t + ∇ ⋅ j = 0. {displaystyle {frac {partial ho }{partial t}}+ abla cdot mathbf {j} =0.}

Электромагнетизм

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

div ⁡ j + ∂ ρ ∂ t = 0. {displaystyle operatorname {div} mathbf {j} +{partial ho over partial t}=0.}

Вывод

Закон Ампера гласит:

rot ⁡ H = j + ∂ D ∂ t . {displaystyle operatorname {rot} mathbf {H} =mathbf {j} +{frac {partial mathbf {D} }{partial t}}.}

Взяв дивергенцию от обеих частей выражения, получим

div ⁡ rot ⁡ H = div ⁡ j + ∂ ∂ t div ⁡ D , {displaystyle operatorname {div} operatorname {rot} mathbf {H} =operatorname {div} mathbf {j} +{frac {partial }{partial t}}operatorname {div} mathbf {D} ,}

но дивергенция ротора равняется нулю, таким образом

div ⁡ j + ∂ ∂ t div ⁡ D = 0. {displaystyle operatorname {div} mathbf {j} +{frac {partial }{partial t}}operatorname {div} mathbf {D} =0.}

По теореме Гаусса,

div ⁡ D = ρ . {displaystyle operatorname {div} mathbf {D} = ho .}

Подставляя это выражение в предыдущее уравнение, получаем искомое уравнение непрерывности.

Интерпретация

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

Теория волн

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

div ⁡ j + ∂ w ∂ t = 0 , {displaystyle operatorname {div} mathbf {j} +{frac {partial w}{partial t}}=0,}

где j = j ( x , y , z , t ) {displaystyle mathbf {j} =mathbf {j} (x,y,z,t)} — вектор плотности потока энергии в точке с координатами ( x , y , z ) {displaystyle left(x,y,z ight)} в момент времени t {displaystyle t} , w = w ( x , y , z , t ) {displaystyle w=w(x,y,z,t)} — плотность энергии.

Вывод

По определению, вектор плотности потока энергии — это вектор, модуль которого равен энергии, переносимой через единичную площадку, перпендикулярную направлению переноса энергии, за единицу времени, то есть j = d W d t d S ⊥ {displaystyle j={frac {dW}{dtdS_{ot }}}} , а направление его совпадает с направлением переноса энергии. Тогда энергия, вытекающая в единицу времени из некоторого макроскопического объёма V,

∮ S j d S = d W out d t . {displaystyle oint limits _{S}mathbf {j} ,dmathbf {S} ={frac {dW_{ ext{out}}}{dt}}.}

По закону сохранения энергии, d W out d t = − d W in d t {displaystyle {frac {dW_{ ext{out}}}{dt}}=-{frac {dW_{ ext{in}}}{dt}}} , где W in {displaystyle W_{ ext{in}}} — энергия, находящаяся в объёме V. По определению, плотность энергии — энергия единицы объёма, тогда полная энергия, заключённая в данном объёме, равна

W in = ∫ V w d V . {displaystyle W_{ ext{in}}=int limits _{V}w,dV.}

Тогда выражение для потока энергии примет вид

∮ S j d S = − d d t ∫ V w d V = − ∫ V ∂ w ∂ t d V . {displaystyle oint limits _{S}mathbf {j} ,dmathbf {S} =-{frac {d}{dt}}int limits _{V}w,dV=-int limits _{V}{frac {partial w}{partial t}},dV.}

Применяя формулу Гаусса — Остроградского к левой части выражения, получим

∫ V div ⁡ j d V = − ∫ V ∂ w ∂ t d V . {displaystyle int limits _{V}operatorname {div} mathbf {j} ,dV=-int limits _{V}{frac {partial w}{partial t}},dV.}

В силу произвольности выбранного объёма заключаем, что подынтегральные выражения равны, откуда и получаем дифференциальную форму уравнения непрерывности.

Гидродинамика и механика деформируемого твёрдого тела

Варианты названия

В гидродинамической литературе, например в работах Жуковского, Чаплыгина, Кочина, Лойцянского, уравнение, выражающее закон сохранения массы, называют уравнением неразрывности (условием неразрывности), тогда как в физической литературе, например в курсе Ландау и Лифшица, Зельдовича и Райзера, русском переводе курса Фейнмана, используется термин уравнение непрерывности. В старой литературе встречалось также название уравнение сплошности. Все три названия являются различными вариантами перевода введённого Эйлером названия уравнения в западноевропейских языках (англ. continuity equation, фр. équation de continuité и подобн.).

Различные формы записи

Уравнение выражает собой закон сохранения массы в элементарном объёме, то есть связь пространственного изменения потока массы жидкости или газа и скорости изменения плотности со временем. Его дифференциальная форма

∂ ρ ∂ t + div ⁡ ρ v = ∂ ρ ∂ t + ρ div ⁡ v + v grad ⁡ ρ = 0 , {displaystyle {frac {partial ho }{partial t}}+operatorname {div} ho mathbf {v} ={frac {partial ho }{partial t}}+ ho operatorname {div} mathbf {v} +mathbf {v} operatorname {grad} ho =0,}

где ρ = ρ ( x , y , z , t ) {displaystyle ho = ho (x,y,z,t)} — плотность жидкости (или газа), v = v ( x , y , z , t ) {displaystyle mathbf {v} =mathbf {v} (x,y,z,t)} — вектор скорости жидкости (или газа) в точке с координатами ( x , y , z ) {displaystyle (x,y,z)} в момент времени t {displaystyle t} .

Вектор j = ρ v {displaystyle mathbf {j} = ho mathbf {v} } называют плотностью потока жидкости. Его направление совпадает с направлением течения жидкости, а абсолютная величина определяет количество вещества, протекающего в единицу времени через единицу площади, расположенную перпендикулярно вектору скорости.

Для однородных несжимаемых жидкостей ρ = const {displaystyle ho ={ ext{const}}} . Поэтому уравнение принимает вид

div ⁡ v = 0 , {displaystyle operatorname {div} mathbf {v} =0,}

из чего следует соленоидальность поля скорости.

Для течений в каналах (течения в трубах, кровеносных сосудах и т. п.) уравнение неразрывности может быть записано в терминах средних значений по поперечному сечению канала. Например, для течения в канале с известной зависимостью площади поперечного сечения S {displaystyle S} от координаты x {displaystyle x} вдоль канала, S = S ( x ) {displaystyle S=S(x)} , (приближенное) уравнение неразрывности имеет вид

S ∂ ρ ∂ t + ∂ ∂ x ( ρ v S ) = 0 , {displaystyle S{frac {partial ho }{partial t}}+{frac {partial }{partial x}}( ho vS)=0,}

где ρ = ρ ( x , t ) {displaystyle ho = ho (x,t)} и v = v ( x , t ) {displaystyle v=v(x,t)} суть средние значения плотности и осевой проекции скорости по поперечному сечению. Здесь предполагается, что площадь поперечного сечения канала изменяется достаточно медленно (так называемое гидравлическое приближение), что позволяет при выводе уравнения заменять среднее значение от произведения на произведение от средних. В частном случае стационарного течения отсюда получается уравнение неразрывности в виде

ρ v S = const , {displaystyle ho vS={ ext{const}},}

имеющее очевидный физический смысл постоянства потока массы, а в случае среды с постоянной плотностью — уравнение

v S = const , {displaystyle vS={ ext{const}},}

выражающее постоянство объёмного расхода.

Аналогичную структуру имеет уравнение неразрывности для течений в каналах со свободной поверхностью, которое широко используется в гидравлике для описания русловых потоков (течения в реках, каналах и проч., движение селей, лавин и т. д.), для описания течений в плёнках и т. п. В простейшем случае течения жидкости с постоянной плотностью в канале с прямоугольным поперечным сечением точное уравнение неразрывности (иногда называемое уравнением Сен-Венана) имеет вид

∂ h ∂ t + ∂ ∂ x ( v h ) = 0 , {displaystyle {frac {partial h}{partial t}}+{frac {partial }{partial x}}(vh)=0,}

где h = h ( x , t ) {displaystyle h=h(x,t)} — глубина жидкости, v = v ( x , t ) {displaystyle v=v(x,t)} — средняя скорость жидкости по поперечному сечению.

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

ρ = ρ 0 ( 1 − div ⁡ w ) , {displaystyle ho = ho _{0}(1-operatorname {div} mathbf {w} ),}

где ρ 0 {displaystyle ho _{0}} , ρ {displaystyle ho } — соответственно начальная и конечная плотности материальной частицы, w {displaystyle mathbf {w} } — вектор перемещения (в случае малых перемещений и деформаций с одинаковой степенью точности можно брать дивергенцию как по эйлеровым переменным, так и по лагранжевым).

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

Историческая справка

В частных случаях, например для осесимметрических течений несжимаемой жидкости, уравнение неразрывности (в виде дифференциального уравнения в частных производных) было впервые получено Д’Аламбером, в общем виде — Эйлером в 1750-х годах. В форме алгебраического соотношения, выражающего (для случая несжимаемой жидкости) постоянство объёмного расхода вдоль трубки тока, уравнение неразрывности было впервые опубликовано Кастелли в первой половине XVII века.

Квантовая механика

В нерелятивистской квантовой механике сохранение вероятности также приводит к уравнению непрерывности. Пусть P ( x , t ) {displaystyle P(x,t)} — плотность вероятности, тогда уравнение запишется в виде

div ⁡ j + ∂ ∂ t P ( x , t ) = 0 , {displaystyle operatorname {div} mathbf {j} +{frac {partial }{partial t}}P(x,t)=0,}

где j {displaystyle j} — ток вероятности.