[物理学与PDEs]第2章第2节 粘性流体力学方程组 2.3 广义 Newton 法则---本构方程

1.  ${\bf P}=(p_{ij})$, 而 $$\bex p_{ij}=-p\delta_{ij}+\tau_{ij}, \eex$$ 其中 $\tau_{ij}$ 对应于摩擦切应力.

 

 

2.  由于内摩擦力只与相对运动有关, 而 $\tau_{ij}$ 与速度无关, 而只与速度梯度有关, 且为线性的 (实验已很好的证实): $$\bex \tau_{ij}=c_{ijkl}\cfrac{\p u_k}{\p x_l}. \eex$$ 由于 $(\tau_{ij})$ 和 $\sex{\cfrac{\p u_k}{\p x_l}}$ 均为二阶张量, 而由张量识别定理, $(c_{ijkl})$ 为四阶张量. 又由 $p_{ij}$ 而 $\tau_{ij}$ 对称知 $$\bex c_{ijkl}=c_{jikl}. \eex$$

 

 

3.  设流体各向同性, 则 $c_{ijkl}$ 为各向同性张量, 有形式 $$\bex c_{ijkl}=\lm \delta_{ij}\delta_{kl} +\alpha \delta_{ik}\delta_{jl} +\beta\delta_{il}\delta_{jk}. \eex$$ 令 $\alpha=\mu+\nu$, $\beta=\mu-\nu$, 则由 $c_{ijkl}=c_{jikl}$ 知 $$\bex c_{ijkl}=\lm \delta_{ij}\delta_{kl} +\mu\sex{\delta_{ik}\delta_{jl}+\delta_{il}\delta_{jk}}. \eex$$ 由此, $c_{ijkl}=c_{ijlk}$, $$\bee\label{2_2_3_tau} \tau_{ij}=\lm \Div{\bf u}\delta_{ij}+2\mu s_{ij},\quad s_{ij}=\cfrac{1}{2}\sex{\cfrac{\p u_i}{\p x_j}+\cfrac{\p u_j}{\p x_i}}, \eee$$$$\bee\label{2_2_3_p} p_{ij}=(-p+\lm \Div{\bf u})\delta_{ij}+2\mu s_{ij}, \eee$$$$\bex {\bf P}=(-p+\lm \Div{\bf u}){\bf I}+2\mu {\bf S}. \eex$$

 

 

4.  $\lm$, $\mu$ 的物理意义

 

 

(1)  考虑沿 $x_1$ 方向的剪切运动 $$\bex u_1=u_1(x_3),\quad u_2=u_3=0.  \eex$$ 则由 \eqref{2_2_3_p}, $$\bex p_{13}=\mu \cfrac{\p u_1}{\p x_3}. \eex$$ 这就是 Newton 法则. 称 $\mu$ 为第一粘性系数 (动力学粘性系数).

 

 

(2)  由 \eqref{2_2_3_tau}, $$\bex \cfrac{1}{3}\sum_{i=1}^3 \tau_{ii} =\sex{\lm+\cfrac{2}{3}\mu }\Div{\bf u} =\sex{\lm+\cfrac{2}{3}\mu}\cfrac{1}{\tau}\cfrac{\rd \tau}{\rd t}\quad\sex{\tau=\cfrac{1}{\rho}:\mbox{ 比容}}. \eex$$ 记 $$\bex \mu'=\lm+\cfrac{2}{3}\mu, \eex$$ 则其为平均摩擦正应力与体积变化率之比, 描述流体运动过程中由膨胀或收缩引起的平均摩擦正应力的变换; 称为第二粘性系数 (膨胀粘性系数).

 

 

5.  总结:

 

 

(1)  应力张量---本构方程: $$\bex p_{ij}=-p\delta_{ij}+ 2\mu\sex{s_{ij}-\cfrac{1}{3}\Div{\bf u}\delta_{ij}} +\mu'\Div{\bf u} \delta_{ij}. \eex$$

 

 

(2)  广义Newton 法则: $$\bex \tau_{ij}=2\mu\sex{s_{ij}-\cfrac{1}{3}\Div{\bf u}\delta_{ij}} +\mu'\Div{\bf u} \delta_{ij}, \eex$$ 其中 $\mu>0$, $\mu'\geq 0$.

 

时间: 2024-10-08 16:55:48

[物理学与PDEs]第2章第2节 粘性流体力学方程组 2.3 广义 Newton 法则---本构方程的相关文章

[物理学与PDEs]第2章第2节 粘性流体力学方程组 2.1 引言

1.  实际的流体与理想流体的主要区别在于: 前者具有粘性 (内摩擦) 和热传导.     2.  内摩擦     (1)  当两层流体有相对运动时, 方有摩擦力; 它是一种内力; 单位面积上所受的内力称为应力; 而它通常与表面相切, 而称为切应力.     (2)  Newton 假设摩擦力与速度梯度成正比; 满足此假设的称为 Newton 流体; 而不满足的称为非 Newton 流体.     3.  热传导     (1)  Fourier 实验定律: $$\bex \rd q=-\kap

[物理学与PDEs]第2章第2节 粘性流体力学方程组 2.4 粘性热传导流体动力学方程组

粘性热传导流体动力学方程组: $$\beex \bea \cfrac{\p \rho}{\p t}+\Div(\rho{\bf u})&=0,\\ \rho \cfrac{\rd {\bf u}}{\rd t} +\n p -\n\sez{ \sex{\mu'-\cfrac{2}{3}\mu}\Div{\bf u} } -2\Div(\mu {\bf S})&=\rho {\bf F},\\ \rho\cfrac{\rd e}{\rd t} +p\Div{\bf u} -\mu\sum_{

[物理学与PDEs]第2章第2节 粘性流体力学方程组 2.6 一维粘性热传导流体动力学方程组

一维粘性热传导流体动力学方程组: $$\beex \bea \cfrac{\p\rho}{\p t}+\cfrac{\p }{\p x}(\rho u)&=0,\\ \cfrac{\p u}{\p t}+u\cfrac{\p u}{\p x} +\cfrac{1}{\rho}\cfrac{\p p}{\p x} -\cfrac{1}{\rho}\cfrac{\p }{\p x}\sez{\sex{\cfrac{4\mu}{3}+\mu'}\cfrac{\p u}{\p x}}&=F,\\

[物理学与PDEs]第2章第1节 理想流体力学方程组 1.2 理想流体力学方程组

1.  质量守恒定律: 连续性方程 $$\bee\label{2_1_2_zl} \cfrac{\p\rho}{\p t}+\Div(\rho{\bf u})=0.  \eee$$   2.  动量守恒定律: $$\bee\label{2_1_2_dl} \cfrac{\p}{\p t}(\rho{\bf u})+\Div(\rho{\bf u}\otimes {\bf u}+p{\bf I})=\rho{\bf F}. \eee$$ 用 \eqref{2_1_2_zl} 可化简 \eqref{

[物理学与PDEs]第2章第1节 理想流体力学方程组 1.3 理想流体力学方程组的数学结构

1.  局部音速 $c$: $c^2=\cfrac{\p p}{\p \rho}>0$.     2.  将理想流体力学方程组 $$\beex \bea \rho\cfrac{\p {\bf u}}{\p t} +(\rho {\bf u}\cdot\n){\bf u}+\n p&=\rho{\bf F},\\ \cfrac{1}{\rho c^2}\cfrac{\p p}{\p t} +\n\cdot{\bf u}+\cfrac{1}{\rho c^2}({\bf u}\cdot\n)p&

[物理学与PDEs]第2章第5节 一维流体力学方程组的 Lagrange 形式 5.4 一维粘性热传导流体力学方程组的 Lagrange 形式

1. 一维粘性热传导流体力学方程组 $$\beex \bea \cfrac{\p\rho}{\p t}+\cfrac{\p }{\p x}(\rho u)&=0,\\ \cfrac{\p u}{\p t}+u\cfrac{\p u}{\p x} +\cfrac{1}{\rho}\cfrac{\p p}{\p x} -\cfrac{1}{\rho}\cfrac{\p }{\p x}\sez{\sex{\cfrac{4\mu}{3}+\mu'}\cfrac{\p u}{\p x}}&=F,\\

[物理学与PDEs]第2章第1节 理想流体力学方程组 1.1 预备知识

1.  理想流体: 指忽略粘性及热传导的流体.   2.  流体的状态 (运动状态及热力学状态) 的描述   (1)   速度向量 $\bbu=(u_1,u_2,u_3)$: 流体微元的宏观运动速度.   (2)   质量密度 $\rho$: 单位体积流体的质量. a.  质量流向量 (动量密度向量) $\rho\bbu$; b.  动量流张量 $\rho \bbu\otimes \bbu$; c.  比容 $\tau=\cfrac{1}{\rho}$: 单位质量流体的体积.   (3)   压

[物理学与PDEs]第2章第1节 理想流体力学方程组 1.4 一维理想流体力学方程组

1.  一维理想流体力学方程组 $$\beex \bea \cfrac{\p\rho}{\p t}+\cfrac{\p}{\p x}(\rho u)&=0,\\ \cfrac{\p}{\p t}(\rho u) +\cfrac{\p}{\p x}(\rho u^2+p)&=\rho F,\\ \cfrac{\p}{\p t}\sex{\rho e+\cfrac{1}{2}\rho u^2} +\cfrac{\p}{\p x}\sez{\sex{ \rho e+\cfrac{1}{2}\rh

[物理学与PDEs]第4章第2节 反应流体力学方程组 2.3 混合气体状态方程

1.  记号与假设   (1)  已燃气体的化学能为 $0$.   (2)  单位质量的未燃气体的化学能为 $g_0>0$.     2.  对多方气体 (理想气体当 $T$ 不高时可近似认为), $$\bex p=(\gamma-1)e^\frac{S-S_0}{c_V}\rho^\gamma,\quad e=e^\frac{S-S_0}{c_V}\rho^{\gamma-1}\ra p=(\gamma-1)\rho e =(\gamma-1)\rho (E-Zg_0). \eex$$