@coolhhh 在 Wray-Agarwal湍流模型 中说:
按照这种方式自定义湍流模型,会提示 Duplicate entry 的警告
是的,这种加 libs 的方法就经常出现这种情况。直接忽略即可。
@coolhhh 在 Wray-Agarwal湍流模型 中说:
不要有这种警告,就要复制很多文件
是的,如果不想要这种警告,依赖这个模型的文件就都要修改。不难,但比较麻烦。
@coolhhh 在 Wray-Agarwal湍流模型 中说:
按照这种方式自定义湍流模型,会提示 Duplicate entry 的警告
是的,这种加 libs 的方法就经常出现这种情况。直接忽略即可。
@coolhhh 在 Wray-Agarwal湍流模型 中说:
不要有这种警告,就要复制很多文件
是的,如果不想要这种警告,依赖这个模型的文件就都要修改。不难,但比较麻烦。
@先先贝 那瞧不出来啥问题,得慢慢debug了。
@先先贝 看上去是#4 sqrt 开方函数有问题,是不是出现负数了?
@李东岳 CFD中文网陪伴了我整个读博生涯,感谢给我答疑解惑的大佬们。感谢李老师!
@ann ECN网站现在可以进去了。
@钟余 如果 $\nu$ 和 $\varepsilon$ 边界条件都没有设置为 0,那就要试试一个核跑去慢慢 debug 了。
@钟余 边界条件里湍动能 $k$ 给一个很小的值,不要设置为0。EDC模型里需要把 $k$ 作为分母。
@李东岳 哈哈,不是大佬。向东岳老师学习!
@wangfei9088 追踪了cpc的文章,看到名字,第一想到了你
@尚善若水 哈哈
@wangfei9088 再多问大佬一句,代码是open source的吗?我看到了CPC Library link to program files。我们用的求解器是基于OpenFOAM 12的内部求解器和v2312俄罗斯的一个求解器,都没考虑混合平均多组分扩散。感觉这对算缓燃还是挺重要的。
@尚善若水 是open source的,校对完正式发布就可以看到完整代码,提供了三个版本,可以满足大多数用户的需求。
@尚善若水 哈哈,是的,刚online的。消息这么灵通啊
@BlancheDH 上面的文件有好几个打不开。可以仿照这个帖子的#15楼附件代码写,用的就是OpenFOAM-10。
https://cfd-china.com/topic/6744/wray-agarwal湍流模型/15?_=1732935616393
@zju_klh cSpace和cTime需要用户给定,对应文献 (https://doi.org/10.1016/j.compfluid.2015.01.018) 公式里的$C_1$和$C_2$。
$$
\begin{equation}
P_{ij}
=\left(\frac {r_i + r_j} {\max (r_i+r_j, \mu_0)} \right)^{C_1}
\exp \left( -C_2 \frac{|\alpha_0 - \beta_0|}{\Delta t} \right)
\end{equation}
$$
其中,模型常数 $C_1 = 1$, $ C_2 = 0-0.3$。
一般就取1和0.3。
trajectoryCoeffs
{
cSpace 1;
cTime 0.3;
coalescence on;
}
@李东岳 没有下载权限,看不了
@yjc95 在 OpenFoam不可压+瞬态 对流项隐性离散推导看不懂 中说:
为啥是二重积分?
分别对体积 $V$ 和时间 $t$ 积分。
@yjc95 在 OpenFoam不可压+瞬态 对流项隐性离散推导看不懂 中说:
等式左边是如何转化到右边的呢?
体积分转为面积分,用的高数里的高斯定理。李老师的《笔记》里有推导的。
@HIT_lly 试试这个:
mpirun --use-hwthread-cpus -np 32 compressibleInterIsoFoam -parallel
@tiankai U 没问题,试试将fields后面的大括号{}改成小括号(),最后面再加一个分号;,如下。
fields
//{
(
U
{
mean on;
prime2Mean off;
base time;
}
//}
);
@13615617682 我不清楚fluent里具体用的单位,可能得查手册。但我见过的文献指前因子A的单位都是 $\rm m^3/(mol \cdot s)$。
@13615617682 啊,我不知道fluent中的指前因子是什么样的,它的单位是什么?能看一下不?
@13615617682 看 $A_i$ 的单位 $\rm kg/(m^2sPa)$,反应速率 $R_i$(单位:kg/s)应该是这样的:
$$
R_i =A_i \exp \left (-\frac{E_i}{RT} \right) S_i P
$$
其中,$S_i$是面积($\rm m^2$),$P$是压力($\rm Pa$)。
@hehaohao2013 报错提示,湍流模型只能选择kEpsilon或kOmegaSST,没有kineticTheory这个模型。