CFD中文网

@东岳 也是一个比较强大的软件，不过主要在化工方面用的多，像一些化工单元流程，如精馏，吸收，蒸发，萃取什么的

这个不是用Aspen比较好吗？

有的，可以参考毕超的《计算流体力学有限元方法及其编程详解 》，上面有详细的利用FEM法对NS方程离散的过程，还有MATLAB程序

@东岳 是的，我看的文献大多数是关注速度，温度和涡量场变化，几乎没有涉及压力场变化，于是我编程计算的时候也没有考虑，结果老师问我为什么，我都不知道如何回答了

是有什么原因吗？

看了许多文献，在方腔自然对流中大家都分析速度，温度，怒塞尔数的变化，为什么没有考虑压力的变化？

是的，我编程做一个三维后台阶的例子，如果按照h , L, k求得话，那这个h对流换热系数如何确定?

@东岳 在资料上看的有平均努塞尔数和局部努塞尔数，在求解的时候有用微分方程表达（但都是无量纲形式），如果是原始变量的方式应该如何计算？

也就是说如何根据原始变量方程得出来的解求得努塞尔数

我是用MATLAB编程的

努赛尔数为无量纲的一个准则数，为壁面法向无量纲过余温度梯度，最后的定义我们可以看到是这个形式，是用无量纲形式表达；
但我在编程计算的时候采用的是原始变量，没有采用无量纲方程，这样的话如何计算努塞尔数？

@李东岳 谢谢李老师

@李东岳 温度和热流密度应该是标量吧？那热流密度为0就是对称边界了

对称边界条件是法向梯度为0，热流密度是法向为q，如果q=0，不是和对称一样吗？

我想问一下，如果热流密度为0，是不是就变成对称边界条件了？

@李东岳 哦，明白，谢谢李老师

我要模拟一个简单的二维强迫对流换热，边界条件为速度入口，压力出口0，入口温度为0，上下壁面温度为1，出口温度不知道，在压力出口时该如何设置，出口温度应该是热流密度为0啊！
他这个回流温度是什么？

毕超的《计算流体力学有限元方法及其编程详解》这本书里面详细介绍了方程各项的离散，不过使用MATLAB编写的，我现在主要就是参考这本书的

时间项的处理基本上都是差分发，然后引入sita因子

@李东岳 对于时间项的处理也是利用差分格式比如速度时间项的偏导：(u^n+1-u^n)/delta t