@李东岳 感谢李老师的回复，根据查阅的资料：this->nuEff() = nu() + nut() interFoam求解器使用的this->nuEff()应该是laminarModel下的Stokes.C这个模型，因为在求解器计算时输出器选用的stress model 是Stokes model，那divDevTau()函数下的this->nuEff()就如Stokes.C内所描述如下，就是指向nu()并不是为this->nut() + this->nu()，这里就很疑惑，这里将nu()换成自己设定的Nu(T)，就得不到合理的计算结果（流场分布）。此外，我选用的是层流在laminarModel.C内也将nut初始化为0,如下代码，所以不知道interFoam这里调用了哪里的的nuEff()，可能不是目前我认知的Stokes.C内部的:141:

template<class BasicMomentumTransportModel> tmp<volScalarField> Stokes<BasicMomentumTransportModel>::nuEff() const { return volScalarField::New ( IOobject::groupName("nuEff", this->alphaRhoPhi_.group()), this->nu() ); } template<class BasicMomentumTransportModel> Foam::tmp<Foam::volScalarField> Foam::laminarModel<BasicMomentumTransportModel>::nut() const { return volScalarField::New ( IOobject::groupName("nut", this->alphaRhoPhi_.group()), this->mesh_, dimensionedScalar(dimViscosity, 0) ); }

@samuel-tu
你好这个问题解决了吗？我也是遇到了同样的问题。

@李东岳
李老师好！抱歉回复晚了。
算例文件在这，https://www.jianguoyun.com/p/DejQRSIQ-4XWCRiVid0EIAA，稍微有些大，压缩包16.8MB

@zy-gong 好像确实可以，太惭愧了，总是问这种低级问题，学习之路任重而道远
谢谢李老师了:xiexie:

@李东岳 不好意思老师回复晚了，是openfoam，老师

@bestucan 感谢前辈

blastFoam 好象有个公司专门做这个。所以，看到效果好的基本都是公司放出来的，开源版本很优秀的话，公司就挣不到钱了:chigua2:

pointMotionU和pointDisplacement这些控制动网格节点的边界条件文件的数据类型是pointVectorField，对应的patch类是Foam::pointPatch而不是Foam::fvPatch，这应该就是报错的根本原因。

我的思路如下：

// 获取pointPatch的ID，它和相同边界上的fvPatch的ID是一致的 label patchIndex = patch().index(); // 在objectRegistry随意获取一个存在的物理场 const volScalarField & p ( this->db().objectRegistry::lookupObject<volScalarField>("p"); ); // 从上述物理场中访问fvMesh中的体心/面心场 const vectorField& Cfp = p.mesh().Cf().boundaryField()[patchIndex];

如题，字数补丁

https://www.openfoam.com/documentation/guides/latest/doc/guide-fos-field-wallHeatFlux.html
https://www.cfd-online.com/Forums/openfoam-post-processing/91755-calculating-heat-flux.html
https://www.cfd-online.com/Forums/openfoam-post-processing/237465-how-get-wall-heat-flux-using-bouyantbossinesqsimplefoam-openfoam-v2006.html

@李东岳 李老师，你好！我是兰大的，我们正常开学的，除了疫情严重的省份地区不允许返校，其他的正常。

fvc::grad(T).component(2);

@李东岳 多谢李老师！

一般放外面就行，如果要放pimple.loop()里的话，标量方程的时间步长需要除以nOuterCorrector。

yum install centos-release-scl -y yum install devtoolset-9 -y source /opt/rh/devtoolset-9/enable

OpenMPI aclocal

cd /usr/local wget https://ftp.gnu.org/gnu/automake/automake-1.15.tar.gz --no-check-certificate tar -xzf automake-1.15.tar.gz cd aclocal-1.15 ./configure --prefix=/usr/local/automake make make install export PATH=/usr/local/automake/bin:$PATH cd openmpi-4.1.2 autogen.pl ./makeOPENMPI yum install blas-devel lapack-devel wget https://ftp.mcs.anl.gov/pub/petsc/release-snapshots/petsc-lite-3.17.2.tar.gz --no-check-certificate ./makePETSC -hypre=http://dyfluid.com/hyper-2.24.0.tar.gz

@李东岳 谢谢李老师，问题已经解决了，对于DNS的槽道流需要用统计数据，流场计算时间不够就很难对得上，还需要等到湍流完全发展后再开始计算平均场。

@李东岳 我找到如何输出endTime方法了，用this->owner().db().time().endTime().value()就可以。但是放在判断语句中出现如下错误：

note: ‘Foam::word’ is not derived from ‘const std::fpos<_StateT>’ if ((this->owner().db().time().timeName()) == (this->owner().db().time().endTime().value())) /usr/include/c++/5/bits/stl_iterator.h:349:5: note: candidate: template<class _IteratorL, class _IteratorR> bool std::operator==(const std::reverse_iterator<_Iterator>&, const std::reverse_iterator<_IteratorR>&) operator==(const reverse_iterator<_IteratorL>& __x, ^ /usr/include/c++/5/bits/stl_iterator.h:349:5: note: template argument deduction/substitution failed:

李老师能指点一下么？

感谢分享！！

@winsway_zero 嗯嗯，我懂了，谢谢，renumberMesh是对网格的操作，setFields也要在renumberMesh之后操作，不然也会乱掉:146:

你这个哪个是进出口，重力是向下的么，

alpha场也发一下