记录一次失败的UDF经历

基于布辛尼斯克假设的湍动能浮升力生成项表达式:
$$G_b=\beta g_i\frac{{\mu }_t}{{\mathrm{Pr}}_t}\frac{\partial T}{\partial x_i}$$
由于热膨胀系数beta不能从mem.h储存文件中直接获得，所以我最开始尝试使用以下公式计算：
$$\beta =-\frac{1}{\rho }{\left(\frac{\partial \rho }{\partial T}\right)}_p$$
但是此项也是无法获得：
$$\frac{\partial \rho }{\partial T}$$
因此我尝试用以下的方式简化处理：
$$\frac{\partial \rho /\partial x}{\partial T/\partial x}$$
但是这种方式在处理的时候面临一个很大的问题，在热流密度几乎为0的区域，会计算得到非物理的解，主要是温度的偏导数会是一个极小值，这样会过度放大浮升力的影响。