CFD中文网

@hangsz
高手！

各位专家好：
在多孔介质模拟的过程中，我发现了一个问题，虽然在发表文章中省略了。但是确实是一个直观的问题，我的模型，RANS方法中，不管采用哪一种，在多孔介质的近壁面处的拟合的速度与实测值总是具备一定的差距。这种差距我分析认为是在较高阻力降的多孔介质环境下，气流在流经多孔介质表面时，一部分流体会因为高阻降而产生流向变化，形成比较复杂的近壁面流场。而在这个近壁面流场中，因为Fluent默认是采用的多孔介质域或者是多孔跳跃面，应该是对这层近壁面没有处理的，依然按照充分发展层处理，造成这层近壁面的流速（甚至是流向）与实测值都有较大区别。
虽然能够较好的拟合实时压降（Endo based equations），但是如果在近壁面流场不符合实际的话，那么很显然对于最后的表面覆尘的预算难以符合实际情况。不知道大家对于这种多孔介质的近壁面流是否有过研究，或者能够有比较好的解决思路？欢迎大家提出宝贵的想法和建议，感谢，再次感谢！
下图为提出新模型对于滤管外壁的模拟和实测结果对比：
0_1544421734755_compare.png

请问openfoam输出槽道底面摩擦系数的参考长度和参考面积如何设置。

带摩擦的运输方程是 u_t = lambda * u_x = eps * u_xx, 通过变化, 我们得到了, 0_1545168405877_3bd01f52-6748-4ece-bac9-e55eeef96e15-image.png ,其中D是导数矩阵, D2 是二阶导数矩阵
然后, 我们对M进行求特征值(EW), 得到了这种图片, 我想知道这是为什么呢?
0_1545168538892_dba5eac1-7474-4ce6-b089-77d3d8903c03-image.png

@fireztw 嗯，那看起来还需要两相流求解器。
不清楚，如果哪位大佬有类似的sci可以发来瞅瞅

请问像这种狭窄处要划分边界层，怎样处理能保证网格高质量呢

我导收了一个湘潭大学的搅拌器的英文sci。实验+cfd，用openfoam搞的。我针对算例设置这一段提出了1个不痛不痒的问题，不能一个不提啊。目前已返回我导。估计明后天就返回编辑吧？哈哈。我导没有抵制，不过提了2页问题...anyway，问题不大。:sunglasses:

@东岳 感谢东岳老师 目前我还在完善过程中 这个模型目前还存在相当多的问题 这个模型主要是集中在工程应用上 他的大部分数据是利用单管试验来获取的 所以是一个靶向非常集中的宏观模型 对其他工艺的适用性可能非常一般 再次感谢李老师的提议 我想在更进一步完善这个模型之后 :xinxin: :xinxin: 对一些研究多孔过滤的人才能真正起到一定的意义

很遗憾，木有啊:sad:

DEM加速比很高
CFD好像一般

就是我实际网格模型是左边的哪一个，是右边的模型简化的结果，
现在想要任意找一个一个过对称轴的平面，不知道如何选取在这个二维fluent中如何选取

fluent中的压力出口边界条件，资料上写着

压力出口边界条件在流场出口边界上定义静压，而静压的值仅在流场为亚音速时使用。如果出口边界上流场达到超音速，那么边界上的压力将从流场内部通过插值得到。
其他流场变量均从流场内部插值获得。

这里说的插值是怎样的插值，导数为零？还是二阶导数为零？

@李东岳 感谢李老师，我试试看能不能画出来

谢谢老师，我找出来看一下

你不觉得大材小用了，浪费那么多的时间搞科研，不说为人类的进步献一份力把，至少让其学有所用啊，他的那些东西中学根本用不上。 人生价值何在啊。何必当年读博士。

目前这个求解器在OpenFOAM里面已经有了：shallowWaterFoam

TIM图片20190509144012.png TIM图片20190509144004.jpg
第一幅图是edem中的传送带例子，第二幅图现在算出的静态（支撑颗粒的板不动）结果，其实我们模型需要板动起来，类似第一幅图的传送带，最终得到颗粒水分、温度分布，以及后期还会加上化学反应。希望大家给点意见，可不可行，然后有没有类似的案例。先谢！

@l-j刘侃 好的，谢谢刘老师！

@ibelief 在 颗粒表面化学反应的文献 中说：

基于简单碰撞理论煤粉燃烧动力学模型的研究

非常感谢

1.我猜是symmetry这个边界条件，试一下改成wall，然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查，外部软件的网格节点编号和blockMesh的逻辑不一样。
2.可能是前后的empty边界问题，你可以在starccm里生成3D网格后，重新在openfoam里extrude来生成前后empty边界
类似这样
constructFrom patch;
sourceCase "../w3-d_hc1-3";
sourcePatches (symFront);

// If construct from patch: patch to use for back (can be same as sourcePatch)
exposedPatchName symBack;

// Flip surface normals before usage. Valid only for extrude from surface or
// patch.
flipNormals false;

//- Linear extrusion in point-normal direction
extrudeModel linearNormal;

nLayers 1;

expansionRatio 1.0;

linearNormalCoeffs
{
thickness 1;
}

// Do front and back need to be merged? Usually only makes sense for 360
// degree wedges.
mergeFaces false; //true;

// Merge small edges. Fraction of bounding box.
mergeTol 0;