CFD中文网

请教大家一个问题，笛卡尔网格，已知网格8个顶点的坐标，以及网格中心的物理量的值，但网格排序是杂乱无章的，有什么好方法进行后处理呀？

@浪迹天大 在mathjax里面要设置跳过某些关键词，比如skipTags: ["<pre>","<span>","<code>"]

这样啊谢谢，我看论文中测试的时间都没有这么久，就看着能加速反应速率，时间加快些不

icoFoam里面压力是除掉密度的压力

很久之前写过。具体忘了。但 return 返回的才是source

ds[eqn]是在定义源的表达式是ax还是ax+b，udf 手册里有对ds[eqn]的详细解释。

是 return absorption_rate ?我记得之前写的好像是这样的套路：ds[eqn]=0 然后 retrun source

可以找udf手册对比着写

@东岳 惭愧啊。整个实验室就自己在零零碎碎做一些流体力学模拟。大多数人都偏催化、化学实验类的。一般上论坛来也就看看帖子。都忘记cookie过期没登陆了。。
我们这边还是很认可AIChE J的。学校认定上还是top级别的。

我模拟的对象是反应堆内的稳态流动传热，由于堆芯内组件数量庞大，精细建模复杂，因此采用多孔介质模型对堆芯区域进行简化。

为了验证这种等效的合理性和误差，我先对单个组件栅元进行多孔介质简化，并将多孔介质模拟结果与未简化的模型结果对比验证多孔介质模型的合理性。

单个组件栅元的截面图如下：
图片1.png

其二维轴对称模型如下图：
f510f5de-0d5e-4389-b577-ad2d3defe29c-image.png

模拟结果
(1)未采用多孔介质模型
我用了二维轴对称模型来模拟，设置固体流体内热源，入口质量流量和出口压力边界条件，流体和固体分别加载体积热源P1和P2（W/m3），流体和固体的导热系数均为常数，得到的稳态温度场和压损分布如下：
bca065f1-7cc4-4385-8113-d7e0acd6d23e-image.png
流体出口温度700℃，固体由于带内热源温度高于流体，固体最高温度为799℃，固体平均温度为756℃，组件压损为2400Pa

(2)采用多孔介质模型模拟
由于固体带内热源，流固存在较大温差，采用非热平衡模型
根据栅元几何结构指定：孔隙率γ=0.1 ，界面换热系数hfs=5000W/M2K（与未采用多孔介质模型算得的换热系数一致），界面面积密度Afs=6.71m-1

通过流速-压力实验关系确定粘性损失项系数和惯性损失项系数，只考虑x方向的流动损失，y方向系数放大100倍

别给多孔介质流体和多孔介质固体设置了内热源，分别为P1×孔隙率和P2×（1-孔隙率），初始化后检查内热源加载总量无误；流体固体的热导率输入值和上个模型一致

得到的模拟结果：
af1efa10-a95f-4e49-a69c-63b4b4357ca0-image.png
采用多孔介质非热平衡模型模拟得到的流体温度是正确的，压损也是正确的，但固体温度分布明显不对，固体域温度热点仅为725℃，平均温度为711℃，均远小于实际模型算得的固体温度值。

我自己的猜测：
多孔介质模型将单组件变为流固均匀介质后，不能模拟固体区域沿导热路径带来的固体温升？刚好我这个栅元模型中固体由于导热产生的温升非常高，如果是这样，说明多孔介质非热平衡模型不适用于我所研究的问题，因为简化后会显著低估固体温度热点。

请问我的理解是正确的么

@一颗李子
@一颗栗子
@一颗橙子
@一颗鸭蛋

:ok: :ok: :ok:

@benf 不客气不客气，欢迎引用我们组的论文~这个模型有外国组已经验证过了，所以如果你有兴趣编应该可以原样复现

@jinlinna 感谢！

各位前辈，大佬们，最近在学习中遇到了几个问题希望各位前辈指教一下：
目前在大多数文献中我看到气液两相流中的曳力模型系数都是以下形式：
$k=\frac{3}{4}\alpha_d\rho_c\frac{C_D}{d_d}\left|U_c-U_d\right|$
但是最近学习文献的时候在核口的一些CFD代码文章比如SIMMER-III，AFDM，还有一些其他程序里出现了另外一种形式：
$k=A+B*\left|U_c-U_d\right|$
B项一般是跟$C_d$相关，在SIMMER-III中其形式为：
$B=\frac{1}{8}a_i\rho_cC_D,a_i是交界面面积$
而A项给出的解释一般是comes from Stokes law for viscous flow，形式一般是：
$A=\frac{3\alpha_c\mu_c}{d_d}$
比较好奇这种形式的曳力系数是怎么得到的，尤其是A项，在文献里也没看到相关的引用。
是否是由于这些代码一般是用来计算液态金属-气体两相流，因此比较特殊呢？
希望各位大佬不吝赐教，非常感谢！:xiexie:

@紫色星光 在设置节点数的小窗口，往下拉，有个copy parameter，不勾选。如果不勾选还会联动，那就是切块的问题了。

替代文字 我模拟的是一个矩形泡沫铜浸润石蜡的传热传质过程，现在模型验证一直没弄好。首先是编写了一个简单的UDF，采用热平衡模型，导入fluent后提醒我定义比热的那一块有错误，求各位前辈帮我看一下；其次，我添加了Boussinesq项后，优化了一点点液相云图中熔融线一直平行的问题，但是与论文中的还是有差距，且我的熔融线上下厚度保持一致（实际上应该下面厚很多），UDF中还需要添加什么项来解决这个问题呢？以下三张图分别是我的UDF、我跑出来的液相云图以及论文的液相云图，求各位帮我解惑，不甚感激。2.png 3.pn替代文字g 4.png

@李东岳 好的，谢谢老师！！

需要处理滑移交界面，metis可以对分块非结构网格进行划分吗？交界面处的网格分区会是什么情况？比如局部交界面两侧的网格会划分为同一并行区域，还是两侧网格属于不同并行区域？

本人按在实际出口下方5cm处建立一个出口边界，实际出口则设为内部边界，最后调整参数，得到了与文献相似的结果。感谢各位的帮助！！

@yfclark

实际工程上说的压力是什么压，压力表测到的是什么压，还有就是数值模拟计算得到的压力结果我们应该关注什么压呢？或者说从动压静压总压分布图上我们能看出什么东西。谢谢

@TINGHAIK 谢谢您，问题已解决

不清楚商软那面怎么回事，本身我也没太关注这个，我查了一下，OpenFOAM本身这面的streamline看起来是你说的patiline：

This function object generates streamline data by sampling a set of user-specified fields along a particle track, transported by a user-specified velocity field.

后处理软件paraview的streamline是如何定义的还没找到。