CFD中文网

@东岳 李老师，我没推出来

@东岳 在 ANSYS Fluent较大体积力模拟问题 中说：

PRESTO!

Fluent中的PRESTO!压力插值格式就是用交错网格吗？把压力存储在网格中心，把速度分量存储在网格界面上？这种PRESTO!插值格式FLUENT理论指导也不说清楚，就一段话就解释完了，真是操蛋。

根据这个公式也推导不出公式3啊

李老师，公式3是怎么来的啊？我自己推了半天也没推出公式3是怎么回事，那个转置怎么处理的？怎么最后给处理没了呢？

非牛顿问题中的交叉扩散项.pdf
各位前辈，我最近遇到一个问题，困扰良久，非牛顿流体由于黏度是变化的，所以交叉扩散项不能省略，所以我尝试着在结构化网格中进行离散，但是有两项我实在不知道怎么处理，如图所示，希望各位前辈能够指点一下，特别是fluent中关于这个交叉扩散项又是怎么处理的？如果图片看不清楚，可以下载PDF文件，我PDF文件也上传了的。

@东岳 好的，李老师，我先按照《computational methods for fluid dynamics》这本书上说的那种把这一项也处理成梯度的形式试着去推导一下，看能否推导出来

@东岳 李老师，我看了你这个方程推导的，我只是比较疑惑的是，如果是以比焓作为变量来求解，那么扩散项的温度是不是要先变成比焓的形式才能进行离散？

目前有限体积法多用于牛顿流体，当用于处理非牛顿流体时，暂且不考虑非牛顿流体的黏弹性问题，这里的非牛顿流体处理为广义牛顿流体，也即满足牛顿流体的本构方程，但是即便是这样简化，仍然存在黏度随温度和剪切速率的变化需要考虑，因此当黏度不是定值时，图片中方程标黄的部分就不能省略，我查阅了很多的CFD教材，所有的算法都是黏度为定值下开发的，那么问题是当黏度不为定值时，标黄的那一项应该怎么处理，或者说怎么离散？虽然我知道fluent可以计算这种情况，但是我不知道fluent是怎么处理标黄的那一项的。
目前我只在《computational methods for fluid dynamics》一书中提到过这一项在黏度变化的情况下不能省略，具体怎么处理也没详细的说明，这一项的处理实在是很困惑。

Fluent中的能量方程中既存在比焓变量（非稳态项和对流项）又存在温度变量（扩散项），这跟CFD相关教材中所述的广义标量方程不一致，那么Fluent中的能量方程在求解的时候到底是以温度作为变量求解还是以比焓作为变量求解呢？

@东岳 课题组之前的师兄是用matlab做的，但是主要是求解热传导方程的，跟流动已经不沾边了

Fluent的theory guide中chapter21:solver theory一章中在广义标量方程求解小节中写的是扩散项的离散中关于控制体界面处的梯度值用的是中心差分，而在空间离散小节中写的是用Green-Gauss Cell-Based,Green-Gauss Node-Based和Least Squares Cell-Based三种方法去计算扩散项中的控制体界面处的梯度值，那么有三个问题：1、前后两种说法是矛盾的，fluent如果是用的中心差分去计算界面处的梯度，那么在非结构化网格中，这个中心差分的表达式是什么？2、如果是用另外三种用于计算控制体中心的梯度的方法去计算控制体界面处的梯度值，那么如何通过控制体中心的梯度值来计算控制体界面的梯度值？3、最终极的问题，fluent关于扩散项在界面的梯度值的计算究竟采用的是什么方法？是中心差分还是空间离散小节中描述的三种方法？

PISO算法是不需要对压力进行亚松弛处理的，可是为什么fluent中PISO算法中默认情况下对压力是要进行亚松弛处理的，且松弛因子默认值为0.3。在实际计算过程中，我尝试过在PISO算法下将压力松弛因子修改为1，但是计算没有默认的0.3稳定，容易发散，实在搞不懂是怎么回事。

第一次在CFD中文网中发帖，比较激动，还在摸索发帖的方式之中。
最近我们碰到一个棘手的问题
如图所示，先往一个模具型腔里面注入一种聚合物熔体（红色，称之为皮层熔体），之后再注入另外一种皮层熔体（芯层熔体），麻烦的地方在于这个芯层熔体里面填充了大量的固体颗粒，这些固体颗粒在整个材料里面起着功能化的作用，比如导热导电传感等等，在实验过程中，我们发现固体颗粒是倾向于朝着皮-芯熔体界面运动的。前期探索性地用VOF+DPM模拟了这些颗粒在这个过程中的运动情况，但与实验观测差异较大，并且由于芯层熔体里面固体颗粒含量特别高，体积分数达到70-80%左右，用DPM并不合适，所以我们想用欧拉-欧拉模型来说明这个固体颗粒在这个过程中的运动情况。
现在我们面临的问题是，究竟能不能用欧拉-欧拉模型来说明这个体系下固体颗粒的运动情况，如果没有前面的皮层熔体，我们能够确定一定可以用欧拉-欧拉模型来说明，但是问题是前面有皮层熔体，并且这个皮层熔体和芯层熔体不互溶，界面非常清晰，似乎得用VOF+欧拉-欧拉才能描述这一个过程。
我们对多相流中的数学模型理解并不是很透彻，所以上面的描述可能存在一定的错误，并且时间很紧张，只能考虑用fluent实现这一过程，用openfoam时间来不及。
所以请相关的大神们讨论一下这个体系，包括如何描述这一个非牛顿多相非等温的流动与传热的过程，以及能不能在fluent中模拟，如果能，应该怎样实现。

对于CFD和NHT的模拟来说，大家都是用Fluent或OPENFOAM吗？有不有人用MATLAB从头到尾模拟的呢？