Загальна транспортна теорема
Нехай $\mathbf F$ буде гладким векторним або скалярним полем в області $\mathcal R(t)$ чиї межі це $\mathcal S(t),$ і нехай $\mathbf W$ буде полем швидкостей часозалежного руху $\mathcal S(t).$ Тоді $$\frac{d}{dt}\int_{\mathcal R(t)}\mathbf F(x, t)\ dV = \int_{\mathcal R(t)}\frac{\partial \mathbf F}{\partial t}\ dV + \int_{\mathcal S(t)}\mathbf F\mathbf W \cdot \mathbf n\ dA.$$Транспортна теорема Рейнольдса
Нехай $\mathbf\Phi$ буде гладким векторним або скалярним полем, і припустимо, що $\mathcal R(t)$ це пакунок плину з поверхнею $\mathcal S(t),$ що мандрує зі швидкістю потоку $\mathbf U.$ Тоді $$\frac{D}{Dt}\int_{\mathcal R(t)}\mathbf\Phi\ dV =\int_{\mathcal R(t)}\frac{\partial\mathbf\Phi}{\partial t}\ dV +\int_{\mathcal S(t)}\mathbf\Phi \mathbf U \cdot \mathbf n\ dA.$$Виведення з Загальної транспортної теореми.
Якщо наша область рухається зі швидкістю потоку, то $\mathbf W = \mathbf U$ і загальна похідна перетворюється на матеріальну.1Рівняння неперервності для контрольного об'єму
Почнемо з розглядання фіксованої маси плину $m,$ що міститься в довільній області $\mathcal R(t).$ $$m=\int_{\mathcal R(t)}\rho\ dV.$$ $\mathcal R(t), m$ можуть змінюватись з часом, але маса повинна залишатись незмінною. Наприклад, можна розглянути повітряну кульку з теплим повітрям у холодному середовищі. $$\frac{dm}{dt}=\frac{d}{dt}\int_{\mathcal R(t)}\rho\ dV.$$ Застосувавши загальну транспортну теорему маємо: $$\int_{\mathcal R(t)}\frac{\partial \rho}{\partial t}\ dV + \int_{\mathcal S(t)}\rho \mathbf W\cdot \mathbf n\ dA = 0.$$ У плинових системах зазвичай зручно вважати, що поле швидкостей $\mathbf W$ було полем швидкостей плину, що відповідає розгляданню на локальному рівні нашої області як пакунок плину. $$\int_{\mathcal R(t)}\frac{\partial \rho}{\partial t}\ dV + \int_{\mathcal S(t)}\rho \mathbf U\cdot \mathbf n\ dA = 0.$$ Використавши формула Остроградського-Гауса ми позбуваємось двох типів інтегралів залишаючи лише інтеграл по об'єму: $$\int_{\mathcal R(t)}\frac{\partial \rho}{\partial t} + \nabla \cdot \rho \mathbf U\ dV = 0.$$ З того, що область $\mathcal R(t)$ було обрана довільно ми можемо вважати, що це діє і для наскільки завгодно малої області (в рамках нашої гіпотези неперервності), тобто інтеграл мусить бути нулем усюди в $\mathcal R(t)$. Якщо це не так ми можемо розділити $\mathcal R(t)$ на дві підобласті додатну і від'ємну і показити цим порушення рівності нулю. Отже, $$\begin{equation}\frac{\partial \rho}{\partial t} + \nabla \cdot \rho \mathbf U=0\label{eq:cont_eq_df}\end{equation}.$$ Це диференціальна форма рівняння неперервності.Якщо розписати оператор дивергенції, то маємо $$\begin{equation} \rho_t + (\rho u)_t + (\rho v)_t + (\rho w)_t = 0\label{eq:cont_eq_div_decomp}\end{equation}.$$ Тепер можна використати добутку для похідних $$\frac{D\rho}{Dt} + \rho\nabla\cdot \mathbf U = 0\label.$$ У випадку нестисного потоку, $\rho - $ стале, маємо $$\begin{equation}u_x + v_y + w_z = 0\label{eq:incompres_uvw}\end{equation}.$$
Аналіз у випадку контрольного об'єму
Оскільки $\ref{eq:cont_eq_df}$ виконується для будь-якого об'єму, то ми можемо записати його і для пакунку плину $$\int_{\mathcal R(t)} \frac{\partial \rho}{\partial t} + \nabla \cdot \rho \mathbf U\ dV = 0.$$ Тепер ми використаємо формулу Остроградського-Гауса в іншому напрямку $$\int_{\mathcal R(t)} \frac{\partial \rho}{\partial t}\ dV + \int_{\mathcal S(t)} \rho \mathbf U \cdot \mathbf n = 0.$$ Тепер Загальну транспортну теорему $$\frac{d}{dt} \int_{\mathcal R(t)}\rho\ dV + \int_{\mathcal S(t)}\rho(\mathbf U - \mathbf W)\cdot \mathbf n \ dA = 0.$$ Якщо $\mathcal S_e(t)$ - зона входів і виходів через які плин може втікати і витікати, то $$\frac{d}{dt} \int_{\mathcal R(t)}\rho\ dV + \int_{\mathcal S_e(t)}\rho(\mathbf U - \mathbf W)\cdot \mathbf n \ dA = 0.$$ Це рівняння називається рівняння неперервності для контрольного об'єму і виражає очевидний факт, що$\{$ | швидкість збільшення маси контрольного об'єму |
$\}$ | = | $\{$ | сумарний потік крізь поверхню |
$\}$ |
Спрощення
Якщо густина незмінна з часом, а контрольний об'єм нерухомий, тобто $\mathbf W = 0,$ тоді можна записати $$\begin{equation}\int_{\mathcal S_e(t)}\rho\mathbf U \cdot \mathbf n\ dA = 0\label{eq:cv_simple}\end{equation}.$$ Рівняння $\ref{eq:cv_simple}$ вимірюється в одиницях $\rho U A,$ тобто $\frac{M}{L^3}\frac{L}{T}L^2 = \frac{M}{T}.$ Ми називаємо це масовою ви́тратою $\dot m = \rho U A.$ У випадку часозалежної поверхні маємо $$\dot m = \int_{\mathcal S_e(t)} \rho(\mathbf U - \mathbf W)\cdot n\ dA.$$Рівняння імпульсу потоку плину
Тут нам знадобиться закон Ньютона, який ми запишемо в альтернативній формі $$F/V = \frac{d}{dt}(\rho u).$$ На відміну від випадку з рівнянням неперевності наша область одразу буде пакунком плину. Якщо ми використовуємо Ейлерову точку зору, то швидкість зміни вектора імпульсу $\frac{D}{Dt}\int_{\mathcal R(t)}\rho \mathbf U\ dV.$ Маємо основні два типи сил, які треба врахувати:- масова сила - $\int_{\mathcal R(t)}\mathbf F_B\ dV,$
- поверхнева сила - $\int_{\mathcal S(t)}\mathbf F_S\ dA.$
Ну що ж, нам залишились два кроки: перевести рівняння $\ref{eq:momentum_under_int}$ до диференціальної форми і розібратись із поверхневими силами. Зауважимо, що $\mathbf F_S$ - це вектор і це говорить, що має існувати матриця $\Lambda,$ така що $\mathbf F_S = \mathbf{\Lambda} \cdot \mathbf n.$ З того, що $\mathbf n$ - це чисто геометрична величина, нормаль до поверхні, то фізична сутність поверхневих сил мусить зберігатись в матриці $\mathbf{\Lambda}.$
Тепер ми готові записати базову форму диференціального рівняння імпульсу: $$\int_{\mathcal R(t)}\rho \frac{D\mathbf U}{Dt}\ dV = \int_{\mathcal R(t)}\mathbf F_B\ dV + \int_{\mathcal S(t)}\mathbf{\Lambda}\cdot\mathbf n\ dA,$$ Тому що $\mathcal R(t)$ - це довільний елемент плину і, отже, ми можемо обрати його як завгодно малим, то підінтегральний вираз також дорівнює нулю: $$\begin{equation}\rho \frac{D\mathbf U}{Dt} - \mathbf F_B - \nabla\cdot\mathbf{\Lambda} = 0\label{eq:momentum_dif_form}\end{equation}.$$ Це дуже загальний баланс імпульсу, що виконується для всіх точок будь-якого потоку плину (в межах гіпотези неперервності).Трактування поверхневих сил
Поверхневе напруження складніша величина. Причиною цьому є те, що ми не можемо говорити про напруження в точці не визначивши поверхню на яку це напруження діє. $$\Lambda = \begin{bmatrix} -p + \tau_{xx} & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & -p + \tau_{yy} & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & -p + \tau_{zz} \\ \end{bmatrix} = -p\mathbf I + \mathbf\tau.$$ $\mathbf\tau$ відомий як в'язкісний тензор напружень.Рівняння Нав'є-Стокса
$$ \begin{align*} \nabla\cdot\Lambda &= \big(-\frac{\partial p}{\partial x} + \frac{\partial\tau_{xx}}{\partial x} + \frac{\partial\tau_{yx}}{\partial y} + \frac{\partial\tau_{zx}}{\partial z}\big)\mathbf e_1 +\\ &= \big(-\frac{\partial p}{\partial y} + \frac{\partial\tau_{xy}}{\partial x} + \frac{\partial\tau_{yy}}{\partial y} + \frac{\partial\tau_{zy}}{\partial z}\big)\mathbf e_2 +\\ &= \big(-\frac{\partial p}{\partial z} + \frac{\partial\tau_{xz}}{\partial x} + \frac{\partial\tau_{yz}}{\partial y} + \frac{\partial\tau_{zz}}{\partial z}\big)\mathbf e_3. \end{align*} $$ Тепер ми можемо переписати $\ref{eq:momentum_dif_form}$ у вигляді $$\begin{equation}\rho\frac{Du}{Dt} = -\frac{\partial p}{\partial x} + \frac{\partial\tau_{xx}}{\partial x} + \frac{\partial\tau_{yx}}{\partial y} + \frac{\partial\tau_{zx}}{\partial z} + F_{B,x}\label{eq:momentum_df_tau}\end{equation}$$ Згадаємо багатовимірний закон в'язкості Ньютона, який надає такі формули для компонентів напруження в'язкості: $$\tau_{xx} = 2\mu u_x, \tau_yx = \mu(u_y+v_x), \tau_zx = \mu=\mu(u_z+w_x).$$ Підставляючи це у другий, третій і четвертий доданки $\ref{eq:momentum_df_tau}$ маємо $$2\mu u_{xx}+\mu(u_y+v_x)_y+\mu(u_z+w_x)_z.$$ Ми можемо переформатувати це так: $$\mu (u_{xx} + u_{yy} + u_{zz}) + \mu (u_x + v_y + w_z)_x,$$ але ми маємо справу з нестисним потоком, тому згідно з $\ref{eq:incompres_uvw}$ дивергенція вектора швидкості дорівнює нулю і залишається лише перший доданок. $$\rho\frac{Du}{Dt}=-p_x+\mu(u_{xx}+u_{yy}+u_{zz})+F_{B,x}.$$ Зазвичай це рівняння ділять на $\rho$ і представляти його як $$u_t + uu_x+vu_y+wu_z = -\frac{1}{\rho}p_x+\nu\Delta u+\frac{1}{\rho}F_{B,x}.$$ Тут, $\nu$ - це кінематична в'язкість.Рівняння для всіх трьох координат називаються рівняннями Нав'є-Стокса, і вони надають опис в кожній точці для по суті будь-якого нестисного потоку.
Аналіз рівнянь Нав'є-Стокса
$$\overbrace{\underbrace{u_t}_{\substack{локальне \\ прискорення}} + \underbrace{uu_x+vu_y+wu_z}_{\substack{конвекційне\\ прискорення}}}^{повне\ прискорення} = \overbrace{-\frac{1}{\rho}p_x}^{сили\ тиску}+ \overbrace{\nu\Delta u}^{сили\ в'язкості}+ \overbrace{\frac{1}{\rho}F_{B,x}}^{масові\ сили}.$$ Доданки ліворуч в контексті рівнянь Н.-С. часто називають інерційними.Сили в'язкості
$$\nu\Delta u = \nu(u_{xx}+u_{yy}+u_{zz}).$$ Загалом, доданки у вигляді другої похідної у диференційному рівнянні зазвичай пов'язують із дифузією, і в в обох фізичному і математичному контекстах представляють розмазування, згладжування або змішування.Ех, зібратись би колись і закінчити це
Примітки
- 1↩Stephen Whitaker. Fundamental Principles of Heat Transfer, стор. 197