@mania 在 请教各位前辈OF中湍流入口DFSEM的设置问题 中说：

U、R、L和point

U、R、L和point可以通过RAS计算得到

@mania 你用RAS计算得到U、R、L和points，直接放到主程序计算就可以。of2006有自带的算例

@tens 好的谢谢:146: 我试一试

@jqbu 基于OpenFOAM的甲烷湍流燃烧LES数值研究_龙文增;不知道这篇文章能不能帮助到你

最近几天调程序老是对不上，今天才发现只有标准的结构性网格才是这样，如果是非结构网格，或者SHM生成的网格在过渡层处这个数量不是固定的
还是要细心点:136:
d0883aaa-8a7f-4939-9a9f-cf30552dfc9a-image.png

额 好像发现个问题，质量守恒方程应该是这个
7973d26e-efef-4af5-9e55-686a447c4996-image.png

其中𝑆(𝜌)与 𝑆(𝛾)形式一样，只不过自变量由相体积分数变成了密度，不知道这样对不对？

感觉有点弄乱了，求大佬指点，不胜感激:135:

@小考拉 自问自答来解决这个问题。

首先定义一个全为0的流场命名为c_flow

然后利用下面的命令

c_flow.boundaryField()[patchID] = c_flow_origin.boundaryField()[patchID]

最后输出c_flow流场即可:chitang:

@鲸落 请问你解决了吗~动态读取cellZone~

众所周知，当我们在进行并行计算的时候，需要把网格进行剖分，OF中可以使用第三方库进行区域分解，可是我们知道，OF中的区域分解是不带重叠层的（或者说没有halo单元），也就是，同一个网格仅可能属于一个处理器。
那么问题来了，我现在想实现带有halo层的重叠型的区域分解，我该怎么实现呢？或者说，有没有那位大佬做过类似的工作或者相应的资料吗？跪谢了！

@hongjiewang 在 温度边界条件 中说：

1632279639(1).png
气液界面是不断变化的，请问红色方框内的边界条件如何添加~

可以识别出气液界面，但是想知道如何在同一个边界里使用第二类和第三类边界条件

@prometheus10 然后，文中提到的另一个链接:papa:
https://www.cnblogs.com/liusuanyatong/p/11259653.html

请教各位同学和老师，我在尝试用MPPIC在一个存在斜面的三维复杂网格中模拟固液两相流时，计算一直没办法收敛。尝试不同的边界条件和不同的k和epsilon值，都没有任何改善。以下是stack trace：

Lookup solver for epsilon.waterFinal [16] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [16] #1 Foam::sigFpe::sigHandler(int) at ??:? [16] #2 ? in "/lib64/libc.so.6" [16] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? [16] #4 Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > const&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:? [16] #5 Foam::RASModels::kEpsilon<Foam::PhaseIncompressibleTurbulenceModel<Foam::singlePhaseTransportModel> >::correct() at ??:? [16] #6 ? at ??:? [16] #7 __libc_start_main in "/lib64/libc.so.6" [16] #8 ? at ??:? [cas019:410749:0:410749] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2)) ==== backtrace ==== 0 /lustre/opt/cascadelake/linux-centos7-skylake_avx512/gcc-8.3.0/ucx-1.6.0-dacxhp7gepme2j7jmvpcmbm2qrptc6k7/lib/libucs.so.0(+0x1bf44) [0x2b968ef64f44] 1 /lustre/opt/cascadelake/linux-centos7-skylake_avx512/gcc-8.3.0/ucx-1.6.0-dacxhp7gepme2j7jmvpcmbm2qrptc6k7/lib/libucs.so.0(+0x1c2ca) [0x2b968ef652ca] 2 /lib64/libc.so.6(gsignal+0x37) [0x2b9689945337] 3 /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6_+0xd5) [0x2b9688a0dca5] 4 /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4FoammlINS_7volMeshEEENS_3tmpINS_16DimensionedFieldIdT_EEEERKS6_RKS5_+0x164) [0x2b9683ba5d74] 5 /lustre/home/acct-naozzn/naozzn-user1/openfoam/OpenFOAM-5.0/platforms/linux64GccDPInt32Opt/lib/libDPMTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_34PhaseIncompressibleTurbulenceModelINS_25singlePhaseTransportModelEEEE7correctEv+0x547) [0x2b9687558727] 6 MPPICFoam() [0x433d76] 7 /lib64/libc.so.6(__libc_start_main+0xf5) [0x2b9689931505] 8 MPPICFoam() [0x4352e0]

边界条件如下表

f3ef028e-2057-43c0-8e87-fe42f54246e7-image.png

感觉应该是颗粒太大或者是网格划的不太好，有没有大佬能帮我看看怎么调整更准确一些

各位老师，大家好，想请教一下，在拥有管理员账号的前提下，怎样在联想集群的LSF作业调度系统中添加新的普通账户？
之前是按照这样的方式，也忘记了是在哪里找到的了：:135:

//在管理节点上添加新用户(root)： useradd -d /gpfsdata/home/cfdem -m cfdem //在管理节点上给新用户添加密码(root) passwd cfdem //添加到用户组 vim /gpfsdata/apps/lsf/conf/lsbatch/cluster1/configdir/lsb.users //添加完后运行 badmin ckconfig //在计算节点上添加新用户 ssh c01n01 useradd -m cfdem passwd cfdem并输入密码 //同理在c01n02~c01n03同样上述操作 //ssh到管理节点mgt01并切换到cfdem ssh mgt01 su cfdem ssh-keygen//然后一直回车

但是这种方式感觉有点繁琐。:136:
后来了解到PBS管理系统里面可以直接通过命令

clusconf -au 账号名

就可以完成普通账户的添加。
现在想问一下，LSF管理系统添加账户的时候，有没有类似于PBS管理系统的这种一气呵成的简单命令？:146:

@冠竹 这个是$\overline{u'u'}$, $\overline{u'v'}$等一共6个，就是解析的雷诺应力除以密度

@五好青年 你可以参考一下Holzmann CFD官网的算例, 关于Magnus Effect的算例中有加入圆柱转动函数条件的相关设置.

目前我只会做到在边界某段添加固定温度，整个边界都使用热流密度条件也是可以的，请问如何在边界的某段添加热流密度边界~比如下面这段0.000333<=x<=0.000666使用热流边界 应该如何修改~

OUTLET { type codedFixedValue; name heated; value uniform 300; codeInclude #{ #include "fvCFD.H" #}; codeOptions #{ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I$(LIB_SRC)/meshTools/lnInclude #}; //libs needed to visualize BC in paraview codeLibs #{ -lmeshTools \ -lfiniteVolume #}; code #{ /*const IOdictionary& d = static_cast<const IOdictionary&> ( dict.parent().parent() ); const fvMesh& mesh = refCast<const fvMesh>(d.db()); const label id = mesh.boundary().findPatchID("OUTLET"); const fvPatch& patch = mesh.boundary()[id];*/ scalarField T(patch().size(), 300); forAll(T, i) { const scalar x = patch().Cf()[i][0]; //const scalar y = patch.Cf()[i][1]; //const scalar z = patch.Cf()[i][2]; if(x>=0.000333||x<=0.000666) { T[i] = 400; } else { T[i] = 300; } } (*this)==T;//T.writeEntry("", os); #}; }

@队长别开枪 谢谢老师指点，我好好研究一下

@七辆战车
你好！ 请问一下，利用rigidbody求解器，有实现双浮体吗，是直接在dynamicdict文件下给两个浮体的参数吗？ 对于外部导入的水平圆柱浮体，刚体类型可以给rigidbody吗？ 期待你的回答，不胜感激。

@队长别开枪 谢谢老师您的指点，已经安装上了。但我是用hypremesh画的斜坡网格，是否可以不使用stl文件来进行相初始化。

@tens 谢谢啦