Классический теоретический вопрос: в распределении скоростей молекул идеального газа Максвелла возникают отклонения при сильных градиентах температуры. Какие поправки к классической теории необходимы и как они проявятся в переноса импульса и тепла
Коротко: при сильных градиентах температуры классическая Максвеллова форма распределения f(0)f^{(0)}f(0) нарушается — нужно учитывать нетепловые поправки, вытекающие из кинетического уравнения Больцмана. Практически это делается через расширение Чепмена–Энского или через метод Грэда (моменты). Последствия — переход от линейных законов Навье–Стокса/Фурье к нелинейным и ненлокальным соотношениям (Burnett, super‑Burnett, моменты), появление нормальных растяжений, нелинейного теплового потока, тепловой тяги (thermal creep) и слоев Кнудсена у границ. Ключевые формулы (схема): - Больцман и разложение ∂tf+v⋅∇f+F⋅∂f∂p=C[f],f=f(0)+f(1)+f(2)+… \partial_t f + \mathbf v\cdot\nabla f + \mathbf F\cdot\frac{\partial f}{\partial\mathbf p}=C[f],\qquad f=f^{(0)}+f^{(1)}+f^{(2)}+\dots ∂tf+v⋅∇f+F⋅∂p∂f=C[f],f=f(0)+f(1)+f(2)+…
- Локальный максвелловский нулевой порядок f(0)=n(m2πkBT)3/2exp (−mc22kBT),c=v−u. f^{(0)}=n\Big(\frac{m}{2\pi k_B T}\Big)^{3/2}\exp\!\Big(-\frac{m c^2}{2k_B T}\Big),\quad \mathbf c=\mathbf v-\mathbf u. f(0)=n(2πkBTm)3/2exp(−2kBTmc2),c=v−u.
- Первый порядок (Navier–Stokes–Fourier), линейные законы: Πij(1)=−η(∂iuj+∂jui−23δij∇⋅u),q(1)=−κ∇T. \Pi_{ij}^{(1)}=-\eta\Big(\partial_i u_j+\partial_j u_i-\tfrac{2}{3}\delta_{ij}\nabla\cdot\mathbf u\Big),\qquad \mathbf q^{(1)}=-\kappa\nabla T. Πij(1)=−η(∂iuj+∂jui−32δij∇⋅u),q(1)=−κ∇T.
Эти поправки предполагают малый Кнудсен Kn=ℓ/L≪1 \mathrm{Kn}=\ell/L\ll1 Kn=ℓ/L≪1 и малые градиенты. При сильных градиентах (Kn≪̸1 \mathrm{Kn}\not\ll1 Kn≪1 или большие относительные градиенты ℓ∣∇lnT∣∼1 \ell|\nabla\ln T|\sim1ℓ∣∇lnT∣∼1) необходимы высшие порядки: - Burnett / super‑Burnett: вторые и высшие по Kn \mathrm{Kn} Kn члены. Схематично Π=Π(1)+Π(2)+…,Π(2)∼α1ℓ2∇∇u+α2ℓ2(∇u)2+α3ℓ2∇∇T, \Pi=\Pi^{(1)}+\Pi^{(2)}+\dots,\qquad \Pi^{(2)}\sim \alpha_1 \ell^2 \nabla\nabla\mathbf u+\alpha_2 \ell^2 (\nabla\mathbf u)^2+\alpha_3\ell^2\nabla\nabla T, Π=Π(1)+Π(2)+…,Π(2)∼α1ℓ2∇∇u+α2ℓ2(∇u)2+α3ℓ2∇∇T,q=q(1)+q(2)+…,q(2)∼β1ℓ2∇∇T+β2ℓ2(∇u)⋅∇T+β3ℓ2∇(∇⋅u). \mathbf q=\mathbf q^{(1)}+\mathbf q^{(2)}+\dots,\qquad \mathbf q^{(2)}\sim \beta_1\ell^2\nabla\nabla T+\beta_2\ell^2(\nabla\mathbf u)\cdot\nabla T+\beta_3\ell^2\nabla(\nabla\cdot\mathbf u). q=q(1)+q(2)+…,q(2)∼β1ℓ2∇∇T+β2ℓ2(∇u)⋅∇T+β3ℓ2∇(∇⋅u).
Коэффициенты αi,βi\alpha_i,\beta_iαi,βi зависят от межмолекулярного потенциала и вычисляются из коллизионного оператора. Практические следствия: - Нелинейная зависимость вязкости и теплопроводности от градиентов (эффективные коэффициенты зависят от ∇T,∇u \nabla T,\nabla u∇T,∇u). - Появление нормальных растяжений (различные диагональные компоненты тензора напряжений). - Тепловой поток не только −κ∇T-\kappa\nabla T−κ∇T, но и вклад от градиентов скорости (взаимные эффекты), и второпорядковые пространственные производные (неместные эффекты). - Вблизи стен — слои Кнудсена с температурными скачками и упруго-термической тягой (thermal creep). - Для сильных неравновесий распределение становится анизотропным и приобретает высокоэнергетические отклонения, что требует решения полной кинетической задачи или применения многомоментных модельных уравнений (например, Грэд‑13), либо численных методов (DSMC, решение Больцмана). Итого: надо перейти от первого порядка (NSF) к высшим порядкам Чепмена–Энского (Burnett и далее) или к многомоментным/кинетическим решениям; в переносе импульса и тепла это проявится как дополнительные (второго и высших порядков по Kn \mathrm{Kn}Kn) члены, нелинейность и неместность в соотношениях для тензора напряжений и теплового потока.
Ключевые формулы (схема):
- Больцман и разложение
∂tf+v⋅∇f+F⋅∂f∂p=C[f],f=f(0)+f(1)+f(2)+… \partial_t f + \mathbf v\cdot\nabla f + \mathbf F\cdot\frac{\partial f}{\partial\mathbf p}=C[f],\qquad
f=f^{(0)}+f^{(1)}+f^{(2)}+\dots ∂t f+v⋅∇f+F⋅∂p∂f =C[f],f=f(0)+f(1)+f(2)+… - Локальный максвелловский нулевой порядок
f(0)=n(m2πkBT)3/2exp (−mc22kBT),c=v−u. f^{(0)}=n\Big(\frac{m}{2\pi k_B T}\Big)^{3/2}\exp\!\Big(-\frac{m c^2}{2k_B T}\Big),\quad \mathbf c=\mathbf v-\mathbf u. f(0)=n(2πkB Tm )3/2exp(−2kB Tmc2 ),c=v−u. - Первый порядок (Navier–Stokes–Fourier), линейные законы:
Πij(1)=−η(∂iuj+∂jui−23δij∇⋅u),q(1)=−κ∇T. \Pi_{ij}^{(1)}=-\eta\Big(\partial_i u_j+\partial_j u_i-\tfrac{2}{3}\delta_{ij}\nabla\cdot\mathbf u\Big),\qquad
\mathbf q^{(1)}=-\kappa\nabla T. Πij(1) =−η(∂i uj +∂j ui −32 δij ∇⋅u),q(1)=−κ∇T. Эти поправки предполагают малый Кнудсен Kn=ℓ/L≪1 \mathrm{Kn}=\ell/L\ll1 Kn=ℓ/L≪1 и малые градиенты.
При сильных градиентах (Kn≪̸1 \mathrm{Kn}\not\ll1 Kn≪1 или большие относительные градиенты ℓ∣∇lnT∣∼1 \ell|\nabla\ln T|\sim1ℓ∣∇lnT∣∼1) необходимы высшие порядки:
- Burnett / super‑Burnett: вторые и высшие по Kn \mathrm{Kn} Kn члены. Схематично
Π=Π(1)+Π(2)+…,Π(2)∼α1ℓ2∇∇u+α2ℓ2(∇u)2+α3ℓ2∇∇T, \Pi=\Pi^{(1)}+\Pi^{(2)}+\dots,\qquad
\Pi^{(2)}\sim \alpha_1 \ell^2 \nabla\nabla\mathbf u+\alpha_2 \ell^2 (\nabla\mathbf u)^2+\alpha_3\ell^2\nabla\nabla T, Π=Π(1)+Π(2)+…,Π(2)∼α1 ℓ2∇∇u+α2 ℓ2(∇u)2+α3 ℓ2∇∇T, q=q(1)+q(2)+…,q(2)∼β1ℓ2∇∇T+β2ℓ2(∇u)⋅∇T+β3ℓ2∇(∇⋅u). \mathbf q=\mathbf q^{(1)}+\mathbf q^{(2)}+\dots,\qquad
\mathbf q^{(2)}\sim \beta_1\ell^2\nabla\nabla T+\beta_2\ell^2(\nabla\mathbf u)\cdot\nabla T+\beta_3\ell^2\nabla(\nabla\cdot\mathbf u). q=q(1)+q(2)+…,q(2)∼β1 ℓ2∇∇T+β2 ℓ2(∇u)⋅∇T+β3 ℓ2∇(∇⋅u). Коэффициенты αi,βi\alpha_i,\beta_iαi ,βi зависят от межмолекулярного потенциала и вычисляются из коллизионного оператора.
Практические следствия:
- Нелинейная зависимость вязкости и теплопроводности от градиентов (эффективные коэффициенты зависят от ∇T,∇u \nabla T,\nabla u∇T,∇u).
- Появление нормальных растяжений (различные диагональные компоненты тензора напряжений).
- Тепловой поток не только −κ∇T-\kappa\nabla T−κ∇T, но и вклад от градиентов скорости (взаимные эффекты), и второпорядковые пространственные производные (неместные эффекты).
- Вблизи стен — слои Кнудсена с температурными скачками и упруго-термической тягой (thermal creep).
- Для сильных неравновесий распределение становится анизотропным и приобретает высокоэнергетические отклонения, что требует решения полной кинетической задачи или применения многомоментных модельных уравнений (например, Грэд‑13), либо численных методов (DSMC, решение Больцмана).
Итого: надо перейти от первого порядка (NSF) к высшим порядкам Чепмена–Энского (Burnett и далее) или к многомоментным/кинетическим решениям; в переносе импульса и тепла это проявится как дополнительные (второго и высших порядков по Kn \mathrm{Kn}Kn) члены, нелинейность и неместность в соотношениях для тензора напряжений и теплового потока.