CFD中文网

@李东岳 李博士，

请问东岳流体官网的塔式工作站（DM0812， DM3.1）cpu 具体参数，像型号，主频，睿频等可否简要分享下。

谢谢！

@是嘿嘿牛啊

这个系泊是自己开发的吗？
看着不像悬链线，实际实验也是这样系泊的吗？如果不一致，这可能是你数值有偏差的原因。

@blookcfd 感谢您的回复，
请问这种语法比较常见吗？
这个objectRegistry::template是指代的 Type类型吗？

在fvPatchFvMeshTemplates.C 中发现一个很怪异的函数调用，

boundaryMesh().mesh().objectRegistry::template

找了半天没找到代码定义在哪里，请教各位大佬指点。具体代码如下：

   30 template<class GeometricField, class Type>
   31 const typename GeometricField::Patch& Foam::fvPatch::lookupPatchField
   32 (
   33     const word& name,
   34     const GeometricField*,
   35     const Type*
   36 ) const
   37 {
   38     return patchField<GeometricField, Type>
   39     (
   40         boundaryMesh().mesh().objectRegistry::template
   41             lookupObject<GeometricField>(name)
   42     );
   43 }

可以参考邱小平大佬的博客
链接文本

@Zhy2022 因为tmp<volScalarField> 是个指针，而你的pc是个volScalarField 变量.

@李东岳 在 solve函数为啥不需要用对象去访问调用？ 中说：

应该用(Ueqn==fvc::grad(p)).ref()->solve()，这是一个tmp类型

@wangfei9088 大佬太厉害了

这里.ref() 是取引用吧，所以应该是

(Ueqn==fvc::grad(p)).ref().solve() 或者

(Ueqn==fvc::grad(p))->solve() ？

@Micro 大佬试过李东岳老师修改后的代码了吗？

列表在小册子第37页，公式（3.33）上面那句话：
“则pEqn.flux()即为方程(3.30)中的下划线部分”

根据上下文语义，应改为：
“则pEqn.flux()即为方程(3.32)中的下划线部分”

@东岳 谢谢李老师，

您可以提示一下，这个“==”会更新边界条件的定义是在哪块代码可以找到的吗？

今天看到李东岳老师发的朋友圈，想起昨天学习interFoam.C的时候,发现这样的代码：

 p == p_rgh + rho*gh;

请问此处“==”是否表示赋值？可以把“==” 换为“=”？
如果“==”不表示判断，是否还有别的用途？

• IncompressibleTurbulenceModel.C 中的代码：
  首先给出devReff() 和devRhoReff()的代码，如下：

template<class TransportModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::IncompressibleTurbulenceModel<TransportModel>::devReff() const
{
    return devRhoReff();
}

...

template<class TransportModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::IncompressibleTurbulenceModel<TransportModel>::
devRhoReff() const
{
    NotImplemented;

    return devReff();
}

可以看到，devReff()调用了 devRhoReff()。然而在IncompressibleTurbulenceModel.C中，并没有找到这两个函数更进一步的细节

• 接着，在linearViscousStress.C 中的的代码中，
  发现devRhoReff()的定义

template<class BasicTurbulenceModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::linearViscousStress<BasicTurbulenceModel>::devRhoReff() const
{
    return tmp<volSymmTensorField>
    (
        new volSymmTensorField
        (
            IOobject
            (
                IOobject::groupName("devRhoReff", this->U_.group()),
                this->runTime_.timeName(),
                this->mesh_,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            (-(this->alpha_*this->rho_*this->nuEff()))
           *dev(twoSymm(fvc::grad(this->U_)))
        )
    );
}

• 先将问题总结如下：
  1）NotImplemented是什么意思？为什么不可压中devRhoReff()返回的是devReff()？
  是否可以理解为，不可压缩流里面的devReff()执行的是其它代码（linearViscousStress.C）中的devRhoReff()？

template<class TransportModel>
Foam::tmp<Foam::volSymmTensorField>
Foam::IncompressibleTurbulenceModel<TransportModel>::
devRhoReff() const
{
    NotImplemented;

    return devReff();
}

2）若devReff() 实际上执行的是devRhoReff()（linearViscousStress.C）

并且在CFDonline一个帖子里，链接文本，有人说devReff()实际上是devRhoReff()中rho为1时的情况
''it seems like it's best to simply search for "devRhoReff()", because "devReff()" is just when "rho" is simply 1.0.''。
若是如此，这里面rho_=1是在哪里定义的，我一直没找到相应代码

望各位不吝赐教 ：）！

@东岳 ，李老师您好。
今天学习了icoFoam里面的解析，发现两个不理解的地方，列出如下：
1）由公式（27）和（18）相加得到公式（28），在公式（28）下面提到： p*=p(r)+p'; 但是在公式（18）中应为p(n). 不知是否可以理解为p*=p(n)+p'。

2）公式（32）中的（1/A（P，f））是否应改为1/A（P）（和公式（26）对应）

@李东岳， 李老师今天我在最新的用户手册里发现一个小错误。在PDF26页，2.1.1.4矢量图小节下，第一段第一句话“在我们绘制矢量图之前,我们有必有**屏蔽其它模块,...” 应该改成必要吧。

我遇到的问题是case在粗网格下单核以及并行没有问题，但是加密网格之后只能在单核下跑，一并行就出问题，第一个时间步还没迭代完库朗数就一下很大了，然后程序就崩了。

我试过下面，都没有成功：
减小时间步，在同学建议下调到了10e-6，（原来是dt=0.005s）
增加PIMPLE中相关的迭代次数
设置了库朗数控制的adjustableRuntime, maxCo=0.5, maxAlphaCo=0.5

求教是不是需要改fvsolution或者fvscheme中的相关内容？或者别的建议？
我用的是FoamExtend 4.0中的fsiFoam.

@东岳 所以要不要买几套房，万一是下一个雄安新区，还可以在边界趁着北京移动的信号。

张家口人表示此贴必火@东岳 ，100块算到天荒地老

@youmengtian 明白，多谢

@东岳 收到，谢谢大佬

各位大佬，小弟openFOAM刚入坑 ，正积极向CFDer的dang组织靠拢。

想请教下这样配置的电脑有什么推荐，如下：

cpu : intel？至强？
主频 : dummy
核数 : 20；
内存 : 32G， //不知道够不够用？
硬盘 : ？？ //想装win10和Ubuntu双系统，不知道要求大不大？
品牌 : Dell //据说戴尔的比较好

欢迎大佬们分享自己电脑的配置