基于布辛尼斯克假设的湍动能浮升力生成项表达式:
G
b
=
β
g
i
μ
t
Pr
t
∂
T
∂
x
i
G_{b}=\beta g_{i} \frac{\mu_{t}}{\operatorname{Pr}_{t}} \frac{\partial T}{\partial x_{i}}
Gb=βgiPrtμt∂xi∂T
由于热膨胀系数beta不能从mem.h储存文件中直接获得,所以我最开始尝试使用以下公式计算:
β
=
−
1
ρ
(
∂
ρ
∂
T
)
p
\beta=-\frac{1}{\rho}\left(\frac{\partial \rho}{\partial T}\right)_{p}
β=−ρ1(∂T∂ρ)p
但是此项也是无法获得:
∂
ρ
∂
T
\frac{\partial \rho}{\partial T}
∂T∂ρ
因此我尝试用以下的方式简化处理:
∂
ρ
/
∂
x
∂
T
/
∂
x
\frac{\partial \rho/\partial x}{\partial T/\partial x}
∂T/∂x∂ρ/∂x
但是这种方式在处理的时候面临一个很大的问题,在热流密度几乎为0的区域,会计算得到非物理的解,主要是温度的偏导数会是一个极小值,这样会过度放大浮升力的影响。
记录一次失败的UDF经历
最新推荐文章于 2021-08-10 00:05:15 发布