CFD中文网

如果你认为动量方程增量很小，可以尝试把动量方程的松弛因子设为1，看看稳定性:duang: 好像问题不是这么简单

@李东岳 噢噢，好的

有哪位大神知道怎么在fluent中使用level set 函数计算物性吗？采用的是vof+coupled level set

@micro 网格数量越多，网格尺度越小，计算肯定更加不稳定，您可以查一下CFL数

@星星星星晴 好的 我搞了一个

对于我目前所计算的模型所涉及是湍流模型以及多相流空化模型，在湍流模型尝试使用RSM（quadratic pressure-strain ）时，总是提示错误，而对于linear pressure strain或者Stress-Omega Model 都能够正常计算，并且每次计算前都使用了RNG k-e计算的结果求得初始值，网格部分没有问题，ICEM结构化网格质量0.5以上，skew在0.65一下，请问各位是怎么原因导致的呢？

您好，我也是做混合物燃烧，想问问您得到等效分子式之后，怎么在fluent中输入的？谢谢

感谢东岳老师，我尝试着改小了，问题得到了部分解决。

forAll循环在openfoam里面特别多，这种小操作不会减慢速度，

https://www.cfd-china.com/topic/5198

参考这个帖子，有人可以做个STAR-CCM、COMSOL、CFX的算例么？

@东岳 在 在OF2.3.0版本如何编译低版本下的湍流模型？ 中说：

模型计算公式拔出来，自

懂了，，，谢谢东岳老师!

下面的图按照初始setField→计算一个时间步排列。
不知道为什么，经过计算后内部场的速度很快发生变化。 而内部场应该是利用势流理论估算的，但是拿N-S算出来的速度不应该跟势流估算的差这么大啊...

初始波峰在左侧
1Abefore.png
1Before.0000.png
1AAfter.png
1After.0001.png

波峰传递过来时
0A.0000.png
0.0000.png
0A.0001.png
0.0001.png

波峰在右侧
2.5A.0000.png
2.5.0000.png
2.5A.0001.png
2.5.0001.png

波峰离开
-5A.0000.png
-5.0000.png
-5A.0001.png
-5.0001.png

@一朵苍凉 我也不太懂，可以看看控制领域的相频特性，也许有帮助:zoule:

@李东岳 哦哦，没问题了。:142:

SV_locate_point (real [], CX_Cell_Id *);
你可以用这个函数来获取这个点的所在网格的信息 要自己去看下头文件
如果这是你的要求的话
而如果不是的话 那么只是一个纯解析几何的问题

@李东岳 谢谢老师解答，之前一直出现这个问题导致每次画完网格都得要重启，我重新装下系统在试试

各位大佬好~向大家请教一个问题。

式（1）是大涡模拟法夫尔滤波操作以后的连续性方程
\begin{equation}
\frac{\partial{\bar{\rho}}}{t}+\frac{\partial{\bar{\rho}\tilde{u}_i}}{x_i}=0
\end{equation}

如果对$\tilde{u}_i$取时均值，即$<\tilde{u}_i>$, 这样得到的是雷诺时均值的速度吗？如果想要取法夫尔时均值，是不是需要做$<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>$的处理？P.S(这里用<φ>表示对φ取时间平均)

第二个问题与第一个问题相似：基于法夫尔时均的可解尺度的湍动能的求解是否可以按照式（2）不可压缩流体的求解方式求解呢？我感觉式（2）的解更像是雷诺时均的解，因为可解尺度的湍流里面也包含了密度脉动。
\begin{equation}
k_{resolved}=\frac{1}{2}(<\tilde{u}_i\tilde{u}_i>-<\tilde{u}_i>^2)
\end{equation}

还是说需要用式（3）继续用法夫尔时均进行计算呢？
\begin{equation}
k_{resolved}=\frac{1}{2}\big(<\bar{\rho}\tilde{u}_i\tilde{u}_i>-<\bar{\rho}>(<\bar{\rho} \tilde{u}_i>/<\bar{\rho}>)^2\big)/<\bar{\rho}>
\end{equation}

@zhangxc0223
Hi
因为并没有做deterministic的模拟，所以code也基本没看过。。
stochastic collision 在计算collision probablity的时候是有限制的，具体你可以看一下O‘roucke的论文。具体是怎么定义的，怎么得到一个non-dimensional collision probablity的，还有如何决定这两个parcel发生了碰撞。

虽然碰撞在这个instantaneously 随机的，但是在统计学上，比如你有1M次碰撞的话，总的碰撞的结果就是相似的。
这也就是为什么stochastic碰撞会比deterministic相对来说’便宜'的原因，通过统计学的方法达到相似的碰撞结果，但是不需要追踪所有的parcel。

lagnrangian maxCo 存在的意义，个人认为，在一定程度上可以理解为再一次离散了每个euler time step，所以maxCo一定程度上决定了每个lagrangian time step的大小，和Euler 的Co 定义是相似的。

碰撞是不会被错过的，因为deterministic理论上来说是追踪所有的parcel，所以下一步这个parcel会在哪，走的路径应该也是计算的，所以不会存在错过的，stochastic collision 会计算两个parcel的碰撞概率，当概率很高的时候，也是不会错过的，而且就算是这次错过了，在茫茫parcel海中，总会碰到的。。统计学上是不会错过的

以上仅为个人理解

TIM图片20190509144012.png TIM图片20190509144004.jpg
第一幅图是edem中的传送带例子，第二幅图现在算出的静态（支撑颗粒的板不动）结果，其实我们模型需要板动起来，类似第一幅图的传送带，最终得到颗粒水分、温度分布，以及后期还会加上化学反应。希望大家给点意见，可不可行，然后有没有类似的案例。先谢！

对于给定的NDF，划分为$i$个$N_{pp}$，对每个$i$上定义$k$阶矩$m_k^i$，给定$m_k^i$，可以计算第$i$区间的节点$d^i_0,d^i_1$以及权重$w^i_0,w^i_1$：
\begin{equation}
\begin{split}
w^i_0&=w^i_1=0.5
\\
d^i_0&=m_1^i-\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\\
d^i_1&=m_1^i+\frac{1}{\sqrt{3}}\sqrt{\frac{m_3^i}{m_1^i}-{m_1^i}^2}
\end{split}
\end{equation}
对于仅考虑破碎的PBE：
\begin{equation}\label{pbe}
\frac{\p n(d)}{\p t}=\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'-g(d)n(d)
\end{equation}
对方程\eqref{pbe}在$i$上取$k$阶矩：
\begin{equation}\label{m}
\frac{\p m_k^i}{\p t}=\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')d^k\beta(d|d')n(d')\rd d'\rd d-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}
\begin{equation}
\begin{split}
\int_{d_{i-1/2}}^{d_{i+1/2}}\int_d^{d_{max}}g(d')\beta(d|d')n(d')\rd d'\rd d&=
\int_{d_{i-1/2}}^{d_{max}}g(d')n(d')\left(\int_{d_{i-1/2}}^{d'}\beta(d|d')\rd d\right)\rd d'
\\&=
\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)
\end{split}
\end{equation}
Therefore
\begin{equation}
\frac{\p m_k^i}{\p t}=\sum_{m=i}^{N}\sum_{j=0}^2g(d_j^m)w_j^m\left(\int_{d_{i-1/2}}^{d_j^m}d^k\beta(d|d_j^m)\rd d\right)-\sum^2_{j=0} g(d_j^i)w_j^i(d_j^i)^k
\end{equation}