文章目录
问题的由来
在splishsplash论文教程当中,在讲到粘性的时候,有这样一页PPT
这是流体力学的基础知识,即利用牛顿本构模型来构建的NS方程。
我们先看左边这个公式
被称为柯西动量方程。这个方程对流体力学和固体力学全都适用。因此可以视为NS方程更基础的存在。它其实就是一个线动量方程。左边的加速度,右边是内力(由应力表示,为面力)和外力(一般是体力)。
其中T是应力张量。
我们再看右边这个公式。
这个公式就是牛顿本构方程。针对于牛顿流体所构建出来的方程。如果不是牛顿流体,就不适用。所谓的本构方程,就是找到应力与应变的关系。实际上,按照固体力学的说法,第二个方程是一个几何方程。第一个方程是物理方程,或者叫本构方程。我们这里暂且不区分,都叫他本构方程。
但是,我们的重点并不在于方程本身。我们的目的是让读者熟悉这些数学记号的含义。也就是说,我们从张量数学的角度来看看这些量写成我们熟悉的矩阵形式是什么样的?
其他的都没有什么难度。唯一我们不熟悉的就是速度梯度了。即
∇
v
\nabla \mathbf{v}
∇v
我们之前写过一篇关于张量的通俗入门,如下:
https://blog.csdn.net/weixin_43940314/article/details/123559800
速度梯度
速度显然是个向量
v
=
(
u
,
v
,
w
)
T
\mathbf{v} = (u, v, w)^T
v=(u,v,w)T
而nabla算子也只是个向量
∇
=
(
∂
∂
x
,
∂
∂
y
,
∂
∂
z
)
T
\nabla = (\frac{\partial}{\partial x}, \frac{\partial}{\partial y}, \frac{\partial}{\partial z})^T
∇=(∂x∂,∂y∂,∂z∂)T
向量与向量写到一块,中间没有任何符号,实际上是并乘(也叫并矢)。有时候我们也主动的添加一个符号,写成 ⨂ \bigotimes ⨂。
其实,梯度就是nabla算子与物理量的并乘。
并乘遵循什么规则呢?
并乘和点乘、叉乘不同,它什么也不做!只是单纯的把矢量并排着写在一起。
此外,记住:并乘是升阶的,结果的阶数就是把原来的阶数加起来!
下面我们把并乘写成熟悉的矩阵形式。
我们一行一行地写
∇
⊗
v
=
(
∂
∂
x
∂
∂
y
∂
∂
z
)
⊗
(
u
v
w
)
\nabla \otimes \mathbf{v} = \begin{pmatrix} \frac{\partial}{\partial x} \\ \frac{\partial}{\partial y}\\ \frac{\partial}{\partial z} \end{pmatrix} \otimes \begin{pmatrix} u \\ v\\ w \end{pmatrix}
∇⊗v=⎝⎛∂x∂∂y∂∂z∂⎠⎞⊗⎝⎛uvw⎠⎞
1 第一行
对应的项简单地并排写在一起
1.1 分量和形式
写成分量和的形式
∂
u
∂
x
i
i
+
∂
v
∂
x
i
j
+
∂
w
∂
x
i
k
\frac{\partial u}{\partial x} \mathbf{ii} + \frac{\partial v}{\partial x }\mathbf{ij} + \frac{\partial w}{\partial x} \mathbf{ik}
∂x∂uii+∂x∂vij+∂x∂wik
1.2 矩阵形式
或者写成矩阵的形式
也就是
( ∂ u ∂ x ∂ v ∂ x ∂ w ∂ x ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ) \begin{pmatrix} \frac{\partial u}{\partial x} & \frac{\partial v}{\partial x } & \frac{\partial w}{\partial x} \\ \cdots & \cdots & \cdots\\ \cdots & \cdots & \cdots\\ \end{pmatrix} ⎝⎛∂x∂u⋯⋯∂x∂v⋯⋯∂x∂w⋯⋯⎠⎞
2 第二行
同理
2.1 分量和形式
写成分量和的形式
∂
u
∂
y
i
i
+
∂
v
∂
y
i
j
+
∂
w
∂
y
i
k
\frac{\partial u}{\partial y} \mathbf{ii} + \frac{\partial v}{\partial y }\mathbf{ij} + \frac{\partial w}{\partial y} \mathbf{ik}
∂y∂uii+∂y∂vij+∂y∂wik
2.2 矩阵形式
或者写成矩阵的形式
也就是
( ⋯ ⋯ ⋯ ∂ u ∂ y ∂ v ∂ y ∂ w ∂ y ⋯ ⋯ ⋯ ) \begin{pmatrix} \cdots & \cdots & \cdots\\ \frac{\partial u}{\partial y} & \frac{\partial v}{\partial y} & \frac{\partial w}{\partial y}\\ \cdots & \cdots & \cdots\\ \end{pmatrix} ⎝⎛⋯∂y∂u⋯⋯∂y∂v⋯⋯∂y∂w⋯⎠⎞
3 第三行
同理
3.1 分量和形式
写成分量和的形式
∂
u
∂
z
i
i
+
∂
v
∂
z
i
j
+
∂
w
∂
z
i
k
\frac{\partial u}{\partial z} \mathbf{ii} + \frac{\partial v}{\partial z }\mathbf{ij} + \frac{\partial w}{\partial z} \mathbf{ik}
∂z∂uii+∂z∂vij+∂z∂wik
3.2 矩阵形式
或者写成矩阵的形式
也就是
( ⋯ ⋯ ⋯ ⋯ ⋯ ⋯ ∂ u ∂ z ∂ v ∂ z ∂ w ∂ z ) \begin{pmatrix} \cdots & \cdots & \cdots\\ \cdots & \cdots & \cdots\\ \frac{\partial u}{\partial z} & \frac{\partial v}{\partial z} & \frac{\partial w}{\partial z}\\ \end{pmatrix} ⎝⎛⋯⋯∂z∂u⋯⋯∂z∂v⋯⋯∂z∂w⎠⎞
4 合在一起
分量和形式
∇
v
=
∂
u
∂
x
i
i
+
∂
v
∂
x
i
j
+
∂
w
∂
x
i
k
+
∂
u
∂
y
i
i
+
∂
v
∂
y
i
j
+
∂
w
∂
y
i
k
+
∂
u
∂
z
i
i
+
∂
v
∂
z
i
j
+
∂
w
∂
z
i
k
\nabla \mathbf{v} = \\ \frac{\partial u}{\partial x} \mathbf{ii} + \frac{\partial v}{\partial x }\mathbf{ij} + \frac{\partial w}{\partial x} \mathbf{ik} + \\ \frac{\partial u}{\partial y} \mathbf{ii} + \frac{\partial v}{\partial y }\mathbf{ij} + \frac{\partial w}{\partial y} \mathbf{ik} + \\ \frac{\partial u}{\partial z} \mathbf{ii} + \frac{\partial v}{\partial z }\mathbf{ij} + \frac{\partial w}{\partial z} \mathbf{ik}
∇v=∂x∂uii+∂x∂vij+∂x∂wik+∂y∂uii+∂y∂vij+∂y∂wik+∂z∂uii+∂z∂vij+∂z∂wik
矩阵形式
( ∂ u ∂ x ∂ v ∂ x ∂ w ∂ x ∂ u ∂ y ∂ v ∂ y ∂ w ∂ y ∂ u ∂ z ∂ v ∂ z ∂ w ∂ z ) \begin{pmatrix} \frac{\partial u}{\partial x} & \frac{\partial v}{\partial x } & \frac{\partial w}{\partial x} \\ \frac{\partial u}{\partial y} & \frac{\partial v}{\partial y} & \frac{\partial w}{\partial y}\\ \frac{\partial u}{\partial z} & \frac{\partial v}{\partial z} & \frac{\partial w}{\partial z}\\ \end{pmatrix} ⎝⎛∂x∂u∂y∂u∂z∂u∂x∂v∂y∂v∂z∂v∂x∂w∂y∂w∂z∂w⎠⎞
对称张量
我们写出了速度梯度,现在来看一下对称张量
∇
v
+
(
∇
v
)
T
\nabla \mathbf{v} + (\nabla \mathbf{v})^T
∇v+(∇v)T
不过就是速度梯度加上它自己的转置而已。
转置就是行写成列,列写成行。
我们先写出转置
(
∇
v
)
T
=
(
∂
u
∂
x
∂
u
∂
y
∂
u
∂
z
∂
v
∂
x
∂
v
∂
y
∂
v
∂
z
∂
w
∂
x
∂
w
∂
y
∂
w
∂
z
)
(\nabla \mathbf{v})^T= \begin{pmatrix} \frac{\partial u}{\partial x} & \frac{\partial u}{\partial y } & \frac{\partial u}{\partial z} \\ \frac{\partial v}{\partial x} & \frac{\partial v}{\partial y} & \frac{\partial v}{\partial z}\\ \frac{\partial w}{\partial x} & \frac{\partial w}{\partial y} & \frac{\partial w}{\partial z}\\ \end{pmatrix}
(∇v)T=⎝⎜⎛∂x∂u∂x∂v∂x∂w∂y∂u∂y∂v∂y∂w∂z∂u∂z∂v∂z∂w⎠⎟⎞
然后再把它们加起来
( ∂ u ∂ x ∂ v ∂ x ∂ w ∂ x ∂ u ∂ y ∂ v ∂ y ∂ w ∂ y ∂ u ∂ z ∂ v ∂ z ∂ w ∂ z ) + ( ∂ u ∂ x ∂ u ∂ y ∂ u ∂ z ∂ v ∂ x ∂ v ∂ y ∂ v ∂ z ∂ w ∂ x ∂ w ∂ y ∂ w ∂ z ) = ( 2 ∂ u ∂ x ∂ v ∂ x + ∂ u ∂ y ∂ w ∂ x + ∂ u ∂ z ∂ u ∂ y + ∂ v ∂ x 2 ∂ v ∂ y ∂ w ∂ y + ∂ v ∂ z ∂ u ∂ z + ∂ w ∂ x ∂ v ∂ z + ∂ w ∂ y 2 ∂ w ∂ z ) \begin{pmatrix} \frac{\partial u}{\partial x} & \frac{\partial v}{\partial x } & \frac{\partial w}{\partial x} \\ \frac{\partial u}{\partial y} & \frac{\partial v}{\partial y} & \frac{\partial w}{\partial y}\\ \frac{\partial u}{\partial z} & \frac{\partial v}{\partial z} & \frac{\partial w}{\partial z}\\ \end{pmatrix} + \begin{pmatrix} \frac{\partial u}{\partial x} & \frac{\partial u}{\partial y } & \frac{\partial u}{\partial z} \\ \frac{\partial v}{\partial x} & \frac{\partial v}{\partial y} & \frac{\partial v}{\partial z}\\ \frac{\partial w}{\partial x} & \frac{\partial w}{\partial y} & \frac{\partial w}{\partial z}\\ \end{pmatrix} =\\ \begin{pmatrix} 2\frac{\partial u}{\partial x} &\frac{\partial v}{\partial x }+\frac{\partial u}{\partial y } & \frac{\partial w}{\partial x}+\frac{\partial u}{\partial z} \\ \frac{\partial u}{\partial y}+\frac{\partial v}{\partial x} & 2\frac{\partial v}{\partial y} & \frac{\partial w}{\partial y}+\frac{\partial v}{\partial z}\\ \frac{\partial u}{\partial z}+\frac{\partial w}{\partial x} & \frac{\partial v}{\partial z}+\frac{\partial w}{\partial y} & 2 \frac{\partial w}{\partial z}\\ \end{pmatrix} ⎝⎛∂x∂u∂y∂u∂z∂u∂x∂v∂y∂v∂z∂v∂x∂w∂y∂w∂z∂w⎠⎞+⎝⎜⎛∂x∂u∂x∂v∂x∂w∂y∂u∂y∂v∂y∂w∂z∂u∂z∂v∂z∂w⎠⎟⎞=⎝⎜⎛2∂x∂u∂y∂u+∂x∂v∂z∂u+∂x∂w∂x∂v+∂y∂u2∂y∂v∂z∂v+∂y∂w∂x∂w+∂z∂u∂y∂w+∂z∂v2∂z∂w⎠⎟⎞
我们观察这个矩阵
∇
v
+
(
∇
v
)
T
=
(
2
∂
u
∂
x
∂
v
∂
x
+
∂
u
∂
y
∂
w
∂
x
+
∂
u
∂
z
∂
u
∂
y
+
∂
v
∂
x
2
∂
v
∂
y
∂
w
∂
y
+
∂
v
∂
z
∂
u
∂
z
+
∂
w
∂
x
∂
v
∂
z
+
∂
w
∂
y
2
∂
w
∂
z
)
\nabla \mathbf{v} + (\nabla \mathbf{v})^T= \begin{pmatrix} 2\frac{\partial u}{\partial x} &\frac{\partial v}{\partial x }+\frac{\partial u}{\partial y } & \frac{\partial w}{\partial x}+\frac{\partial u}{\partial z} \\ \frac{\partial u}{\partial y}+\frac{\partial v}{\partial x} & 2\frac{\partial v}{\partial y} & \frac{\partial w}{\partial y}+\frac{\partial v}{\partial z}\\ \frac{\partial u}{\partial z}+\frac{\partial w}{\partial x} & \frac{\partial v}{\partial z}+\frac{\partial w}{\partial y} & 2 \frac{\partial w}{\partial z}\\ \end{pmatrix}
∇v+(∇v)T=⎝⎜⎛2∂x∂u∂y∂u+∂x∂v∂z∂u+∂x∂w∂x∂v+∂y∂u2∂y∂v∂z∂v+∂y∂w∂x∂w+∂z∂u∂y∂w+∂z∂v2∂z∂w⎠⎟⎞
可以很容易的发现两个特点:
- 它是对称的
- 对角线元素有个系数2
实际上,任何矩阵加上它的转置都是对称的。这点很容易就能想通。
所以,在流体力学中,我们就管它叫对称张量好了。(狂野但简单粗暴的名字)
顺带一提,在OpenFOAM当中,它的代码表示就是
twoSymm(gradU)
为啥加个two呢?因为它是2倍的应变张量。我们其实可以这样记忆:它对角元有个系数2。
应变张量
刚才说了,应变张量其实就是1/2的对称张量。
E
=
1
2
(
∇
v
+
(
∇
v
)
T
)
E = \frac{1}{2}(\nabla \mathbf{v} + (\nabla \mathbf{v})^T)
E=21(∇v+(∇v)T)
即
E
=
1
2
(
2
∂
u
∂
x
∂
v
∂
x
+
∂
u
∂
y
∂
w
∂
x
+
∂
u
∂
z
∂
u
∂
y
+
∂
v
∂
x
2
∂
v
∂
y
∂
w
∂
y
+
∂
v
∂
z
∂
u
∂
z
+
∂
w
∂
x
∂
v
∂
z
+
∂
w
∂
y
2
∂
w
∂
z
)
\mathbf{E}= \frac{1}{2}\begin{pmatrix} 2\frac{\partial u}{\partial x} &\frac{\partial v}{\partial x }+\frac{\partial u}{\partial y } & \frac{\partial w}{\partial x}+\frac{\partial u}{\partial z} \\ \frac{\partial u}{\partial y}+\frac{\partial v}{\partial x} & 2\frac{\partial v}{\partial y} & \frac{\partial w}{\partial y}+\frac{\partial v}{\partial z}\\ \frac{\partial u}{\partial z}+\frac{\partial w}{\partial x} & \frac{\partial v}{\partial z}+\frac{\partial w}{\partial y} & 2 \frac{\partial w}{\partial z}\\ \end{pmatrix}
E=21⎝⎜⎛2∂x∂u∂y∂u+∂x∂v∂z∂u+∂x∂w∂x∂v+∂y∂u2∂y∂v∂z∂v+∂y∂w∂x∂w+∂z∂u∂y∂w+∂z∂v2∂z∂w⎠⎟⎞
简单明了。
应力张量
到上面为止,我们都没用到牛顿本构。也就是说,上面的内容对任何流体都是适用的。
这时候再套入牛顿本构方程。
假如为牛顿流体
T
=
−
p
1
+
2
μ
E
\mathbf{T}=-p \mathbb{1}+2 \mu \mathbf{E}
T=−p1+2μE
这里的1是单位阵。
所以无非就是乘以了几个系数。
T = − p 1 + 2 μ E = ( 2 μ ∂ u ∂ x − p μ ( ∂ v ∂ x + ∂ u ∂ y ) μ ( ∂ w ∂ x + ∂ u ∂ z ) μ ( ∂ u ∂ y + ∂ v ∂ x ) 2 μ ∂ v ∂ y − p μ ( ∂ w ∂ y + ∂ v ∂ z ) μ ( ∂ u ∂ z + ∂ w ∂ x ) μ ( ∂ v ∂ z + ∂ w ∂ y ) 2 μ ∂ w ∂ z − p ) \mathbf{T}=-p \mathbb{1}+2 \mu \mathbf{E}=\\ \begin{pmatrix} 2\mu \frac{\partial u}{\partial x} -p & \mu ( \frac{\partial v}{\partial x }+\frac{\partial u}{\partial y }) & \mu ( \frac{\partial w}{\partial x}+\frac{\partial u}{\partial z}) \\ \mu ( \frac{\partial u}{\partial y}+ \frac{\partial v}{\partial x}) & 2\mu \frac{\partial v}{\partial y}-p & \mu ( \frac{\partial w}{\partial y}+\frac{\partial v}{\partial z})\\ \mu ( \frac{\partial u}{\partial z}+ \frac{\partial w}{\partial x}) & \mu ( \frac{\partial v}{\partial z}+\frac{\partial w}{\partial y}) & 2\mu \frac{\partial w}{\partial z}-p\\ \end{pmatrix} T=−p1+2μE=⎝⎜⎛2μ∂x∂u−pμ(∂y∂u+∂x∂v)μ(∂z∂u+∂x∂w)μ(∂x∂v+∂y∂u)2μ∂y∂v−pμ(∂z∂v+∂y∂w)μ(∂x∂w+∂z∂u)μ(∂y∂w+∂z∂v)2μ∂z∂w−p⎠⎟⎞
结束
2022-6-10