CFD中文网

或者是每个cell的最大压力值也可以

@李东岳 李老师有没有兴趣看下风工程这块的有个商业软件WT，里面是速度-压力耦合求解器，真的非常快非常快。我看了下，软件底层是从PHOENICS来的，然后搞了一套叫MIGAL的算法。

现在CFD用的SIMPLE这种segregated solver确实太慢了，城市这种尺度下的流动本身大概率又比较简单，不可压，也没有化学反应之类的，可以暂时不考虑温度？可能会有颗粒流的问题，但是不是也可以先只考虑流动对颗粒的影响。这种情况下感觉流场本身的求解其实就看求解器的效率了。

我搜了下，OpenFOAM这边其实已经有一些coupled solver 的工作，不知道成不成熟。

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

我认为这个是和几何无关的，衡量网格质量的参数中有一个叫aspect ratio，这个值为1时最好，像那种类似薄片或者细长的网格，这个值会很大，这种网格一般质量很差。

如果我没记错的话，流体中文网上的ICEM版块里有人画过4119桨的结构网格，也传了文件在上面，你可以去上面找一下看还在不在。（耐心往后翻一下）

这个放在他们搞涡那面，估计能讨论一上午 :duang:

我个人也不太熟悉。在我的文章里面，也不太关心涡的情况。

这个其实是个好话题，很多大佬在sci里面也明确说了，LES之类的模拟，如果对于这种类似的圆管，施加cyclic边界条件，会对湍流产生抑制。但是我也没见过大佬们进行对比。

另一方面，感觉边界条件，直觉感觉，单单的一个patch引起的差异，不至于那么大。LES模型参数的影响更大。如果谁有类似的文章，可以发来看看。

各位前辈好！
我对fluent的export ASCII file 文件有些疑惑如下：
该输出数据可以选择节点或者是网格，我输出两个点在同一个网格内，按网格输出是两点的数据是一样的，但是按节点输出是不一样的。（我现在对节点的认识将其认为是网格的端点，不知这种认识是否正确。）按我现在的认识的话，两个点的位置是不在任一节点上，但又因为两个点的输出数据不一样，所以应该是有数值差分做了处理，那么fluent是如何进行差分？fluent的理论指南和用户指南，只粗略介绍是按节点输出，不是很能理解。还望大佬们指点一二，感谢不尽！！

@程迪 666啊老铁，听君一席话，胜读百篇论文:146: :146: :146: :146: :146:

https://pubs.cstam.org.cn/data/article/mie/preview/pdf/lxysj2011-045.pdf
感觉就是附近涡造成的。我在哈尔滨上学的时候也见过树周围的这种现象，确实都是出现在风比较大的地方。

对于二维模型采用边界框方式进行划分，网格无法在Mesh中完成更新（只显示几何，网格生成不出来），而在三维模型采用边界框方式进行划分后可以在Mesh中完成网格的更新，想问一下大家用没有碰到过这种问题，是我的操作有什么问题吗？
1.二维边界框划分网格导入更新
cd591abd-856b-4041-9be9-f0e15412d5a1-image.png
2.三维边界框划分网格导入更新
80fc451e-60da-4712-984d-548cded523cd-image.png

VOF+动网格是interDyMFoam，用这个求解器就行

网上搜了一下，基本没有人讨论这个。记录一下，方便其他人能搜索到。

https://www.service-public.fr/particuliers/vosdroits/F2718

法国目前给未成年人颁发的居留卡是这个DECM。虽然DECM赋予未成年人从非申根地区，返回申根地区的权利。但实际操作过程中，很有可能无法入境除了法国之外的申根地区。个人经历：2024年从洛杉矶直飞西班牙，持有DECM，不给发登机牌。说不认识这个DECM。

个人第二个经历：从马尔代夫去巴塞罗那，购票后跟航空公司合适，航空公司说不认识法语的DECM，让我们跟海关确认。海关说这个需要跟航司确认。所以基本无解。有很大风险。

群友经历：

捕获s.JPG

希望以后某个人google搜索的时候，能搜到。

AI-CFD有时候可以实现CFD完全没法实现的问题。这样就比较好。但是有很多AI-CFD跟经典CFD是完全重复的，在这种情况下，目前来看还有很大提升。比如PINN的训练效率问题。但如果用PINN做隐藏流体力学，这就好多了。改天给你发个公众号问问。哈哈。看看评论区怎么说。

好的感谢李老师

@bart 慢慢增大进出口压差计算，在看看结果

是用八叉树算法吗？设置part上的最大尺寸也需要是2的幂吗？对这方面有了解吗？我也有这方面的一个疑问呢，之前在论坛发帖子还没人回。帮助文档里面也没写清楚，网上查了一些也没查到

@Tens 恭喜恭喜！！！！隔壁师兄又要嫉妒死了！:mihu: :mihu: :mihu: