CFD中文网

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

有的时候，某些内部的情况，没有几何，但是需要对网格进行一些处理，就需要手动的在不能关联的情况下，切一下。 某些需要投影的情况，比如地面上悬空的一个圆筒。在地面上只有圆筒的投影，也不存在几何关系。需要把edge投影过去，就需要spline方法。

后来这个怎么样了

@五好青年 大佬，后来者大胆问下这个有解决的好办法或者好思路了么

参考文献
Demirdzic, I.; Peric, M., Space conservation law in finite volume calculations of fluid flow.

@doubletrump 您好，我想向您请教，投影面积指的是？如果是一个圆柱，风垂直吹过来，指的是直径乘以高还是圆的面积？

OpenFOAM中使用multiDimAMR编译和运行正常，使用reconstructParMesh也正常，但是之后再使用reconstructParMesh会报错。
第一个时刻还是正常的，第二个时刻发生报错，还有一个算例是第5个时刻开始报错，报错内容一样。使用paraview中的decomposePar也能查看结果，为什么会出现这个问题呢~
Reconstructing fields for mesh region0
Time = 0.01
Reconstructing FV fields
Reconstructing volScalarFields
alpha.water
alpha.water_0
p
p_rgh
rAU
Reconstructing volVectorFields
U
Reconstructing surfaceScalarFields
alphaPhi0.water
phi
Reconstructing surfaceVectorFields
Uf
Reconstructing point fields
No point fields
No lagrangian fields
No FA fields

Time = 0.02

--> FOAM Warning :
From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 210
Number of patches has changed. This may have unexpected consequences. Proceed with care.
--> FOAM Warning :
From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 210
Number of patches has changed. This may have unexpected consequences. Proceed with care.

**--> FOAM FATAL ERROR: (openfoam-2012)
Processor 2 has a different polyMesh at time 0.02 compared to any previous processors.
Please check time 0.02 directories on all processors for consistent mesh files.

From Foam::polyMesh::readUpdateState Foam::processorMeshes::readUpdate() in file processorMeshes.C at line 192.**

FOAM exiting

块就是将你的整个计算区域分割成几个区域，每个区域有各自的网格分布规律；
模型应该是你计算的几何外形吧

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

@百善孝为先0626 您好这是用的我们课题组的超算账号，天河那边的

最近要做一些跨临界的问题，但是在跨临界中，如下图，由于伪沸腾的影响，流体的热力学性质会在伪沸腾区呈现剧烈的变化，形成类似于激波的区域。但是在亚临界区域，压力对密度的变化极为敏感，所以传统的全守恒格式会引起压力震荡。很多研究者针对这个问题做了很多修正。比如：求解压力传播方程并加入人工扩散、计算左右通量求黎曼解并对最后的结果熵修正。但是目前的算法总体来说很繁琐，且由于通过压力反推能量会引起能量不守恒。主要有两个问题，非常感谢您的解答：

求解压力传播方程的那个算法很类似simple算法的pressure based方法，为什么现在几乎所有的计算空气动力学论文里面几乎没有使用过simple算法，更倾向于TVD RK这种显示迭代？ 显式格式很适合GPU加速计算，想开发一套GPU的代码用来计算。但是GPU对双精度计算的支持很差。当在激波附近使用高阶离散格式的时候如果使用单精度浮点数计算，截断误差是否会对计算结果产生比较大的影响？
3.png

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

@bestucan 是的 老师 就是直接采用C_T进行读取多孔介质年内的气体温度，然后通过努塞尔数获得对流换热系数，然后可以求解出多孔介质骨架的温度

@余正东 这个我关注了，未来可期，目前只是初步移植了一些

最近老师让做一个垂直轴风力机的模拟，但是没做过类似的，想请问一下各位前辈，怎么在运动中设置旋转之后，增加一个旋转时的阻尼用来模拟有负载下的风机转动,麻烦前辈们了

@firejet 补一下图9a7cf44e-8063-4afa-ac14-ffff490811e2-Screenshot from 2022-04-26 23-16-44.png be16cb5a-0889-4fc1-b536-d22654dae7d1-Screenshot from 2022-04-26 23-17-57.png 2e8ac1c2-54b9-4b90-84f2-b1ae0e2c515d-Screenshot from 2022-04-26 23-18-35.png 8c2ea078-a7aa-43b1-880e-149f6a54fe15-Screenshot from 2022-04-26 23-33-06.png

$\tau$这个雷诺应力，对角线的是normal stress，非对角线的是shear stress。不过你说的更详细，感谢分享！

@hy1112006 您好，如果方便的话可以有您的邮箱或者联系方式交流一下使用pisoFOAMopenfast吗？

往一个体积为30平方米的水池里喷射带有泥沙的水流，开始阶段水池内没有水。计算域中包含
gas（空气）
liquid （水）
solid （泥沙颗粒）
目的：观察泥沙颗粒对水池底部的冲蚀程度。

应该采用 Mixture+DPM模型
还是采用 Euler+DPM模型 ？

我觉得这个是一个自圆其说的东西。如果看sci，肯定有人对比LES和DNS的结果，证明他们的LES或者DNS的结果更好。

但我认为，本身欧拉拉格朗日里面存在了很多模化，在粒子受力的角度，存在更多的模化，尤其是湍流分散力那面，基本都是参数。本身DPMFoam和DNS就是冲突的东西。所以如果做这方面研究，自圆其说是最重要的。