CFD中文网

@zzcfd 改fvSolution，把p改成这个solver PCG; preconditioner DIC;试试

@cfd-china 这个问题，有意思。去知乎上搜一下

不匿名，中航工业624所搞航空发动机的，工作9年，本科生，明年评高工，普通职工，不是领导干部。
税后年收入十万上下，在成都市郊区，房子单价五千左右。住140㎡的房子，开十万左右的车。自己觉得还行，就是每隔一段时间工作压力略大。

9年年收入10万年。
看来还是得CS

@XieXiaoyang 好文章，我去研究研究

@李东岳 李老师您好，之前确实没考虑到压力波。不过我这里的泄压过程是缓慢过程，是通过密闭空间周围的多孔介质渗流漏出去的，给的流速很低，大概在10e-4m/s量级，也会有这种问题么？然后考虑到空间内的热对流，我可能更倾向于使用buoyantPimpleFoam求解器？

你加的那些私有成员，声明了么

@李东岳 我在KinematicParcel.H文件，以及这个文件的前面，我都对我新加变量声明了，不然不声明这些东西源代码都无法进行编译，会有很多error。（目前只不过编译时候有一些warning，连着源代码以前自带的一些东西都在warning）。
目前这个问题是我开始编译求解器时候，求解器识别源代码时候出现的情况。

我再看看我代码声明问题吧，我也不知道和这个后没有关系

另：之前我代码是这样

template<class ParcelType> inline Foam::KinematicParcel<ParcelType>::KinematicParcel ( const polyMesh& owner, const vector& position, const label cellI, const label tetFaceI, const label tetPtI ) : ParcelType(owner, position, cellI, tetFaceI, tetPtI), active_(true), typeId_(-1), nParticle_(0), d_(0.0), dTarget_(0.0), U_(vector::zero), rho_(0.0), age_(0.0), tTurb_(0.0), dRt_(0.0), //- myadd Ro_(0,0), //- myadd pgo_(0.0), //- myadd Pc_(0.0)， //- myadd UTurb_(vector::zero), rhoc_(0.0), Uc_(vector::zero), muc_(0.0) {}

他求解器会报错

/home/zly/OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/intermediate/lnInclude/KinematicParcelI.H:105:12: error: expression list treated as compound expression in mem-initializer [-fpermissive] muc_(0.0)

看起来只报错最后一个地方，我也不知道为什么。这些代码的声明和修改，一部分是我看网上的一篇资料修改，另一部分是我参考着原始代码格式，进行改写。

你看看是不是有U=0的点，specify时强制要求选中全部zone，可能是固体边界也被计入，导致无法计算。

@sxz0823 你好，请问有关concentration和parcelconcentration您是怎么计算的，求指导。

请问像这种狭窄处要划分边界层，怎样处理能保证网格高质量呢

@程迪 在 GPU VS CPU 中说：

隐式的得回头修改这个自由度那个自由度

你好，请问您所说的自由度是什么意思，能够稍微解释一二，不胜感激，谢谢！

icem网格划分导入fluent时检查失败，显示Error: interface zone 3 has two adjacent cell zones.
4f68eb65-132e-4314-b0f6-6964dc466673-image.png

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

@J 我们课题还是偏这个软件，谢谢您的回复

你是什么情况下遇到的这个问题呢？

说一点我的看法，不一定对，仅供参考：

non-orthogonality

dca809a8-9cd5-4f74-9329-6d8459966264-image.png

non-orthogonality由snappyHexMeshDict.meshQualityControls.maxNonOrtho直接控制，你的snappyHexMeshDict里这个值是45，这个值一般不需要控制的这么小，65是更常见的配置值。

另外，你的snappyHexMeshDict里，relaxed.maxNonOrtho给的是75，这个值通常是和maxNonOrtho = 65相配合的。如果你真的要把maxNonOrtho设为45，按逻辑来说relaxed.maxNonOrtho也应该相应的减少一点。

网格纵横比

这里截取你的背景网格配置：

vertices ( (-340 -302 -3) // Slightly smaller than the STL model bounds (332 -302 -3) (332 341 -3) (-340 341 -3) (-340 -302 95) // Slightly larger than the STL model bounds (332 -302 95) (332 341 95) (-340 341 95) ); blocks ( hex (0 1 2 3 4 5 6 7) (100 100 100) simpleGrading (1 1 1) // Adjust mesh density as needed );

如上面B老师所说，你网格的大纵横比就是这里导致的。如果不是有意要在Z方向上做加密的话没必要这么分块，我口算一下，大概(70 65 10)就可以。

castellate

从上面可以看到，你背景网格的量级是100 * 100 * 100 = 1,000,000；而你SHM的maxGlobalCells给了2,000,000。这看起来不太对，这几乎没给SHM进一步细分的空间，从结果来看也是这样，网格划分的不够细。

另外，截取一段level的配置：

features ( { file "building.eMesh"; level 3; } ...... ); refinementSurfaces { building { level (3 4); } ...... }

我理解features的level应该不低于Surface的。

胡乱写了一些个人看法，希望能有所帮助。

JCP
让大家膜拜一下。:sunglasses:

0_1468926164745_QQ图片20160719190226.png

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

wedge几何只能用wedge，wedge不能用symmetry。wedge与symmetry的流动特征不一样

https://www.wikiwand.com/en/Centroid

While in geometry the word barycenter is a synonym for centroid, in astrophysics and astronomy, the barycenter is the center of mass of two or more bodies that orbit each other. In physics, the center of mass is the arithmetic mean of all points weighted by the local density or specific weight. If a physical object has uniform density, its center of mass is the same as the centroid of its shape.