Skip to content
  • CFD中文网 定位与规范

    已固定 公告
    101
    0 赞同
    101 帖子
    606k 浏览
    H

    我写 固我在,嘿嘿

  • 这种尖角梯形网格画出来应该是啥样的

    Meshy
    3
    0 赞同
    3 帖子
    3k 浏览
    S

    @李东岳 尖角部位,加一个竖直的翼墙,这样平分后角度好点,但是附近网格不再是水平的,这样对interFoam计算有影响吗?

  • SCI文章从投稿到录用11天接收

    CFD彩虹条
    4
    0 赞同
    4 帖子
    6k 浏览
    C

    @小八 ICER,文章已经找不到了。化工方向的:+1:

  • 红皮书出中文版了!

    OpenFOAM
    3
    0 赞同
    3 帖子
    4k 浏览
    C

    @liujuncfd 我买了,翻译上还是可以的,但是不知道是出版社的原因,还是本身译者矫正工作做的不到位,这本书有一些符号和标志性的错误,比如有的函数是S(x,y,z)书里给写成S(x,x,z)这种的小错误,其他的问题倒是没什么。

  • 关于两相流二维涡旋的问题

    Fluent
    11
    0 赞同
    11 帖子
    15k 浏览
    S

    @l-j刘侃 目前这个现象只是刚刚发现,暂未做深入研究。李老师对流动本身就是三维的论述点醒了我,流动本身就可以分解为平均速度和脉动速度,所以,二维发文章还是慎重点好。

  • 0 赞同
    4 帖子
    5k 浏览

    @如果随风 pointswise可以画

  • 北师大学生性行为调查报告2016版

    CFD彩虹条
    3
    0 赞同
    3 帖子
    6k 浏览
    C

    @Wayne

    我已联系你:sunglasses:

  • 关于旋转圆筒的压力边界条件问题

    OpenFOAM
    1
    0 赞同
    1 帖子
    2k 浏览

    冒昧打扰各位,我最近正在做一个封闭旋转圆筒半浸没气液两相流动的case,速度边界对于侧面和两端面选择的RotationalWallVelocity,压力边界条件,听从Openfoam官方文件的建议

    In the above example, it can be seen that all the wall boundaries use a boundary condition named fixedFluxPressure. This boundary condition is used for pressure in situations where zeroGradient is generally used, but where body forces such as gravity and surface tension are present in the solution equations. The condition adjusts the gradient accordingly).

    将侧面和两端面设置成fixedFluxPressure固定值为0的格式,结果报错为:Unable to set reference cell for field p,Please supply either pRefCell or pRefPoint。结合问题提示和网上查找解决方案,在fvsolution中的pimple加入参考网格或者参考点的压力值,因此我在pimple的子目录下添加最后两行如下所示:

    PIMPLE { nCorrectors 2; nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 0; }

    结果报错如下:
    98414f11-7705-4f4a-a503-bc89ed313159-image.png
    于是选择第二种解决方法,设置参考点的压力,如下所示:

    PIMPLE { nCorrectors 2; nNonOrthogonalCorrectors 0; pRefPoint (0 0 0); pRefValue 0; }

    结果又是报错,问题如下:
    d88d6776-a2f0-4df3-ba30-a1fb49516c6f-image.png
    谷歌搜索解决方案,很多人说需要设置一个压力场的fixedvalue才可以,我于是尝试分两次将两端面或者侧面的其中一个压力边界条件设置成fixedvalue 0,跑是能跑了,但是0.几秒就发散了,而且本人也觉得这种压力边界条件设置方法很不合理且不符合实际情况(明明都是封闭区域的墙啊)。实在找不到可能的解决措施了,麻烦各位指教一下,谢谢啦。

  • 关于用ICEM画圆柱网格平滑过度问题

    Meshy
    7
    0 赞同
    7 帖子
    9k 浏览
    K

    好的好的谢谢大家的意见~!!!!非常感谢大家,根据大家的意见已经调整好了,最近再准备毕业事宜所以没有及时感谢大家!

  • OpenFOAM 模型编译资料求助

    OpenFOAM
    6
    0 赞同
    6 帖子
    7k 浏览
    李东岳

    啊 我以为 @Zhoujq 老铁太忙了 也没好意思问 妥了 !

  • 0 赞同
    3 帖子
    5k 浏览

    @tens 好的,谢谢你,我试一下

  • ANSYS FLUENT UDF学习

    Fluent
    12
    0 赞同
    12 帖子
    18k 浏览
    L

    @l-j刘侃 刘老师这个成本很高吧:chigua:

  • 论坛有人搞矩方法么?

    Algorithm
    22
    0 赞同
    22 帖子
    27k 浏览
    bestucanB

    @东岳 Yu老师的工作重心最近两年主要放在把TEMOM应用到更多的场景,就是推广它。深挖的工作,我们这边团队小,没博士,老师越来越忙……老师是有意向的,只是条件不允许。对于其他矩方法的研究目的更多是为了对比验证,增强说服力,然后可以更好的推广。当然这些都是我猜的,毕竟老师不跟我汇报工作:chigua2:

  • 关于fluent中的DEM模型请教

    Fluent
    3
    0 赞同
    3 帖子
    6k 浏览
    L

    一般人还是不太行的 整这个玩意 其实DEM非常简单了 我原来用matlab直接开发过一个做沉积的
    fluent也有自带的 但是fluent的基本属于MP-PIC这种 针对DEM是coarsed的那种方式了 也挺工程学 未必你是想要的
    你想要的话 还得自己把颗粒整成1/parcel
    fluent可以调用其他lib的 但是有时候不是稳定
    目前DEM的方式很多 伊朗人16年出的CFD-DEM coupling 那本是是CFD-DEM的比较好的专著吧
    DEM挺简单的 要有自学 但是目前这股热是退了 因为工程价值比较低

  • 0 赞同
    4 帖子
    4k 浏览
    李东岳

    主任招博士后不,我给您发个简历

  • 以读取系数的方式修改湍流模型

    OpenFOAM
    1
    0 赞同
    1 帖子
    2k 浏览
    F

    请教各位老师~
    ****目标:***希望给雷诺应力增加高阶项来修改湍流模型,高阶项的形式是固定的(标量对称张量)
    ***想法:***因为雷诺应力本身就是对称张量,其中的标量是随时变化的,机器学习生成的。如果用include **.H的方式来修改湍流模型,则每次需要编译,比较耗时,(因为增加项会不断变化)。所以想通过在constant文件夹中以.H形式存储标量表达式,从而在turbulenceProperties中读取文件,以系数的形式修改表达式,由于系数是实时读取的,因此只需要编译一次。.H文件中就是一个标量表达式,命名为zetai.H
    ****问题:****如果仅有zetai.H,内容如下(I1 I2是已经在湍流模型中计算好了的体标量场):
    8b12a85a-ca5c-48ca-8d0a-0e643d6d5861-image.png
    并在湍流模型中读取:
    fb6a553a-b074-470e-97ef-4ccc9cc304ec-image.png
    如何能将zetai中的内容转换为体标量场进行计算?

  • 编译新边界时报错

    OpenFOAM
    2
    0 赞同
    2 帖子
    3k 浏览
    GarethG

    我已经解决了这个编译错误,主要是自己粗心大意导致的:135: :136: 。
    在***FvPatchFields.C文件中增加了#include surfaceFields.H即能编译通过

  • shm生成边界层没有完全覆盖

    Meshy
    4
    0 赞同
    4 帖子
    6k 浏览
    S

    @samuel-tu
    前辈好,请问交界面网格畸形的问题您是否解决了?我也遇到了类似的问题:135:

  • 两相流模型库朗数问题(二)

    OpenFOAM
    3
    0 赞同
    3 帖子
    4k 浏览

    @bestucan 好的,我去按照这个思路研究研究,非常感谢!

  • 说说关于dpm系统的半透膜的问题

    Fluent
    3
    0 赞同
    3 帖子
    5k 浏览
    L

    因为时间比较紧,叙述不是很明了,现在按照李老师的建议,加上标点。

    一般dpm的碰撞判断基本是BC宏,这个不多说。

    碰撞算法fluent未透露,但是在ANSYS博客中间有一部分提及了,是采用了颗粒轨迹线与面相交的办法。这个办法是属于碰撞检测算法 line-surface。 这个方法非常的复杂。同时在CFD-ONLINE中间,也有人讨论过这个问题。 从我观点,这点Fluent的保密算法不可能是采用最原始的intersection算法;如果是的话,那么巨量的颗粒沉降计算很明显没有fluent这样的高效。

    好了,这是说了题外话,但是也跟接下来的相关。

    也就是当你不利用BC宏的时候,你需要写这样一个算法(line-surface)来判断 dpm 的截留。而这显然是在半透处理上是较为简单的。因为这样的话,你只要一个在某一个面上,甚至是一个解析几何的构成虚拟面上,就可以构造你的半透面。只要颗粒符合你的条件。 当然这种复杂在于你的颗粒将在第几步碰撞? 他能够在途中的每个网格的步数是多少?(假设颗粒在第一次进入一个cell内被捕获,这一般是现实的,因为你基本上只要处理碰撞表面网格的数据)。 那么又有一个很关键的问题: 你需要重新编制一个颗粒更新步的计算。 假设你是最简单的欧拉蛙跳,也显然你需要了解一个parcel的合外力,所以你必须了解整个dpm.h和其他几个.h的头文件。 那么,你将可以直接解决这个半透的问题。

    而如果你不写trap算法。那么很明显,你需要使用fluent所谓的DPM体系的BC宏。这个宏在颗粒碰撞你挂载的壁面的时候执行,记住这个宏执行将中断颗粒时间步。 而他本身是不出网格的。 因为出网格需要赋值新的P0值。(P0属于particle结构体而p属于tracked particle )。也就是说,当颗粒碰到这个面的时候,碰撞数据被计算,碰撞点被获取,颗粒必须更新他的P_v0值,而particle依然是在当前cell,且颗粒的fate只有active, end, abort, escape 四种。 你的trap可以使用abort,但是他的显示数值决于于p->gvtp, 如果要他在head中显示正常, 你还得做一番处理。

    最要的问题是,颗粒将在当前cell被分割,但是这个这个新的P0值依然是在当前cell内。 那么就是说,这些颗粒继续会被这个BC执行。 最终的结果是所有颗粒都被在封闭在这个cell内来回碰壁, 截留, 最后造成完全的incomplete。
    那么如何半透? 现在已经说明了在fluent实现这个半透可能的方法。

    设想一个现实的条件。颗粒一部分会被反弹(他的vo重新被设置) 一部分被捕获(path_abort) 一部分透过(记住你的时间步已经被分割, 而透过是不分割时间步的)。但就是这样最简单的日常物理现象,因为时间步被分割,碰撞算法被黑盒子,依然是十分困难的。(并且fluent的更新过程并不是属于经典的更新过程,时间步的更新会被更小的切分,根据你的设定,并且每一个切分步,你的UPDATA_SCALAR宏都会被执行。也就是说一个理论上的碰撞算法,很可能并不具备真正的可实现性。)