CFD中文网

@东岳 如此说来，除了加密网格，没有好的解决方案啦:zoule: :zoule:

论坛最近迁移服务器，才发现搜索功能丢了，我去找找，铁铁
先用 @random_ran 的方法

得看人是不是明白人。明白一口口饭都是从哪来的，明白jio底下站的地为什么这么结实，明白前方的路为什么这么清晰。

碰上明白人，给高了会用富裕的资源提升自己（给多少钱干多少活，否则钱拿不稳）；给低了也不担心，毕竟瞎还能活着的概率不是很大，找个不瞎没有特别的困难。

依靠某个人、某个环境下对自己的态度、待遇来判断自己的价值。这显然不是个明白人……给高了骄傲，给低了还怨恨呢。不能自知，怎么都是一件很危险的事情。有核心能力的也有不明白事的：仰天大笑出门去，归来依杖自叹息:136:

而能碰上多少明白人又取决于所在行业对社会精英的吸引能力和自己在行业的层次。总之不努力就整天和糊涂蛋扯吧……:135:

还有这种组织心态，把员工当小孩、追胡萝卜的驴子（有些事情你不需要知道，我在掌着。我要塑造你的认知、影响你的判断）。当然在这个繁复多样的世界，这种方式无法完全避免。但是显然有更好的方式，只是对组织中的每个人要求比较高。

@李东岳 找到了 https://vimeo.com/148670249

您好 您后来展向平均怎么做的呢

@jdzhang 你的颗粒会动么？不会动，为什么颗粒间还有距离？会动，那对着“多相流”查比对着“多孔介质”查更容易找到解决方案些。

这有一个画网格的：
https://www.bilibili.com/video/BV1eS4y1F7Wm/?vd_source=c01e8d4fff0db460c75b32e0d0191b24

这有些相关的：
https://www.researchgate.net/publication/308996206_Pore-scale_numerical_simulation_of_mud_erosion_in_the_subsea_sand-mud_alternate_layer_using_lattice_Boltzmann_method

https://ascelibrary.org/doi/10.1061/(ASCE)HY.1943-7900.0001489

@刀尔東 在 向各位大佬请教一下网格划分的基础问题 中说：

感谢各位大佬不吝赐教。

链接: https://pan.baidu.com/s/11z9O_mzDLuqEp30mWRr5Ow?pwd=1234 提取码: 1234 复制这段内容后打开百度网盘手机App，操作更方便哦
--来自百度网盘超级会员v6的分享

以上是我的两个文件

网格文件为ENSIGHT GOLD FORMAT

压差是分析很多工程设备的关键参数，在之前，主要都是针对多孔介质的Ergun模式来处理这一块。
也跟试验对比很好。但是，最近新开发的一个模型，虽然与试验的结果对比比较好，但感觉总是有问题。
于是不知道CFD中静压值到底是怎么得出来的？
对于一块多孔介质好像比较好处理，但是如果这个时候多孔介质表面沉积了灰，这个时候这块的整个压降会升高，然后灰越堆越多整个表面的压降都会升高。然而，这个局部压降和表面灰的分布十分相关，造成整个表面的这个灰分的压降相差一千多Pa。而最后一般的检测压降点，都是压差变松器的两个接口，一般都是设置在容器外壁的。但是如果我模拟检测点和试验一样的时候，这个时候的压差比我的表面计算压差平均值要小得多。（压差是1200多Pa，但是我的多孔介质有800多Pa，加上平均粉尘层的压降有640多Pa，合起来应该有1400多Pa，而不是1200多Pa.）
这就产生了一个问题，这个通过N-S方程求解的压力项到底是怎么处理的?
当他通过一层多孔介质的时候，多孔介质内部有压力源项，而在外部没有，那么当他气体从多孔介质外向内流通的时候，他的这层压力和速度是怎么计算的？如果这个时候我再添加一个源项在多孔介质层，而这个源项在多孔介质外是没有的，那么这层额外的源项是不是也会通过增加流体的静压来表达呢？（通过伯努利方程似乎很好理解，但是N-S是怎么实现的，他的这个\Delta P 是怎么获得的？）

@xuanze 3374769460@qq.com

@guanghuilineu 那看起来发散是个很正常的结果，如果相变区速度不为0时，也变速也变密度。

如果只变速，和周围空间交换能量，即动能转换为其他形式的能量，没问题。
如果速度和密度一起变，动能是依赖于质量的。这时候，质量也变了，不守恒了，就崩了。

消失的那部分质量，带有的动能如何传递出去。

压力就是传递能量的桥梁。在传递的时候一头突变了，就断了。可不可以先使相变区域停止流动再变密度。如果这样，不会出现动能突变的状况。或者在变密度操作后立刻更新一下速度；变前计算动能；变后用动能，根据变后的质量，计算新的速度。

或者最粗鲁的方法，用极小的时间步，这样突变也可以稳住。

各位老师，晚上好
我参考了$InterAdsFoam: An Open-Source CFD Model for Granular Media–Adsorption Systems with Dynamic Reaction Zones Subject to Uncontrolled Urban Water Fluxes$和$phaseScalarTransport.C$在compressibleInterFoam中添加了组分方程，个人感觉这个组分方程形式也不是非常复杂，但是在求解过程中会出现局部质量分数增加的情况，如下图所示，想请教一下大家有没有什么建议？

volScalarField Deff("Deff",turbulence.muEff()); volScalarField alpharho ( "alpharho", alpha2*rho2 ); surfaceScalarField alphaPhi ( "alphaPhi", linearInterpolate(alpha2*rho2*U) & mesh.Sf() ); fvScalarMatrix alpha2YEqn ( fvm::ddt(alpharho, Yi) + fvm::div(alphaPhi, Yi) - fvm::laplacian ( fvc::interpolate(Deff)*fvc::interpolate(alpha2), Yi ) == fvOptions(alpha2, rho2, Yi) - fvm::ddt(residualAlpha_*rho2, Yi) + fvc::ddt(residualAlpha_*rho2,Yi) ); alpha2YEqn.relax(); fvOptions.constrain(alpha2YEqn); alpha2YEqn.solve(mesh.solver("Yi")); fvOptions.correct(Yi); Yi.max(0.0); //update the results Yi = Yi*alpha2; // should use this guys as Yi Yt += Yi;

O2:
ede448a4-2b25-4d6a-99de-8f78100ea78b-O2.jpeg
alpha:
820eb593-7b64-4668-8589-a37ace97e464-alpha.jpeg
T:
24cd8ab9-c4ea-463c-b7c5-061e76f80f51-T.jpeg
U:ceee3886-0e70-4e37-9112-801c5f0fc52f-u.jpeg

自己写的代码终于可以出结果了，但是用cgns存坐标数据的时候发现只能存一半，很奇怪，代码举例如下
三维的面网格，节点数2x11，因此是1x10个网格, 坐标储存在x,y,z数组中

real, dimesion(2,11) :: x,y,z CALL CG_BASE_WRITE_F(INDEX_FILE, BASENAME, 2, 3, INDEX_BASE, IER) ISIZE(1,1) = 2 ISIZE(2,1) = 11 ISIZE(1,2) = 1 ISIZE(2,2) = 10 ISIZE(:,3) = 0 CALL CG_ZONE_WRITE_F(INDEX_FILE, INDEX_BASE, ZONENAME, ISIZE, Structured, INDEX_ZONE, IER) CALL CG_COORD_WRITE_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, RealSingle, & 'CoordinateX', x, INDEX_COORD, IER) YZ坐标同上

读取时代码为

real :: node(2,11) node = 0.0 call CG_COORD_READ_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, 'CoordinateX', RealSingle, (/1,1/),(/2,11/),node,IER ) write(*,*) node

node的读取结果为

x(1,1) x(2,1) ... x(1,6) 0.0 0.0 ... 0.0

请问有没有大佬能指点一下

@ZY-GONG 谢谢，已经解决的，之前的错误是因为错误声明了压力，重新在终端中运行setFields命令就可以了

@东岳 谢谢！

好久没来了，过来充充电。另外问下东岳老师，之前注册的邮件可以更改吗？

多谢，已更新

报的错都没贴出来，人眼debug啊:135:

从m_lg=0开始试，一项一项往里加东西，看加到哪一项出的问题

@李东岳 在 momentumLookupTableInjection这个粒子注入模型一直无法正确注入到相应的位置 中说：

那个帖子很老了，我看是用在openfoam2.0.0的，相当老了
我还没测试过你这么玩，感觉你可以把算例发到论坛试试。不过我马上要做国际航班了。够呛能给你debug了近期。

谢谢！我已经解决了，自己编了一个injection模型。

@李东岳 谢谢李老师。这个问题我当时解决了，忘记回论坛答复了。我的处理办法是，在并行的时候，分节点node_n把所有壁面信息如壁面摩擦速度、壁面摩擦速度传递给node0，由node0统一处理把全部壁面信息传递给所有分节点node_n，然后分节点分别计算壁距离、y+。思路是这样，写在这里供大家参考吧