粘度模型，求解随温度变化的粘度？

[小龙]

我的模型要考虑一下粘度随温度的变化，在OF里我看到了这样的一个模型：ArrheniusViscosity，头文件里给出的解释是这样的

Description
    The Arrhenius temperature-dependent viscosity model multiplies the viscosity
    of a base model by an Arrhenius-type temperature function:

        mu = mu0*exp(k1*(1/(T + k2) - 1/(Tref + k2)))

    Where:
        mu0 is the base-model viscosity
        k1 and k2 are Arrhenius coefficients
        T is the local temperature
        Tref is the reference temperature

这个模型是考虑了温度对粘度的影响的模型吗？里面的k1和k2解释是 Arrhenius系数，这个东西是什么？这个模型怎么用啊，没有找到OF里的自带算例。

[小龙]

还有一个粘度模型，power-law non-Newtonian viscosity model，

    // Private data

        dictionary powerLawCoeffs_;

        dimensionedScalar k_;
        dimensionedScalar n_;
        dimensionedScalar nuMin_;
        dimensionedScalar nuMax_;

        volScalarField nu_;

这里的n，nuMin_和nuMax_是什么？

[东岳]

这都是非牛顿流体的本构方程，实际上非牛顿流体的书上都有的，我建议你参考Non-Newtonian Fluids: An Introduction。有问题再回帖 :big_mouth:

[小龙]

@李东岳 东岳前辈，弱弱的问一下，考虑水的粘度随温度变化，还是牛顿流体吗？:big_mouth:

[东岳]

牛顿流体定义为剪切力和形变率为线性关系。如果粘度只随着温度变化，剪切力和形变率依然为线性关系，但却不是点对点的对应关系。从概念上我不确定是不是牛顿流体。但如果你要编程植入的话，这个非常简单。吧粘度设置为温度的参数就可以了。可以基于非牛顿流体class进行修改。

[yuanlee2011]

@东岳 看过一个chalmers的例子，是基于powerLaw类进行的修改，我想知道可否基于Newtonian进行修改？

[东岳]

是基于powerLaw类进行的修改，我想知道可否基于Newtonian进行修改？

这俩个类大同小异，个人觉得没什么区别。

[yuanlee2011]

@东岳 基于powerLaw类，自定义了一个粘度模型，能正常编译，但是运用至求解器中后会出现这样的错误，请问是什么原因？

[东岳]

看你的TypeName("tempdeppowerLaw");有没有写

[yuanlee2011]

@yuanlee2011 您指的是这里吗？（基于of231）

[yuanlee2011]

@东岳 您指的是这里吗？（基于of231），

[yuanlee2011]

@东岳 新的粘度模型编译成库(libusertempdeppowerLaw.so)，通过controDict中添加“libs (“libusertempdeppowerLaw.so”);”调用，当使用串行计算时，case正常运行，但是采用并行计算，在decomposePar的过程中会出现该warning(如图所示)，请问该怎样解决？

[东岳]

decomposePar只要分解了就可以，这个可以忽略

[yuanlee2011]

@东岳 谢谢！