CFD中文网

Why the two-fluid model fails to predict the bed expansion characteristics of Geldart A particles in gas-fluidized beds: A tentative answer

http://www.sciencedirect.com/science/article/pii/S0009250908005137

最近用pointwise中的非结构网格生成了一个风力机旋转域中的block，检查质量发现最小角度0.14，网格质量特别差，但目前只会调整domain的网格质量，block的没看到相关教程，有大神可以指导一下吗？谢谢！

我觉的画边界层pointwise在这方面或许会提高网格质量

如果你的d是个dimensionedScalar，那We是个dimensionedScalar

更新科氏力源项：$\mathbf{f}\times\bfU=2\Omega\sin \lambda \times\bfU$，其中$\mathbf{f}$表示科氏力参数，一些文献会直接给出$\mathbf{f}$的值。如果绕z旋转，那么$\mathbf{f}$在$z$方向存在值。在这种情况下，$\mathbf{f}\times\bfU$的值为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
f_2u_3-f_3u_2\\
f_3u_1-f_1u_3\\
f_1u_2-f_2u_1\
\end{matrix}
\right]=\left[
\begin{matrix}
-f_3u_2\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
如果进一步考虑速度方向是x方向的，那么变为：
\begin{equation}
\mathbf{f}\times\bfU=\left[
\begin{matrix}
0\\
f_3u_1\\
0\
\end{matrix}
\right]
\end{equation}
可以看出如果地球绕z旋转，对于x方向的速度，科氏力导致y方向的偏离。

$\Omega$表示地球旋转速率，单位为[弧度$/\mathrm{s}$]=$1/\mathrm{s}$，$\lambda$表示维度[无单位]。因此$f$的单位为[$1/\mathrm{s}$]

在国外的时候我每天都去办公室，最近搬家之后，每天都做公交车去。

因为我们所离市中心很远，我住的介于市中心和所之间，这个公交车的终点站就是我们所，因此只要我上车遇到的人，不管周几，大部分都是去上班的。

工作日的时候，人还是很多的，当然不到北京那个程度，但是有时候也没有座位。除了工作日，我做的公交车里面有一半是中国人。

今天德国万圣节，这面一旦过节，就是个死城。这个节和周六周日连起来，连着4天放假。这4天，我每天早上来上班，有一半是中国人..

中国人都太勤奋了。按照这个发展速度，中国以后绝对世界第一。祝愿！！！

http://advances.sciencemag.org/content/2/6/e1501869

@赵一铭

@xjwang

CFDRFG的附录code还有个问题，求解$K$，不是按照论文中的原理计算的。他先用$PQ$没有代入$S(f)$的值，算了一组$K$，然后所有坐标点都用了同一组$K$。而根据理论，应该是要根据每个坐标点的非均匀$PQ$，求对应该点的$K$。DSRFG如果生成非均匀，理论上也要这么计算。
739858e5-ac4f-4736-94bd-9f4ec9f07a9d-image.png

关于CDRFG方法，他们的code里是对每一个点都生成一组PQK

最近看到一篇JFM文章《Guo, H., et al., An efficient and low-divergence method for generating inhomogeneous and anisotropic turbulence with arbitrary spectra. Journal of Fluid Mechanics, 2023. 970: p. A2.》，用的就是三维能谱，不是频谱，生成非均匀湍流场也是low-divergence。

通过vorticity求解涡量； 写个程序求解$\nabla w$
2.1 在这个程序里面，$\nabla w$的边界场$\cdot\bfn$再乘以粘性系数

说简单也简单，只不过繁琐了一些，好像没有几个命令就能出来的东西

太厉害了。。。。。。又做实验又写代码， 全才啊

射流 假设你是超音速流体 那么这是一个非常复杂的问题 这个问题是直接和初始条件相关的 对于瞬态 就是每一步的上部计算结尾的状态 因为它是对于初始条件非常敏感的 所以一个微小的积累可以导致很大不同的后续结果 实际上就是一句话 超音速这块实际上还是很难控制的 最好还是辅以实验 一般优化效果需要达到15%以上 才能说是设计真的有优势 如果要做严肃精确研究 fluent目前对于这种还不是太适合 但是fluent的优势是可以综合多场一起做出一个工程上可以大致量化的结论

最近在看Turbulence Modeling for CFD，里面有个CD，谁有CD的内容呀！

@李东岳 +1

感谢，有一个问题，为什么不把入口直接设置为湍流入口，是因为计算的雷诺数不高嘛

pointwise中，三组connector，AB、AC可以生成domain，BC生成一个空domain是什么情况？

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

连续方程 如果你是在模拟一个工程实际设备 这种不能达到-4的情况是比较常见 通常在-2左右 这代表在你的设备设计中间 真正的运行的时候 设备并不能保证一个进气和出气均衡的稳态（比如内部有一部分涡旋致使一部分进气永远在这里堆积） 而其他方面的原因上述几位已经讲解过了 对于工程模拟来说 连续性不达标的问题是比较常见的 你可以看看你的设备内部是不是存在涡旋流 或者是因为本身你的边界条件的设定是与实际不符合的 比如你在一个明显不均匀的出口面取了一个平均压力出口值 这样有的时候你需要扩充出口 直到真正符合工况的外界条件

这已经是比较成熟的方法了 应该aixprocess在15年前就已经成型了 做木材洗涤工艺的除尘优化 如果说无脑的办法 还是直接写碰撞列表做加减法 也可以直接设置某一方质量半径都为0 然后另外一方做简单的质量累加 但是实际上真实的情况比这复杂得多 液固的碰撞的结局以及最后的补集情况是多样性的 这有很多相关文献 但是对于UDF编写来说 因为是不同injection 还是有一定难度的

@李东岳 您好，抱歉，很久没来了。当时在学习FVM，最近在尝试利用openfoam解决NS方程。

vof模型只有一个混合相的温度 所以你这个方式其实难以实现 现在新出来了一个欧拉-欧拉-vof 可以尝试一下 实际上更好的办法是直接采用欧拉-欧拉同样算一遍 这样你可以方便的获取出口温度 实际上这样的方法也不会有太大的误差