最近在做一个污染物浓度扩散的算例，需要选定一条线并将其设置为污染物浓度扩散的线源，但在setFields中好像没有选定线的设定，所以哪位大佬要是知道请告知，万分感谢。

你说的把Q的插值变成比值，应该是Ω方法吧？？？
Liutex-Ω的话好像比较麻烦

@西湖冷月 你具体编程过吗？

@yfleung @荷兰后腰阿贾克斯 @波流力
大佬朋友，请问你们处理400万以上网格，速度很慢，无法出结果的问题，解决了吗？
感谢

@队长别开枪 明白了，谢谢老师！

@oitocfd 好的，我尝试tecplot打开了其他格式是binary的算例，发现可以打开，那么就说明是之前的算例的问题，我再修改试试。谢谢您。

@星星星星晴 谢谢，我再看看

@西极若木
这个算例中在Stokes2里的angle是不是可以设置波浪来的角度？也就是模拟斜浪

@tens 谢谢您的热心指点，虽然还有点晕，我后面仔细体会体会:quwan: ，再次感谢:xinxin:

@sungda 这是二维的自适应网格吗？是怎么实现的

单核计算：
ff8a928a-d108-454d-8548-a662114c6a34-image.png
并行计算：
2674cfa8-2f86-410a-9ebc-144e98f89800-image.png

@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
气液界面是不断变化的，请问红色方框内的边界条件如何添加~

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

请教各位同学和老师，我在尝试用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

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