CFD中文网

86f6fd36-007a-42a5-b105-de5d4fcefffa-image.png
cited:http://openfoamwiki.net/index.php/See_the_MRF_development

图为在旋转坐标系内解绝对速度的控制方程，其中 \vec \omega \times \vec U_i 这项称为科氏力，旋转坐标系内独有的力。他的效果如下面这个视频所示https://haokan.baidu.com/v?vid=4227630901810560602&pd=bjh&fr=bjhauthor&type=video
我的理解是可以将上面的式子理解为绝对速度产生组成的动量的输运方程，其中的通量是相对通量（\vec U_R \vec U_I）,同时由于动量是一个矢量，因此在旋转坐标系中还受到科氏力的作用。等式坐标是压力的贡献，和由于运输的物理量U_I的扩散贡献。

因此，我在思考是不是诸如湍动能，湍流耗散率这些标量，也可以认为在旋转坐标系内，运输物理量为湍动能，湍流耗散率这些标量，通量为相对速度的通量，但是由于这些量都是标量，因此没有多的项需要添加。另外比如SA模型中的发生项需要计算到voricity的大小，此处的voricity是不是也应该用相对速度？

2422ba2b-f9b8-402e-936e-8075325aae10-image.png 例如我生成一个part名字叫多“cavity1”，但是导入到fluent里面会多出“cavity1：007，cavity1：008，cavity1：008shadow，”这几个边界fb6ad2d6-0c7d-4266-a7b4-fa752482f0a6-image.png 。相应的网格也被分裂出好多细块。微信图片_20190126110946.jpg 微信图片_20190126110950.jpg

@xiaofei6538567 谢谢，确实没有找到办法直接编辑，最后非常窄的地方都单独分出来不加边界层了

可能北京今天发誓的人有点多，老天看不过去了

@blgy 谢谢您的回复，
1.第一个问题明白，我继续试一下。
2.对的，这个模型原来的确是采用混合网格，部分采用四面体，部分采用六面体，但是现在想整体都用六面体尝试，所以难度确实比较大。硬件资源应该问题不大，但是只是想用六面体完成，有这个要求
3.周期性网格我还了解的不太多，包括您后面说到的一些，愿闻其详
4.较早前曾使用gambit进行网格划分，但是网格量比较大，之前的计算网格量大概是是过亿了，gambit似乎当初处理不了这么多的网格量，所以才放弃的。
5.确实目前因为icem的网格思路与操作还比较陌生，整体还是比较困难

CFDer的旗帜 4600！

@qjh888 h好的，谢谢啦

adept

@两月三年 good ，好主意。我试试哈。

近壁面网格长宽比本来就很大，如果用RANS计算，有的模型对近壁面y+要求很小，比如SST k-omega要求y+≈1，但是对x+或者z+基本没有明确的限制，为了减少网格数量，一般流向和展向的网格无量纲参数可以在百千量级，这样的话近壁面网格长宽比就会很大，也是可以正常计算的。

请教各位，该怎么对边界附近的单元赋值，比如下面这个UDF.我想将一个区域的 C_UDMI(c, tt, 0)全部设置为0，然后边界附近的 C_UDMI(c, tt, 0)设置为1，该怎么操作？谢谢！

void set_zone(Thread* t) { cell_t c; real vol; Thread *tt = THREAD_SUPER_THREAD(t); begin_c_loop(c, t) { if (THREAD_ID(t) == 1047) { C_UDMI(c, tt, 0) = 1; } end_c_loop(c, lct) } void set_cells(Thread *t) { Thread *tt = THREAD_SUPER_THREAD(t); face_t f; Thread *t_nbr; cell_t c0, c1; Domain *mix_domain = Get_Domain(1); for (int i = 0; i < num_boundaries; i++) { t_nbr = Lookup_Thread(mix_domain, 1); begin_f_loop(f, t_nbr){ c0 = F_C0(f, t_nbr); c1 = F_C1(f, t_nbr); if (c0 != -1 && THREAD_ID(THREAD_T0(t_nbr)) == 1047) { C_UDMI(c0, tt, 1) = 0.0; } if (c1 != -1 && THREAD_ID(THREAD_T1(t_nbr)) == 1047) { C_UDMI(c1, tt, 1) = 0.0; } } end_f_loop(f, t_nbr) } } DEFINE_ADJUST(aaa, mix_domain) { Thread* mix_thread; int phase_domain_index = 0; thread_loop_c(mix_thread, mix_domain) { Thread* t = THREAD_SUB_THREAD(mix_thread, phase_domain_index); if (t != NULL) { set_capillary_pressure_sources(t); set_boundary_cells(t); } } }

可以参考用户指南4.4，里面还是有不错的概括的。

@shuchang769 一样的思路啊

对于给定的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}

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

是OpenFOAM吗？
如果是的话参考
https://www.cfd-online.com/Forums/openfoam-meshing/61596-cyclic-bcs-pointwiseopenfoam-export.html

我设置了interface，check mesh的时候还是显示Unassigned interface zone detected for这是为什么呢？

@winsway_zero 多谢大神指点

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

@金石为开 Foam::fvPatchField::evaluate(const Pstream::commsType) 函数里面有调用 updateCoeffs()