CFD中文网

注意格式还有先调试一下吧 这个问题其实最难的是最后的各个方向上的分量处理
而如果是平面各向同性的话 难点在于处理最后源项中针对的方向长度问题
注意源项的单位 你好像并没有关心

放点图有利于讨论..

哈哈哈哈有点搞笑

@东岳 https://www.linkedin.com/jobs/view/postdoc-in-cfd-of-turbulent-environmental-flows-at-university-of-houston-1628249421
我离青千条件远着呢，先在这边混着，等我老婆博士毕业了再说:chigua:

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

顺便测试一下侧边栏回复

@李东岳 感谢李老师的回复，一下子就通透了

@东岳 内网格，尝试着接着划吧

@chapaofan 没啥没啥，下次如果自己解决了也可以自己贴出来奥

见过一些，目前也不好解释。:141:

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

@ljpgbl
这个是全貌
0_1539505430619_Clip_20181014_162251.png

@Hungryandfool
Operator 翻译过来叫做 ”算子“，也就是说，微分方程的空间项可以看作是不同的算子，比如 “对流”（散度），“扩散”（梯度）；
Splitting 的意义在于不同项可以相对独立的去求解，甚至采用完全不同的求解器，这样有利于提高计算效率，减小计算复杂性

@sxz0823 楼上大佬做湍流干了1篇JCP 1篇POF 相当流弊

@东岳 请问东岳老师有关于双欧拉解决相变问题的例子吗，附带算例内没有找到有关相变的设置

@李东岳 他应该测试
关闭超线程，核心选48
打开超线程，核心选48

关闭超线程，核心选96
打开超线程，核心选96

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

@李东岳 我在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)

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

对对对，这样是对的，我之前把应变率和湍动能搞混了，谢谢李老师了！

@子仲无未 大佬你是出现了ill defined primitiveentry starting at keyword actions on line 19吗？怎么解决的啊？求教

TEMOM只展开$(v+v_1)^{1/2}$，任意阶矩展开任意阶都是泰勒的截断误差，是个原式的高阶无穷小，原式也就是$1/2$次方。
但DEMM把$(v+v_1)^k-v^k-v_1^k$也带上展开了，这个截断误差是变的，k越大，误差越大。越高阶矩，误差的量级越大，虽然是都是无穷小。

所以我觉得越高阶这个方法的优势越不明显。当然都是我猜的:135:
在提出DEMM的那篇文章里的数据也有这样的苗头，低阶矩符合很好，高阶的出点问题又修复了。
我把相关资料发给老师。