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

18 Ноя в 10:18
5 +5
0
Ответы
1
Коротко: при сильных градиентах температуры классическая Максвеллова форма распределения 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
t f+vf+Fpf =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=vu. - Первый порядок (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=/L1 и малые градиенты.
При сильных градиентах (Kn≪̸1 \mathrm{Kn}\not\ll1 Kn1 или большие относительные градиенты ℓ∣∇ln⁡T∣∼1 \ell|\nabla\ln T|\sim1∣∇lnT1) необходимы высшие порядки:
- 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 uT,u).
- Появление нормальных растяжений (различные диагональные компоненты тензора напряжений).
- Тепловой поток не только −κ∇T-\kappa\nabla TκT, но и вклад от градиентов скорости (взаимные эффекты), и второпорядковые пространственные производные (неместные эффекты).
- Вблизи стен — слои Кнудсена с температурными скачками и упруго-термической тягой (thermal creep).
- Для сильных неравновесий распределение становится анизотропным и приобретает высокоэнергетические отклонения, что требует решения полной кинетической задачи или применения многомоментных модельных уравнений (например, Грэд‑13), либо численных методов (DSMC, решение Больцмана).
Итого: надо перейти от первого порядка (NSF) к высшим порядкам Чепмена–Энского (Burnett и далее) или к многомоментным/кинетическим решениям; в переносе импульса и тепла это проявится как дополнительные (второго и высших порядков по Kn \mathrm{Kn}Kn) члены, нелинейность и неместность в соотношениях для тензора напряжений и теплового потока.
18 Ноя в 10:50
Не можешь разобраться в этой теме?
Обратись за помощью к экспертам
Гарантированные бесплатные доработки в течение 1 года
Быстрое выполнение от 2 часов
Проверка работы на плагиат
Поможем написать учебную работу
Прямой эфир