CFD中文网

@bestucan 谢谢老师，我试了下，还是不行~

他是由于none量 他把它删了，看来我的问题就是那个omega，只要用

const volScalarField& omega_ = ( this->mesh_.objectRegistry::template lookupObject<volScalarField>("omega") );

就会找不到，估计是结构问题，难道只能重写了。

@Chen_hao
当这个模型用于多相流时，IOobject::groupName("omega", alphaRhoPhi.group()) 会变成类似 omega.air 这样的名字。

你可以直接采用湍流模型来算试一下。用湍流模型计算层流，只要湍流粘度足够小，影响就很小。

openFoam可以无量纲的，应该在etc/controlDict这个文件里里吧dimensionset 1改为0

@chen_hao
解决了，改成了无量纲，就避免出问题：

volSymmTensorField Sij(twoSymm(tgradU())/omega_); volTensorField Wij(2.0*skew(tgradU())/omega_);

@李东岳
针对 ncompressible::turbulenceModel 模板特化后的 ShihQuadraticKE 模型 ，仔细研究了一下其怎么实现非线性项的实现，发现其是在 NonlinearEddyViscosity.H 中定义 nonlinearStress_ 项，然后用 ShihQuadraticKE 作为子类继承父类 NonlinearEddyViscosity 以便在 ShihQuadraticKE.C 中直接引用omega_。

但由于我修改的是SSTkomegabase.H( 继承于TurbulenceModel类 )与 linearViscousStrss.H ( 继承于BasicturbulenceModel原则类 ) 没有特例化，不能简单类似上面的ShihQuadraticKE 模型。需要从新仿造框架，是一项大工程！

我的策略是：在现有的文件上直接硬植入，简单很多，毕竟涉及到湍流模型的类错综复杂，针对找不到omega，按照我的分析应该是对的，运行求解器时linearViscousStrss.C先运行此时并未运行SSTkomegabase.C 所以没有omega场注册，找不到他，所以我在creatfield.H中IOobject一下omega场，目前可以算了，结果还未检查，不知道我的理解有没有问题。

@李东岳
谢谢老师，大概理解this用法，这个是指向name()，还以为这个命名是起到什么作用，既然一般情况下没区别，那我就不深究了。

@chen_hao

const volScalarField& omega_ = this->omega_;

这一种

@李东岳
谢谢老师，虽然只是一个小建议和思路，却给我节约了大量的时间，目前已经在OF7中调通了mutiphaseInterFoam求解器。:xiexie:

可以。我记得湍流方程右边都是只考虑网格点的作用，不需要考虑边界。因此你的G声明了内部场，没有问题。