目录
本文将简要介绍常见的三种Boussinesq approximation。
法国数学家和物理学家 Joseph Valentin Boussinesq对于流体力学的研究十分深入,曾提出过多种用于简化建模的假说。这些假设都有着共同的名字,即Boussinesq approximation。上述的近似主要组成部分有以下三种:
-
buoyancy(浮力)
-
water waves(水波)
-
Turbulence modeling and eddy viscosity(湍流建模和涡粘)
Boussinesq近似(浮力)
事实上,对于以下可压缩流体的N-S方程组来说,几乎每一项中都存在密度。方程组高度非线性,同时求解代价非常大。所以,我们会考虑在一定情况下进行简化。
特别的,以浮力为主要驱动的流动,我们知道他的密度变化往往不会太大,此时就会考虑能否进行某种简化。
该近似多用于自然对流场景,即参与换热的流体由于各部分温度不均匀而形成密度差,从而在重力场或其他力场中产生浮升力所引起的对流换热现象。常见于自然通风、建筑环境等场景。对于许多这样的流动来说,Boussinesq近似值非常准确,能够大大提高模拟效率。值得注意的是,该近似显然不适用于那些高压缩程度的流体。
Boussinesq近似定义如下:
除了动量方程中的浮力项中的密度是温度的函数外,其他所有求解方程中的密度均认为常数。
能进行该近似的原因是,浮力驱动项是非常大的,相比于对流和时间中的密度变化的影响,要大得多。
列出连续性方程:
根据Boussinesq近似,此处认为密度是常数,可简化为下式
动量方程如下,除等式最右边的F之外,任何密度也都是常数:
F为单位体积流体所受体积力,假设该项只有重力
代入上式,F=ρg,动量方程如下:
一部分以重力压力的形式放进压力项,另一部分留在后面。
最后,给出该近似的应用判定范围如下,变化的密度要求远小于初始密度
Boussinesq近似(湍流)
Boussinesq 假设具体如下:
Reynolds应力项实质上体现了湍动产生的湍动应力,而湍动应力与粘性产生的分子粘性应力具有比拟性,见下式。事实上,由于分子运动与湍流脉动之间存在本质上的差别,上述简单地袭用分子运动来比拟湍流脉动在理论上是不严谨的。
Boussinesq近似(水波)
Boussinesq 近似的基本思想是从流动方程中消除垂直坐标,同时保留水波动对垂直结构的一些影响。
Boussinesq 近似中的步骤为:
-
围绕特定高度的水平和垂直流速进行泰勒展开
-
将泰勒展开截断为有限数量的项
-
使用不可压缩流的质量守恒(连续性方程)和无旋转流的零卷曲条件,将泰勒展开中量的垂直偏导数替换为水平偏导数。
此后,将 Boussinesq 近似应用于其余的流动方程,以消除对垂直坐标的依赖性。因此,得到的偏微分方程是根据水平坐标(和时间)的函数得出的。
考虑 (x,z) 平面中水平床上的势流,其中 x 为水平坐标,z 为垂直坐标。对z=0 处,我们对速度势 φ(x,z,t) 进行泰勒展开
在上式的基础上,进一步对x展开。之后截断为有限数量的项