CFD中文网

@xpqiu 感谢老师详细指点，试了下这个方法确实简洁

@浪迹天大 感谢回复~如果是二维的可以用白色隐藏，但三维的应该不行了

我刚试出一个方法，先通过Calculator计算出速度的大小作为标量结果，因为Threshold貌似只能用标量分割，然后使用Threshold根据上述结果选择区间，这样可实现。如有更简单的方法还请各位老师指点~

@shepherd 感谢回复~ 将下限提高后，速度小的网格部分没有隐藏掉，而是和区间内最小速度一个颜色

各位老师好，paraview能否显示一定速度区间的网格部分？比如一个算例速度较大的部分很少，想单独显示这部分，而隐藏速度小的网格部分，不知在paraview能否实现？

@东岳 多谢老师！

@东岳 谢谢东岳老师，那就是说这个值设置0还是1没影响么，我还以为是和alpha有关系

各位老师好，关于interFoam的算例capillaryRise 0/alpha.water 中 constantAlphaContactAngle这个边界条件是怎么工作的呢？里面“value”这个关键词是什么意义？网上资料没有找到，而代码没看懂，请老师指点

 walls
   {
       type           constantAlphaContactAngle;
       theta0         45;
       limit          gradient;
       value          uniform 0;
   }

各位老师好，我想在OF里植入一个与流动直径有关的边界条件。通过sHM后得到的网格，其各流动通道直径不一且比较复杂，能否获取各网格所在流动通道的当量直径呢？在网上查阅wallDist(mesh).y()可以计算网格到壁面的最近的距离，但不知道能否用来计算出直径。

@bestucan 非常感谢老师的详细建议。因为集群之前就已经安装过OF，我是在个人账户里安装OF7版本的，所以可能就差这一个包。我已经使用了跑了计算，目前来看没什么问题

@xpqiu @bestucan 非常感谢两位老师的建议！
我的目的是：在虚拟机内自定义修改了一个求解器，想放在集群里编译使用。试了两个方法似乎都编译成功了，我把我的步骤简要说一下，如果有什么错误请各位老师指正：
1、用集群的OpenFOAM编译。之前不敢用集群的OpenFOAM编译自己个人文件夹的求解器，因为我不知道 $FOAM_USER开头的这些环境变量具体的是什么目录，在终端键入env后就可以看到这些环境变量了，其实用集群的个人账户登录，这些目录都在自己的个人文件夹里，所以接下来和在虚拟机操作类似了。
2、在集群个人文件夹内安装OpenFOAM（如果想用其他版本的OF的话）。开始报错说我缺少flex，我找集群管理员要了root账号，因为我们集群的系统是Centos，所以我键入yum -y install flex 来安装flex。flex安装成功后，登陆个人账户，在终端键入env后可以看到 $FOAM_INST_DIR是在我个人文件夹内的目录。编译OF，经过了很长时间编译成功，接下来对自定义的求解器编译，也类似于在虚拟机操作的步骤了。如果是不同版本，编译自定义的求解器会有代码语句上的一些不同提示报错，需要修改下。

@bestucan 感谢您的回复，我没有超级用户权限。那如果不在个人用户文件夹安装OF，将在虚拟机修改的求解器拷贝到个人文件夹，使用集群的OF能否编译我个人用户文件夹的求解器？这样的话环境变量需要怎么修改呢？

各位老师好，我在虚拟机内自定义修改了一个求解器，想放在集群里使用。因为在虚拟机里自定义的文件是放在OpenFOAM的run里的，所以想在集群的个人文件夹内安装一个OpenFOAM。不知能否在集群的个人文件夹内安装OF。我将OpenFOAM-7和ThirdParty-7下载并解压，配置环境变量，然后编译OpenFOAM7，报错如下：

flex wmkdep.l; \
gcc -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -O3   -fPIC lex.yy.c -o /gpfs/xxx/OpenFOAM/OpenFOAM-7/wmake/platforms/linux64Gcc/wmkdep
/bin/sh: flex: command not found
gcc: error: lex.yy.c: No such file or directory
gcc: fatal error: no input files
compilation terminated.

请问是哪儿出了问题？我看了一下集群的gcc版本是5.3.0

@东岳 谢谢东岳老师！

@东岳 东岳老师您好，我按照上述步骤运行，会报错Cannot find patchField entry for fixedWalls_xxx （xxx是我stl中solid的命名）。

我把上述decomposePar改为decomposePar -copyZero，运行会报错Cannot find patchField entry for procBoundary6to0等等。后来我在0文件夹的各初始条件文件中添加了

"proc*.*"
 {
    type processor;
 }

成功运行了。

请问老师这是不是正确的方法？网上也有写加上value uniform的，比如下面，两者有区别吗？

"proc.*"
{
     type processor;
     value uniform 0;
}

@东岳 好的，谢谢东岳老师~

各位老师好，使用sHM并行后要先reconstructParMesh再使用求解器计算吗？还是先使用求解器计算再reconstructParMesh和reconstructPar？
举个栗子，如下两个步骤哪个正确？

blockMesh 
decomposePar 
mpirun  -np 24 snappyHexMesh -parallel -overwrite 
reconstructParMesh  -mergeTol 1e-6 -constant
mpirun  -np 24 simpleFoam -parallel 
reconstructPar  

blockMesh 
decomposePar 
mpirun  -np 24 snappyHexMesh -parallel -overwrite 
mpirun  -np 24 simpleFoam -parallel 
reconstructParMesh  -mergeTol 1e-6 -constant
reconstructPar  

@东岳 老师，应该放在什么位置呢，我放在了crateFields.H 里面新定义的量下面

Info<< "Reading field nuk\n" << endl;
volScalarField nuk
(
    IOobject
    (
        "nuk",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    ),
    mesh
);
  
   const fvPatchList& patches = nuk.mesh().boundary();
   forAll(patches, patchi)
   {
        scalarField one(3.0, nuk.mesh().boundary()[patchi].size()); 
        nuk.boundaryFieldRef()[patchi].patchInternalField() = one;
   }

还是这种问题，这是为什么呢？

@东岳 老师您好，我尝试了您给出的这行代码

 scalarField one(1.0, nu_.mesh().boundary()[patchi].size());
 nu_.boundaryFieldRef()[patchi].patchInternalField() = one;

成功编译，但是运行时报错

Attempted assignment to a const reference to an object of type N4Foam5FieldIdEE

    From function void Foam::tmp<T>::operator=(const Foam::tmp<T>&) [with T = Foam::Field<double>]

是不是不能这样赋值给patchInternalField()？应该怎么修改呢?

@东岳 老师，这个公式（1）是他们参考另一篇文章
Wu, K., Chen, Z., Li, J., Li, X., Xu, J., & Dong, X. (2017). Wettability effect on nanoconfined water flow. Proceedings of the National Academy of Sciences, 114(13), 3358-3363. https://www.pnas.org/content/114/13/3358.short
我的理解是因为在纳米管内靠近壁面的位置，粘度变化没有统一的解析解，所以看到一些文章就这么分成“多层”粘度，算是粗略简化的方法吧