CFD中文网

@yingqing 不知道为啥，有一次成功了，在弄就不能成功了。节点的位置进行了设置，但网格的节点与我设置的节点不重合。就像这样
177924ad-d438-468f-8cf5-a612fe049465-image.png

已经解决。
（在一个视频中提到这一点）网格统计问题有三种情况：
第一种，是显示问题的网格数目很少，就几个，占总体的0.0002%。对于这种情况直接更改限定值。
第二种，是显示问题的网格数目很多，占总体的超过10%，甚至更多，对于这种情况需要更改MeshData的参数设置，重新对网格进行绘制。
必须保证，不能出现负网格。

自己写的代码终于可以出结果了，但是用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

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

@vonboett

if I may implement two kinematic clouds in DPMFoam, one that accounts for particle collisions (larger stones but smaller than gridsize) and one that only accounts particle-wall collisions (many gravel particles). This way, I could model erosion..."

Dear Albrecht,

Thanks for your comments. The E-L algorithmn is expensive since it needs to track the lagrangian particles. So it depends on your particles' number. From what you've written, it looks you are on your right way doing this.

D

@羽之下 谢谢

@hxy970122 您好，请问这个问题您解决了吗？

有没有人用OpenFOAM里面的滑移网格AMI功能的，想了解一下AMI的并行做到什么程度了？对网格面的分区划分有限制吗？可以支持任意网格并行分区吗

@星星星星晴 明白了.我这个方法确实直接套用不上.并且碰撞模型从同一个cell扩展到周围cells之间的碰撞,还会有其他的问题.我再想想,谢谢了.以后常联络!

@李东岳 谢谢:threaten: :qinqin:

@东岳 :huahua: :huahua: :huahua: :huahua: :huahua: :huahua:

WeChat8db3756fb66113a48305c51893ed1215.png

模型还是很简单的，就是一个柱体+方形，但是我用icem画必须要用两个O grid，网格过渡很不好（至少我看起来很难受），不知有什么其他思路？

@haining-luo
我现在用ansa处理极复杂三维模型，保存.mesh格式面网格，用fluent meshing或star ccm生成多面体网格或切割体网格，挺方便。
如果模型太复杂，搞科研可以用icem，工程上不太建议用ICEM生成非结构网格。

各位前辈好，我最近准备从其他版本移植一个动网格功能时发现新旧版本之间motionSolver类在调用declareRunTimeSelectionTable函数时存在不同，旧版本中argList使用了Istream，新版本改为了Iodictionary。这导致旧版本的派生类移植到新版本遇到问题，如果直接改新版本的话又会导致其他派生类出现问题。因此可否在基类motionSolver中调用两次declareRunTimeSelectionTable，argList分别使用Istream和Iodictionary？
事实上我在polyPatch.H中看到了调用两个declareRunTimeSelectionTable的例子：

TypeName("patch"); //- Debug switch to disallow the use of genericPolyPatch static int disallowGenericPolyPatch; // Declare run-time constructor selection tables declareRunTimeSelectionTable ( autoPtr, polyPatch, word, ( const word& name, const label size, const label start, const label index, const polyBoundaryMesh& bm, const word& patchType ), (name, size, start, index, bm, patchType) ); declareRunTimeSelectionTable ( autoPtr, polyPatch, dictionary, ( const word& name, const dictionary& dict, const label index, const polyBoundaryMesh& bm, const word& patchType ), (name, dict, index, bm, patchType) );

但是不清楚是否是为了实现上述的功能？

@李东岳 好的，谢谢李老师耐心解答

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到OF里面是会报错的，错误是：--> FOAM FATAL ERROR:
face 11528151 does not belong to cell 210100. Face owner: -1 neighbour: 210099

From Foam::cellShape Foam::extrudedQuadCellShape(Foam::label, const labelList&, const faceList&, const labelList&, const labelList&, Foam::label, Foam::faceList&) in file extrudedQuadCellShape.C at line 105.

FOAM aborting
感觉好像是在z方向拉伸出来的面没有识别到具体的cell中

使用UDS计算电势和流体所受到的电场力，但是在边界定义标量数值的时候默认单位是1吗？怎么弄成电势单位？

现在的个性签名是胖猫说的，当时的背景是我创立各种网站，翻译OpenFOAM用户指南，做很多人觉得毫无意义的事情。于是胖猫说

你去了那么多地方，走了那么多的路，你要学会的就是包容和接受这个世界的一切观念

博士联合培养期间我去了意大利，博士后我去了德国，中途为了见各种CFD大牛一直在欧洲奔波。这句话确实很适合我。我曾经跨过山和大海，也穿过人山人海，去了很多地方，应该包容和接收其他人的一切观念。

今天，或许应该改一改个性签名咯？:baobao:

有些站友提出问题之后，其他站友们给了追问或者讨论之后，就没有消息了.....
方案能不能用，倒是回复一下啊
完全不能理解这种人提出问题就再也不管的人，不靠谱

@李东岳 这点我清楚，实际上这么多参数里也就总摩擦系数对不上，其他的基本上全对上了。LES的话得在文献里是说在入口用vortex method施加扰动，而且模拟的是有肋的，所以基本上不担心涡的问题。流向周期性的模型也可以在RANS模拟之后用TUI命令施加扰动，之前试过了涡能长期存在，不过那个case没算完结果对不对就不清楚了。