在编译求解器时引用弹性项

[冠竹]

想在求解器中直接使用粘弹性流体模型，但是发现使用divTau引用弹性项的时候产生了报错，说明在of8中的粘弹性流体模块里弹性项不是这么引用的，所以想问问各位大佬，of8中弹性项应该如何引用呢？

[sungda]

您好！想问下您是在做涉及粘弹的两相流吗？

[冠竹]

目前还在编译求解器，但我目前应该是只涉及到粘弹性流体，不涉及到两相流

[sungda]

@冠竹 单相流的话，foam-extend4.0中好像有一个自带的求解器viscoelasticFluidFoam

[冠竹]

@sungda 十分感谢您

[guangzhou]

@冠竹 你好，请问你说的这个求解器算的准吗

[冠竹]

@guangzhou 这个求解器在高维森贝格数下会产生稳定性问题，所以一般采用对数重构法去解决这个问题。论文可以参考张红娜老师在17年的一篇论文《Numerical Simulation of Heat Transfer Process of Viscoelastic Fluid Flow at High Weissenberg Number by Log-Conformation Reformulation》，或者可以看看rheoTool工具箱相关介绍

[guangzhou]

好的，万分感谢！
其实我viscoelasticFluidFoam和rheotool都用了，但结果发现viscoelasticfluidfoam结果更准确，rheotool则完全不对，因为我模拟的条件不能忽略雷诺数， Re和Wi差不多大，两个的值都在3-50左右，惯性的影响不能忽略。rheotool是不是只适用于惯性被忽略的情况啊，因为我看里面的案例都是不考虑Re的影响

[冠竹]

@guangzhou 这个不影响吧，我用rheoTool的目的是他的粘弹性流体库。他的求解器对我来说限制有点多，所以我就直接把他当库文件用自己新写了一个DNS求解器。如果觉得他的不适合的话就引用库自己改一个就好了啊

[guangzhou]

@冠竹 好的，谢谢回复，我刚学Openfoam没多久，还是个菜鸟，
自己改求解器还不会，我再好好学学吧