剪切应力张量

最新推荐文章于 2024-03-12 09:00:00 发布

寒冰彻骨

最新推荐文章于 2024-03-12 09:00:00 发布

阅读量4.5k

点赞数 3

分类专栏：控制方程文章标签： openfoam

hty

本文链接：https://blog.csdn.net/hanbingchegu/article/details/107429620

版权

控制方程专栏收录该内容

13 篇文章 11 订阅

订阅专栏

文章目录

剪切应力方程

对于剪切应力张量 $\vec \tau$ ，其九个分量定义如下：
在这里插入图片描述
其中，

$\kappa$ -体积粘度，这种粘度对于稠密的气体和液体并不重要，在许多文献中被忽略了;
$\mu$ -动力粘度;
$\lambda$ - 膨胀粘度(dilatation term)，大小为 $\frac23\mu$ .

根据N-S方程，有：
$\frac{\partial( \rho \mathbf U)}{\partial t}+\nabla \cdot {(\rho\mathbf U \mathbf U ) }=-\nabla \cdot \mathbf \tau -\nabla p+\rho \mathbf g$
对于各向同性流体，形变率 $\mathbf D$ :
在这里插入图片描述
则 $\tau$ （RHS的第一项的负号与剪切应力张量前面的负号相互抵消，令其变为正值）：

在这里插入图片描述
此时， $\lambda=-\frac23\mu$ 。这里，我们令 $\kappa$ =0,则：

在这里插入图片描述
若将压力梯度项放入剪切应力张量中，可得到柯西应力张量（Cauchy stress tensor） $\sigma$ :

在这里插入图片描述

可压缩流

对于可压缩流，连续性方程为：
$\frac{\partial \rho}{\partial t}+\nabla \cdot(\rho \mathbf U)=0$
拆分为：
在这里插入图片描述
将剪切应力 $\tau$ 引入，则：
$\tau=2\mu \mathbf D-\frac23\mu(\nabla \cdot \mathbf U)\mathbf I\\[2ex] \qquad \qquad\qquad\quad=2\mu \mathbf D \underbrace{ -\frac23\mu \left(-\frac 1\rho[\frac{\partial \rho}{\partial t}+\mathbf U\cdot \nabla\rho]\right) \mathbf I }_{`膨胀粘度项`}$
该项主要与体积的膨胀和压缩有关。

不可压缩流

对于不可压缩流体来讲， $\rho=constant$ ，膨胀粘度项可忽略。方程两边同时除以密度，有：

$\frac{\partial \mathbf U}{\partial t}+\nabla \cdot {(\mathbf U \mathbf U ) }=\nabla \cdot \frac { \tau}{ \rho}-\nabla \frac p{\rho}+ \mathbf g$
令 $p=\frac p \rho$ （运动压力）， $\tau=\frac{\tau}{\rho}$ ,
$\frac{\partial \mathbf U}{\partial t}+\nabla \cdot {(\mathbf U \mathbf U ) }=\nabla\cdot {\mathbf \tau} -\nabla p+ \mathbf g$

此时，
$\tau=2\;\nu \mathbf D$

其散度项

$\nabla\cdot\tau=\nabla\cdot\left[\nu\underline{(\nabla\mathbf U+(\nabla\mathbf U)^T)}\right]=\nabla\cdot(\nu\nabla\mathbf U)+\nabla \mathbf U\cdot\nabla \nu$