CFD中文网

请问大家，如果a/b 是一条曲线， b有的时候为0，如何在编程中很好的处理？我一般采用 a/(0.000000000001+b)来实现。感觉并不是很好，

@bestucan 真的非常感谢您的建议。对于二维的case尝试了一下，我设置一个不变形的界面，进入filter区域的粒子数量大约是等同于流通量。但我手上的计算资源和代码的优化方面并不乐观，耗时久。如果界面变成可变形的情况，计算时间应该会加倍，时间上等不起。

大家好，我在课题中遇到大量颗粒穿过一个2D可变形界面的case。我需要计算穿过界面的法向速度，来计算流通量。界面通过immersed boundary method 来表达，我的想法是：穿过界面上相邻lagrangian点做矢量为基准求每个颗粒的相应的法向速度。但总的颗粒量500万左右，这个方法凭空想想，需要的耗费很多资源，时间上不允许。
可否有大神能提供一些思路，或者建议，万分感谢！

@zhouxu 感谢您的回复，我尝试了一下，的却可以做出来。这部分本就有代码在计算。

大家好，我尝试复现组里前辈的论文结果。有一张图我不太明白怎么做出来，希望能得到大家的一些帮助和提示。万分感谢！

尝试了一下各个选项，发现这个地方改变就会解决这个问题，那么这个是干嘛用的。

自己琢磨、研究挺久的了，不知道有没有大神可以指点一下

准确来讲，不需要在圆周上产生part中定义的prism，只要网格数量增加即可，还尝试过edit mesh中的Adjust mesh density。总觉得不够好。希望大家指点一下。

几何结构非常简单，几乎等同于三个圆柱的绕流。
要求对于圆的内部和边界层进行加密。
我尝试使用非结构part mesh setup改变指定圆内部的maximum size
和curve mesh setup的圆周上的节点数来实现加密

单个圆柱能够很好的实现，我所要求的加密。内部的网格数量可以达到10000+。
但三个圆周按照我的方法无法实现加密，圆内部的网格数量只有区区2000-4000.不满足要求。

希望得到指点，谢谢。
附上三个圆和单个圆的网格图。

@ibelief 非常感谢

DEFINE_SOURCE(up_y,cell,water_thread,dS,eqn)
  {real source;
   if(C_UDMI(cell,water_thread,4)==-2)
    {source=-C_UDMI(cell,water_thread,1)/C_VOLUME(cell,water_thread);
    }
    else
     source=0;
    return source;
  }

其中C_UDMI(cell,water_thread,1)是剪切应力y方向的分量。

C_UDMI(cell,water_thread,4)==-2这个在前面的各个函数体内都没有出现过，后续会define-init中会用到。实在看不懂了。同时什么作用也不清楚

这个函数体是传质过程的y方向的源项定义宏

最近在学习fluent udf方面的知识。看到一句话“CFD问题只有三维”感觉熟悉又陌生。
尤其是我在fluent中打开二维网格，check的时候显示会是体积单位。非常困惑
同时，我看到早期的udf无论是二维还是三维问题，都会选择float三维数组引入坐标。
据说如果二维问题只float二维数组会出现错误，有没有理解比较深刻的大神指导一下。